Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 14767.
Published online 2017 March 21. doi:  10.1038/ncomms14767
PMCID: PMC5364416

Periodic potentials in hybrid van der Waals heterostructures formed by supramolecular lattices on graphene


The rise of 2D materials made it possible to form heterostructures held together by weak interplanar van der Waals interactions. Within such van der Waals heterostructures, the occurrence of 2D periodic potentials significantly modifies the electronic structure of single sheets within the stack, therefore modulating the material properties. However, these periodic potentials are determined by the mechanical alignment of adjacent 2D materials, which is cumbersome and time-consuming. Here we show that programmable 1D periodic potentials extending over areas exceeding 104 nm2 and stable at ambient conditions arise when graphene is covered by a self-assembled supramolecular lattice. The amplitude and sign of the potential can be modified without altering its periodicity by employing photoreactive molecules or their reaction products. In this regard, the supramolecular lattice/graphene bilayer represents the hybrid analogue of fully inorganic van der Waals heterostructures, highlighting the rich prospects that molecular design offers to create ad hoc materials.

Van der Waals (vdW) heterostructures, in which different two-dimensional (2D) materials are superimposed in a stacked configuration, represent a versatile experimental platform to study fundamental properties for device physics and materials science1,2,3. Based on these mechanically assembled stacks, atomically thin novel devices have been demonstrated, such as (tunnelling) transistors4, p-n and tunnelling diodes5,6, photovoltaic elements7 and light emitting diodes (LED)8.

From a more fundamental point of view, the mutual interaction between vertically-stacked 2D materials generates electronic properties that are different from those of the isolated materials, as experimentally demonstrated for graphene on boron nitride9,10,11,12,13. In the latter case, a Moiré pattern develops due to the lattice mismatch between the two materials14, resulting in a periodic potential (superlattice) with hexagonal geometry capable of modifying profoundly the band structure of graphene. Similarly, other geometries of periodic potentials are predicted to affect the electronic properties of graphene in different ways. For example, a one-dimensional (1D) Kronig–Penney periodic potential with nanoscale periodicity is predicted to create an anisotropic propagation of charge carriers along the different directions of the potential15. Hitherto, 1D graphene superlattices that could be pre-programmed with atomic precision have not been demonstrated.

In this context, the use of molecules offers two attractive features. First, organic molecules can prompt doping effects in 2D materials16,17,18,19, causing local modifications of their surface potential20,21,22. Second, these molecules form ordered 2D crystalline structures when physisorbed on graphene and other 2D materials23,24,25,26,27,28,29,30,31,32,33,34,35. Given the immense opportunities offered by hybrid organic–inorganic vdW heterostructures in modifying the fundamental electronic properties of the pristine materials, the field is still widely unexplored3. There are only a few reports connecting doping effects with the specific position of functional groups26,36,37,38,39, and combining molecules and 2D materials for the technological need of forming p–n junctions40,41.

In this work, we demonstrate that a tunable periodic potential with 1D geometry on graphene can be realized at ambient conditions in a hybrid vdW heterostructure composed of a 2D supramolecular lattice (SL) self-assembled on graphene, with single domains extending over areas exceeding 250 × 250 nm2. The amplitude of the periodic potential is mediated by the specific interaction between graphene and molecular dipoles, while the nanoscale periodicity is determined by the geometry of the self-assembly. In particular, by using a photoreactive organic molecule, we show how a subtle photo-induced change in the chemical structure of the starting molecule leads to a different amplitude of the potential, while leaving the periodicity unaltered.


Molecular approach to the formation of periodic potentials

The analogy between hybrid and fully inorganic vdW heterostructures goes beyond their layered nature; in analogy to stacks of inorganic 2D materials, vdW forces drive the interaction between graphene and the self-assembled molecular layer, and similar inter-layer interactions can be envisaged. However, the paradigm leading to the formation of such heterostructure is substantially different. While in inorganic vdW heterostructures the different monolayers are mechanically superimposed with an empirical (optical) alignment42, in the molecular case a careful choice of the molecular unit allows the spontaneous creation of self-assembled and ordered layers with predictable geometry and atomic precision. Hence, single molecules can be regarded as molecular building blocks (MBBs) that determine both geometry and functionality of the SL.

Our approach is schematically illustrated in Fig. 1a. We design an MBB composed of a head, bearing a light-reactive moiety, and a tail, consisting of a long aliphatic chain. This linear tail drives the molecular self-assembly43,44 and acts as a spacer between adjacent functional heads, forming a lamellar assembly which determines the 1D periodicity of the potential (see Fig. 1a). Electric fields generated by molecular dipoles within the functional heads are responsible for the modulation of the surface potential of graphene, introducing an electric-field effect analogous to that of a constant external gate, and ultimately determining the amplitude and sign of the periodic potentials. Moreover, the functional headgroup of the photoreactive MBB can be modified before deposition on graphene by simple irradiation in different solvents, yielding new MBBs. Specifically, we used 3-trifluoromethyl-3-(3-octadecyloxyphenyl)diazirine (MBB-1), shown in Fig. 1b, as photoreactive MBB. The diazirine headgroup consists of a carbon bound to two nitrogen atoms, and was selected because it can undergo photolysis very efficiently, forming a reactive carbene and nitrogen gas under irradiation with ultraviolet light (Fig. 1b)45. If irradiated in chloroform, the in situ generated carbene binds to the solvent molecules in its proximity, generating modified MBBs in a mixture which we call MBB-2. As detailed in the Supplementary Figs 1–3, we found that the main reaction product in the mixture has the structure shown in Fig. 1b, in which a Cl atom substitutes the diazirine moiety. This evidence is further corroborated by a joint optical and electrical characterization (see Supplementary Fig. 4) performed on ad hoc synthesized molecules (called P-1 and P-2). Absorption spectra allow to follow the change in the molecular head of MBB-1, as shown in Fig. 1c. After 30 min of irradiation at λ=365 nm (areal power density 1.7 mW cm−2), the characteristic absorption at λ=365 nm of MBB-1 disappears, revealing that MBB-1 was efficiently photolyzed.

Figure 1
Cartoon of our approach.

Nanoscale characterization of the supramolecular lattices

The periodicity and geometry of the SL-induced potential are determined by the nanoscale molecular arrangement which is investigated through scanning tunneling microscope (STM) imaging. We performed STM imaging in ambient conditions on dry SLs, formed by simply spin-coating a solution of either MBB-1 or MBB-2 (as displayed in Fig. 2a–c or 2d–f, respectively). In view of its atomic flatness over hundreds-of-micrometer-sized terraces, the SL can be mapped with increased spatial resolution on highly ordered pyrolitic graphite (HOPG) substrates (Fig. 2a,b,d,e), yet its motif is identical to that obtained when using graphene grown by chemical vapour deposition (CVD) supported on SiO2 as substrate (Fig. 2c,f).

Figure 2
Scanning tunnelling microscopy images of supramolecular lattices.

The self-assembly is driven by the interplay of molecule–substrate and molecule–molecule interactions; in particular vdW interactions occur between the different alkyl chains, which adsorb flat on graphite/graphene, generating well-defined, crystalline lamellar architectures.

In the case of MBB-1, the functional heads lie close to each other, with the alkyl chains forming an interdigitated structure, exhibiting a unit cell where a=(3.8±0.1) nm, b=(0.9±0.1)  nm and α=(84±2)° therefore leading to an area A=(3.4±0.1) nm2, with each unit cell containing two molecules (Fig. 2a). Similar SL structures are also observed for MBB-2 (Fig. 2d). Again, the alkyl chains lie flat on the surface by aligning in interdigitated crystalline structures with a similar unit cell where a=(3.8±0.2) nm, b=(0.9±0.1) nm, α=(84±2)° giving rise to an area A=(3.5±0.2) nm2. The nearly identical unit cells are justified if one considers that on ultraviolet irradiation only a few atoms in the molecular head are replaced. Instead, the long alkyl chains are unaltered in both cases. Hence, MBB-1 and MBB-2 form SLs that possess the same lattice parameters but differ in the head groups, resulting in different electronic interactions with graphene.

For both MBB-1 and MBB-2, crystalline domains with different orientations of the lamellae may form on the graphene surface. This evidence reflects the three-fold symmetry of the substrate and may give rise to misorientation in the direction of the 1D potential. Therefore the size of single-crystalline domains is a key parameter for the relevance of this study on the device physics. To gain detailed insight into this issue, we recorded survey STM images mapping the surface over a scale of a few hundred nanometres. Figure 2b,e show large-area STM images of MBB-1 and MBB-2, respectively. In both cases, an ordered array of lamellae with only one orientation extends over the whole image. Domains with different orientations could be observed in other large-area STM images. After careful analysis of 15 large-area STM images for each MBB, we can conclude that the typical single-domain size is ~300 nm × 300 nm, as shown in Supplementary Fig. 5 and in Supplementary Table 1. Compared with the length scale relevant for the charge transport in graphene, these single domains are typically two orders of magnitude larger than the typical size of the potential puddles that are responsible for limiting the mobility at the graphene/SiO2 interfaces46. Moreover, such area is still suitable for nano-fabrication and optical/spectral inspection of the physical properties generated by the presence of the SLs. The assembly of MBB-1 and MBB-2 was investigated also on CVD graphene on SiO2 to consider a situation closer to that of the actual devices, as shown in the Fig. 2c,f. In this case, the image resolution is lower, due to the intrinsic roughness of the underlying SiO2 substrate. However, even on CVD graphene, it is possible to clearly observe the existence of ordered lamellar structures whose size is in good agreement with that monitored on HOPG. This finding confirms that the MBBs assemble in the same way on both HOPG and graphene/SiO2. Large-size STM images show that the crystalline order is maintained over wide areas even on CVD graphene on SiO2, as displayed in the Supplementary Fig. 5 for MBB-2. We stress that the STM images in Fig. 2 were measured in air at room temperature and that crystalline domain size and orientation were found to be unchanged even after several hours of continuous acquisition. Moreover, STM images revealed that the unit cell of the molecular assemblies recorded within a few minutes after the SL formation would be identical (within experimental error) to those recorded a few days later. Indeed, the formation of a tight molecular packing along with the strong interaction with graphene promoted by the long alkyl chains would hinder molecular diffusion and thus stabilize the assembly.

Electrical characterization of devices

To determine whether the molecules are effectively introducing a potential, we study the doping caused by the different SLs on graphene-based field-effect devices (see the ‘Methods' section for details on the device fabrication). The devices were measured before and after the formation of the MBB-1 and MBB-2 SLs, to evaluate their effect on the electrical characteristics of graphene. The measurements were reproduced on three different devices for each MBB, and the results were in qualitative and quantitative agreement (see the ‘Methods' section for details). Representative dependence of the drain current IDS on the gate voltage VGS is displayed in Fig. 3. In the case of MBB-1 (Fig. 3a), the charge neutrality point was not shifted significantly by the presence of the SL, indicating minor doping effects. Interestingly, the hole mobility increased after the formation of the SL (from 2,200 to 2,650 cm2 V−1 s−1 for the device in Fig. 3a). The increase in mobility was directly related to the presence of the SL, since the initial electrical characteristics of the clean graphene were recovered after the molecules were washed away by rinsing the sample with CHCl3, and the mobility decreased to its initial value (see the Fig. 3a).

Figure 3
Effect of supramolecular lattices onto the electrical characteristics of graphene.

In the case of MBB-2, after the formation of the SL the charge neutrality point was shifted towards positive voltages (ΔV>25 V), as shown in Fig. 3b. This effect corresponds to hole accumulation in the graphene channel, with an induced charge density of Δp>5 × 1012 cm−2, according to the widely used parallel-plane capacitor model36. The electrical characterization provides unambiguous evidence that the effect of the MBB-2 SL–averaged over several randomly oriented domains–is analogous to that of a (fixed) top gate. Moreover, the electrical characterization demonstrates that the doping effect is solely determined by the head group, which is the only part of the MBB modified by the ultraviolet irradiation.

Generally, our devices possess a channel length of a few micrometres that is wider than the typical size of a single domain. As a consequence, our measurements probe not a single but rather a few domains with different random orientations, which impedes measuring the anisotropy in electrical conductance possibly induced by the SL. However, the information extracted by the electrical characterization provides very useful insight on effects, averaged over several randomly oriented domains, associated with the formation of the hybrid vdW heterostructure. While the STM nanoscale characterization provides information about the geometry and periodicity of the potential, the information extracted through device characterization is related to the amplitude of the periodic potential. The demonstration of anisotropy in the conductance of graphene would require a miniaturization of devices that is beyond the scope of this work. To better quantify the effects on the nanoscale electrical properties in presence of the SL, we have performed conductive atomic force microscopy (C-AFM) measurements on a contacted graphene flake, acting as the device active layer, before and after formation of a MBB-1-SL. After formation of the SL, some inhomogeneities in the conductance are introduced with a length scale of a few hundred nanometres, which is the typical size of the single domains imaged by STM. Although the lamellas cannot be visualized by C-AFM, this finding supports the presence of anisotropic conductance within graphene covered by a single-domain SL, as detailed in Supplementary Fig. 6.

Nanoscale origin of the doping effect

To gain more insight into the origin of the doping effects, the interactions of the molecular assembly with graphene were elucidated through molecular mechanics/molecular dynamics (MM/MD) simulations. For the case of MBB-2, we assumed that all the molecules at the surface possess the structure shown in Fig. 1b, following the findings described in Supplementary Figs 1–4. For both the assembly of MBB-1 and MBB-2, (MM/MD) simulations show that the alkyl chains are interdigitated and the carbon backbones are lying parallel to the graphene surface. The lattice parameters of the calculated unit cell in both cases, that is MBB-1 and MBB-2, are the same, in agreement with the experimental data (see Supplementary Fig. 7). From a closer look at the functional head groups of the calculated assemblies, a major difference appears between the positioning of the molecular head in MBB-1 versus MBB-2 (see Fig. 3c,d). In the former, the diazirine moiety is lying flat on graphene, so that the CF3 group is pointing parallel to graphene, whereas the axis of the N=N double bond is normal to the graphene plane (Fig. 3c). Conversely, the CF3 unit of MBB-2 is pointing in the z-direction normal to the graphene surface (Fig. 3d). Interestingly, such a difference in the orientation of the functional head groups has very limited impact on the overall organization of the molecular adlayer assembly, but it drastically affects the intrinsic dipole moment of the SLs in the z-direction. Indeed, in the MBB-1 case the molecules are almost parallel to the graphene surface, so that a small net dipole moment per molecule is expected perpendicular to the surface (Fig. 3c). Instead, in the MBB-2 case, a larger vertical dipole arises from the CF3 groups, which are found to be almost perpendicular to the graphene layer (Fig. 3d). Based on this observation, one can qualitatively understand the p-type doping measured in graphene with the MBB-2 adlayer. The electric-field generated by the vertical dipoles acts as a top gate and increases the work function of graphene, effectively inducing p-type doping. At a more quantitative level, density functional theory (DFT) calculations made it possible to estimate the doping induced by the two SLs on the basis of their calculated arrangement, as detailed in Supplementary Methods and Supplementary Figs 8–11. More specifically, molecular doping of graphene is usually sourced by two effects17,47, namely charge transfer from the molecules, and the presence of molecular dipoles exerting a local gating. DFT calculations allow to disentangle these two contributions for each SL and demonstrate that the charge transfer contribution is almost identical in both cases. Instead, in the case of MBB-1 the calculated vertical molecular dipole is rather low (μz=−0.17 D mol−1), while in the case of MBB-2 it is significantly higher (μz=−0.82 D mol−1). The doping effect in the MBB-2 case, therefore, results from the higher vertical dipole moment of the SLs, as expected from the geometry of the assembly. By accounting for molecular motion at room temperature from MD simulations, we predict a negligible and a significant work function shift for MBB-1 and MBB-2, respectively, as experimentally verified by photoemission spectroscopy in air (see Supplementary Fig. 12). The work function shift is accompanied by an induced p-type doping for MBB-2 (on the order of Δp=5 × 1012 cm−2) and a minor p-type doping for MBB-1 (Δp<5 × 1011 cm−2), in excellent agreement with the device experiments. Moreover, our theoretical analysis confirms that the origin of the measured effects can be entirely ascribed to the geometry of the head groups, while the assembly of the alkyl chains, similar for MBB-1 and MBB-2, is not a source of doping.

Spatially resolved electrostatic potential

Proceeding with the analysis, the spatially resolved electrostatic potential at the graphene/molecule interface was calculated by means of classical microelectrostatic calculations (see Supplementary Methods for details).

In Fig. 4a,b we show the differential potential V (z=8 Å)–V (z=0 Å) for both MBB-1 and MBB-2, where V (z=8 Å) is the potential calculated on the parallel plane just above the SLs, and V (z=0 Å) is the potential at the graphene surface. Such differential potential probes the dipolar field associated with the out-of-plane component of the molecular dipoles and has an effect similar to that of a fixed externally applied gate voltage. In both assemblies, such induced gating potential is characterized by a 1D-modulation, with higher negative values localized at the molecular heads separated by inert alkyl chains. The amplitude of the induced gating effect increases by a factor of 4 from MBB-1 to MBB-2, as expected on the basis of the different vertical dipoles. Note that image charge and depolarization effects associated with intermolecular interactions have been estimated using DFT electronic structure calculations implementing periodic boundary conditions. These were found to reduce the electric dipole per molecule by ~30% for MMB-1 and 20% for MBB-2. Thus, image effects are found to be fairly limited in our case and the graphene layer acts primarily as a template breaking up symmetry in picking selected dipolar conformations for the MBBs at the surface, which then translates into the observed electrostatic potential shown in Fig. 4.

Figure 4
Periodic potentials introduced by the supramolecular lattices.

This different amplitude induced by the MBB-1 and MBB-2 SL demonstrates the ability of introducing significant variations in the SL-induced periodic grating with a subtle light-induced change in the molecular structure. In both cases, the in-plane distribution of the potential can be considered as a nanoscopic experimental realization of a Kronig–Penney potential. We highlight that not only its geometry, but also its amplitude is within the same order of magnitude of that considered in the initial prediction of anisotropic behaviour of charge carriers in graphene superlattices15.

Modification of the periodic potential

Finally, we show that periodic potentials with the same geometry but different intensity can be achieved by preparing SLs of MBB-1 after ultraviolet irradiation in different solvents. To prove further this concept, we present here a detailed nanoscale and electrical characterization of MBB-1 after ultraviolet irradiation in diethylamine, which results in a mixture hereafter referred to as MBB-3. The assembly of MBB-3 is shown in Fig. 5a,b on HOPG and CVD graphene, respectively. Similarly to MBB-1 and MBB-2, MBB-3 self-assembles forming ordered lamellae in which the functional heads lie close to each other, with the alkyl chains forming an interdigitated structure. The unit cell parameters were a=(3.8±0.2) nm, b=(0.9±0.1) nm and α=(84±2)° which lead to an area A=(3.5±0.2) nm2, analogous to that of MBB-1 and MBB-2. Again, the atoms incorporated through photolysis do not perturb the assembly at the nanoscale, which is determined and imposed by the alkyl chains. Instead, the intensity of the induced potential is determined by the interaction with the head groups. The overall effect can be measured in a three-terminal device, by covering a graphene device with a MBB-3 SL. In this case (Fig. 5c), the charge neutrality point shifts towards negative values, corresponding to electron accumulation (n-doping) in the graphene channel. A shift in the charge neutrality point ΔV=−11 V corresponds to an induced electron density Δn=2.6 × 1012 cm−2. Even in this case, the electrical properties of the clean graphene are recovered when the molecules are washed away from the graphene surface with a simple rinsing step in chloroform.

Figure 5
Modification of the periodic potential.

Interestingly, the effect on the IDSVGS traces of graphene is opposite to that of MBB-2. Since the induced periodic potential is primarily determined by the orientation of local molecular dipoles, one is thus led to conclude that the sign of such a potential for MBB-3 is opposite to that of MBB-2.


The concept that interplanar interactions are capable of modifying fundamental electronic properties of pristine materials lies at the heart of the research field of vdW heterostructures. In this context, the possibility of introducing tunable 1D potentials through organic SLs is particularly appealing. For instance, one could combine graphene covered by a SL with other inorganic 2D materials to fabricate more complex multi-layered hybrid vdW heterostructures. In such systems, an ordered lamellar assembly modifying the momentum dispersion of graphene would induce anisotropy not only on the in-plane (graphene) but also on the inter-layer charge transport, for which momentum conservation plays an important role48. Moreover, a fully anisotropic heterostructure might be obtained by combining our SL-covered graphene with other 2D materials intrinsically possessing 1D anisotropy (such as ReS2 (ref. 49) and black phosphorus50) to demonstrate novel device architectures actively exploiting charge transport along a preferential direction. Finally, we point out that bottom–up supramolecular approaches for the realization of periodic potentials is not limited to graphene, since analogous molecules form ordered 1D assemblies also on the surface of other 2D materials51, widening the horizons of the present study. For instance, a SL similar to those described in this study assembled on semiconducting transition metal dichalcogenides (such as WSe2 and MoS2) might introduce a 1D periodic potential capable of locally varying the position of the Fermi level, resulting in alternating hole- or electron- rich regions. In this way, one would obtain a series of consecutive nanoscale p–n junctions with distinctive charge transport properties and optical response.

In conclusion, our study shows that organic SLs represent an ideal system for generating periodic potentials whose periodicity, amplitude and sign can be pre-programmed by careful molecular design. In particular, here we have shown that a periodic potential with 1D geometry can be generated at the graphene surface and manipulated by making use of molecular photo-reactivity. In perspective, a 1D potential will offer the opportunity to introduce anisotropy in otherwise isotropic materials, hence paving the way for design and implementation of novel functional vdW heterostructures, for example, for 1D charge transport. The control over periodic potentials can hardly be achieved in vdW heterostructures based solely on inorganic 2D materials, while their combination with SLs allows expanding substantially the 2D material library. Chemistry can offer an almost unlimited choice of SL that can form hybrid vdW heterostructures with controllable structural and electronic properties, exerting periodic potentials with tunable amplitude and periodicity. In perspective, a great deal of novel electrical, magnetic, piezoelectric and optical functionalities arising from such hybrid vdW heterostructures are expected by taking full advantage of the infinite degrees of freedom offered by the design of the MBBs.


Sample preparation

Full details regarding the characterization and synthesis of the MBBs are given in Supplementary Figs 13–19 and in Supplementary Note 1.

Ultraviolet–visible spectroscopy

Absorption spectra were recorded at room temperature (~25 °C) with a JASCO V-670 spectrophotometer and all solutions were examined in quartz cells with 2-mm pathlength (HELLMA) with a concentration of 2.2 × 10–3 M. Photochemical reaction of MBB-1 was performed in air-equilibrated CHCl3 (Uvasol Merck-Millipore) or in degassed diethyalmine (Sigma-Aldrich) by using an Ultraviolet lamp (UV-6 L/M Herolab) with λirr=365 (±5) nm and a power density=1.7 mW cm−2.

Device fabrication and characterization

Back gated devices were fabricated on Si/SiO2 (90 nm) substrates with a Microtech laser writer, exposing with a λ=405 nm laser a standard photoresist (AZ1505, Microchemicals). Gold (without adhesion layer) was thermally evaporated onto the patterned photoresist and lift-off was carried out in warm acetone (40 °C). After fabrication, the devices were immersed in warm NMP (40 °C) overnight, rinsed with chloroform, acetone and isopropanol. The devices were kept in a nitrogen-filled glove box in which they could be measured in a probe station connected to a Keithley 2636. By employing this procedure, the graphene device showed almost-ideal, symmetric and stable current–voltage IDSVGS characteristics, with the charge neutrality point close to a gate voltage VGS=0 V, and hole mobility reproducibly above 1,800 cm2 V−1 s−1. The carrier mobility was determined by the conventional parallel-plane capacitor model:

An external file that holds a picture, illustration, etc.
Object name is ncomms14767-m1.jpg

where L and W are the graphene channel length and width, t is the thickness of gate oxide, epsilonr is the relative dielectric permittivity of SiO2. The mobility was extracted at the carrier concentration value of |n|=1 × 1012 cm−2, by measuring dIDS/dVGS at 4.5 V away from the charge neutrality point. In all the electrical measurements, the potential applied between the graphene electrodes was V=10 mV. The reproducibility of the effect of each SL on the electrical characteristics of devices was tested by repeating the measurement on three different devices for each MBB. For every test, a fresh solution of MBB was employed. The doping effect is reproducible: MBB-1 was found to induce small p-type doping (ΔV<5 V in the three tested devices); MBB-2 significantly higher p-type doping (ΔV>25 V in the three tested devices); and MBB-3 n-type doping (ΔV in the range between −7 and −12 V in the three tested devices). Moreover, the experiments were repeated several times on the same device, by washing the MBBs away and re-forming the SL. In the cases of MBB-1 and MBB-3, the increase in mobility after the formation of the SL was also found in the three tested devices.

Supramolecular lattice formation

To form the SLs, molecules were spin-coated onto either HOPG or CVD graphene on SiO2 and onto devices from chloroform solutions (1 mg ml−1). For the sake of consistency, the solutions were spin-coated onto the graphene devices by employing the same parameters used for the STM imaging. In the case of irradiation in diethylamine, the molecules were dried after irradiation and re-solubilized in chloroform before spin-coating. In the case of devices, the solution were spin-coated in situ in the nitrogen-filled glovebox. The SLs could be washed away by rinsing the substrates or the devices with CHCl3.

Scanning tunneling microscopy

STM measurements were carried out by using a Veeco scanning tunneling microscope (multimode Nanoscope III, Veeco) operating with an A piezoelectric scanner which allowed the mapping of a maximum area of 1 μm × 1 μm. As substrates, we used highly oriented pyrolytic graphite and commercial CVD graphene supported on Si/SiO2 (300 nm) purchased from Graphenea. The graphene sample has been used as received, without any additional cleaning step, and has been stored in air for 4 months before the STM experiments. The substrates were glued onto a magnetic disk and an electric contact was made with (conductive) silver paint (Aldrich Chemicals).

The STM tips were mechanically cut from a Pt/Ir wire (90/10, diameter 0.25 mm). The images were obtained in air at room temperature. The raw STM data were processed through the application of background flattening, and in the case of the HOPG substrates in Figs 2a,d and and5a5a the drift of the piezo was corrected using the underlying graphite lattice as a reference. The lattice of the underlying substrate was visualized by lowering the bias voltage Vt to 10 mV and setting the average tunneling current It=60 pA. Tip height and current were measured for all STM images.

MD/MM and DFT calculations

Full details regarding the MD/MM and DFT calculations are given in Supplementary Methods and in Supplementary Table 2.

Data availability

The data that support the findings of this study are available from the corresponding authors on request.

Additional information

How to cite this article: Gobbi, M. et al. Periodic potentials in hybrid van der Waals heterostructures formed by supramolecular lattices on graphene. Nat. Commun. 8, 14767 doi: 10.1038/ncomms14767 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:

Supplementary Figures, Supplementary Tables, Supplementary Methods, Supplementary Note and Supplementary References.


We thank Wen Zhang and Joachim Räder for their help in the ASAP MS analysis. We acknowledge funding from the European Commission through the FET project UPGRADE (GA-309056), the Graphene Flagship (GA-696656), the Marie Sklodowska-Curie projects IEF GALACTIC (PIEF-GA-2014-628563) and ITN SYNCHRONICS (GA 643238), the Agence Nationale de la Recherche through the Labex projects CSC (ANR-10-LABX-0026 CSC) and Nanostructures in Interaction with their Environment (ANR-11-LABX-0058 NIE) within the Investissement d'Avenir programme (ANR-10-120 IDEX-0002-02), and the International Center for Frontier Research in Chemistry (icFRC). The work in Mons is supported by BELSPO through the PAI P6/27 Functional Supramolecular Systems project and the National Science Foundation for Research, FNRS. D.B. is an FNRS research director.


The authors declare no competing financial interests.

Author contributions The experiment was conceived by M.G., S.B., P.S. and E.O.; M.G. performed the STM characterization and the device study; S.B. carried out the optical characterization and M.-A.S. the PESA measurements; M.A.S. and M.G. performed the C-AFM experiments.; J.X.L. performed the MD/MM/DFT and G.D. the microelectrostatic calculations under the supervision of Y.O. and D.B.; Y.L. and X.-Y.W. synthesized the molecules under the supervision of A.N., K.M. and X.F.; P.S. and E.O. supervised the entire study. M.G. wrote the manuscript with comments and suggestions from all co-authors.


  • Geim A. K. & Grigorieva I. V. Van der Waals heterostructures. Nature 499, 419–425 (2013). [PubMed]
  • Novoselov K. S., Mishchenko A., Carvalho A. & Castro Neto A. H. 2D materials and van der Waals heterostructures. Science 353, aac9439 (2016). [PubMed]
  • Jariwala D., Marks T. J. & Hersam M. C. Mixed-dimensional van der Waals heterostructures. Nat. Mater. 16, 170–181 (2017). [PubMed]
  • Britnell L. et al. . Field-effect tunneling transistor based on vertical graphene heterostructures. Science 335, 947–950 (2012). [PubMed]
  • Lee C.-H. et al. . Atomically thin p–n junctions with van der Waals heterointerfaces. Nat. Nanotechnol. 9, 676–681 (2014). [PubMed]
  • Britnell L. et al. . Resonant tunnelling and negative differential conductance in graphene transistors. Nat. Commun. 4, 1794 (2013). [PMC free article] [PubMed]
  • Britnell L. et al. . Strong light-matter interactions in heterostructures of atomically thin films. Science 340, 1311–1314 (2013). [PubMed]
  • Withers F. et al. . Light-emitting diodes by band-structure engineering in van der Waals heterostructures. Nat. Mater. 14, 301–306 (2015). [PubMed]
  • Ponomarenko L. A. et al. . Cloning of Dirac fermions in graphene superlattices. Nature 497, 594–597 (2013). [PubMed]
  • Dean C. R. et al. . Hofstadter's butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598–602 (2013). [PubMed]
  • Hunt B. et al. . Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure. Science 340, 1427–1430 (2013). [PubMed]
  • Gorbachev R. V. et al. . Detecting topological currents in graphene superlattices. Science 346, 448–451 (2014). [PubMed]
  • Wang L. et al. . Evidence for a fractional fractal quantum Hall effect in graphene superlattices. Science 350, 1231–1234 (2015). [PubMed]
  • Yankowitz M. et al. . Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382–386 (2012).
  • Park C.-H., Yang L., Son Y.-W., Cohen M. L. & Louie S. G. Anisotropic behaviours of massless Dirac fermions in graphene under periodic potentials. Nat. Phys. 4, 213–217 (2008).
  • Schedin F. et al. . Detection of individual gas molecules adsorbed on graphene. Nat. Mater. 6, 652–655 (2007). [PubMed]
  • Wehling T. O. et al. . Molecular doping of graphene. Nano Lett. 8, 173–177 (2008). [PubMed]
  • Liu H., Liu Y. & Zhu D. Chemical doping of graphene. J. Mater. Chem. 21, 3335–3345 (2011).
  • Samuels A. J. & Carey J. D. Molecular doping and band-gap opening of bilayer graphene. ACS Nano 7, 2790–2799 (2013). [PubMed]
  • Cervenka J. et al. . Graphene field effect transistor as a probe of electronic structure and charge transfer at organic molecule-graphene interfaces. Nanoscale 7, 1471–1478 (2014). [PubMed]
  • Hapala P. et al. . Mapping the electrostatic force field of single molecules from high-resolution scanning probe images. Nat. Commun. 7, 11560 (2016). [PMC free article] [PubMed]
  • Wickenburg S. et al. . Tuning charge and correlation effects for a single molecule on a graphene device. Nat. Commun. 7, 13553 (2016). [PMC free article] [PubMed]
  • Prado M. C. et al. . Two-dimensional molecular crystals of phosphonic acids on graphene. ACS Nano 5, 394–398 (2011). [PubMed]
  • Wang Q. H. & Hersam M. C. Room-temperature molecular-resolution characterization of self-assembled organic monolayers on epitaxial graphene. Nat. Chem. 1, 206–211 (2009). [PubMed]
  • Deshpande A. et al. . Self-assembly and photopolymerization of sub-2 nm one-dimensional organic nanostructures on graphene. J. Am. Chem. Soc. 134, 16759–16764 (2012). [PubMed]
  • Zhang T. et al. . Self-assembled 1-octadecanethiol monolayers on graphene for mercury detection. Nano Lett. 10, 4738–4741 (2010). [PubMed]
  • Järvinen P. et al. . Molecular self-assembly on graphene on SiO2 and h-BN substrates. Nano Lett. 13, 3199–3204 (2013). [PubMed]
  • Riss A. et al. . Imaging and tuning molecular levels at the surface of a gated graphene device. ACS Nano 8, 5395–5401 (2014). [PMC free article] [PubMed]
  • Tsai H.-Z. et al. . Molecular self-assembly in a poorly screened environment: F4TCNQ on graphene/BN. ACS Nano 9, 12168–12173 (2015). [PMC free article] [PubMed]
  • Mali K. S., Greenwood J., Adisoejoso J., Phillipson R. & De Feyter S. Nanostructuring graphene for controlled and reproducible functionalization. Nanoscale 7, 1566–1585 (2015). [PubMed]
  • Alaboson J. M. P. et al. . Templating sub-10 nm atomic layer deposited oxide nanostructures on graphene via one-dimensional organic self-assembled monolayers. Nano Lett. 13, 5763–5770 (2013). [PubMed]
  • Macleod J. M. & Rosei F. Molecular self-assembly on graphene. Small 10, 1038–1049 (2014). [PubMed]
  • Zeng M. et al. . Self-assembly of graphene single crystals with uniform size and orientation: the first 2D super-ordered structure. J. Am. Chem. Soc. 138, 7812–7815 (2016). [PubMed]
  • Huang H., Chen S., Gao X., Chen W. & Wee A. T. S. Structural and electronic properties of PTCDA thin films on epitaxial graphene. ACS Nano 3, 3431–3436 (2009). [PubMed]
  • Zhou H. T. et al. . Direct imaging of intrinsic molecular orbitals using two-dimensional, epitaxially-grown, nanostructured graphene for study of single molecule and interactions. Appl. Phys. Lett. 99, 2009–2012 (2011).
  • Li B. et al. . Toward tunable doping in graphene FETs by molecular self-assembled monolayers. Nanoscale 5, 9640–9644 (2013). [PubMed]
  • Zhang X. et al. . Supramolecular chemistry on graphene field-effect transistors. Small 10, 1735–1740 (2014). [PubMed]
  • Yu Y. et al. . Epitaxially self-assembled alkane layers for graphene electronics. Adv. Mater. 29, 1603925 (2017). [PubMed]
  • Phillipson R. et al. . Tunable doping of graphene by using physisorbed self-assembled networks. Nanoscale 8, 20017–20026 (2016). [PubMed]
  • Vélez S. et al. . Gate-tunable diode and photovoltaic effect in an organic–2D layered material p–n junction. Nanoscale 7, 15442–15449 (2015). [PubMed]
  • Jariwala D. et al. . Hybrid, Gate-tunable, van der Waals p−n heterojunctions from pentacene and MoS2. Nano Lett. 16, 497–503 (2016). [PubMed]
  • Castellanos-Gomez A. et al. . Deterministic transfer of two-dimensional materials by all-dry viscoelastic stamping. 2D Mater. 1, 011002 (2014).
  • De Feyter S. & De Schryver F. C. Self-assembly at the liquid/solid interface: STM reveals. J. Phys. Chem. B 109, 4290–4302 (2005). [PubMed]
  • De Feyter S. & De Schryver F. C. Two-dimensional supramolecular self-assembly probed by scanning tunneling microscopy. Chem. Soc. Rev. 32, 139–150 (2003). [PubMed]
  • Lawrence E. J. et al. . 3-Aryl-3-(trifluoromethyl)diazirines as versatile photoactivated ‘linker' molecules for the improved covalent modification of graphitic and carbon nanotube surfaces. Chem. Mater. 23, 3740–3751 (2011).
  • Zhang Y., Brar V. W., Girit C., Zettl A. & Crommie M. F. Origin of spatial charge inhomogeneity in graphene. Nat. Phys. 5, 722–726 (2009).
  • Wehling T. O., Lichtenstein A. I. & Katsnelson M. I. First-principles studies of water adsorption on graphene: the role of the substrate. Appl. Phys. Lett. 93, 202110 (2008).
  • Massicotte M. et al. . Photo-thermionic effect in vertical graphene heterostructures. Nat. Commun. 7, 12174 (2016). [PMC free article] [PubMed]
  • Lin Y. C. et al. . Single-layer ReS2: two-dimensional semiconductor with tunable in-plane anisotropy. ACS Nano 9, 11249–11257 (2015). [PubMed]
  • Xia F., Wang H. & Jia Y. Rediscovering black phosphorus as an anisotropic layered material for optoelectronics and electronics. Nat. Commun. 5, 4458 (2014). [PubMed]
  • Cincotti S. & Rabe J. P. Self-assembled alkane monolayers on MoSe2 and MoS2. Appl. Phys. Lett. 62, 3531 (1993).

Articles from Nature Communications are provided here courtesy of Nature Publishing Group