Abstract
Electron-electron interactions in materials lead to exotic many-body quantum phenomena, including Mott metal-insulator transitions (MITs), magnetism, quantum spin liquids, and superconductivity. These phases depend on electronic band occupation and can be controlled via the chemical potential. Flat bands in two-dimensional (2D) and layered materials with a kagome lattice enhance electronic correlations. Although theoretically predicted, correlated-electron Mott insulating phases in monolayer 2D metal-organic frameworks (MOFs) with a kagome structure have not yet been realised experimentally. Here, we synthesise a 2D kagome MOF on a 2D insulator. Scanning tunnelling microscopy (STM) and spectroscopy reveal a MOF electronic energy gap of ∼200 meV, consistent with dynamical mean-field theory predictions of a Mott insulator. Combining template-induced (via work function variations of the substrate) and STM probe-induced gating, we locally tune the electron population of the MOF kagome bands and induce Mott MITs. These findings enable technologies based on electrostatic control of many-body quantum phases in 2D MOFs.
Similar content being viewed by others
Introduction
Strong electronic correlations arise in a material at specific electron fillings of its bands, provided that the on-site Coulomb repulsion (characterised by the Hubbard energy, U) is of the order of, or larger than, the bandwidth, W. These electronic correlations can result in a wide range of exotic many-body quantum phases. Examples include correlated insulating phases, quantum spin liquids, correlated magnetism, and superconductivity—phenomena that have been realised in monolayer transition metal-dichalcogenides1,2,3,4,5,6, twisted few-layer graphene7,8, inorganic kagome crystals9,10,11, and organic charge transfer salts12,13,14.
Tuning of the chemical potential via electrostatic gating can allow for control over such band electron filling, enabling reversible switching between correlated phases7. This makes these systems amenable to integration as active materials in voltage-controlled devices, offering enticing prospects for applications in electronics, spintronics, and information processing and storage12,15.
Two-dimensional (2D) materials have emerged as particularly promising candidates for realising strongly correlated phenomena as the absence of interlayer hopping and screening can contribute to decreasing W and increasing U4. Additionally, some 2D crystal geometries—such as the kagome structure—give rise to intrinsic flat electronic bands16,17. When these extremely narrow bands are half-filled, even weak Coulomb repulsion can open an energy gap and give rise to a Mott insulating phase12. Away from half-filling, the gap closes, and the system becomes metallic.
Metal-organic frameworks (MOFs) are a broad class of materials whose properties are highly tunable through careful selection of constituent organic molecules and metal atoms18. There has been growing interest in 2D MOFs for their electronic properties19,20,21. In particular, layered 2D MOF structures have been recently shown to host strongly correlated superconductivity22. Monolayer 2D MOFs have attracted attention for their magnetism23,24,25, including ferromagnetism resulting from exchange interaction between unpaired metal centre electrons26. However, despite theoretical predictions27,28, correlated Mott phases have not yet been realised experimentally in monolayer 2D MOFs.
Here, we demonstrate local electrostatic control over a Mott metal-insulator-transition (MIT) in a single-layer 2D kagome MOF, in excellent agreement with theoretical predictions.
Results
A monolayer 2D MOF on an atomically thin insulator
We synthesised the monolayer MOF—consisting of 9,10-dicyanoanthracene (DCA) molecules coordinated to copper (Cu) atoms—on monolayer hexagonal boron nitride (hBN) on Cu(111) (see Methods for sample preparation). A scanning tunnelling microscopy (STM) image of a crystalline single-layer MOF domain grown seamlessly across the hBN/Cu(111) substrate is shown in Fig. 1a. We observe some defects within the MOF domain, as well as some DCA-only regions (discussed in Supplementary Note 11). The long-range modulation of the MOF STM apparent height follows the hBN/Cu(111) moiré pattern, which arises due to a mismatch between the hBN and Cu(111) lattices (giving rise to pore, P, and wire, W, regions—see upper inset)29,30,31. This moiré pattern has been shown to affect the electronic properties of adsorbates31, including one previous example of a MOF32.
The MOF is characterised by a hexagonal lattice (lattice constant: 2.01 ± 0.06 nm), with a unit cell including two Cu atoms (honeycomb arrangement; bright protrusions in Fig. 1b) and three DCA molecules (kagome arrangement, with protrusions at both ends of anthracene backbone in STM image in Fig. 1b), similar to previous reports (including on a decoupling graphene monolayer)25,33,34,35.
We calculated the band structure of this monolayer DCA3Cu2 MOF on hBN/Cu(111) by density functional theory (DFT; with U = 0); Fig. 1d. Projection of the Kohn–Sham wavefunctions onto MOF states shows the prototypical kagome energy dispersion with two Dirac bands and a flat band, consistent with prior theoretical calculations for the freestanding MOF19,27,36. This near-Fermi band structure has predominantly molecular DCA character, and is well described by a nearest-neighbour tight-binding (TB) model (see the corresponding density of electronic states, DOS, as a function of energy E in Fig. 1e)36. The hBN monolayer, a 2D insulator with a bandgap >5 eV31, prevents electronic hybridisation between the underlying Cu(111) surface and the 2D MOF31. This allows the MOF to preserve its intrinsic electronic properties, in contrast to previous findings on metal surfaces25,35,36,37. These U = 0 calculations predict that the MOF on hBN/Cu(111) is metallic, with some electron transfer from substrate-to-MOF leading to the chemical potential lying above the Dirac point, close to half-filling of the three kagome bands (Fig. 1d).
Strong electronic interactions have been theoretically predicted in DCA3Cu227, and a signature was recently detected experimentally25. We, therefore, calculated the many-body spectral function A(E)—analogous to the DOS in the non-interacting regime—of the freestanding MOF via dynamical mean-field theory (DMFT). In contrast to the TB model or DFT, DMFT explicitly captures local electronic correlations caused by the Hubbard energy U (see Methods)38,39,40,41. In Fig. 1e, A(E) is shown for U = 0.65 eV (consistent with previous experimental estimates33) and for a chemical potential that matches the DFT-predicted occupation of the kagome bands for the MOF on hBN/Cu(111) (see Supplementary Note 2 and Supplementary Note 23 for further DMFT calculations, including temperature dependence). We observe two broad peaks (lower and upper Hubbard bands) separated by an energy gap of ∼200 meV, dramatically different from the non-interacting kagome DOS, and indicative of a Mott insulating phase27,42.
Observation of ∼200 meV Mott energy gap
To experimentally probe the electronic properties of DCA3Cu2/hBN/Cu(111), we conducted differential conductance (dI/dV) scanning tunnelling spectroscopy (STS); dI/dV is an approximation of the local DOS [A(E)] in the non-interacting (interacting, respectively) picture. We performed STS at the ends of the DCA anthracene moiety and at the Cu sites of the MOF—locations where we expect the strongest signature of the kagome bands based on the spatial distribution of the orbitals that give rise to these bands25,33,34. These spectra (Fig. 2a), taken at a pore region of the hBN/Cu(111) moiré pattern, both show broad peaks at bias voltages Vb ≈ −0.2 and 0.2 V. In a bias voltage window of ∼0.2 V around the chemical potential EF (Vb = 0), the dI/dV signal is low, significantly smaller than that for bare hBN/Cu(111).
STM images acquired within this low-dI/dV bias voltage window (Fig. 2b, c) show mainly the topography of the MOF, with the molecules appearing as ellipses of uniform intensity and the Cu atoms as weak protrusions. Outside the low-dI/dV bias voltage window (∣Vb∣ > 200 mV), Cu sites and the ends of the DCA anthracene moieties appear as bright protrusions (Fig. 2d, e), similar to the spatial distribution of the electronic orbitals of the DCA3Cu2 MOF associated with the near-Fermi kagome bands (right inset of Fig. 2d; see Supplementary Figs. 11 and 12 for more Vb-dependent STM images and dI/dV maps)25,33,34.
This suggests that the dI/dV peaks at ∣Vb∣ ≈ 0.2 V in Fig. 2a are related to intrinsic MOF electronic states near EF, with the low-dI/dV bias voltage window of ∼0.2 V around EF representing an energy gap, Eg, between these states. This is consistent with in-gap topographic STM imaging in Fig. 2b, c43. These dI/dV peaks cannot be attributed to inelastic tunnelling (e.g., MOF vibrational modes) as they are not always symmetric about EF (see Fig. 3). The gap is much larger than that predicted from spin-orbit coupling in such a system19, and the monocrystalline growth of the MOF domains (Fig. 1a) makes a large disorder-related gap unlikely. Furthermore, the gap is inconsistent with DFT and TB calculations (Fig. 1d, e).
The MOF spectra in Fig. 2a strongly resemble the DMFT-calculated spectral function A(E) in Fig. 1e, including an energy gap of the same magnitude between two similar peaks (Supplementary Fig. 9). This suggests that these dI/dV spectra are hallmarks of a Mott insulator.
Template-assisted Mott metal-insulator transition (MIT)
We further measured dI/dV spectra at Cu sites of the DCA3Cu2 MOF across the hBN/Cu(111) moiré pattern (Fig. 3a, b). In Fig. 3b, Eg is centred symmetrically about EF for spectra taken in the middle of a pore region, while those taken closer to the wire region show the Hubbard bands shifting upwards in energy (lowering the barrier to creation of a hole). At the centre of the wire region, the gap at the Fermi level vanishes with a clear increase in Fermi level dI/dV signal (Fig. 3e). The same behaviour was observed for DCA lobe sites of the MOF (see Supplementary Note 14).
The hBN/Cu(111) moiré pattern consists of a modulation of the local work function Φ (with little structural corrugation), where the quantity \(\Delta \Phi={\Phi }_{{{{{{{{\rm{wire}}}}}}}}}-{\Phi }_{{{{{{{{\rm{pore}}}}}}}}}\) depends on the period of the moiré superstructure, λ29,30,31. For the hBN/Cu(111) domain in Fig. 3a with λ ≈ 12.5 nm, ΔΦ ≈ 0.2 eV (Fig. 3c; see Supplementary Note 15 for moiré domains with different periods)30. Due to energy level alignment44,45,46, this corrugation of Φ affects substrate-to-MOF electron transfer and hence the effective electron filling of the MOF bands, with this filling smaller at wire than pore regions44,47. This is consistent with the effective reduction of the hole-creation barrier at the wire relative to the pore in Fig. 3b.
To capture the effect of this moiré-induced modulation of Φ on the MOF electronic properties, we conducted further DMFT calculations. Using U = 0.65 eV (the same as Fig. 1e), we calculated A(E) for a range of EF assuming a uniform system. We considered a sinusoidal variation of EF from a minimum value corresponding to half-filling of the kagome MOF bands, with an amplitude of 0.2 eV, to match the experimental ΔΦ for this specific hBN/Cu(111) moiré domain30,47 (Fig. 3c; see Methods). The obtained A(E) (Fig. 3d) reproduce the experimental spectral features in Fig. 3b, including the shifting of the lower Hubbard band maximum (LHBM) and upper Hubbard band minimum (UHBM) (Fig. 3f), and the vanishing of the gap and increase of the spectral function at the Fermi level for the wire region (Fig. 3e, g).
The DMFT-calculated spectral functions A(E) with no gap at the Fermi level at the smallest electron filling (ΔΦDMFT = 0.2 eV) show peaks near EF (Fig. 3d), however, which were not observed in the experimental spectra at the wire region (Fig. 3b). These peaks are indicative of coherent quasiparticles38, with their width associated with the quasiparticle lifetime and quasiparticle mean free path ℓ. Via our DMFT and TB calculations (Supplementary Note 3), we estimate ℓ ≈ 10 nm, much larger than the wire region width of ∼4 nm. We hypothesise that the coherence peaks are suppressed in the experiment as quasiparticles are strongly scattered by the pore regions (where the MOF remains insulating). This is a key difference between Fig. 3b, d: the experimental measurements represent changes in electron population at finite size MOF regions due to a locally varying work function, whereas each theoretical spectrum corresponds to an infinite uniform system with a constant chemical potential.
Furthermore, some of the experimental spectra in the vicinity of the wire region in Fig. 3b feature narrow peaks at energies of ∼0.4 eV, not present in the theoretical A(E) in Fig. 3d. We claim that these peaks do not represent intrinsic electronic states, but are instead related to charging phenomena not captured by the DMFT calculations (see Fig. 4 and Supplementary Note 17).
Despite these discrepancies with experiment, our DMFT calculations capture the fundamental electronic properties of the 2D DCA3Cu2 MOF, hosting a Mott insulating phase with Eg ≈ 200 meV, and a metal-like phase (with no gap at the Fermi level) at the wire region (for the specific tip-sample distance considered in Fig. 3). DMFT is a well-established method for understanding the Mott insulator and Mott MITs38,39,40,41. This claim of strong electronic correlations and of a Mott insulating phase for the DCA3Cu2 kagome MOF is consistent with previous literature25,27,36,42.
Tip-assisted Mott MIT at moiré wire region
To explore the nature of the metal-like phase observed at the wire region, we conducted further dI/dV STS of the MOF as a function of tip-sample distance Δz + z0 (where z0 is set by tunnelling parameters), at a DCA lobe site within the wire region (Fig. 4d). For large Δz + z0, these spectra feature an energy gap Eg, with a small dI/dV signal at the Fermi level (Fig. 4f), similar to spectra in the pore regions (Figs. 2a, 3b; Supplementary Note 20), with a sharp peak at positive Vb (purple circles in Fig. 4d) and a subtler band edge (red squares; similar to band features in Fig. 3b) at negative Vb. As Δz decreases, the energy position of the sharp peak decreases linearly (Fig. 4e). Conversely, the energy position of the subtler band edge increases with decreasing Δz, non-linearly and at a lower rate (Fig. 4e). These features cross the Fermi level at intermediate Δz, and the spectrum becomes gapless (Eg = 0), with a significant increase in dI/dV signal at the Fermi level (Fig. 4f). Note that an intermediate Δz was also used for all spectra in Fig. 3b where similar metal-like signatures were observed at the wire region. We found a similar Δz-dependent trend for wire region Cu sites (Supplementary Note 21). Notably, moiré pore regions remain gapped for all Δz values (for both Cu and DCA lobe sites; Supplementary Note 20 and Supplementary Note 21).
As Vb is applied between the tip and Cu substrate, the STM double-barrier tunnel junction (DBTJ)—where the vacuum between tip and MOF is a first tunnel barrier and the insulating hBN is a second one—causes a voltage drop at the MOF location. This can lead to energy shifts of MOF states and/or charging of such states when they become resonant with the Cu(111) Fermi level (Fig. 4a)33,46,48. These phenomena have previously been observed for the DCA3Cu2 MOF on a decoupling graphene surface33. In this scenario, the bias voltages corresponding to an intrinsic electronic state, Vstate, and to charging of such a state, Vcharge, vary as a function of Δz as:
where deff is the effective width of the hBN tunnel barrier, V∞ is the bias voltage corresponding to the electronic state as Δz → ∞, and ΔΦts is the difference between tip and sample work functions48.
We fit the Δz-dependent bias voltage associated with the subtle band edge (red squares) in Fig. 4d with Eq. (1), and the bias voltage associated with the sharp peak (purple circles) with Eq. (2) (Fig. 4e). The agreement between experimental data and fits indicate that the subtle spectral band edge (red) represents an intrinsic MOF electronic state, with its energy shifting as Δz varies, and with the sharp peak (purple) corresponding to charging of such a state.
Discussion
We interpret these results as follows. For large Δz (Fig. 4c), the MOF electronic states are strongly pinned to the substrate, and the MOF near-Fermi kagome bands are approximately half-filled. Here, the dI/dV spectra feature an energy gap at the Fermi level, both at moiré pore (Figs. 2a, 3b; Supplementary Note 20 and Supplementary Note 21) and wire regions (bottom spectra of Fig. 4d; Supplementary Note 21): the entire monolayer 2D MOF is intrinsically a Mott insulator featuring localised electrons. This is consistent with DMFT calculations, which show that the system can remain Mott insulating even with a 0.2 eV modulation of EF (corresponding to the experimental moiré work function corrugation; Supplementary Fig. 4).
Due to the hBN/Cu(111) moiré work function corrugation, the LHBM at a wire region is very close to the Fermi level (in comparison to the pore region). As Δz is reduced, the MOF states become less pinned to the Cu(111) substrate and more pinned to the tip via the DBTJ effect (Fig. 4; Supplementary Note 19, Supplementary Note 20, Supplementary Note 21). Given that \({\Phi }_{{{{{{{{\rm{tip}}}}}}}}} \, > \, {\Phi }_{{{{{{{{\rm{wire}}}}}}}}}\) (see Supplementary Note 18), this leads to an energy upshift of the MOF states with respect to the Cu(111) Fermi level (with the LHBM susceptible to charging as Vb becomes more positive). At intermediate Δz, this energy upshift depopulates the LHB and leads to the transition from the Mott insulating phase (only existing at half-filling of the three near-Fermi kagome bands; Fig. 1) to the metallic phase. This is concomitant with a dramatic change in the dI/dV spectra, including the vanishing of the energy gap and the increase in dI/dV signal at the Fermi level. As Δz is further reduced (Fig. 4a), the near-Fermi kagome bands become fully depopulated (with the bottom of these bands susceptible to charging as Vb becomes more negative; top spectra of Fig. 4d), and the MOF becomes a trivial insulator (with a gap between these near-Fermi bands and lower energy bands; Supplementary Fig. 1). The DBTJ effect also manifests itself at the pore region, with energy shifts of the LHBM and UHBM as Δz varies (Supplementary Note 20, Supplementary Note 21). Yet, due to the smaller \({\Phi }_{{{{{{{{\rm{pore}}}}}}}}}\) (Fig. 3c), the Fermi level lies close to the centre of the energy gap (Fig. 3b), making the Mott insulating phase robust for the considered Δz range, consistent with DMFT (Fig. 3d, Supplementary Figs. 2b, 4). In the wire region, the STM tip, via the DBTJ effect, in combination with the large \({\Phi }_{{{{{{{{\rm{wire}}}}}}}}}\), acts as a local electrostatic gate, switching the 2D MOF from Mott insulator to metal (Supplementary Fig. 26).
This tip-induced gating at the wire region is inherently local, occurring within the cross-section of the DBTJ (typically with a diameter of ∼10 nm given by the tip radius of curvature). This locality could lead to dI/dV features (e.g., due to electronic confinement) not captured by DMFT. Whether these local changes produce truly delocalised metallic states across an extended area of the sample remains an open question, beyond the scope of this investigation. Future work could investigate the DCA3Cu2 MOF on single-crystal exfoliated hBN, within a gated heterostructure, allowing for uniform electrostatic control and bulk (e.g., transport) measurements.
Our assertion that the MOF is a Mott insulator at the moiré pore regions for all tip-sample distances, and at the wire regions for large tip-sample distances, is well supported by DMFT calculations, and by previous literature25,27,36,42. The DMFT spectral functions demonstrate excellent agreement with dI/dV spectra at the pore regions, including the ∼200 meV gap at the Fermi level, and the energy modulation of the LHB and UHB due to variations in electrostatic potential. Our assertion of a metallic phase for the MOF at the wire region for intermediate tip-sample distances is well supported by the DBTJ model, and by the qualitative agreement with DMFT, including the absence of an energy gap at the Fermi level and an increase in Fermi level dI/dV signal and spectral function.
Mott insulating—with intrinsically localised electronic states—and metallic phases have been observed in TaS2 and TaSe2 monolayers, which feature frustrated lattice geometries and lattice constants (due to charge density wave distortions) similar to our work3,4,49,50,51, for finite domain sizes as small as ~10 × 10 nm2 49. Also, local topological phase transitions induced by an STM tip have been demonstrated52. Our interpretation of a local population-induced Mott metal-insulator transition is consistent with these findings.
Monolayer DCA-based MOFs have been studied on other substrates25,26,33,34,35,37,53,54,55, without observing a Mott phase. In our case, the combination of the wide bandgap hBN as a template (allowing the MOF to retain its intrinsic electronic properties), and of the adequate energy level alignment given by the hBN/Cu(111) substrate (resulting in half-filling of kagome bands; Fig. 1d), plays a key role in the realisation of the correlated-electron Mott phase.
Note that the vanishing of the energy gap at the Fermi level for the MOF metallic phase at the moiré wire region (Fig. 3b) can be reminiscent of spectral features for the pseudogap phase in cuprates and other transition metal oxides56,57,58. Future work on controllable Mott MITs in MOFs might shed light on the general understanding of doped Mott insulators.
We have demonstrated that single-layer DCA3Cu2 not only hosts a robust Mott insulating phase (with Eg ≫ kBT at T = 300 K), but also that Mott MITs can be achieved via the combination of template- (Fig. 3) and tip- (Fig. 4) induced gating, consistent with DMFT and the DBTJ model. This shows that such phase transitions can be controlled in monolayer MOFs via electrostatic tuning of the chemical potential.
Our findings represent a promising step towards incorporation of 2D MOFs as active materials in device-like architectures (e.g., van der Waals heterostructures based on 2D materials), benefiting from efficient synthesis approaches and versatility offered by MOFs, and allowing for access and control of correlated-electron phases therein via electrostatic gating59. Our work establishes single-layer 2D MOFs—with crystal geometries allowing for flat bands—as promising platforms for controllable switching between diverse many-body quantum phenomena, potentially including correlated magnetism, superconductivity, and quantum spin liquids.
Methods
Sample preparation
The monolayer (ML) DCA3Cu2 kagome MOF was synthesised on hBN/Cu(111) in UHV (base pressure ∼2 × 10−10 mbar). The Cu(111) surface was first cleaned via 2–3 cycles of sputtering with Ar+ ions and subsequent annealing at ∼770 K. A hBN ML was synthesised on Cu(111) via the thermal decomposition of borazine31. We dosed a partial pressure of borazine of ∼9 × 10−7 mbar for 45 minutes with the Cu(111) sample maintained at 1140 K. We kept the Cu(111) sample at this temperature for a further 20 mins to ensure a complete reaction. We then cooled the sample to room temperature and deposited the DCA molecules via sublimation at 390 K, corresponding to a deposition rate of 0.007 ML/sec. In our experiments we considered DCA coverages of ∼0.4–0.6 ML. We then further cooled the sample to ∼77 K before depositing Cu via sublimation at 1250 K (Cu deposition rate: ∼0.002 ML/sec; typical Cu coverages in our experiments: ∼0.05 ML). Finally, the sample was annealed to ∼200 K for 15 minutes. Further details are in Supplementary Note 6.
The DCA3Cu2 MOF crystalline structure was found to be commensurate with the hBN lattice but incommensurate with the long-range hBN/Cu(111) moiré patterns of different sizes, across which the MOF grows without disruption (see Supplementary Fig. 8). The hBN/Cu(111) moiré pattern is clearly visible in large-scale STM images of the MOF (Fig. 1a). This is consistent with the modulation of the MOF’s electronic properties illustrated in Fig. 3 (also see dI/dV maps in Supplementary Fig. 12).
STM and STS measurements
All STM and dI/dV STS measurements were performed at 4.5 K (except measurements in Supplementary Note 23 performed at 77 K), at a base pressure <1 × 10−10 mbar, with a hand-cut Pt/Ir tip. All STM images were acquired in constant-current mode with tunnelling parameters as reported in the text (bias voltage applied to sample). All dI/dV spectra were obtained by acquiring I(V) at a constant tip-sample distance (stabilised by a specified setpoint tunnelling current and bias voltage), and by then numerically differentiating I(V) to obtain dI/dV as a function of bias voltage. Tips were characterised on regions of bare hBN/Cu(111) prior to spectroscopy measurements, where the Shockley surface state of Cu(111) could be observed (grey curve in Fig. 2a, onset shifted upwards in energy due to confinement by hBN monolayer)29.
DFT calculations
We calculated the non-spin-polarised band structure of DCA3Cu2 on hBN on Cu(111) via DFT (Fig. 1d), using the Vienna Ab-Initio Simulation Package60 with the Perdew-Burke-Ernzerhof functional under the generalised gradient approximation61. We used projector augmented wave pseudopotentials62,63 to describe core electrons, and the semi-empirical potential DFT-D364 to describe van der Waals forces.
The substrate was modelled as a slab three Cu atoms thick, with the bottom layer fixed at the bulk lattice constant65. A layer of passivating hydrogen atoms was applied to the bottom face to terminate dangling bonds.
A 400 eV cut-off was used for the plane wave basis set. The threshold for energy convergence was 10−4 eV. The atomic positions of the DCA3Cu2/hBN/Cu(111) were relaxed until Hellmann-Feynman forces were <0.01 eV/Å, using a 3 × 3 × 1 k-point grid for sampling the Brillouin zone and 1st order Methfessel-Paxton smearing of 0.2 eV. The charge density for the relaxed structure was calculated using an 11 × 11 × 1 k-point grid, Blöchl tetrahedron interpolation, and dipole corrections. The band structure was determined non-self-consistently from the charge density.
Note that small (<1 Å) perturbations in the height of DCA3Cu2 above the hBN/Cu(111) do not appreciably affect the calculated band structure. As such, small perturbations in height related to the hBN/Cu(111) moiré pattern (of at most 0.7 Å66) were not captured by these calculations36.
DMFT calculations
We performed DMFT calculations on the freestanding DCA3Cu2 kagome MOF. We used the Hubbard model for a kagome lattice with nearest-neighbour hopping,
where the first term is the TB Hamiltonian with nearest-neighbour hopping energy t, ∑〈i, j〉 is a sum over nearest-neighbour sites, and the second term is the interaction Hamiltonian with on-site Coulomb repulsion U. The operator \({\hat{c}}_{i,\sigma }^{{{{\dagger}}} }\) (\({\hat{c}}_{i,\sigma }\)) creates (annihilates) an electron at site i with spin σ; \({\hat{n}}_{i,\sigma }={\hat{c}}_{i,\sigma }^{{{{\dagger}}} }{\hat{c}}_{i,\sigma }\) is the density operator. We take t = 0.05 eV to match prior DFT calculations of DCA3Cu219,27,36.
We first calculated the non-interacting (U = 0) density of states (DOS; blue curve in Fig. 1e) by numerically integrating over all momenta in the first Brillouin zone. The chemical potential EF in Fig. 1e was chosen to be consistent with the electron filling predicted by DFT (Fig. 1d). We applied a thermal broadening (kBT = 2.5 meV) to this non-interacting TB DOS to make it consistent with the thermal broadening of the DMFT-generated (see below) spectral function A(E) in Fig. 1e.
To account for electronic correlations, we then implemented the DMFT formalism38,67,68 using the Toolbox for Researching Interacting Quantum Systems (TRIQS)69, with the continuous-time hybridisation expansion solver (CTHYB)70,71 at a temperature of ∼29 K (kBT ≈ 0.05t, unless specified otherwise; see Supplementary Note 23 for temperature-dependent calculations), with U = 0.65 eV. To use a single-site DMFT formalism38 with the kagome band structure, the non-interacting DOS of the three kagome bands were combined into a single function for use as the input into the DMFT procedure.
We calculated the many-body spectral functions A(E) (analogous to the DOS, but in the interacting regime; Fig. 3d) via analytic continuation using the maximum entropy method (MaxEnt) as implemented72 in TRIQS. The meta-parameter, α, was determined from the maximum curvature of the distance between the MaxEnt fit and data, χ2, as a function of α73.
Each DMFT calculation assumed a spatially uniform work function Φ; long-range modulation of Φ is beyond the capabilities of DMFT. As such, the spatially varying sample work function Φ resulting from the experimental hBN/Cu(111) moiré pattern was not explicitly captured in the individual A(E) spectra in Fig. 3d. This spatial variation of Φ was approximated by varying the uniform EF of the system for each individual A(E) spectrum. Each of these A(E) spectra was then associated to a specific location of the hBN/Cu(111) moiré pattern (and hence to a specific experimental dI/dV curve) based on how this variation of EF would translate to a local Φ. These calculations assume that the theoretical A(E) spectra, calculated with a uniform EF (and hence uniform Φ), are reasonable representations of the locally acquired experimental dI/dV curves, which are affected by a spatially varying Φ. This assumption is reasonable for the insulating phase. Indeed, in the Mott insulating phase, electronic states are localised at the kagome sites, confined within areas that are small in length74 compared to the distance between nearest-neighbour kagome sites (∼1 nm) and to the periodicity λ of the hBN/Cu(111) moiré domains considered in our experiments (λ > 5 nm). DMFT indicates that this Mott insulating phase is robust to variations in chemical potential EF larger than 0.2 eV, with the spectral function A(E) shifting in energy as EF is varied within this range, without other significant qualitative changes (see Supplementary Fig. 2b for EF = 0.25 to 0.5 eV). In our experiments, the hBN/Cu(111) moiré pattern imposes a periodic modulation of the local work function Φ, with a peak-to-peak modulation amplitude of ∼0.2 eV and a modulation periodicity λ ≈ 12.5 nm (Fig. 3c; this amplitude becomes smaller with decreasing λ, see Supplementary Fig. 22). That is, Φ varies slowly across the MOF kagome lattice. The effect of such long-range modulation of Φ is to shift the energy of the MOF localised states accordingly. As long as the Φ modulation period is larger than the distance between nearest-neighbour kagome sites and the Φ modulation amplitude is smaller than a critical value inducing the transition to the metallic phase, there is no other dramatic qualitative effect on these localised electronic states. This explains the excellent agreement between experimental dI/dV spectra for the MOF at the hBN/Cu(111) moiré pore region and DMFT-calculated spectral functions A(E) for the system in the Mott insulating phase (Fig. 3b, d), with LHB and UHB modulated in energy following the variation in electrostatic potential. For the MOF in the metallic phase at the moiré wire region (Figs. 3, 4), discrepancies between theory and experiment can be explained (as discussed in the main text) by long electronic coherence lengths (Supplementary Note 3), and by effects of long-range moiré Φ modulation and of finite DBTJ cross section on the potentially delocalised metallic MOF states.
DBTJ model
In Fig. 4d, e, the bias voltage, Vstate, corresponds to the energy level of an intrinsic MOF frontier electronic state (which is susceptible to charging), at a band edge in either the Mott or trivial insulator regime. Red square markers indicating Vstate were assigned based on the method outlined in Supplementary Note 13. The bias voltage, Vcharge, corresponds to the peak associated with charging of such a MOF state. Purple circle markers indicating Vcharge were assigned by finding a local maximum in dI/dV.
In Fig. 4e, we considered four experimental datasets for fitting with Eqs. (1) and (2): Vstate(Δz) (red squares) and Vcharge(Δz) (purple circles) for small values of Δz (trivial insulator phase), and Vstate(Δz) and Vcharge(Δz) for large values of Δz (Mott insulator phase). Given the Mott MIT, we considered two different intrinsic MOF band edges susceptible to charging, embodied in two different values of V∞: one for the trivial insulator phase (small values of Δz) and one for the Mott insulator phase (large Δz). This phase transition is evident from the offset in Vstate (red squares in Fig. 4e) observed when Δz varies from small to large (through the metallic phase at intermediate Δz). Accordingly, we used a global fitting approach to obtain the same fitting parameters deff, z0, and ΔΦts (characteristic of the DBTJ and the acquisition location) for these four experimental datasets, and a separate V∞ value for each regime.
Schematics in Fig. 4a–c represent cartoon illustrations of our proposed physical mechanism for tip-induced gating.
It is important to note that the DBTJ inherently affects all dI/dV measurements in this work, including those in Figs. 2 and 3, for both pore and wire regions of the hBN/Cu(111) moiré pattern. The DBTJ effect does not cause phase transitions at the pore regions, however (see Supplementary Note 20, Supplementary Note 21). The measurements in Fig. 3 were performed with an intermediate tip-sample distance—which is why metallic properties were observed at the wire region.
Data availability
The data supporting the findings of this study is available from the authors upon request.
Code availability
All codes relating to DMFT and MaxEnt are available at https://doi.org/10.5281/zenodo.7439858. The authors can provide code for data analysis and theoretical calculations upon request.
References
Nakata, Y. et al. Monolayer 1T-NbSe2 as a Mott insulator. NPG Asia Mater. 8, e321 (2016).
Liu, M. et al. Monolayer 1T-NbSe2 as a 2D-correlated magnetic insulator. Sci. Adv. 7, eabi6339 (2021).
Lin, H. et al. Scanning tunneling spectroscopic study of monolayer 1T-TaS2 and 1T-TaSe2. Nano Res. 13, 133–137 (2020).
Chen, Y. et al. Strong correlations and orbital texture in single-layer 1T-TaSe2. Nat. Phys. 16, 218–224 (2020).
Ruan, W. et al. Evidence for quantum spin liquid behaviour in single-layer 1T-TaSe2 from scanning tunnelling microscopy. Nat. Phys. 17, 1154–1161 (2021).
Chen, Y. et al. Evidence for a spinon Kondo effect in cobalt atoms on single-layer 1T-TaSe2. Nat. Phys. 18, 1335–1340 (2022).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).
Park, J. M., Cao, Y., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Tunable strongly coupled superconductivity in magic-angle twisted trilayer graphene. Nature 590, 249–255 (2021).
Zheng, L. et al. Emergent charge order in pressurized kagome superconductor CsV3Sb5. Nature 611, 682–687 (2022).
Ye, L. et al. Massive Dirac fermions in a ferromagnetic kagome metal. Nature 555, 638–642 (2018).
Yin, J.-X. et al. Quantum-limit Chern topological magnetism in TbMn6Sn6. Nature 583, 533–536 (2020).
Powell, B. J. & McKenzie, R. H. Quantum frustration in organic Mott insulators: from spin liquids to unconventional superconductors. Rep. Prog. Phys. 74, 056501 (2011).
Kanoda, K. & Kato, R. Mott physics in organic conductors with triangular lattices. Annu. Rev. Condens. Matter Phys. 2, 167–188 (2011).
Miksch, B. et al. Gapped magnetic ground state in quantum spin liquid candidate κ-(BEDT-TTF)2Cu2(CN)3. Science 372, 276–279 (2021).
Shao, Z., Cao, X., Luo, H. & Jin, P. Recent progress in the phase-transition mechanism and modulation of vanadium dioxide materials. NPG Asia Mater. 10, 581–605 (2018).
Balents, L., Dean, C. R., Efetov, D. K. & Young, A. F. Superconductivity and strong correlations in moiré flat bands. Nat. Phys. 16, 725–733 (2020).
Neupert, T., Denner, M. M., Yin, J.-X., Thomale, R. & Hasan, M. Z. Charge order and superconductivity in kagome materials. Nat. Phys. 18, 137–143 (2022).
Cui, Y. et al. Metal-organic frameworks as platforms for functional materials. Acc. Chem. Res. 49, 483–493 (2016).
Zhang, L. Z. et al. Intrinsic two-dimensional organic topological insulators in metal-dicyanoanthracene lattices. Nano Lett. 16, 2072–2075 (2016).
Jin, Y. et al. Large-gap quantum anomalous Hall phase in hexagonal organometallic frameworks. Phys. Rev. B 98, 245127 (2018).
Jiang, W., Ni, X. & Liu, F. Exotic topological bands and quantum states in metal-organic and covalent-organic frameworks. Acc. Chem. Res. 54, 416–426 (2021).
Takenaka, T. et al. Strongly correlated superconductivity in a copper-based metal-organic framework with a perfect kagome lattice. Sci. Adv. 7, eabf3996 (2021).
Yamada, M. G., Fujita, H. & Oshikawa, M. Designing Kitaev spin liquids in metal-organic frameworks. Phys. Rev. Lett. 119, 057202 (2017).
Zhang, L.-C. et al. Two-dimensional magnetic metal-organic frameworks with the Shastry-Sutherland lattice. Chem. Sci. 10, 10381–10387 (2019).
Kumar, D. et al. Manifestation of strongly correlated electrons in a 2D kagome metal-organic framework. Adv. Funct. Mater. 31, 2106474 (2021).
Lobo-Checa, J. et al. Ferromagnetism on an atom-thick & extended 2D metal-organic coordination network. Nat. Commun. 15, 1858 (2024).
Fuchs, M. et al. Kagome metal-organic frameworks as a platform for strongly correlated electrons. J. Phys. Mater. 3, 025001 (2020).
Nourse, H. L., McKenzie, R. H. & Powell, B. J. Multiple insulating states due to the interplay of strong correlations and lattice geometry in a single-orbital Hubbard model. Phys. Rev. B 103, L081114 (2021).
Joshi, S. et al. Boron nitride on Cu(111): an electronically corrugated monolayer. Nano Lett. 12, 5821–5828 (2012).
Zhang, Q. et al. Tuning band gap and work function modulations in monolayer hBN/Cu(111) heterostructures with moiré patterns. ACS Nano 12, 9355–9362 (2018).
Auwärter, W. Hexagonal boron nitride monolayers on metal supports: versatile templates for atoms, molecules and nanostructures. Surf. Sci. Rep. 74, 1–95 (2019).
Urgel, J. I. et al. Controlling coordination reactions and assembly on a Cu(111) supported boron nitride monolayer. J. Am. Chem. Soc. 137, 2420–2423 (2015).
Yan, L. et al. Synthesis and local probe gating of a monolayer metal-organic framework. Adv. Funct. Mater. 31, 2100519 (2021).
Yan, L. et al. Two-dimensional metal-organic framework on superconducting NbSe2. ACS Nano 15, 17813–17819 (2021).
Hernández-López, L. et al. Searching for kagome multi-bands and edge states in a predicted organic topological insulator. Nanoscale 13, 5216–5223 (2021).
Field, B., Schiffrin, A. & Medhekar, N. V. Correlation-induced magnetism in substrate-supported 2D metal-organic frameworks. npj Comput Mater. 8, 1–10 (2022).
Pawin, G. et al. A surface coordination network based on substrate-derived metal adatoms with local charge excess. Angew. Chem. Int. Ed. 47, 8442–8445 (2008).
Georges, A., Kotliar, G., Krauth, W. & Rozenberg, M. J. Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions. Rev. Mod. Phys. 68, 13–125 (1996).
Kotliar, G. & Vollhardt, D. Strongly correlated materials: insights from dynamical mean-field theory. Phys. Today 57, 53–59 (2004).
Vollhardt, D. Dynamical mean-field theory for correlated electrons. Ann. Phys. 524, 1–19 (2012).
Adler, R., Kang, C.-J., Yee, C.-H. & Kotliar, G. Correlated materials design: prospects and challenges. Rep. Prog. Phys. 82, 012504 (2018).
Ohashi, T., Kawakami, N. & Tsunetsugu, H. Mott transition in kagome lattice hubbard model. Phys. Rev. Lett. 97, 066401 (2006).
Repp, J., Meyer, G., Stojković, S. M., Gourdon, A. & Joachim, C. Molecules on insulating films: scanning-tunneling microscopy imaging of individual molecular orbitals. Phys. Rev. Lett. 94, 026803 (2005).
Joshi, S. et al. Control of molecular organization and energy level alignment by an electronically nanopatterned boron nitride template. ACS Nano 8, 430–442 (2014).
Zimmermann, D. M. et al. Self-assembly and spectroscopic fingerprints of photoactive pyrenyl tectons on hBN/Cu(111). Beilstein J. Nanotechnol. 11, 1470–1483 (2020).
Pörtner, M. et al. Charge state control of F16CoPc on h-BN/Cu(111). Adv. Mater. Interfaces 7, 2000080 (2020).
Kumar, D., Hellerstedt, J., Lowe, B. & Schiffrin, A. Mesoscopic 2D molecular self-assembly on an insulator. Nanotechnology 34, 205601 (2023).
Kumar, D., Krull, C., Yin, Y., Medhekar, N. V. & Schiffrin, A. Electric field control of molecular charge state in a single-component 2D organic nanoarray. ACS Nano 13, 11882–11890 (2019).
Vaňo, V. et al. Artificial heavy fermions in a van der Waals heterostructure. Nature 599, 582–586 (2021).
Bu, K. et al. Possible strain induced Mott gap collapse in 1T-TaS2. Commun. Phys. 2, 1–7 (2019).
Fei, Y., Wu, Z., Zhang, W. & Yin, Y. Understanding the Mott insulating state in 1T-TaS2 and 1T-TaSe2. AAPPS Bull. 32, 20 (2022).
Collins, J. L. et al. Electric-field-tuned topological phase transition in ultrathin Na3Bi. Nature 564, 390–394 (2018).
Kumar, A., Banerjee, K., Foster, A. S. & Liljeroth, P. Two-dimensional band structure in honeycomb metal-organic frameworks. Nano Lett. 18, 5596–5602 (2018).
Yan, L., Pohjavirta, I., Alldritt, B. & Liljeroth, P. On-surface assembly of Au-dicyanoanthracene coordination structures on Au(111). ChemPhysChem 20, 2297–2300 (2019).
Vaňo, V. et al. Emergence of exotic spin texture in supramolecular metal complexes on a 2d superconductor. https://arxiv.org/abs/2309.02537v1 (2023).
Cai, P. et al. Visualizing the evolution from the Mott insulator to a charge-ordered insulator in lightly doped cuprates. Nat. Phys. 12, 1047–1051 (2016).
Battisti, I. et al. Universality of pseudogap and emergent order in lightly doped Mott insulators. Nat. Phys. 13, 21–25 (2017).
Zhong, Y. et al. Direct visualization of ambipolar Mott transition in cuprate CuO2 planes. Phys. Rev. Lett. 125, 077002 (2020).
Riss, A. et al. Imaging and tuning molecular levels at the surface of a gated graphene device. ACS Nano 8, 5395–5401 (2014).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
Blöchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758–1775 (1999).
Grimme, S., Antony, J., Ehrlich, S. & Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 132, 154104 (2010).
Rumble, J. R. CRC handbook of chemistry and physics, 100 edn (CRC Press/Taylor & Francis, Boca Raton, FL, 2019).
Schwarz, M. et al. Corrugation in the weakly interacting hexagonal-BN/Cu(111) system: structure determination by combining noncontact atomic force microscopy and X-ray standing waves. ACS Nano 11, 9151–9161 (2017).
Georges, A. & Kotliar, G. Hubbard model in infinite dimensions. Phys. Rev. B 45, 6479–6483 (1992).
Georges, A. Strongly correlated electron materials: dynamical mean-field theory and electronic structure. AIP Conf. Proc. 715, 3–74 (2004).
Parcollet, O. et al. TRIQS: a toolbox for research on interacting quantum systems. Comput. Phys. Commun. 196, 398–415 (2015).
Seth, P., Krivenko, I., Ferrero, M. & Parcollet, O. TRIQS/CTHYB: a continuous-time quantum Monte Carlo hybridisation expansion solver for quantum impurity problems. Comput. Phys. Commun. 200, 274–284 (2016).
Gull, E. et al. Continuous-time Monte Carlo methods for quantum impurity models. Rev. Mod. Phys. 83, 349–404 (2011).
Kraberger, G. J., Triebl, R., Zingl, M. & Aichhorn, M. Maximum entropy formalism for the analytic continuation of matrix-valued Green’s functions. Phys. Rev. B 96, 155128 (2017).
Bergeron, D. & Tremblay, A.-M. S. Algorithms for optimized maximum entropy and diagnostic tools for analytic continuation. Phys. Rev. E 94, 023303 (2016).
Fazekas, P. Lecture notes on electron correlation and magnetism, https://www.worldscientific.com/doi/abs/10.1142/2945 (WORLD SCIENTIFIC, 1999).
Acknowledgements
A.S. acknowledges funding support from the ARC Future Fellowship scheme (FT150100426). B.J.P. acknowledges funding support from the ARC Discovery Project scheme (DP180101483). H.L.N. acknowledges funding support from the MEXT Quantum Leap Flagship Programme (JPMXS0118069605). B.L., J.H., J.C., and N.V.M. acknowledge funding support from the Australian Research Council (ARC) Centre of Excellence in Future Low-Energy Electronics Technologies (CE170100039). B.L., B.F., and J.C. are supported through Australian Government Research Training Programme (RTP) Scholarships. B.F. and N.V.M. gratefully acknowledge the computational support from the National Computing Infrastructure and Pawsey Supercomputing Facility. The authors also thank Prof. Michael S. Fuhrer, Prof. Jaime Merino Troncoso, and Dr. Daniel Moreno Cerrada for their valuable discussions.
Author information
Authors and Affiliations
Contributions
B.L., J.H., and A.S. conceived and designed the experiments. B.L., J.H., and J.C. performed the experiments. B.F. and H.L.N. performed the theoretical calculations with guidance from B.J.P. and N.V.M. All authors contributed to writing the manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Ignacio Piquero-Zulaica, Amadeo L. Vázquez de Parga, and the other, anonymous, reviewer for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lowe, B., Field, B., Hellerstedt, J. et al. Local gate control of Mott metal-insulator transition in a 2D metal-organic framework. Nat Commun 15, 3559 (2024). https://doi.org/10.1038/s41467-024-47766-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-47766-8