Orbital degree of freedom in artificial electron lattices on metal surface

Orbital degree of freedom in artificial electron lattices on metal surface


Orbital degree of freedom plays a fundamental role in condensed matter physics. Recently, a new kind of artificial electron lattice has been realized in experiments by confining the metal surface electrons with adsorbed molecular lattice. A most recent example is the Lieb lattice realized by CO adsorption on Cu(111) surface [M. R. Slot, et al., Nat. Phys. 13, 672 (2017)]. The Lieb lattice is of special interest due to its flat band physics. Here, by first-principles calculations, muffin-tin potential model and tight binding model, we demonstrate that, the high energy states observed in the experiment actually correspond to the artificial -orbitals of the electron lattice. Our numerical results, together with the experimental observation, show that artificial -orbital fermionic lattice has already been realized in solid state system. This opens a new avenue to investigate the orbital degree of freedom in a controllable way.


The orbital degree of freedom, which refers to the orbital degeneracy and orientational anisotropy, is a fundamental attribute of the Bloch electrons in crystal, in addition to their charge and spin. When coupled to the charge and spin, orbital degree of freedom of electrons can give rise to many important phenomena, such as metal-insulator transition, superconductivity, and colossal magneto-resistanceTokura and Nagaosa (2000); Hotta (2006). However, understanding the orbital physics in real materials is still a big challenge due to some realistic reasons, e.g, the constraints of materials, the coupling among multiple degrees of freedom and the lack of controllability.

Utilizing artificial quantum systems is a promising way to study the orbital physics. The most successful example is the cold atoms in optical lattices, where atoms can be excited into higher orbital states of the optical lattice. In last decade, -orbital related novel quantum states in optical lattices have been intensively studied (see, for example, Ref. Li and Liu, 2016 and the references therein). For example, for -band bosons, an unconventional BEC has been realized in experimentMüller et al. (2007); Wirth et al. (2011); Wu (2009). However, an experimental realization of -band fermions has not been reported in any artificial quantum lattices so far.

Recently, a new kind of artificial two dimensional electron lattice has been realized in experiments. The first example is the molecular grapheneGomes et al. (2012); Wang et al. (2014), where the Cu surface electrons, a perfect two dimensional electron gas (2DEG) with dispersion, is confined in a honeycomb lattice by lateral periodic potential induced by adsorbed molecules. Even more exciting is the realization of artificial Lieb latticeQiu et al. (2016); Slot et al. (2017); Drost et al. (2017), that has not been found in natural materials and is of special interest due to its flat band physics. These pioneering works demonstrate that artificial electron lattice on metal surface could be a promising solid-state quantum simulation platform.

An interesting issue in the recent experimentSlot et al. (2017) is that some high-energy electron states with complex local density of states (LDOS) pattern are observed in the scanning tunnelling microscope (STM), which are drastically different from the low energy states. Here, combining first-principles simulation based on Density Functional Theory (DFT), muffin-tin potential model and tight-binding calculations, we demonstrate that these high energy states originate from artificial -orbitals of the electron lattice. Through carefully comparison of our theoretical calculation with the experimental observation, we show that the experiment in Ref. Slot et al., 2017 has actually realized an artificial -orbital square (and Lieb) electron lattice in solid state system. It may give the first example of the two dimensional -orbital fermion lattice in artificial quantum systems, and offers an ideal solid state platform to study orbital physics in a controllable way.

Figure 1: (Color online) (a) Schematic of the artificial electron square lattice on Cu(111) surface. The red discs are the -orbitals and the green lobes denote the -orbitals. The lattice vectors are nm and nm, respectively. (b),(c),(d) are the DFT simulated STM images, where the energy region, with respect to the Fermi level, are eV, eV and eV, respectively. (e) Schematic of the -orbital tight binding model on square lattice. (f),(g),(h) are the LDOS calculated from the muffin-tin potential model, where the energy are eV, eV and eV, respectively. nm, nm, eV, nm. Scale bars, 1 nm.

Let us start with the square electron lattice on metal surface. The structure is shown in Fig. 1 (a), where the CO molecules (black balls) are arranged on the Cu(111) surface to form a square lattice. Note that, due to the geometry of the Cu(111) surface, the lattice constant (CO molecule lattice) in direction is slightly different from that in direction . For simplicity, we ignore this difference in the muffin-tin potential model unless specified otherwise. Here, each CO molecule exerts a repulsive potential on the Cu surface electrons, which can be described by a muffin-tin potential [see Fig. 1 (a), inside the blue circles and zero elsewhere, is the diameter of circular potential]. The Hamiltonian of the Cu surface states now is


where is the effective mass of the Cu surface electrons, and eV, nm are reasonable values for the muffin-tin potential of CO/Cu(111) systemLi et al. (2016); Ropo et al. (2014). It is actually an anti-dot lattice, since the surface electrons under the CO molecules are depleted by the repulsive potential . Consequently, surface electrons are forced into a square lattice, where the electron sites are in the center of the squares formed by four adjacent CO molecules. Theoretically, the corresponding band structure can be obtained by solving the Hamiltonian of the muffin-tin potential model with plane wave basis. Here, the LDOS of surface electrons is the most important quantity, which actually represents the electron wave functions and can be directly measured by STM. Theoretically, the LDOS can be calculated from


where is the wave function of surface electron, , and are the indices of band, momentum, and spin, respectively.

In addition to the muffin-tin potential model, we perform DFT calculations to direct simulate the experimentally realized systemRopo et al. (2014). We use the Vienna - Simulation Package (VASP), which is based on the projector-augmented wave method and a plane wave basis setKresse and Furthmüller (1996). We choose the Perdew-Burke-Ernzerhof (PBE) version of the generalized gradient approximationPerdew et al. (1996), and the energy cutoff is 400 eV. A four-layer slab is used to model the Cu(111) surface and a 10 Å vacuum region separates slabs between nearest supercells to avoid their interaction. The positions of the CO molecules and the top layer Cu atoms are optimized in a smaller supercell and used in present calculation. We can directly simulate the STM image (LDOS of surface electron states) using the Tersoff-Hamann schemeTersoff and Hamann (1983) from the DFT electronic structure.

Figure 2: (Color online) (a) Energy bands of the electron square lattice, where green lines are the -bands, blue lines are the -bands. eV, nm, nm. (b) Fitting the -bands in (a) (blue, solid) with a -orbital tight binding model on square lattice (red, dashed). PW means the muffin-tin model with plane wave basis. []=[0.55, -0.0275, 0.0131, 0.0963] eV. eV. (c) Energy bands of the electron square lattice in the atomic limit. nm, eV and nm. (d) Fitting the -bands in (c) (blue, solid) with the tight binding model (red, dashed). []=[0.21, -0.0055, 0.0026, 0.0014] eV, eV. (e) LDOS of the -bands in atomic limit, eV. (f) LDOS of the -bands in atomic limit at eV. Scale bars, 1 nm.

First, we show that the LDOS of this artificial square lattice observed in experiment can be well reproduced by the both muffin-tin potential model and the first-principles calculation. The numerical results are summarized in Fig.1, where Fig. 1 (b),(c),(d) are results of first-principles calculation, and Fig. 1 (f),(g),(h) are those from the muffin-tin potential model. Basically, three typical surface electron LDOS patterns at different energies have been observed in experiment: (1) at low energy [Fig. 1 (b),(f)], the electron states are localized around the lattice sites, so we get a square lattice pattern; (2) for a high energy state [Fig. 1 (c),(g)], the electron states are now located in between the lattice sites; (3) increasing the energy further, though the electron states are still mainly in between the lattice sites, some fine structures of the LDOS pattern appear, i.e., we find a LDOS node between two adjacent lattice sites [Fig. 1 (d),(h)]. Both DFT and the muffin-tin potential model nicely reproduce the experimental observations [Fig.  (e), (f), (g) in Ref. Slot et al., 2017].

The main focus of this work is to understand the physical origin of these complex electron LDOS (wave functions) pattern. Our viewpoint here is that, the low energy states are from the -band of the electron square lattice, while the high energy states are from the -bands. In other words, even for this anti-dot lattice, the orbital degree of freedom is still a valid concept. Here, each electron lattice site could be considered as an artificial atom in two dimensions with various atomic orbitals, such as -and -orbitals (only and here since it is a two dimensional system). These artificial orbitals form the corresponding energy bands due to hopping between them. Thus, Fig. 1 (b-d) and (f-h) actually show the real space distribution of these orbital bands. Importantly, in Fig. 1 (c), we see that the -orbitals of two neighboring sites form a bond, which is localized in between the sites. Thus the STM measurements actually have observed the artificial -orbital bonds. Meanwhile, at higher energy, the -orbitals form anti-bonding bonds, which have function nodes in between the artificial atoms. This is the origin of the LDOS nodes observed in the STM measurement.

Given the success of our muffin-tin potential model, we further use it to calculate the energy bands of the square lattice in order to prove the above physical picture. We show that the calculated band structure can be well understood from the orbital point of view. The calculated bands are shown in Fig. 2 (a), where we set nm and eV, nm are the parameters for CO/Cu(111) system. Note that is the distance between molecules which is tunable in experiment. Different gives different hopping amplitude, but the corresponding band shape is similar. We see that there is one band with the lowest energy (green solid line) gapped from the others. This is the -band, because in square lattice each site contributes one -orbital and these -orbitals form one -band with lowest energy. The -orbitals are around the lattice sites, and thus should give an electron square lattice. In our calculation, in the energy interval of the -bands, the corresponding LDOS are all like Fig. 1 (b), which show a clear square lattice LDOS pattern and are consistent with the experimental observation. Note that the -band from the muffin-tin potential calculation can be well described by a single band tight-binding model on square lattice. Similarly, the upper two bands above the -band (blue solid lines) should be the -bands, because each 2D atom has two -orbitals, i.e. and . We use a two-band (, ) tight binding model on square lattice, which is illustrated in Fig. 1 (e), to fit the -bands. The tight binding Hamiltonian is




Here, . The nearest neighbour hopping , and the next nearest neighbour hopping , are considered, which are illustrated in Fig. 1 (e). The fitting is quite well as shown in Fig.2 (b), except that the tight binding model can not reproduce the gap at the point. We attribute this difference to the band overlap with the -bands. As can be seen in Fig. 2 (a), there is no gap between the - and -bands. This means that there should be some overlap between them, which is not included in the tight binding model.

Figure 3: (Color online) (a) Schematic of Lieb lattice. (b),(c),(d) are the DFT simulated STM image (LDOS), where the energy regions are eV, eV and eV, respectively. nm, nm. (e) Schematic of the -orbital tight binding model on Lieb lattice. The -orbitals with different energy are denoted with different colors. (f),(g),(h) are the calculated LDOS with muffin-tin potential model, the energy of which are eV, eV, and eV, respectively. nm, nm, eV, nm. Scale bars, 1 nm.

In order to make the -orbital picture more clear, we further consider an extreme case in the muffin-tin potential model, where we set nm, and use a fictional muffin-tin potential eV and nm. With given and very large values of and (), the hopping between adjacent sites is greatly suppressed, so that the artificial atoms are nearly isolated. Thus, it is just the atomic limit. The band structure in this atomic limit is given in Fig. 2 (c). It is now obvious that the two -bands are separated from the higher -bands, i.e. there is no band overlap any more. As a consequence, the gap of the -bands at the point disappears . The fitting of -bands is given in Fig. 2 (d). The tight binding model now can give a perfect description of the -bands.

We plot the LDOS pattern in the atomic limit in Fig. 2 (e) and (f). At low energy [Fig. 2 (e)], we see that the isotropic -orbitals are well separated from each other. Increasing the energy to the -band region, we get an LDOS pattern in Fig. 2 (f), resembling the shape of isolated -orbitals. All the discussions above give a clear illustration of the -orbital picture in the square lattice.

Another important issue is the position of Fermi level relative to the -bands, which depends on the square lattice constant Gomes et al. (2012); Li et al. (2016); Qiu et al. (2016). We can give an estimation about the value of to access the -bands in square lattice. In square lattice, the number of electrons in each site is about , where is the electron density of the metal surface. For Cu(111) surface, is about nm. Here, we assume that the adsorbed CO molecules do not modify the surface electron number of Cu. It is reasonable for the experimental situation, since we do not expect that one hundred CO molecules can change the Fermi level of the bulk Cu crystal. To access the -bands, should be in the region nm nm. This estimation is consistent with the experimental observations. In experimentSlot et al. (2017), two square lattices are realized. The of large one is about nm ( nm), two times larger than that of the small one. So, according to the estimation above, the Fermi level of the large square lattice should be at the -bands, while that of the small lattice is at the -bands. In the STM measurement, for the large square lattice, the LDOS patterns of the -band and bonding -band are found below the Fermi level, while that of the anti-bonding -band is above the Fermi level. This is consistent with our estimation. In contrast, for the small square lattice, no matter with positive or negative bias voltage, only the LDOS pattern of the -band is observed , indicating that the Fermi level is at the -band. In order to access the -band, more larger bias voltage is needed for the small square lattice case.

Figure 4: (Color online) (a) Energy bands of the Lieb lattice. nm, nm, eV. (b) Fitting the -bands in (a) (blue, solid) with the -orbital tight binding model (red, dashed) on Lieb lattice. [,,,,]=[0.63,-0.0345,0.011,0.1563,0.096] eV. eV, eV, eV. (c) Energy bands of the Lieb lattice in the atomic limit. nm, nm, eV. (d) Fitting the -bands in (c) (blue, solid) with the -orbital tight binding model (red, dashed) on Lieb lattice. [,,,,]=[0.175,-0.0055,0.0026,0.0016,0.0019] eV. eV, eV, eV. (e)-(h) The LDOS of the Lieb lattice in the atomic limit. The energies are taken at eV, eV, eV, eV, respectively. Scale bars, 1 nm.

Now, we turn to the Lieb lattice. As shown in Fig. 3 (a), there are three lattice sites (A, B,C) in each unit cell of a Lieb lattice. Thus, there are three -bands with lowest energy, of which the middle one is a flat band if only the nearest neighbor hopping is considered. A unique property of this flat band is its special LDOS pattern. The electrons in the flat band are only localized at the B and C sites, while the electrons in other two -bands mainly distribute around the A sites. This unique LDOS pattern actually reflects the flat band localization phenomenon in Lieb lattice, and has been observed in different systemsQiu et al. (2016); Slot et al. (2017); Vicencio et al. (2015); Mukherjee et al. (2015).

We first show that, for the -bands, the special LDOS pattern of the flat band (the flat band localization phenomenon) can be well reproduced by the our calculation. The DFT results are shown in Fig. 3 (b) and (c), while the corresponding muffin-tin results are given in Fig. 3 (f) and (g) as a comparison. In Fig. 3 (b), it is seen that the electron states are absent around the A site, which is just the special LDOS pattern of the Lieb flat band. Increasing the energy to the upper -band, the electrons should distribute mainly around the A site, instead of the B and C sites. This is observed in Fig. 3 (c). The muffin-tin potential model can also reproduce these features, as shown in Fig.3 (f),(g).

The -orbital picture also works well for the high energy states in the Lieb lattice. It is not surprising because a square lattice can be changed into a Lieb lattice by removing one of the four sites. The illustration of the -orbitals in Lieb lattice is given in Fig. 3 (e). Similar to the square lattice, we expect that the -band electrons will mainly distribute in between the sites, due to the formed -orbital bond, as was found in the experiment. The calculated LDOS confirms our expectation. We can see that, the experimental observation, the first-principles calculation [Fig. 3 (d)], and the muffin-tin model [Fig. 3 (h)] all coincide well with the -orbital picture. The main distinction from the square lattice is that the -orbitals at the A, B, C sites are now inequivalent, due to their different geometric environment. They thus have different on-site energies, which are denoted with different colors in Fig. 3 (e).

We plot the band structure of the Lieb lattice calculated from the muffin-tin potential model with plane wave basis in Fig. 4. We first consider a normal situation [Fig. 4 (a)], where nm, eV. In this case, we see that the lowest three bands (green solid lines) are -bands of Lieb lattice. Considering the -orbitals, there should be six -bands above the three -bands. This is shown in Fig. 4 (a) as blue solid lines, separated from both the lower -bands and the upper -bands. These -bands can be qualitatively interpreted by a -orbital tight binding model on Lieb lattice [see in Fig. 3 (e)]. In Fig. 4 (b), we show the tight binding fitting of the band structure. At low energy, the tight binding model works quite well, but it can not describe the top most two -bands. Similar to the square lattice, we attribute this discrepancy to the influence of the upper -bands. The corresponding LDOS of the -bands are like Fig. 3 (d),(h), where the -orbital bonds are shown clearly.

We now consider the atomic limit of the Lieb lattice ( eV, nm, nm) to support our -orbital picture. In the atomic limit, the hopping between adjacent sites is greatly suppressed. The bands of the atomic limit are shown in in Fig. 4 (c). In this case, the -bands are far from the - and -bands. In Fig. 4 (d), we also use the -orbital tight binding model to fit these -bands. Now the agreement is much better. Since the next nearest neighbor hopping becomes very tiny here, some -orbitals on the B and C sites now form dangling bonds, which results in degenerate flat -bands. The LDOS in atomic limit are also given in Fig. 4. At low energy, the -orbitals are around the lattice sites, and form a Lieb lattice [Fig. 4 (e)]. Note that, because the hopping is tiny in the atomic limit, the artificial atoms are nearly isolated. Thus, the unique LDOS pattern of flat band we mentioned above can not be observed here. Continuously increasing the energy, the LDOS can sequentially show different -orbitals as illustrated in Fig. 4 (f), (g) and (h). As we mentioned above, this is because the energy of the -orbitals in Lieb lattice are different.

In summary, we theoretically demonstrate that the high energy states in the artificial electron square (Lieb) lattice, as observed in the recent STM experiment, are from the -bands of the artificial atom confined in the lattice. The orbital degree of freedom is still a valid concept in this artificial electron lattice system. Our results suggest that, the electron lattice realized in Ref. Slot et al., 2017 may be the first artificial -orbital fermionic system in the solid state. Compared with other artificial quantum systems, this kind of electron lattice on metal surface is easy to manipulate, and the electron states can be directly detected. Thus, we believe that it is an ideal solid state platform to study orbital physics. Finally, we comment that, the same physics applies to artificial antidot lattice on 2DEG in semiconductor heterostructuresPark and Louie (2009); Polini et al. (2013); Gibertini et al. (2009).

J.H.G. and J.T.L. are supported by the National Natural Science Foundation of China (Grants No. 11534001, 61371015,11274129).


  1. thanks: These two authors contributed equally
  2. thanks: These two authors contributed equally


  1. Y. Tokura and N. Nagaosa, Science 288, 462 (2000).
  2. T. Hotta, Rep. Prog. Phys. 69, 2061 (2006).
  3. X. Li and W. V. Liu, Rep. Prog. Phys. 79, 116401 (2016).
  4. T. Müller, S. Fölling, A. Widera,  and I. Bloch, Phys. Rev. Lett. 99, 200405 (2007).
  5. G. Wirth, M. Ölschläger,  and A. Hemmerich, Nat. Phys. 7, 147 (2011).
  6. C. Wu, Mod. Phys. Lett. B 23, 1 (2009).
  7. K. K. Gomes, W. Mar, W. Ko, F. Guinea,  and H. C. Manoharan, Nature 483, 306 (2012).
  8. S. Wang, L. Z. Tan, W. Wang, S. G. Louie,  and N. Lin, Phys. Rev. Lett. 113, 196803 (2014).
  9. W.-X. Qiu, S. Li, J.-H. Gao, Y. Zhou,  and F.-C. Zhang, Phys. Rev. B 94, 241409 (2016).
  10. M. R. Slot, T. S. Gardenier, P. H. Jacobse, G. C. P. van Miert, S. N. Kempkes, S. J. M. Zevenhuizen, C. M. Smith, D. Vanmaekelbergh,  and I. Swart, Nat. Phys. 13, 672 (2017).
  11. R. Drost, T. Ojanen, A. Harju,  and P. Liljeroth, Nat. Phys. 13, 668 (2017).
  12. S. Li, W.-X. Qiu,  and J.-H. Gao, Nanoscale 8, 12747 (2016).
  13. M. Ropo, S. Paavilainen, J. Akola,  and E. Räsänen, Phys. Rev. B 90, 241401 (2014).
  14. G. Kresse and J. Furthmüller, Phys. Rev. B 54, 11169 (1996).
  15. J. P. Perdew, K. Burke,  and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
  16. J. Tersoff and D. R. Hamann, Phys. Rev. Lett. 50, 1998 (1983).
  17. R. A. Vicencio, C. Cantillano, L. Morales-Inostroza, B. Real, C. Mejía-Cortés, S. Weimann, A. Szameit,  and M. I. Molina, Phys. Rev. Lett. 114, 245503 (2015).
  18. S. Mukherjee, A. Spracklen, D. Choudhury, N. Goldman, P. Öhberg, E. Andersson,  and R. R. Thomson, Phys. Rev. Lett. 114, 245504 (2015).
  19. C.-H. Park and S. G. Louie, Nano Lett. 9, 1793 (2009).
  20. M. Polini, F. Guinea, M. Lewenstein, H. C. Manoharan,  and V. Pellegrini, Nat. Nanotechnol. 8, 625 (2013).
  21. M. Gibertini, A. Singha, V. Pellegrini, M. Polini, G. Vignale, A. Pinczuk, L. N. Pfeiffer,  and K. W. West, Phys. Rev. B 79, 241406 (2009).
This is a comment super asjknd jkasnjk adsnkj
The feedback cannot be empty
Comments 0
The feedback cannot be empty
Add comment

You’re adding your first comment!
How to quickly get a good reply:
  • Offer a constructive comment on the author work.
  • Add helpful links to code implementation or project page.