Chemical modification of 2D materials using molecules and assemblies of molecules

ABSTRACT In this article, we discuss the recent progress on controlled structural modification of 2D materials by means of molecular functionalization, with a focus on scanning probe microscopy techniques for their characterization. For many practical applications of these novel materials, it is necessary to tune their electronic and optical properties, and molecule-based functionalization is a powerful approach to reach this. We discuss recent covalent and non-covalent approaches, for functionalization of graphene, transition metal dichalcogenides, black phosphorus, and hexagonal boron nitride. Nanostructuring approaches and their impact on 2D materials’ properties are highlighted. Graphical abstract

flexibility lead to potential applications in touch screens, capacitors, spintronic devices, fuel cells, batteries, sensors, transparent conductive films, high-frequency circuits, toxic material removal, and flexible electronics. A limitation of graphene arises from the Dirac-like linear dispersion, which on the one hand limits the range of applications due to the absence of a band gap, but on the other hand, can lead to unexpected physics with a lot of fascinating properties to explore. Thus, it was soon realized that tuning the electronic band structure of this material is an important target. For this purpose different methods to functionalize graphene by chemisorbed or physisorbed molecules [2][3][4][5], substitutional doping [6], intercalation of molecules [7] and metal atoms, and other reactions such as hydrogenation and fluorination were explored.
Modifications of graphene by physisorbed organic molecules or chemisorption of reactive molecules via covalent bonds have resulted in the overall modification of the electronic properties of graphene. Functionalization of graphene by physisorbed molecules (randomly adsorbed or self-assembled) has been shown to induce n-type or p-type doping. Due to the presence of delocalized π electrons, graphene is highly susceptible to adsorbates on the surface. Therefore, organic molecules with the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) at a suitable position with respect to Fermi level of graphene can be used to dope graphene [8,9]. Chemisorption of reactive molecules occurs via formation of covalent bonds between graphene carbon atoms and molecules, resulting in transforming sp 2 carbon atoms of graphene in sp 3 carbon atoms at the site of molecule attachment [9,10]. Chemisorption of molecules is more stable when compared to physisorption but at the same time can cause detrimental effects on the novel properties of graphene. Covalent functionalization of graphene is considered as the best possible way to robustly modify the electronic band structure of graphene thus promoting its application in future nanoelectronic devices. For covalent functionalization of graphene different organic functional groups such as diazonium [10][11][12], and maleimide [13] have been used. Among these, covalent functionalization using diazonium chemistry is most thoroughly studied.
The growth of graphene research has led to exploring other layered 2D inorganic materials such as transition metal dichalcogenides (TMDs), black phosphorus, and hexagonal boron nitride. These can be thinned down to monolayers. The atomically thin nature of these materials results in unique structural and physical properties and provides access to novel fundamental properties which are not present in their bulk counterparts. Synthetic methods such as vapor phase methods, wet chemical synthesis, and liquid phase exfoliation are the most prominent methods for growth of 2D TMDs [14,15]. Bulk TMDs have found applications in energy conversion and storage. The electronic band structure is what makes TMDs interesting, whereby the bandgap of TMDs is highly dependent on the number of layers in a TMD stack and also on the constituent atoms as well as the coordination of these materials [16], which is a basic fundamental property not found in graphene. The combination of atomic scale thickness, layer dependent electronic band structure, strong spin-orbit coupling, and superior mechanical properties makes these materials interesting for applied as well as fundamental research.
However, for developing practical and reliable high performance (opto)electronic devices based on 2D TMDs, a controlled tuning of their electronic properties is necessary. Doping of 2D TMDs is important for minimizing the Schottky barrier at the interface of the TMD and metal contact to reduce the charge carrier recombination rate at the interface. Conventional methods of doping such as ion implantation and substitutional doping are not suitable for 2D TMDs as they can create defects in the layer, having a negative impact on the electronic properties. A number of different approaches have been proposed to overcome this hurdle which are discussed in a recent review on doping engineering of TMDs by different methods [17]. Among these approaches, the most promising are the ones involving physisorption and chemisorption of organic molecules with functional groups having strong interactions with the TMD layers.
Charge transfer interactions and dipolar interactions are two prominent factors that can result in doping of TMDs. Charge transfer doping is mostly dependent on the redox potential of molecules with respect to the Fermi level of TMDs, whereas molecules on TMDs with dipolar interactions create an electric field that modulates the local charge carrier density of TMDs. Creating interfaces of TMDs with molecules can, therefore, improve the device performance characteristics by controlling the doping (n-type or p-type). Self-assembly of molecules with fixed a dipole can create a non-zero field, therefore, creating macroscopic changes to the (opto)-electronic properties of TMDs.
In addition to physisorption, also chemisorption approaches are very promising. A lot of focus has been devoted towards improving the dispersibility of TMDs in solutions by covalent functionalization with molecules to develop methods for large area processing of TMDs. Covalent functionalization of MoS 2 , the most popular TMD, has been carried out by number of approaches such as thiol addition [18][19][20][21][22][23][24], addition of dibenzene thiophene [25], aryl diazonium molecules [26,27] and organohalides [28,29]. Pristine MoS 2 , either CVD grown or chemically exfoliated, always has structural defects in the form of sulfur vacancy sites along with other point defects. These defects act as scattering centers for the majority of charge carriers in MoS 2 . Covalent functionalization by sulfur containing groups or diazonium groups can result in passivating the surface defects by grafting at such defect sites. Additionally, for 2D TMDs covalent functionalization is a facile and robust method of doping compared to traditional ion-implantation methods. Molecular chemistry approaches for covalent and well as non-covalent functionalization of TMDs is comprehensively reviewed by Bertolazzi et al. [30].
In this review, we will focus on recent work on non-covalent and covalent functionalization of substrate-supported graphene and other 2D materials, which result in a modulation of their properties. We mostly focus on recent work that involves the use of scanning probe microscopy (SPM) based techniques to visualize the supramolecular assemblies as well as covalently grafted molecules on the 2D material. SPM techniques provide access to sub-nanometer level resolution, thus, any change in the molecular and 2D material properties upon interaction with the 2D material can be directly visualized. We demonstrate that SPM-based methods, such as Kelvin probe force microscopy (KPFM) and scanning tunneling spectroscopy (STS), can give insight into the charge transfer interactions between molecules and the 2D material. In section 1, the impact and importance of primarily self-assembled molecular layers are discussed. Covalent functionalization is discussed in section 2.

Non-covalent functionalization of 2D materials
Given the ultrathin nature of 2D materials, surface functionalization using molecular species is a promising approach for tuning the properties of these materials. Non-covalent functionalization of 2D materials provides a way to tune the electronic and physical properties without disrupting crystal structure. Physisorbed molecules interact with 2D materials via van der Waals (vdW) interactions or electrostatic forces, which can be used to controllably tune doping levels. Non-covalent functionalization of 2D materials can be carried out by physisorbed molecules, forming selfassembled structures or adsorbing randomly on the top (or bottom) surface of the 2D material.
Self-assembly of organic molecules on highly oriented pyrolytic graphite (HOPG) into self-assembled molecular networks (SAMNs) [31], which can be considered as a model system for molecular assembly on graphene, has been thoroughly studied experimentally and theoretically, resulting in a good understanding of molecular self-assembly on planar surfaces. Similarly, functionalization by SAMNs is being considered a promising approach for other 2D materials. In the case of graphene, a number of reviews published in the past few years have summarized different approaches towards non-covalent functionalization by self-assembled molecules [2][3][4][5].
Thus, for graphene, we will only focus on recent studies on the noncovalent functionalization with molecules leading to a charge transfer. Charge transfer interactions result in a redistribution of electrons between the molecules and the underlying 2D material, thereby tuning the electronic band structure, i.e. doping the 2D material.
Initial work on molecular functionalization of graphene was focused on planar aromatic molecules such as PTCDA deposited via sublimation on epitaxial graphene grown on SiC (EG/SiC) Figure 1(b)). The self-assembly of PTCDA molecules has been studied by different research groups on epitaxial graphene using scanning tunneling microscopy (STM) [32][33][34] in ultra-high vacuum (UHV) conditions at room temperature as well as at 77K. An important aspect of PTCDA SAMN formation on graphene is that the PTCDA domains don't get discontinued at EG-SiC step edges but are in full compliance with graphene topology and propagate over step edges as shown in Figure 1(c). Slight n-type doping of graphene by PTCDA molecules was confirmed by Huang et al. [33] and Lauffer et al. [34] at low temperature in separate studies by scanning tunneling spectroscopy (STS) measurements (Figure 1(f)). However, such n-type doping is not observed upon measurements at room temperature. Recently, H. J. Karmel et al. studied the structural and thermal stability of PTCDI molecules on a monolayer of G/SiC using UHV-STM and temperature-dependent X-ray reflectivity (XRR) [35], showing a high thermal stability compared to PTCDA of up to 311.15 ± 10 K. Recent theoretical studies have suggested strong π-π interactions of these molecules with graphene, which potentially can result in doping as well as the introduction of a small bandgap.
Due to different intermolecular interactions, TCNQ and its fluorinated counterpart F4-TCNQ self-assemble differently on graphene, as observed by Barja et al. and Stradi et al. in separate studies using high-resolution STM measurements [36,37]. With an electron affinity of 5.24eV the LUMO of F4-TCNQ is below the Dirac point of graphene, making it a very strong p-type dopant, as shown by Chen et al. for epitaxial graphene on SiC [38]. Using high-resolution photoemission spectroscopy (PES) a clear evidence of electron transfer from molecules to graphene was obtained.
Metal phthalocyanines (MPc) are another class of molecules with an extended delocalized π-electron system, which have been used to functionalize graphene as well as other 2D materials. They are a class of highly stable molecules consisting of a macrocycle of alternating carbon and nitrogen atoms with a single metal atom in the center. Over the past few years, a number of MPc molecules with different central metal atoms have been studied on graphene grown on various substrates [39][40][41][42]. Thereby, MPc molecules have shown variable adsorption behavior on graphene depending on the level of coverage, peripheral functionalization and central metal atom as revealed by a number of separate STM studies. Mao et al [39]. demonstrated the formation of a Kagome lattice of H 2 Pc, NiPc, and (t-bu)4-ZnPc molecules which follows the moiré pattern of underlying graphene. These Kagome lattices had different spins depending on the metal atom of the MPc molecule, which demonstrates the prospect of using these molecules to tune the spin interaction of graphene. A slight doping of graphene by a single layer of FePc molecules deposited in UHV conditions was revealed by Scardamaglia et al. using photoemission spectroscopy techniques [43].
Additionally, organic molecules with long alkyl chains with functional groups are also used for controlled tuning of the doping levels of graphene. Li et al. showed that SAMNs of oleylamine (OA) molecules on graphene on SiC can result in n-type doping [8] via electron transfer from the NH 2 groups, while the alkyl chains direct the self-assembly of these molecules. Homogeneous SAMNs were observed by STM and AFM measurements in ambient conditions. The n-type doping caused by OA molecules was expected to modify the electronic properties of graphene and confirmed by electric transport measurements carried out on back-gated graphene field effect transistors (FETs), which showed a shift of the charge neutrality point towards more negative values (−16 V and −28 V) after the first and second deposition cycles of OA molecules.
Recently, Phillipson et al. demonstrated the tuning of the Fermi level by varying the chain-length of alkylated amine molecules forming SAMNs on graphene [44]. Control of the surface concentration of amino groups resulted in tunable n-type doping of graphene as shown by Raman spectroscopy measurements and charge transport characteristics measured by FETs ( Figure 2). Raman spectroscopy measurements carried out on octadecylamine (ODA) and nonacosylamine (NCA) functionalized graphene on SiO 2 are shown in Figure 2 e)-h). A higher n-type doping from ODA compared to NCA molecules was observed based on a larger redshift of the G peak position and a higher increase in the I(2D)/I(G) ratio. The charge carrier concentration found was 1.59 times higher with ODA compared to NCA, which is in good agreement with the ratio of alkyl chain length of NCA and ODA (1.6 times longer). The doping mechanism was explained by the formation of a dipole layer between graphene and the NH 2 groups of ODA and NCA molecules. The transport measurements by FET also showed a 1.52 times higher doping for ODA functionalized devices compared to NCA functionalized ones.
Prado et al. showed p-type doping of graphene resulting from noncovalent functionalization via the formation of SAMNs using phosphonic acid molecules with long alkyl chains, namely octadecylphosphonic acid (OPA) and tetradecyl phosphonic acid (TPA), in ambient conditions [45]. The phosphonic acid headgroup is an electron acceptor, leading to p-type doping of graphene. The molecules were deposited on graphene flakes exfoliated on SiO 2 layers. AFM images showed well-organized monolayers consisting of OPA and TPA dimers, respectively. The alkyl chains of the molecules were arranged in a zigzag alignment aligned along a highsymmetry direction of graphene. Raman spectroscopy of the OPA functionalized graphene showed a shift of G peak (ν G = 1586 cm −1 ) towards a higher wavenumber compared to the pristine graphene flake (ν G = 1593 cm −1 ). A decrease in FWHM of the G-peak after functionalization and decrease in the I(2D)/I(G) ratio was also observed, which was attributed to p-type doping of graphene by the OPA molecules. However, no information was given concerning doping levels of graphene functionalized by TPA molecules.
In all of the studies discussed above, tuning of graphene's properties was achieved by varying the type of functional group and density of the molecules. Recently, P. Samorì and coworkers demonstrated the formation of periodic potential on graphene using SAMN formation. The surface potential of graphene was modulated by the change in dipole of the selfassembled molecules [46]. Although this is an interesting aspect for the application of graphene for advanced electronics, it does not offer precise control, which future solid state devices based on graphene will require. One approach towards this goal is to include specific functional groups in self-assembling molecules, which can be modulated by external stimuli, also demonstrated by P. Samorì and coworkers [47]. They showed that the transport properties of graphene and MoS 2 can be modulated by using photochromic molecules. In this study, the authors have used a spiropyran (SP) derivative with an 18 carbon long alkyl chain for self-assembly on 2D materials. SP molecules can be reversibly switched to merocyanine (MC) forms, the latter having a large dipole moment. The conversion from the SP to the MC form takes place upon irradiation with ultraviolet (UV) light as shown by the schematic in Figure 3(a). The reversible reaction can be  either thermally or by irradiation with visible light. Upon irradiation with UV-light, the change in the molecular dipole moments causes a shift in the work function of the 2D material. The isomerization from SP to MC upon UV irradiation was confirmed with STM in ambient conditions, as shown in Figure 3(e-f). The SP isomer is characterized by a lamellar structure in which different rows of alkanes lie flat on the surface and are separated by bright fuzzy regions, while STM images after isomerization show fully immobilized MC groups lying close to each other and the alkyl chains form an interdigitated lamellar structure. Back-gated FET, fabricated using graphene, showed a strong shift in the charge neutrality point upon irradiating SP/graphene with UV light, suggesting a strong n-type doping with an electron density reaching = 4.4 × 10 12 cm −2 . This effect was reverted back upon irradiation with a green light. A similar effect was observed for SP/MoS 2 hybrids. Irradiation at a precise location of the SP covered graphene leads to SP→MC isomerization only at that particular region, resulting in mixed transport characteristics of the devices.

Non-covalent functionalization of TMDs
Following different methods developed for the functionalization and tuning of the electronic properties of graphene by self-assembled molecules, similar approaches for other 2D materials such as TMDs and hBN were also developed. Bulk MoS 2 , one of the most studied material from the TMD family was used a number of times as an alternative substrate to study molecular self-assembly. Here we will discuss some of the studies which dealt with the modulation of properties of TMDs by molecules adsorbed on the surface. The pioneering work on the functionalization of MoS 2 by SAMNs was carried out by M. Hara et. al, who used STM to determine the impact of substrates on the self-assembly of the smectic liquid crystal 4ʹ-n-octyl-4-cyanobiphenyl (8CB) at the air/solid interface [48]. These authors have demonstrated that the molecules assemble in an energetically unfavorable manner due to strong interactions with the underlying substrate suggesting epitaxial growth of molecular layers on the MoS 2 lattice. A year later Cincotti et al. demonstrated SAMN formation of dotriacontane (C 32 H 66 ) molecules at the liquid-MoS 2 (and MoSe 2 ) interface [49]. These initial studies were mostly focused on using TMD as substrates for a better understanding of molecular self-assembly forming SAMNs and the impact of the lattice structure of the underlying surface.
An interesting approach towards tuning the physical and chemical properties of TMD is using self-assembled monolayers (SAMs) [50,51]. In contrast to SAMNs, where the molecules are physisorbed on the surface, typically with their long axis parallel to the surface, SAMs are formed by molecules consisting of (i) a functional anchoring group such as thiol, silane or phosphoric acid to chemisorb molecules on the target substrate, (ii) a molecular backbone, typically alkyl chains, which can stabilize the SAM via intermolecular van der Waals interactions and (iii) a functional head group, which determines the surface energy and dipole moment after SAM formation, thereby allowing to tune the physicochemical properties of the exposed surface. The interactions among the molecular backbone of the SAMs generally result in high crystallinity, which causes molecules to orient edge on with respect to the substrate plane. This makes it possible for single molecule dipoles to sum up and align, resulting in a sufficient electric field to interact with TMDs. In addition to acting as a substrate for the deposition of a SAM [49,52], TMDs can also be deposited on SAMmodified substrates [53]. When TMDs are in contact with a SAM, the dipole moments of the end functional groups affect the Fermi levels of TMDs. An upward or downward shift of the TMD Fermi level may occur depending on the functional group, resulting in doping of TMDs.
Molecular doping of TMDs, such as MoS 2 and WSe 2 , by self-assembled monolayers (SAM) with a silane anchoring group has been demonstrated to result in the tuning of the optoelectronic properties in a number of experimental studies [53][54][55][56][57]. Li et al. studied the doping behavior of different molecules with a silane end group, such as octyltrichlorosilane (C 8 H 17 SiCl 3 )(or OTS); 3-(trimethoxysilyl)-1-propanamine (NH 2 (CH 2 ) 3 Si(OCH 3 ) 3 )(or APTMS); and trichloro(1H,1H,2H,2H perfluorooctyl)silane (CF 3 (CF 2 ) 5 CH 2 CH 2 SiCl 3 ) (or FOTS), with different dipole moments to modify the work function of MoS 2 [53]. These molecules form monolayers on SiO 2 and therefore can act as an interface between MoS 2 and SiO 2 . The self-assembled monolayers of these molecules were characterized using AFM and XPS. KPFM, Raman Spectroscopy and transport characteristics by FET were used to determine the modulation in the doping levels of MoS 2 . Depending on the dipole moment and polarity, these molecules showed different hole and electron doping abilities. MoS 2 on OTS SAMs showed very little change in Raman spectra compared to pristine MoS 2 flakes with the same thickness on SiO 2 , suggesting very little charge transfer between MoS 2 and the octyltrichlorosilane (OTS) molecules. However, for substrates functionalized by APTMS and FOTS, blue shifts of 1.6 cm −1 and 1.7 cm −1 , respectively, were observed. Spatial charge distribution measurements to determine the charge transfer between MoS 2 nanoflakes and SAM-modified SiO 2 surface were carried out using KPFM, showing a shift of the Fermi levels of monolayer MoS 2 by more than 0.45 eV due to the interactions with FOTS and APTMS, respectively. Figure 4 shows a clear contact potential difference (CPD) between different layers of MoS 2 flakes deposited on the SAM-modified SiO 2 surface. A shift in the Fermi level, calculated from CPD measurements, confirms the charge transfer process between the molecules and MoS 2 .
These studies showed the possibility to induce unique characteristics in TMDs by using molecules of interest as SAMs, as recently demonstrated with SAMs composed of thiolated photochromic azobenzene (AZO) molecules (E)-6-(4-(phenyldiazenyl) phenoxy)hexane-1-thiol (HS-C6AZO) and (E)-6-(4-((4-chlorophenyl) diazenyl)phenoxy)hexane-1-thiol (HS-C6AZO-Cl) as electropositive and electronegative SAMs chemically attached to a gold substrate [58] as shown in Figure 5(a)). These molecules were mixed with spacer molecules (6-(2-mercapto)-1-hexanol) to form mixed monolayers. Azo molecules are photochromic and can be switched between a cis and a trans isomer upon irradiation of light at a different wavelengths. In this work, photoswitching of molecules between the cis and trans state resulted in a modification of surface energetics of single MoS 2 layers deposited on top of the SAMs. The change in the surface potential of MoS 2 , resulting from photoswitching between cis and trans state, was demonstrated using KPFM. Photoluminescence (PL) measurements revealed that for n-type phototunable doping the PL intensity between trans-mSAM and cis-mSAM is increased ( Figure 5(c)), while the opposite behavior was recorded for phototunable p-type doping. This Au/ AZO-SAM/MoS 2 stack was also used to demonstrate the application in photo switchable vertical diodes. In these diodes, the current flowing vertically through the Au/SAM/MoS 2 stack is measured by placing a conductive AFM tip in contact with MoS 2 . Current rectification behavior was observed when the AZO-SAM is in the trans state, which is suppressed when the AZO-SAM is isomerized to cis state, thus demonstrating the possibility to reversibly tune the energetics of TMDs using molecular capabilities such as photoswitching.
An n-type doping of MoS 2 nanoflakes on SiO 2 /Si using SAMNs of physisorbed oleylamine (OA) molecules was recently reported as an efficient way to reduce the contact resistance between metal contacts and MoS 2 in FETs [59]. AFM measurements revealed that OA molecules were self-assembled on the MoS 2 surface in the form of lamellar domains as shown in Figure 6(a,d). Raman and PL spectroscopy along with electrical characterization using a four probe back-gated FETs showed clear evidence of n-type doping of MoS 2 resulting from OA molecules shown in Figure 6 (c,f). An increase in I ON by a factor of 1.9 suggested an increase in the carrier concentration arising from doping by amine groups. The increased carrier concentration was also attributed to the negative dipole moment of the amine groups that influence the tunneling of interface-trapped carriers to the conduction band of the MoS 2 film.
Benzyl viologen (BV) is one of the most used molecules to understand and study the impact of doping on properties of MoS 2 . Kiriya et al. have studied confocal absorption spectra for doped and undoped MoS 2 samples [60]. In this study, the authors used BV to induce n-type doping in MoS 2 . Air-stable benzyl viologen was adsorbed on the surface of few-layer MoS 2 , which resulted in a stable charge transfer doping. The charge transfer doping between BV molecules and MoS 2 was probed by using the functionalized MoS 2 in FET devices in back gate configuration and by Raman spectroscopy of the functionalized MoS 2 . After doping with BV molecules, the MoS 2 FET devices did not show gate voltage dependence of the on current (I DS ), which was increased by an order of magnitude (I DS = 1.3 × 10 − 4 A at V GS = 40 V). Poor gate voltage dependence and the increased current density are clear indications that MoS 2 was strongly doped by BV molecules. Devices of doped MoS 2 retained the observed characteristics with minimal changes even after 9 days of exposure to ambient conditions. Reversing the functionalization of MoS 2 by immersing the devices in toluene lead to increasing the gate voltage dependence of the current, similar to the curve obtained for pristine MoS 2 FETs. This shows that BV molecules on MoS 2 form air stable charge transfer complexes, which disintegrate upon immersing the MoS 2 layers in toluene. Raman spectroscopy carried out on BV functionalized MoS 2 showed a shift in the A 1g vibration modes towards lower wavenumbers compared to nonfunctionalized MoS 2 as a result of high electron concentration. These studies suggest that physisorbed molecules can potentially be used for airstable and reversible functionalization of few-layer MoS 2 .
For tuning the properties of MoS 2 , several molecules which have been thoroughly studied on graphene were also used. F4-TCNQ was used as a p-type dopant of MoS 2 in a couple of studies [61,62]. Wang et al. [62]  recently reported that charge transfer from MoS 2 to F4-TCNQ leads to the formation of van der Waals interface when the molecules are deposited on MoS 2 surface in vacuum conditions (Figure 7(a,c)). The molecules lead to p-doping of the MoS 2 layer along with a depletion of residual electrons from the sulfur vacancies, resulting in bringing the onset voltage value of the MoS 2 FET close to zero, while no reduction in the mobility of MoS 2 FETs were observed after the deposition of F4-TCNQ molecules. Scanning Kelvin probe force microscopy studies further verified the charge transfer mechanism (Figure 7(b,d)). The authors have also shown that such interesting charge transfer effects in van der Waals interfaces between organic-inorganic hybrid systems can have potential applications for sensing applications, as illustrated by using this hybrid system in FET configuration for NH 3 sensing. F4-TCNQ was also deposited on MoS 2 flakes by solution-based techniques. Mouri et al. [61]. studied doping characteristics of molecules including F4-TCNQ drop casted on MoS 2 flakes. They have shown that depositing F4-TCNQ molecules on MoS 2 leads to a strong enhancement in PL intensity of the MoS 2 . This enhancement was the result of switching the dominant PL process, i.e. the recombination of the negative trion which usually occurs in as-prepared 1L-MoS 2 , to the recombination of the exciton instead.
Functionalization of MoS 2 based photodetectors by organic dyes, such as perylenediimides (PDIs) and porphyrins (TPP), resulted in an increase of photocurrent and responsivity up to 3 times higher compared to pristine MoS 2 based photodetectors [63]. AFM measurements confirms the formation of 5-36 nm thick layer of TPP molecules on MoS 2 flakes. In another work, perylene-based molecules were used as a seeding layer for atomic layer deposition of Al 2 O 3 gate dielectric.
Non-covalent functionalization of hBN and black phosphorus Its large electronic bandgap (~5.2 eV) makes hBN a good decoupling layer to study the electronic properties of individual or self-assembled molecules without much influence from the underlying substrate. Molecular functionalization is an effective method to tune the electronic, magnetic or metallic properties of hBN. However, here the presence of the large bandgap is a major drawback of hBN. Willi Auwärter discussed the selfassembly of most of the molecular systems on hBN in a comprehensive review [64].
In the case of hBN epitaxially grown on metallic substrates such as Ni (111), Rh(111) and Ru(0001), strong interactions and a lattice mismatch of the hBN with the underlying substrate result in the formation of a nanomesh. This nanomesh can be used for preferential deposition of single molecules at low coverage and hence allows the study of electronic properties of molecules without any external influence, as reported by Schulz et al. [65]. It was shown by STM and STS that CoPc molecules adsorb preferentially in the nanomesh at low coverage, but at high coverage form larger square-shaped islands. STS measurements on the metal central atom of CoPc gave different spectra for molecules adsorbed in the pore and on wire regions of the nanomesh. STM measurements at different biases resulted in different contrasts from the metal center and carbon backbone depending on the applied bias. Based on STS spectra it was concluded that the LUMO of the molecules in the pores are below the Fermi energy compared to a molecule on the wire, suggesting an ordered charging of the molecule in pores as shown in Figure 8(e,f).
The wide bandgap makes it difficult to study the impact of functionalization on hBN properties. However, photoluminescence and Raman spectroscopy were used to study the vacancy defects in the hBN lattice. Nitrogen or boron vacancies in the hBN lattice can create additional states in the electronic bandgap of hBN, which leads to interesting PL behavior. Molecular doping of hBN is another approach towards tuning the electronic bandgap of hBN, which can be studied using PL and Raman spectroscopy.
Functionalization by wet-chemical treatment of liquid phase exfoliated black phosphorus (BP) has been shown to be successful with electron-poor TCNQ as well as a perylene diimide (PDI) derivative by the group of A. Hirsch [66]. In the case of TCNQ, an electron transfer from BP to the dopant took place, while the interactions between BP and PDI were mainly based on van der Waals forces. These strong interactions between the organic molecules and BP supported the liquid phase exfoliation of bulk BP into a few layer BP flakes. Additionally, the functionalization with PDI increased the resistance of the BP flakes against oxygen degradation. Since the lack of stability and degradation of BP in ambient conditions is a major drawback of this material in few of potential applications, the latter finding prompted a systematic study on the oxidation and passivation behavior of mechanically exfoliated BP by the same group [67]. The oxidation of BP was monitored by statistical Raman spectroscopy and various factors affecting the stability of BP flakes were investigated, i.e. layer thickness, lateral dimensions, and exposure to light. Basically, the stability of fewlayer BP flakes increased with thickness and light illumination led to a faster degradation (Figure 9(b,c)). In this study it was also found that treatment with ionic liquids like N-methyl-2-pyrrolidone (NMP), significantly increased the stability of BP by several days (Figure 9(a,d)).
In the context of stabilizing BP under ambient conditions, Korolkov et al. successfully formed SAMNs stabilized by hydrogen bonds by solution deposition to passivate BP and inhibit oxidation under ambient conditions [68]. The molecules used were trimesic acid (TA) and melamine cyanurate (CA.M), which are known to form supramolecular networks on other 2D materials. AFM images with molecular resolution in ambient conditions showing molecular self-assembly could be obtained (Figure 10  and XPS spectra of passivated and non-passivated BP after several days up to three months showed a clear passivation effect of the molecular layers ( Figure 10(e,h)). Even after 1 month in ambient conditions, atomically resolved AFM images of the BP could be obtained.
In addition to studies concerning the stability of BP, non-covalently modified BP has also been used in FETs. Han et al. showed that surface  functionalization of few-layer BP with potassium (K) can increase the performance of FET [69]. Therefore, mechanically exfoliated BP was used for FET fabrication and K was evaporated onto the devices under UHV conditions. Compared to the non-functionalized devices, the electron doping effect of K resulted in significantly improved electron mobility by over one order of magnitude and a clear bandgap reduction. Complementary in-situ UPS and XPS measurements were carried out with K-modified bulk BP and confirmed the interfacial charge transfer taking place at the K/BP interface leading to the electron doping effect.

Covalent functionalization of 2D materials
The functionalization methods discussed so far to tune the properties of 2D materials involved the physisorption of molecules on the basal plane. This approach, however, suffers from drawbacks in terms of robustness of the process and stability of the molecule-2D material interactions. Due to the weak interactions between physisorbed molecules and the 2D material, external stimuli can easily alter the properties of these systems. In this regard, covalent attachment (or grafting) of molecules on 2D materials provides a more robust approach for controllably tuning the properties of 2D materials. Similar to graphene, TMDs, hBN and black phosphorus have also been covalently modified using diazonium chemistry and other chemical approaches. In this section, we will be discussing the most prominent functionalization strategies for covalent functionalization of these 2D materials.

Covalent functionalization of graphene
Covalent functionalization results in rehybridization of sp 2 carbons to sp 3 carbons by forming bonds with reactive molecules on the basal plane of the graphene lattice. This process is similar to creating defects in the graphene lattice and therefore leads to a modification of the electronic band structure, magnetic, optical and mechanical properties of graphene. Additionally, covalent functionalization can lead to better solubility and processibility of graphene, which is an important aspect for more practical applications. Till today, a number of functionalization methods such as addition of diazonium radicals [10,12,[70][71][72][73][74][75][76], Diels-alder reactions [77][78][79][80], cycloaddition reactions [13,81,82], nitrine addition [83,84] and carbene addition [85] have been developed for the covalent grafting of molecules on graphene basal plane. Among these methods, covalent functionalization of graphene using diazonium chemistry is the most efficient and versatile. Diazonium chemistry can also be extended to other 2D materials as evident from the recent studies discussed later in this section. First, we will discuss different methods for covalent functionalization of graphene and other 2D materials.
Functionalization of graphene by diazonium chemistry proceeds via a radical mechanism, where the diazonium salt is reduced to form an aryl radical, which subsequently reacts with the surface as shown in Figure 11(a). Haddon and coworkers were the first to demonstrate covalent functionalization of graphene using diazonium chemistry [70,86]. They showed the functionalization of epitaxial graphene on SiC by spontaneous grafting of aryl groups via reduction of 4-nitrobenzenediazonium (NBD). FTIRspectroscopy and XPS measurements confirmed the presence of nitro groups on the functionalized graphene. The surface coverage was estimated to be 1 × 10 15 molecules/cm 2 based on cyclic voltammetry measurements. Additional work from Haddon and coworkers used spectroscopy methods to study the impact of covalent grafting on the electronic properties of graphene. Using a combination of Raman spectroscopy with ARPES, the authors demonstrated how the covalent functionalization with NBD affects the electronic properties of graphene. Raman measurements carried out on flakes as well as on epitaxial graphene on SiC showed a pronounced D peak of functionalized graphene (Figure 11 (c)). The D peak is associated with defects in the basal plane of graphene, which in this case was caused by the grafted NBD molecules. ARPES measurements showed a slight transformation of linear dispersion of graphene into massive particle dispersion with the Dirac point shifting to 0.36 eV below the Fermi level as shown in Figure 11(b).
Further, SPM-based techniques have been used to determine and control the locations of the molecules covalently bound to graphene. Hossain et al. used UHV-STM to study the morphology of graphene with covalently functionalized 4-NBD molecules [87]. STM images revealed that certain areas of the 4-NBD functionalized graphene did not show any sign of grafting, which was attributed to a blocking effect of oligomers formed during the grafting reaction as shown in Figure 12(a). The aryl oligomers possibly form as a side product of the reaction between nitrophenyl radicals in solution with already grafted molecules. STS spectra revealed no clear modification of graphene's intrinsic properties in blocked areas, whereas STS on functionalized areas showed a high perturbation of electronic band structure of graphene as shown in Figure 12(d,e). Later on, Greenwoood et al. used STM and Raman spectroscopy to study the grafting of different diazonium salts on graphene and HOPG [12]. A comparison between grafting morphology and molecular structure for 4-NBD and bulky 3,5-bis-tert-butylbenzenediazonium (3,5-TBD) was done ( Figure 12(f,g)). It was demonstrated that the density of the grafted molecules on graphene and HOPG can be controlled by using diazonium molecules with appropriate substituents, i.e. bulky substituents help in creating a monolayer of grafted molecules with high density. STM tipinduced lithography showed the potential of this technique to remove the grafted molecules from the substrate and restore the pristine graphene or graphite surface, respectively.
Detailed studies on the electron transfer between grafted molecules and graphene were summarized by Strano and coworkers in comprehensive reviews [74,88]. Wang et al. showed the influence of the substrate underneath graphene on the reactivity of graphene, which can be exploited to control the diazonium chemistry [75]. Localized grafting can be achieved by patterning the substrate supporting graphene. Raman maps of the D/G ratio and AFM measurements confirmed patterned grafting.
Another approach for patterned grafting of diazonium molecules on graphene was demonstrated by Xia et al. by using long alkyl chains attached to diazonium molecules [89]. CVD graphene was functionalized by 4-docosyloxy-benzenediazonium tetrafluoroborate (DBT) molecules. The long alkyl chains promoted self-assembly of DBT molecules, leading to the formation of a SAMN on graphene (Figure 13(a-d)). The diazonium molecules were then activated electrochemically, resulting in patterned grafting. STM, AFM, Raman spectroscopy and conductivity measurements showed significant modifications in the properties of graphene, thus confirming covalent functionalization.
An interesting method for patterning grafted molecules on graphene is using self-assembled molecular monolayers, as recently demonstrated by Tahara et al [90]. In this method, graphene/graphite was covalently functionalized by few nanometer wide rows of 3,5-TBD molecules using SAMNs of long-chain alkanes as a template. This approach relies on phase separated solution double layer as shown in Figure 13(e,f). A first solution layer contains n-alkane molecules which self-assemble on the surface of graphene/graphite. A second solution layer contains diazonium molecules which diffuse through the first layer upon a potential sweep in the electrochemical cell and get reduced to diazonium radicals. These diazonium radicals attack the available graphene/graphite area and form organized covalent bonds with the surface. Using n-alkane networks of different sizes the periodicity of covalently grafted molecules can be controlled. The organic layer plays a critical role in this method by allowing lateral displacement of n-alkane terminal edge around which the alkyl radical attack occurs. Raman spectroscopy and STM confirmed covalent grafting of molecules on HOPG and graphene. STM images clearly showed the patterned grafting of 3,5-TBD molecules as shown in Figure 13(j) Diels-Alder (DA) cycloaddition reactions on graphene were thoroughly studied by Haddon and coworkers [77,79]. Due to the absence of an electronic bandgap, graphene can act as either diene or dienophile for DA reactions. Using STM, Raman spectroscopy and ARPES, the authors have demonstrated the covalent grafting of different molecules such as tetracyanoethylene (TCNE) and maleic anhydride (MA), where graphene acts as diene and with 2,3-dimethoxy-1,3-butadiene (DMBD), where graphene acts as a dienophile. Due to the formation of two sp 3 bonds between the molecules and graphene, this approach is highly promising towards selective functionalization of graphene with the potential for atomically precise patterning for CMOS technology. L. Simon and coworkers [13] showed the modification of the Fermi velocity of electrons upon increasing the grafting density of fluorinated maleimide molecules using ARPES. The authors have proposed a method for determining the location of grafted molecules by systematically analyzing the pattern of standing waves generated by quasiparticle scattering of electrons. Apart from diazonium grafting there are several other approaches towards covalent functionalization of graphene. One such approach was shown by Bueno et al. by functionalizing graphene by p-aminophenol (p-AP), which can covalently bind to graphene on the site of dangling bonds surrounding vacancies in the graphene lattice [91]. UHV-STM images and XPS confirmed the covalent attachment of p-AP molecules on graphene and the impact of grafting on the electronic band structure of graphene was evident by ARPES studies (Figure 14). Cycloaddition reactions of graphene with aromatic molecules can also result from interactions of graphene with the underlying substrate. For graphene on Ir(111), weak covalent interactions between C atoms and Ir atoms arise in the areas of the moiré pattern where C atoms are located on top of Ir atoms due to π-d hybridization. As a result, there can be unsaturated dangling bonds that are located on C atoms in hollow sites of the Ir lattice. Altenburg et al. showed that different MPc molecules at lower coverage adsorb distinctly by binding to the graphene layer with one to three Pc lobes [92]. It was reported that most of the MPc molecules are physisorbed and form dimers or clusters, but some of the MPc molecules show only two or less lobes and are immobile as observed by STM images as shown in Figure 15(a). STS spectra of the four lobes state compared with the three lobes state showed a shift of peaks by 200-300 meV for the three-lobed states, which was further downshifted for two-lobed or one lobed state (Figure 15(b)). DFT calculations confirmed a cycloaddition reaction taking place between MPc molecules and C atoms of graphene. It was concluded that the molecules react on the fcc sites of graphene. Similar studies carried out by Kumar et al. showed the covalent bond formation between graphene and a cyano group of F4-TCNQ molecules [93] by STM/STS and nc-AFM. The molecules showed varying adsorption geometry on different sites of a moiré unit cell of graphene on Ir(111), which resulted in a change of the molecular electronic structure. On fcc sites of the G/Ir(111) moiré, the molecules shows a two-legged structure, whereas on hcp and top sites the molecules were lying flat on the surface (Figure 15(d-f)). Further, the molecular configuration could be manipulated by the STM tip. STS measurements revealed that in two-legged geometry the molecules were charged, whereas in the other geometry the molecules remained uncharged. Two possible reasons were identified for the observed charging: a) a tradeoff between tilt-induced charging and reduced van der Waals interactions and b) partial rehybridization of graphene carbon atoms into sp 3 configuration with enhanced bonding interactions with the cyano groups of F4-TCNQ molecules. The impact of functionalization by aryl diazonium molecules on electronic transport properties of graphene has been explored in a number of different studies [73,94,95]. A. Jacobsen et al. showed that diazonium chemistry has an effect on the transport characteristics of graphene. Transport measurements of conductance (G) as a function of backgate voltage (V BG ) at room temperature showed that after functionalization by NBD, a p-type doping is induced in graphene. A shift in the Dirac point from +9V to +31V were observed with increasing functionalization. H. Zhang et al. performed transport measurements on mechanically exfoliated graphene flakes before and after functionalization with 4-NBD molecules [95]. I-V, as well as conductance measurements for pristine graphene devices, did not show much dependence on temperature characteristics. However, after functionalization the I-V characteristics of the devices showed a linear to non-linear transition upon decreasing temperature from 300 K to 4.2 K. This transport behavior of graphene was in agreement with the variable range hopping (VRH) model, suggesting that diazonium functionalized graphene behaves as granular metal with a bandgap of~100 mV at low temperatures.

Covalent functionalization of TMDs
Several studies have been published in the past few years on covalent functionalization of TMDs by different methods. Most of the studies were focused on covalent functionalization of TMD nanosheets obtained by solution processing via ion intercalation [23,26] or by ultrasonication [18,96]. Following the successful demonstration of TMD nanosheets functionalized with alkyl halides and diazonium chemistry in solution, this approach was also applied for the functionalization of TMD nanosheets supported on different substrates. We have focused mostly on the functionalization of TMDs by diazonium chemistry, reaction of alkyl halides and thiol chemistry.
The pioneering work on covalent functionalization of transition metals in cluster form was reported by Tuxen et al [25]. In this work, the authors have shown the impact of MoS 2 cluster size on the adsorption efficiency of dibenzothiophene (DBT) and methyl substituted DBT (4,6-DMDBT) molecules with high-resolution STM images. MoS 2 clusters of different sizes react differently with the DBT molecules, depending on the cluster edge preference for S vacancy formation.
In 2015, Chhowala and coworkers showed a general route towards covalent functionalization of TMDs by attaching organohalide and aryl diazonium molecules to the basal plane of 2D TMD nanosheets [29]. In this work, the TMD nanosheets were first converted into a metallic 1T/1T' polytype using n-butyllithium, followed by functionalization with 2-iodoacetamide and iodomethane. This approach also allows the functionalization of substrate-supported CVD MoS 2 nanosheets and functionalization with 4-bromobenzene diazonium molecules. The electron rich 1T/ 1T' polytype transfers charge to the organohalide (or diazonium) molecules resulting in around 30% of functionalization as confirmed by XPS studies. Structural characterization by atomically resolved transmission electron microscopy (TEM) measurements suggested that the octahedral coordination of the TMD nanosheets is not affected by the functionalization. Afterwards, studies by Knirsch et al. with a similar approach in liquid phase showed that pre-existing defects are not required for the functionalization of TMDs [26]. To produce MoS 2 nanoflakes without pre-existing defects, chemical exfoliation (CE) by intercalation of n-BuLi was performed in presence of an excess of MoS 2 powder, resulting in an almost 1:1 ratio for 2H and 1T polymorphs as shown by XPS results. Chemically exfoliated MoS 2 has been functionalized with 4-methoxyphenyldiazonium tetrafluoroborate. The covalent functionalization was confirmed by XPS, IR spectroscopy and TGA in combination with mass spectroscopy. After covalent functionalization, the 1T polytype of MoS 2 showed semiconducting behavior. This approach can be used for attaching various organic functionalities on TMDs using diazonium or alkyl halide interface [97,98]. However, while for most of the proposed applications of TMDs it is desirable to study the impact of covalent functionalization on substratesupported TMDs, functionalization of substrate-supported TMDs was not carried out.
Recent studies by Chu et al. have addressed the issue of covalent functionalization on the basal plane of substrate-supported TMDs [27]. The authors have shown by SPM, XPS, Raman and PL spectroscopy, TGA analysis and scanning electron microscopy (SEM) that 4-nitrobenzene diazonium molecules can be grafted on the surface of mechanically exfoliated as well as CVD grown MoS 2 layers. Exfoliated MoS 2 flakes on Si/SiO 2 were subjected to an aqueous solution of 4-NBD molecules. A change in surface topography of MoS 2 flakes was observed by AFM measurements on single layer regions of MoS 2 as shown in Figure 16(b). Theoretical calculations predicted a higher probability of grafting of molecules near the defect sites or close to the already grafted molecules as a result of a higher surface dipole moment in proximity to these sites which is evident by complete coverage of NBD molecules on Ar plasma treated MoS 2 after very short reaction times (5s). A high PL yield was observed for NBD functionalized monolayer MoS 2 when compared to fewlayer MoS 2 due to a direct bandgap in monolayer MoS 2 (Figure 16(c-h)). An increase in C-S peak intensity as a function of increasing functionalization duration on CVD MoS 2 was observed in the C1s XPS spectra, confirming the covalent nature of the bond between 4-NBD molecules and MoS 2 . Further, the authors have shown that diazonium chemistry can be used for biosensing applications. Therefore, 4-carboxylbenzene diazonium molecules were first grafted on MoS 2 samples and subsequently reacted with polyhistidine (His)-tagged green fluorescent protein (GFP) and the red fluorescent protein, mCherry. The surface was probed using confocal fluorescence microscopy, which showing bright field optical images along with fluorescence images in the green and red channels for GFP and mCherry emission, respectively.
Additionally, diazonium chemistry has been shown to be useful in improving the catalytic activity of MoS 2 in hydrogen evolution reactions (HER) [97]. 4-Cyanobenzyl was grafted on MoS 2 nanosheets followed by exposure to Co(dmgBF 2 ) 2 complex, resulting in the coordination of cobalt on 4-cyanobenzyl. Tafel slope for MoS 2 -Co(dmgBF 2 ) 2 was found to be much lower compared to CE-MoS 2 , which suggests that the kinetics of the water molecule dissociation step are efficiently facilitated on the surface of MoS 2 -Co(dmgBF 2 ) 2 . In a similar approach, Chen et al. used the molecular photosensitizer [RuII(bpy) 3 ] 2+ covalently grafted on CE-1T-MoS 2 for light-driven HER [98].
Organic molecules containing a thiol group have also been used for the covalent functionalization of TMDs, particularly MoS 2 . Thiol groups react and form a bond with Mo atoms at the single vacancies present in MoS 2 layers. In almost all of the studies concerning modification of sulfur vacancy (SV) sites using thiol chemistry resulted in an improvement in the charge carrier mobility of the functionalized TMDs substrate. The clear proof of the formation of a covalent bond between thiol molecules and   [24]. Few-layer MoS 2 was functionalized with two thiol derivatives (3-mercaptopropyl)-trimethoxysilane (MPS) and dodecanethiol, which were covalently bound to the SV sites of MoS 2 . The molecules were observed as bright features in STM images. These bright features could be removed by STM measurements carried out at high tunneling current, which authors termed as brushing mode. Additionally, the sulfur vacancy defects, which are usually observed on pristine MoS 2 , also disappeared. Thus, it was concluded that the thiol molecules form bonds with the Mo atoms at the defect sites by transferring electrons from the S lone pair to the Mo 4d orbitals, i.e. a back electron transfer occurs from the Mo 4d orbitals to the antibonding S-C orbital, which is comparatively weak and results in S-C bond cleavage during STM measurements. These studies suggest the possibility for selective chemisorption at pre-determined vacancy sites.
A similar approach for the passivation of sulfur vacancies in monolayer MoS 2 by (3-mercaptopropyl)trimethoxysilane (MPS) molecules was reported by Yu et al [99]. Due to the presence of trimethoxy silane groups, which form SAMs on SiO 2 as discussed in the earlier sections, functionalization of MoS 2 from the bottom was also possible. Transport characteristics from back-gated FET using functionalized MoS 2 as channel samples showed an increase in room temperature field effect mobility up tõ 80 cm 2 V −1 S −1 . Similar studies carried out by Cho et al. on chemisorbed alkanethiol molecules on multilayered MoS 2 showed a shift in Raman spectra, PL spectra and a decrease in channel current in FETs [100]. The observed decrease in the channel current was attributed to the passivation of SV sites by alkanethiol molecules. It was postulated that passivation of SV sites resulted in a decrease of electron hopping, which takes place via SV sites in pristine MoS 2 channels. The channel current was measured for different samples with varying MoS 2 thickness in order to establish a relation between the impact of alkanethiol molecules and surface-volume ratio of MoS 2 . These studies clearly showed the possibility of using thiol molecules to tune the electronic properties of TMDs. In this context, Sim et al. showed controlled n-type or p-type doping of MoS 2 by functionalization with thiol molecules containing donor or acceptor groups [22]. Using thermal annealing, additional SV sites were created in MoS 2 which were then functionalized by mercaptoethylamine (MEA) and 1H,1H,2H,2Hperfluorodecanethiol (FDT) as n-type and p-type dopant, respectively. Covalent functionalization of MoS 2 was confirmed using XPS analysis of C1s and N1s core level spectra. PL measurements showed a redshift for MEA-MoS 2 and blue shift for FDT-MoS 2 , suggesting n-type and p-type doping, respectively (Figure 17(e,f)). XPS results were supported by AFM measurements, which show a clear modification of the MoS 2 flakes after functionalization (Figure 17(a-d)). The doping effect was also observed by transfer characteristics measured by FET devices. However, the doping effect became much more evident after thermal annealing of the samples at 150°C for 10 min, which possibly removed water and oxygen form the surface.

Covalent functionalization of hBN and black phosphorus
Compared to graphene chemistry covalent functionalization of hBN is not extensively studied, due to its perceived inert nature. However, there are still some studies which show covalent chemical functionalization of hBN nanosheets. Lin et al. showed the covalent attachment of octadecylamine (ODA) on the defect sites of hBN nanoflakes [101]. To enhance the grafting density of ODA molecules, additional defects were created in the hBN flakes. TGA measurements confirmed the impact of these additional defect sites on the grafting efficiency on hBN. A signature for amine-boron interactions was observed by NMR spectroscopy. Another approach of covalent functionalization of hBN is solution phase oxygen radical functionalization of boron atoms as demonstrated by Sainsbury et al [102]. In this two-step approach, first the hBN nanosheets are reacted with tert- butoxy radicals which results in the formation of boronate ester bonds C-O-B with boron atoms in the lattice. This forms a chemical handle at the surface which is susceptible to hydrolysis and thus can be used to graft other chemical species.
A more straightforward approach for covalent functionalization of hBN nanosheets was demonstrated by Coleman and coworkers [103]. Covalent functionalization by nitrene groups was carried out by first attaching an organo-azide precursor, 4-methoxybenzyloxycarbonyl azide, which upon thermolysis creates nitrene groups which attach to boron atoms in the hBN lattice.
Kumar et al. have recently studied molecular doping of hBN nanoflakes by TCNQ molecules [104]. UV-vis and PL spectra showed the emergence of new peaks for TCNQ functionalized hBN nanoflakes. On further investigation, it was found that an oxidative decay product of TCNQ, namely α, α-dicyano-p-toluoylcyanide (DCTC), was responsible for new UV-vis peaks and PL peak in the red region. DFT calculations revealed that hybridization between the N atomic states of hBN flakes with boron vacancy (hBN-V B ) and DCTC molecular orbitals resulted in new states at 1.8eV and 1.96eV. Upon excitation, an electron transfer from these states to the LUMO of DCTC molecules leads to emission peaks at~616 nm and~672 nm. These studies show that a charge transfer process between the mid-gap states, resulting from hybridization between molecules, and defect sites of the substrates can also lead to a modification in the optoelectronic properties of 2D materials.
Aryl diazonium based covalent functionalization of black phosphorus was recently reported by Hersam and coworkers [105]. BP flakes prepared by mechanical exfoliation on SiO 2 /Si wafers were immersed in a solution of 4-nitrobenzene-diazonium (4-NBD) or 4-methoxylbenzene-diazonium (4-MBD). The immersion time had a profound impact on the degree of functionalization. Topographical images by AFM showed an increase in the thickness of BP flakes form 7 nm to 8.3 nm after functionalization. Peaks corresponding to P-C bonds were observed by XPS spectroscopy, confirming the covalent functionalization of diazonium molecules. BP flakes functionalized by 4-NBD molecules showed very little or no change in morphology even after 10 days in ambient conditions, confirming the passivating effect of grafted diazonium species as shown in Figure 18(a,b). Transfer characteristics confirmed the electron transfer from BP flakes to NBD molecules resulting in p-type doping of BP flakes as shown in Figure 18(c,f).

Conclusions and outlook
In the first part, highlighting non-covalent functionalization strategies, we focused on supramolecular chemistry-based approaches towards engineering the electronic band structure of 2D materials. Charge transfer doping by self-assembled molecules provides an alternative to substitutional doping towards controlled doping of 2D materials. Organic molecules with electron donating or withdrawing functional groups can be used for inducing n-type or p-type doping in 2D materials. This approach is highly versatile as there is a vast library of functional groups and organic molecules that can be synthesized depending on the nature of the 2D material and the extent of doping required, and can be utilized for the bottom-up synthesis of van der Waals heterostructures with 2D materials. SAMNs can induce a periodic potential modulation which can result in profound modification of electronic band structure. In the case of TMD functionalization, this may also lead to superior transport characteristics and strong photoluminescence yield. Additionally, self-assembled molecules can be used as passivating layers for highly unstable 2D materials like black phosphorus and some TMDs.
Covalent functionalization is an interesting approach which is more robust and leads to a strong modification of the electronic properties of 2D materials. We have discussed the most prominent chemistries for functionalization of graphene and other 2D materials. The strong concerns over the randomness of the covalent chemistry approach have very recently been addressed by combining it with SAMN templating, resulting in patterned grafting on 2D materials. This approach allows a better control over the density and the extent of grafting. Covalent chemistry also allows attaching functional groups to 2D materials, which are sensitive towards different organic and inorganic species, therefore realizing the potential of 2D materials in sensors.
Finally, we have stressed the fact that SPM-based methods are critical in studying the impact of molecular functionalization on the properties of 2D materials. These techniques result in acquiring atomically or sub-molecular -resolved images of the molecule-2D material systems. Beyond providing structural information, scanning probe-based spectroscopy approaches, such as STS and KPFM, are also highly useful in studying the modulation in the electronic band structure of 2D materials interacting with molecules.
Non-covalent and covalent functionalization approaches of 2D materials such as graphene, TMDs, BP and hBN are highly promising methods for the controlled modification of a number of properties such as dispersion in solutions, tuning of the electronic band structure and transport properties. Although the functionalization of graphene and TMDs is well established, the development of functionalization strategies for other 2D materials such as BP and hBN is still ongoing. A number of fundamental challenges such as the development of suitable characterization methods, improving the interactions between molecules and target materials and the reactivity of 2D materials towards various functionalizing compounds remain to be addressed to obtain well-controlled chemical activation strategies for 2D materials. Chemical modification of 2D materials using molecules and assemblies of molecules is not only a promising alternative to other approaches aiming at controlling the (opto)electronic properties of 2D materials, the huge range of compounds and chemistries that can be used for functionalization further offers almost endless opportunities to modify and engineer the properties of 2D materials. New functionalities added in 2D materials using chemical approaches can open up possibilities towards the development of multifunctional materials. As discussed in this review, 2D materials can be functionalized by photo-switchable, electrochemically switchable and magneto-responsive molecules, which can be activated by external stimuli. A controlled activation of molecules at the interface of 2D materials will result in tuning the electronic, chemical and magnetic properties of these materials with much higher control. A systematic control over the degree and site of functionalization would result in better control over the optoelectronic properties of 2D materials. Finally, encapsulation of molecule-2D material hybrid structures must be studied in order to improve the long-term stability of devices based on these hybrids.