Flat Band and Hole-induced Ferromagnetism in a Novel Carbon Monolayer

TL;DR

This study predicts that hole doping in a novel carbon monolayer called cyclicgraphdiyne induces flat bands and ferromagnetism, offering new insights into flat band formation and magnetic manipulation in carbon materials.

Contribution

It introduces cyclicgraphdiyne as a new carbon monolayer with flat bands and demonstrates hole doping induces ferromagnetism, providing an alternative approach to control magnetism in carbon systems.

Findings

Hole doping induces ferromagnetism in cyclicgraphdiyne.

Flat bands are tunable by next-nearest-neighbor hopping.

Cyclicgraphdiyne exhibits high conductivity at certain doping levels.

Abstract

In recent experiments, superconductivity and correlated insulating states were observed in twisted bilayer graphene (TBG) with small magic angles, which highlights the importance of the flat bands near Fermi energy. However, the moir\'e pattern of TBG consists of more than ten thousand carbon atoms that is not easy to handle with conventional methods. By density functional theory calculations, we obtain a flat band at E$_F$ in a novel carbon monolayer coined as cyclicgraphdiyne with the unit cell of eighteen atoms. By doping holes into cyclicgraphdiyne to make the flat band partially occupied, we find that cyclicgraphdiyne with 1/8, 1/4, 3/8 and 1/2 hole doping concentration shows ferromagnetism (half-metal) while the case without doping is nonmagnetic, indicating a hole-induced nonmagnetic-ferromagnetic transition. The calculated conductivity of cyclicgraphdiyne with 1/8, 1/4 and 3/8 hole doping concentration is much higher than that without doping or with 1/2 hole doping. These results make cyclicgraphdiyne really attractive. By studying several carbon monolayers, we find that a perfect flat band may occur in the lattices with both separated or corner-connected triangular motifs with only including nearest-neighboring hopping of electrons, and the dispersion of flat band can be tuned by next-nearest-neighboring hopping. Our results shed insightful light on the formation of flat band in TBG. The present study also poses an alternative way to manipulate magnetism through doping flat band in carbon materials.

Figures

Equations

$$H(k_x,k_y)=\left(\begin{array}{cccc}0& B& C& 0\\ B^{\ast }& A& 0& 0\\ C^{\ast }& 0& 0& B^{\ast }\\ 0& 0& B& A^{\ast }\end{array}\right),$$

$$A=\left(\begin{array}{ccc}0& a& b\\ a^{\ast }& 0& c\\ b^{\ast }& c\ast & 0\end{array}\right),$$

$$H_{NNN}(k_x,k_y)=t_2\left(\begin{array}{cccc}0& D& 0& 0\\ D^{\ast }& 0& 0& 0\\ 0& 0& 0& D^{\ast }\\ 0& 0& D& 0\end{array}\right),$$

$$D=\left(\begin{array}{ccc}0& e& f\\ m& 0& n\\ p& q& 0\end{array}\right),$$

Full text

1 Introduction

A twisted bilayer graphene (TBG) with small angles was predicted to have moiré pattern with a very large unit cell of more than ten thousand carbon atoms. At the magic angles, the two layers become strongly coupled, and a flat band appears 1. A tight-binding model 2 and a continuum model of TBG 3 were thus studied, and the confined states were also proposed in TBG 4. Recently, the superconductivity and correlated insulating states at half-filling were observed in TBG, which highlight the importance of the flat bands near Fermi energy 5, 6. The strong coupling phases in TBG (e.g. 7) were explored. It was shown that the degeneracy of the low-energy flat bands can be lifted by the electric fields 8. The topological properties of electronic band structure 9, 10, ferromagnetic Mott state 11 and metal-insulator transition 12, 13, extended Hubbard model, nematic superconductivity by density wave fluctuations 14, 15, 16 in TBG were addressed.

The study of TBG has inspired extensive studies in other twisted bilayer materials. Moiré excitons, moiré trapped valley excitons, and resonantly hybridized excitons in van der Waals heterostructures were observed in experiments 17, 18, 19, 20. Van der Waals crystals with discretized Eshelby twist in two-dimensional (2D) materials 21, and van der Waals contacts between three-dimensional metals and 2D semiconductors 22 were investigated. Chiral twisted van der Waals nanowires 23, and helical van der Waals crystals 24 were also studied.

In this paper, by means of the first-principles calculations, we find that the flat band near the EF could appear in a carbon monolayer with unit cell of eighteen atoms, which is coined cyclicgraphdiyne. By doping holes into cyclicgraphdiyne to make the flat band partially occupied, cyclicgraphdiyne with 1/8, 1/4, 1/8 and 1/2 hole doping concentration becomes ferromagnetic, while the case without doping holes remains nonmagnetic, showing a hole doping induced nonmagnetic-ferromagnetic transition in carbon materials. After carefully studying the structural stability, we disclose that cyclicgraphdiyne is viable. By investigating several carbon monolayers, we notice that a perfect flat band can occur in the lattices with both separated or corner-connected triangular motifs and nearest-neighboring hopping of electrons, and the dispersion of flat band can be tuned by next nearest-neighboring hopping. Our results shed light on the formation reasoning of flat band in TBG. The present study also provides a simple 2D platform of carbon atoms to probe the flat bands and magnetism in 2D systems.

2 Results

2.1 Structure and Stability

Cyclicgraphyne and cyclicgraphdiyne can be made by inserting an acetylene linkage and two acetylene linkages between two neighboring carbon triangles 25, respectively, as shown in Fig. 1. The crystal structure has the space group $P6/mmm$ (No.191). Our calculations show that the equilibrium lattice constants of cyclicgraphyne and cyclicgraphdiyne are $a_{0}$ = 9.6452 Å and $a_{0}$ =14.0920 Å, respectively. To unveil the stabilities of cyclicgraphyne and cyclicgraphdiyne, the phonon spectra were calculated, as shown in Figs. S1(a) and (b) (supplemental material). It can be seen that there is no negative frequency phonon in the whole Brillouin zone, indicating that cyclicgraphyne and cyclicgraphdiyne are kinetically stable. To further examine the thermal stability, we performed molecular dynamics simulations by considering $3\times 3\times 1$ supercells of cyclicgraphyne and cyclicgraphdiyne with 108 and 162 carbon atoms, respectively. As displayed in Fig. S1(c), after being heated at 500 K, 700 K, and 1000 K for 6 ps with a time step of 3 fs, no structural changes occur, indicating that cyclicgraphyne and cyclicgraphdiyne are also dynamically stable. As shown in Fig. S1(d), we present the free energy of cyclicgraphyne and cyclicgraphdiyne, compared with three experimentally available carbon allotropes: T-carbon 26, graphene 27, and graphdiyne 28. The results suggest that they may be feasible for experimental synthesis.

2.2 Flat Band in Cyclicgraphyne and Cyclicgraphdiyne

By DFT calculations, the band structure of cyclicgraphyne near Fermi energy $E_{F}$ are shown in Fig. 2(a). At 145 meV below Fermi energy $E_{F}$, a flat band with bandwidth of about 2 meV is observed [Fig. 2(b)]. The flat band is mainly coming from $p_{z}$ electrons, as demonstrated in the partial density of state (PDOS) in Fig. 2(a). Our result of flat band in cyclicgraphyne is consistent with previous studies 29, 30. As shown in Fig. 2(c), a very similar electronic band structure is obtained in cyclicgraphdiyne, where a perfect flat band appears nearly at $E_{F}$ (see also Fig. 2(d)).

2.3 Origin of Flat Band in Cyclicgraphdiyne

The total electron charge density of cyclicgraphdiyne is illustrated in Fig. 3(a). It is instructive to note that the peak of total electron charge density occurs at acetylene linkages, while the valley appears at single bonds and triangles. If we fix the energy at the flat band near EF, the obtained electron charge density of the flat band is shown in Fig. 3(b). It has a characteristic of $p_{z}$ orbitals, and is not distributed on the carbon atoms at acetylene linkages in cyclicgraphdiyne. By comparing Figs. 3(a) and (b), we know that the flat band should mainly come from the valley part of the total electron charge density as shown in Fig. 3(c).

2.4 Hole-induced ferromagnetism in Cyclicgraphdiyne

The interplay between flat band and Coulomb interaction may cause the flat band to be magnetic. To investigate the flat band magnetism, we doped a hole in the unit cell of cyclicgraphdiyne to make the flat band near EF half filled. Figure 4(a) shows the band structure of the doped ferromagnetic ground state with an average spin moment about 0.03 $\mu_{B}$ per atom. The spin splitting of the flat band is about 0.28 eV. To study the change of conductivity with different hole doping in cyclicgraphdiyne, we calculated the conductivity divided by the relaxation time [Fig. 4(b)] in a 2$\times$2$\times$1 supercell of cyclicgraphdiyne with 0, 1, 2, 3 and 4 hole dopants, respectively, i.e. the hole doping concentration corresponds to the hole filling of 0, 1/8, 1/4, 3/8 and 1/2, respectively. From Fig. 4(b), we note that the conductivity in half and full fillings of the flat band in cyclicgraphdiyne is much smaller than that in the case of 1/8, 1/4 and 3/8 hole doping. The behavior of conductivity in Fig. 4(b) can be understood in terms of the density of states (DOS). The DOS near Fermi energy with hole filling of 1/8, 1/4, 3/8 and 1/2 is presented in Fig. S2. It is shown that the DOS at Fermi energy becomes finite and appears only for single species electrons with hole hoping, and becomes to nearly zero at half filling. This result reveals that the hole-doped cyclicgraphdiyne is a half-metal, which would be useful for spintronics. The resemblance of conductivity as a function of hole doping between cyclicraphdiyne in Fig. 4(b) and TBG observed in the experiment 5, 6 makes cyclicgraphdiyne really attractive.

2.5 Flat Band and Tight-Binding Model

To understand the origin of the flat band in cyclicgraphdiyne, we construct a specified tight-binding model based on the $p_{z}$ orbitals of cyclicgraphdiyne. A simplified effective model based on $p_{z}$ orbitals of carbon atoms for cyclicgraphdiyne can be written by a 12-bands Hamiltonian, where we only consider the carbon atoms on the triangles and their nearest neighboring carbon atoms, and the interaction between the carbon atoms on the acetylene linkages is summarized to the effective hopping $t^{{}^{\prime\prime}}$ as shown in Fig. 5. So we have

[TABLE]

where

[TABLE]

$a=te^{i\cdot(k_{x}/2+\sqrt{3}k_{y}/2)\cdot l}$, $b=te^{i\cdot k_{x}\cdot l}$, $c=te^{i\cdot(k_{x}/2-\sqrt{3}k_{y}/2)\cdot l}$, $B=Diag[$ $t^{{}^{\prime}}e^{i\cdot 0.5473\cdot(3k_{x}/2+\sqrt{3}k_{y}/2)\cdot l},t^{{}^{\prime}}e^{-i\cdot 0.5473\cdot\sqrt{3}k_{y}\cdot l},$ $t^{{}^{\prime}}e^{-i\cdot 0.5473\cdot(3k_{x}/2-\sqrt{3}k_{y}/2)\cdot l}]$, $C=Diag[t^{{}^{\prime\prime}}e^{i\cdot 1.5476\cdot(3k_{x}/2+\sqrt{3}k_{y}/2)\cdot l},t^{{}^{\prime\prime}}e^{-i\cdot 1.5476\cdot\sqrt{3}k_{y}\cdot l},t^{{}^{\prime\prime}}e^{-i\cdot 1.5476\cdot(3k_{x}/2-\sqrt{3}k_{y}/2)\cdot l}]$, $l$=1.425 Å(the bond length of two carbon atoms at triangles), $t$, $t^{{}^{\prime}}$ and $t^{{}^{\prime\prime}}$ are the nearest-neighbouring (NN) hopping parameters as shown in Fig. 5(a). By DFT calculations 31, these NN hopping parameters were obtained, and the tight-binding (TB) flat band fits well with the DFT bands as shown in Fig. 5(b). If the next-nearest-neighbouring (NNN) hopping (such as $t_{2}$) is included, which is described by

[TABLE]

where

[TABLE]

with $e=e^{i[(0.5473\cdot 3+1)k_{x}/2+(0.5473+1)\sqrt{3}k_{y}/2]\cdot l}$, $f=e^{i[(0.5473\cdot 3+2)k_{x}/2+0.5473\sqrt{3}k_{y}/2)]\cdot l}$, $m=e^{i[-k_{x}/2-(0.5473\cdot 2+1)\sqrt{3}k_{y}/2]\cdot l}$, $n=e^{i[k_{x}/2-(0.5473\cdot 2+1)\sqrt{3}k_{y}/2]\cdot l}$, $p=e^{i[-(0.5473\cdot 3+2)k_{x}/2+0.5473\sqrt{3}k_{y}/2)]\cdot l}$, $q=e^{i[-(0.5473\cdot 3+1)k_{x}/2+(0.5473+1)\sqrt{3}k_{y}/2]\cdot l}$, the obtained perfect flat band can be tuned. A small NNN hopping $t_{2}=-0.05$ can induce a tiny bandwidth as shown in Fig. 5(c).

The above results are usually implemented in a kagome lattice, where a perfect flat band appears with only NN hopping and can be tuned out by NNN hopping. By comparing the results of cyclicgraphdiyne in Fig. 5 and those of kagome lattice in Fig. S3, and studying several materials with triangular and hexagonal structures, we can conclude that a perfect flat band can occur in the lattices with both separated or corner-connected triangular structures with only NN hopping, and can be tuned out by the NNN hopping.

3 Discussions

3.1 Flat Band and Hidden Valley Kagome Lattice in TBG

For TBG with small magic angles, the unit cell can be constructed from moiré pattern, which may be regarded as a huge triangular lattice 5, where the lattice points are comprised of the peaks of local density of states (LDOS) of TBG, while the valleys of LDOS, i.e., the edge center of triangles, form a hidden kagome lattice as shown in Fig. S5. It is known that the kagome lattice with only NN hopping can induce the flat band 32, 33, as we also addressed in the Supplemental Materials. In TBG, the hidden kagome lattice obtained from moiré pattern is quite large, so the NNN hopping should be much smaller than the NN hopping, and a nearly flat band can be obtained from such a hidden kagome moire lattice. Analog to the flat band in cyclicgraphdiyne [Fig. 2(d)], where the flat band mainly comes from the valleys of the total electron charge density, we believe that the flat band in TBG may be attributed to the hidden kagome lattice in moiré pattern.

4 Conclusions

In summary, by density functional theory calculations, we propose a novel carbon monolayer with unit cell of eighteen atoms, coined as cyclicgraphdiyne. A flat band appears at EF in cyclicgraphdiyne. By doping holes to make the flat band in cyclicgraphdiyne partially occupied, we find that cyclicgraphdiyne with 1/8, 1/4, 3/8 and 1/2 hole doping concentration becomes ferromagnetic, and the conductivity of cyclicgraphdiyne with 1/8, 1/4 and 3/8 hole doping concentration is much higher than the case that the flat band in cyclicgraphdiyne is full or half filled. These results make the cyclicgraphdiyne proposed here very attractive. By exploring different carbon monolayers, we show that a perfect flat band could appear in the lattices with both separated or corner-connected triangular structures and only NN hopping, while the flat band can be tuned out by NNN hopping. Owing to the resemblance of the electronic structures between cyclicgraphdiyne and TBG, our results may also shed useful light on the formation reasoning of the flat band in TBG. We believe that the flat band in TBG may be attributed to a hidden valley kagome lattice in moiré pattern. This present study poses a simple carbon monolayer and gives an alternative reasoning to the formation of the flat bands and manipulating magnetism in 2D carbon materials.

5 Computational Method

Our first-principles calculations were based on the density-functional theory (DFT) as implemented in the Vienna ab initio simulation package (VASP) 34, using the projector augmented wave method 35. The generalized gradient approximation with Perdew-Burke-Ernzerhof 36 realization was adopted for the exchange-correlation functional. The plane-wave cutoff energy was set to 550 eV. The Monkhorst-Pack $k$-point mesh 37 of size $11\times 11\times 1$ was used for the BZ sampling. The structure relaxation considering both the atomic positions and lattice vectors was performed by the conjugate gradient (CG) scheme until the maximum force on each atom was less than 0.0001 eV/Å, and the total energy was converged to $10^{-8}$ eV with Gaussian smearing method. To avoid unnecessary interactions between the monolayer and its periodic images, the vacuum layer is set to 15 Å. The crystal structures and charge density were plotted with the software VESTA 38, and the images in the figure legend were created with the software SciDAVis 39. The unit of hole-concentration is dimensionless. The 100$\%$ doping means 2 holes doped in area of 172 Å (the area of unit cell of cyclicgraphdiyne). We used Hubbard $U$ (with values 1 and 4 eV) correction in DFT calculations and the electronic band structure was found unaltered. The conductivity is calculated with the package BoltzTrap 40 based on a smoothed Fourier interpolation of the bands at temperature of 10 K.

6 Acknowledgements

The authors acknowledge Zheng Zhu for valuable discussions. B.G. is supported by the National Natural Science Foundation of China (Grant No. Y81Z01A1A9), the Chinese Academy of Sciences (Grant No. Y929013EA2), the University of Chinese Academy of Sciences (Grant No. 110200M208), and the Beijing Natural Science Foundation (Grant No. Z190011). J.Y.Y. and G.S. are supported in part by the National Key R&D Program of China (Grant No. 2018YFA0305800), the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB28000000), the National Natural Science Foundation of China (Grant No.11834014), and Beijing Municipal Science and Technology Commission (Grant No. Z118100004218001).