Realization of Lieb lattice in covalent-organic frameworks with tunable topology and magnetism
A Lieb-like lattice in a covalent-organic framework and its Stoner ferromagnetism
Wei Jiang, Huaqing Huang & Feng Liu
Coordination assembly of 2D ordered organic metal chalcogenides with widely tunable electronic band gaps
Yanzhou Li, Xiaoming Jiang, … Gang Xu
Atomically precise single-crystal structures of electrically conducting 2D metal–organic frameworks
Jin-Hu Dou, Maxx Q. Arguilla, … Mircea Dincă
From an antiferromagnetic insulator to a strongly correlated metal in square-lattice MCl2(pyrazine)2 coordination solids
Panagiota Perlepe, Itziar Oyarzabal, … Rodolphe Clérac
Implications of bond disorder in a S=1 kagome lattice
Jamie L. Manson, Jamie Brambleby, … Bruce Noll
Controlling pairing of π-conjugated electrons in 2D covalent organic radical frameworks via in-plane strain
Isaac Alcón, Raúl Santiago, … Stefan T. Bromley
A coronene-based semiconducting two-dimensional metal-organic framework with ferromagnetic behavior
Renhao Dong, Zhitao Zhang, … Xinliang Feng
Making large single crystals of 2D MOFs
Renhao Dong & Xinliang Feng
Magnetic ordering through itinerant ferromagnetism in a metal–organic framework
Jesse G. Park, Brianna A. Collins, … Jeffrey R. Long
Bin Cui ORCID: orcid.org/0000-0002-7240-005X1,
Xingwen Zheng1,
Jianfeng Wang2,
Desheng Liu1,
Shijie Xie1 &
Bing Huang ORCID: orcid.org/0000-0001-6735-46372
Nature Communications volume 11, Article number: 66 (2020) Cite this article
Density functional theory
Ferromagnetism
Molecular electronics
An Author Correction to this article was published on 29 January 2020
This article has been updated
Lieb lattice has been predicted to host various exotic electronic properties due to its unusual Dirac-flat band structure. However, the realization of a Lieb lattice in a real material is still unachievable. Based on tight-binding modeling, we find that the lattice distortion can significantly determine the electronic and topological properties of a Lieb lattice. Importantly, based on first-principles calculations, we predict that the two existing covalent organic frameworks (COFs), i.e., sp2C-COF and sp2N-COF, are actually the first two material realizations of organic-ligand-based Lieb lattice. Interestingly, the sp2C-COF can experience the phase transitions from a paramagnetic state to a ferromagnetic one and then to a Néel antiferromagnetic one, as the carrier doping concentration increases. Our findings not only confirm the first material realization of Lieb lattice in COFs, but also offer a possible way to achieve tunable topology and magnetism in organic lattices.
The electronic properties of a crystal are determined by its crystalline lattice symmetry. Lieb lattice, a two-dimensional (2D) edge-depleted square lattice (Fig. 1a), can be regarded as the reduced 3D perovskite lattice1. As one of the most important frustrated lattices, the ideal Lieb lattice can host exotic electronic structures, which is featured by the Dirac cone intersected by a flat band (Dirac-flat bands)2, as shown in Fig. 1c. Interestingly, various physical phenomena, e.g., topological Mott insulator3, superconductivity4, ferromagnetism5,6, and fractional quantum Hall (FQH)7,8,9,10 effects, have been predicted to exist in the ideal Lieb lattice systems. However, until now only a few Lieb lattices have been realized in artificial lattice systems, e.g., molecular patterning lattices on metal substrates11,12, photonic and cold-atom lattices13,14,15,16, rather than a real material system, which significantly prevent the realization of these unusual physical properties of Lieb lattice for practice applications.
Fig. 1: Three-band Lieb lattice models.
a Geometrical structure of a Lieb lattice, which contains one corner site 1 and two edge sites 2 and 3 in the unitcell (marked by the blue dashed-lines). Strength of lattice (geometrical) distortion is indicated by the change of structural angle θ away from 90°. See text for the definitions of other parameters. b Brillouin zone of a Lieb lattice. It notes that M and M1, equal in an ideal Lieb lattice, are inequivalent in a distorted Lieb lattice (θ≠90° or t2≠t3). Projected band structures of an ideal Lieb lattice c without and d with NNN hopping. e–g Projected band structures of distorted Lieb lattices under three representative situations. In (c–g), the colors of the energy bands represent the wavefunction components from different lattice sites shown in (a). Curve width indicates component weight.
Although it is found that the Lieb lattice is extremely difficult to be realized in an inorganic material, it is still possible to realize such a lattice in an organic one, as the structural flexibility allows an organic system exhibiting diverse crystalline symmetries17,18,19,20,21,22,23,24,25. Especially, it is noticed that the covalent-organic frameworks (COFs) can host various 2D lattices18,21,24,26, including square lattice21,26, which makes the discovery of a Lieb lattice in a COF possible. Usually, the molecular orbitals (MOs) play an essential role in determining the electronic properties of a ligand-based organic semiconductor, e.g., a COF. Since the MOs in an organic system can be more easily modulated than the atomic orbitals (AOs) in an inorganic system by the lattice distortions, it is expected that the tunable electronic and magnetic states could be achieved in an organic Lieb lattice.
In this article, using tight-binding (TB) modeling and first-principles density-functional theory (DFT) calculations, we find that the electronic and topological properties of a Lieb lattice can be effectively modulated by its lattice distortions. Interestingly, we discover that sp2C-COF and sp2N-COF, which are synthesized in recent experiments27, are the first two material realizations of organic-ligand-based Lieb lattices. The electronic structures of these two COFs around the band edges can be well characterized by the MOs-based Dirac-flat-band models. Furthermore, we find that the lattice distortions can also dramatically affect the bandwidth of the Dirac-flat bands, which in turn determines its electronic instability against spontaneous spin-polarization during carrier doping. Remarkably, it is found that both ferromagnetic (FM) and Néel antiferromagnetic (AFM) phases can be realized in these distorted Lieb lattices under certain hole doping concentrations (nh). Our discovery not only confirms the first material realization of Lieb lattice in COFs, but also suggests an interesting routine to realize tunable topological and magnetic states in d-(f-) orbital-free organic lattices.
Three-band Lieb lattice models
We begin with the three-band tight-binding Hamiltonian of Lieb lattice [without spin-orbit coupling (SOC) effects], including one corner and two edge sites in a 2D line-centered square lattice unitcell:
$$H_0 = \mathop {\sum}\nolimits_i {E_ic_i^\dagger c_i} + \mathop {\sum}\nolimits_{\left\langle {i,j} \right\rangle } {\left( {t_1c_i^\dagger c_j + h.c.} \right)} \\ + \mathop {\sum}\nolimits_{\left\langle {\left\langle {i,j} \right\rangle } \right\rangle } {\left( {t_{2,3}c_i^\dagger c_j + h.c.} \right)}$$
where ci (\(c_i^\dagger\)) annihilates (creates) an electron with the energy of Ei on the ith site (i = 1, 2, or 3) in the unitcell, and t1 and t2,3 are nearest-neighbor (NN, \(\left\langle {i,j} \right\rangle\)) and next nearest-neighbor (NNN, \(\left\langle {\left\langle {i,j} \right\rangle } \right\rangle\)) hopping integrals, respectively, as shown in Fig. 1a. The details of TB modeling can be found in the Supplementary Notes 1–3.
For an ideal Lieb lattice, the on-site energy of the three sites are identical, i.e., E1 = dE/2 = 0 and E2 = E3 = −dE/2 = 0. Then, we can calculate the (three-band) electronic structure of an ideal Lieb lattice. When the NNN hopping is not included, as shown in Fig. 1c, the calculated band structure is featured by the Dirac cone intersected by a flat band. The wavefunction of flat band is mostly localized at the edge sites, while the wavefunctions of Dirac bands are equally contributed by both corner and edge sites. When NNN hopping is included, e.g., t2 = t3 = 0.2t1, the features of the Dirac-flat bands still maintain besides that the flat band becomes more dispersive around the Γ point in the Brillouin zone (BZ), as shown in Fig. 1d.
Since the ideal Lieb lattice is very rare in solids3,4,5,6,7,8,9,10, it is curious for us to further understand the electronic structures of distorted Lieb lattices. Firstly, we consider the situation of θ ≠ 90°, e.g., θ < 90°, corresponding to the geometrical lattice distortions. In this case, the fourfold rotation symmetry of the system is no longer valid, then t2 ≠ t3. As shown in Fig. 1e, the Dirac-flat bands decompose into two sets of Dirac cones that deviate away from M (M1) to Γ point along the Γ-M (Γ-M1) in the BZ. Secondly, we consider the situation of dE ≠ 0, which corresponding to different atom (ligand) occupations at corner and edge sites in an inorganic (organic) lattice. In this case, the band degeneracy of Dirac-flat bands is broken and a bandgap Δ1 (hereafter, Δ1 is defined as the bandgap between the top Dirac and middle flat bands) is induced, as shown in Fig. 1f. Meanwhile, the flat band and the bottom Dirac band still keep in touch with each other. Interestingly, it is found the wavefunction components of these Dirac-flat bands are dramatically different from the ideal Lieb lattice, i.e., the top Dirac band is contributed by the corner sites while the flat and bottom Dirac bands are contributed by the edge sites. Finally, we consider the distortions with both θ ≠ 90° and dE ≠ 0, as shown in Fig. 1g, which correspond to the realistic situations in many materials, including the COFs focused in our study. In this case, the changes of Dirac-flat bands could be considered as a combined effect of θ ≠ 90° (Fig. 1e) and dE ≠ 0 (Fig. 1f). Interestingly, only one Dirac cone can survive around M1 point along the Γ − M1 line.
Distorted Lieb lattice model with SOC effects
We further consider the SOC effects on the distorted Lieb lattice model, with the total Hamiltonian H = H0 + HSO. HSO is considered as an imaginary hopping between the NNN sites, similar to Kane and Mele's SOC term9,28:
$$H_{{\mathrm{SO}}} = - {\mathrm{i}}{\uplambda} \sum_{\left\langle {\left\langle {i,j} \right\rangle } \right\rangle } \left( {{\mathbf{d}}_{ik} \times {\mathbf{d}}_{kj}} \right) \cdot {\mathbf{s}}_{\alpha \beta }^zc_{i\alpha }^\dagger c_{j\beta }$$
where λ is the SOC constant and dik (dkj) denotes the unit vector from edge (corner) site i (k) to corner (edge) site k (j) (see Supplementary Fig. 1). s is the Pauli matrix representing the electron spin. Generally, we find that SOC effects can lift the degeneracy of the Dirac cone, and the topological properties of these Dirac-flat bands strongly depend on the value of dE (the values of other parameters are the same as those in Fig. 1g).
For dE = 0 (Fig. 2a), the bandgaps of \({\mathrm{\Delta }}_1\)and \({\mathrm{\Delta }}_2\) (hereafter, \({\mathrm{\Delta }}_2\) is defined as bandgap between middle flat and bottom Dirac bands) can be induced around M and M1 points. It is noted that \({\mathrm{\Delta }}_2\) can only be induced by SOC effects. For the case of \(dE \ne 0\), we find that there is a competing mechanism between dE and SOC for the changes of \({\mathrm{\Delta }}_1\). When 0 < dE < dEc (\(dE_{\mathrm{c}} = 2\sqrt {\left( {t_2 - t_3} \right)^2 \, + \, 4{\uplambda}^2}\), see Supplementary Notes 3–5), \({\mathrm{\Delta }}_1\) decreases as dE increases. When dE = dEc, \({\mathrm{\Delta }}_1 = 0\), as shown in Fig. 2b, two Dirac cones can form around M and M1 points. When dE > dEc (Fig. 2c), \({\mathrm{\Delta }}_1\) opens again and increases as dE increases, which is opposite to the case of dE < dEc. Meanwhile, when dE changes through dEc, the wavefunction compositions of top Dirac band and middle flat band are switched, as shown in Fig. 2a, c.
Fig. 2: Distorted Lieb lattice models with SOC effects.
a–c Calculated band structures of distorted Lieb lattices with SOC effects (λ = 0.1t1). TB parameters of distorted Lieb lattice are adopted from Fig. 1g but with variable dE: (a) dE = 0 < dEc, (b) dE = dEc, and (c) dE = 2t1 > dEc. \(dE_{\it{c}} = 2\sqrt {\left( {t_2 - t_3} \right)^2 + \, 4{\uplambda}^2}\). Calculated spin Chern number for each band (\(C_{\mathrm{n}}^{\mathrm{s}}\)) is also labeled. d–f Spin-resolved edge states for each band structure in (a–c), respectively.
The topological properties of each band can be characterized by its spin Chern number, due to the spin degeneracy and conservation of Sz. The spin Chern number for the nth band is defined as \(C_n^{\mathrm{s}} = \left( {C_{n \uparrow } - C_{n \downarrow }} \right)/2\), and \(C_{n{\upsigma}}\) is Chern number for the spin-\({\upsigma}\) (\({\upsigma} = \, \uparrow , \downarrow\)) of the nth band, which can be calculated from the integral of the Berry curvature over the whole BZ29,30 (see "Method" section and supplementary Fig. 2). As shown in Fig. 2, when dE = 0, the calculated \(C_1^{\mathrm{s}} = - 1\) and \(C_3^{\mathrm{s}} = + 1\) for the top and bottom Dirac bands, respectively, and the calculated \(C_2^{\mathrm{s}} = 0\) for middle flat band. Interestingly, when dE > dEc, the top Dirac and middle flat bands can switch their topologies, as shown in Fig. 2c, i.e., \(C_1^{\mathrm{s}} = 0\) (\(C_2^{\mathrm{s}} = - 1\)) for the top Dirac (middle flat) band.
The nontrivial topology of these Dirac-flat bands in Fig. 2a–c can be reflected by edge states calculations, as shown in the according Fig. 2d–f. It can be seen that the helical edge states with opposite spin channels connect the two bands with nonzero \(C_n^{\mathrm{s}}\), indicating the quantum spin Hall effect.
It is emphasized that in a realistic material, only the total spin Chern number Cs, the sum of \(C_n^{\mathrm{s}}\) for all the occupied bands, is associated with the observable quantum conductance of the system. Therefore, the topological properties of a realistic Lieb lattice material sensitively depend on the position of Fermi level, and the charge doping may be needed to achieve a nontrivial state.
Realization of distorted Lieb lattices in COFs
Taking advantage of the structural flexibility of organic materials, numerous COFs with various crystalline symmetries have been synthesized in recent years. Remarkably, after extensive structural analysis, we notice that two COFs with sp2 hybridized backbones fabricated in the very recent experiments, named as sp2C-COF and sp2N-COF, can realize slightly distorted square Lieb lattices27 (θ < 90°), as shown in Fig. 3a, b, respectively. Meanwhile, the ligands at the corner and edge sites are different in monolayer sp2C-COF and sp2N-COF, corresponding to the case of dE ≠ 0. In both COFs, a pyrene surrounded by four phenyl rings is located at the corner site, named as 1,3,6,8-tetrakis(2-phenyl)pyrene (TPPy), while the ligands at the edge sites in these two COFs are slightly different, i.e., 1,4-bis(2-crylonitrile)benzene (BCNB) in sp2C-COF and 1,4-bis(2-aldimine)benzene (BAIB) in sp2N-COF. Meanwhile, crylonitriles (CN) and aldimine (AI) bond to the bis- positions on the phenyl ring of edge sites in sp2C-COF and sp2N-COF, respectively. The enlarged view of connection parts between corner and edge sites are shown in Supplementary Fig. 3.
Fig. 3: Realization of distorted Lieb lattices in COFs.
Top views of (a) sp2C-COF and (b) sp2N-COF in a 2 × 2 supercell. The corner (TPPy) and edge (BCNB or BAIB) ligands are marked by the green-dashed squares and blue-dashed ellipses, respectively. Insets: Enlarged structures of end groups on the edge sites, i.e., crylonitrile for BCNB and aldimine for BAIB. c and d Projected band structures of monolayer sp2C-COF and sp2N-COF, respectively, where the green and blue circles (the sizes represent the weight of MOs) indicate that the states are contributed by corner (TPPy) and edge (BCNB or BAIB) sites, respectively. Insets: Partial charge distributions of top three valence bands. Orange dashed-lines in (c) and (d): TB fittings of DFT band structures based on distorted Lieb lattice model for these top three valence bands. In (c) sp2C-COF: \(t_1 = 0.1\;{\mathrm{eV}},t_2 = 0.4t_1,\;t_3 = 0.45t_2,\;dE = 5.9t_1\); in (d) sp2C-COF: \(t_1 = 0.12\, {\mathrm{eV}},\, t_2 = 0.2t_1,\, t_3 = 0.45t_2,\, dE = 2.5t_1\).
We have calculated the electronic properties of monolayer sp2C-COF and sp2N-COF using DFT-PBE calculations, as shown in Fig. 3c, d, while the results for bulk sp2C-COF and sp2N-COF can be found in Supplementary Fig. 4. Generally, it is found that monolayer and its corresponding bulk COFs have similar electronic properties due to the relatively weak interlayer van der Waals (vdW) interactions. The bandgaps of monolayer sp2C-COF and sp2N-COF are 1.45 and 1.30 eV, respectively, which is smaller than the experimental values (e.g., 1.9 eV for sp2C-COF)27 due to the underestimation of bandgaps in DFT calculations.
Here, we mainly pay our attention to the band edges, especially for the top three valence bands, i.e., VB1–VB3 marked in Fig. 3c, d, in both COFs. In the band structure of sp2C-COF (Fig. 3c), VB1 separates from VB2 and VB3 in energy, opening an energy gap of \(\Delta _1\)~0.4 eV. Notably, VB2 and VB3 touch with each other forming a single Dirac cone near the M1 point on the \({\mathrm{\Gamma }} - {\mathrm{M}}1\) line. The wavefunction distributions of energy bands can be better understood by the ligand-based MOs than AOs in a COF. Generally, the band edges in a ligand lattice are formed by the HOMOs and LUMOs of the ligands, respectively31, and the interaction between the HOMOs and LUMOs are negligible. Therefore, the VB1–VB3 can be considered separately from the CB1–CB3 in a TB modeling based on the Lowdin perturbation theory. As shown in Fig. 3c, the charge densities of VB1–VB3 are projected into the ligands of sp2C-COF. Interestingly, it is found that the top VB1 is solely contributed by the HOMO of TPPy (corner ligands), while the less dispersive VB2 and the bottom VB3 are contributed by the HOMOs of both BCNBs (edge ligands), respectively. Interestingly, it is found that the band dispersions and wavefunction distributions of VB1–VB3 (Fig. 3c) indeed agree well with the three-band distorted Lieb lattice model under the situation of dE ≠ 0 and θ ≠ 90° (Fig. 1g). The band structure of sp2N-COF (Fig. 3d) is similar to that of sp2C-COF. The \(\Delta _1\) in sp2N-COF is ~0.2 eV, which is smaller than that in sp2C-COF. Meanwhile, it is also noted that VB1 in sp2C-COF is much less dispersive than that in sp2N-COF. The calculated SOC effects (\({\mathrm{\Delta }}_2\)~0.03 meV) are negligible in these two COFs, indicating that it could be challenging to observe the nontrivial topological properties (Fig. 2) in these two COFs.
Since VB1–VB3 in these two COFs are similar to that of three-band distorted Lieb lattice model, we can further understand the band dispersions of VB1–VB3 using TB band fitting of DFT calculations. Overall, we find that the differences of VB1–VB3 band dispersions in these two COFs are mainly determined by the difference of dE between the HOMOs of corner and edge sites in these two systems. It is found that dE determines the size of \(\Delta _1\), i.e., the larger the dE, the larger the \(\Delta _1\) will be (without SOC effects). Meanwhile, \(\Delta _1\) can in turn determine the band dispersion of VB1, i.e., the larger the \(\Delta _1\), the flatter the VB1 will be. In addition, the interfacial torsion angle between the corner and edge ligands in the sp2C-COF is larger than that in the sp2N-COF (see Supplementary Fig. 3b, e). The larger torsion angle can reduce the conjugation in a larger degree in the COF and lead to a smaller t1. Therefore, the larger dE and torsion angles in sp2C-COF than in sp2N-COF can give rise to a flatter VB1 in sp2C-COF, agreeing well with the DFT results. On the other hand, our understandings suggest a possible way to design flat bands in organic lattices, as the flat bands could have great importance for realizing many novel physical phenomena.
Similarly, the bottom three conduction bands (CB1–CB3), formed by the LUMOs of the edge and corner ligands, can also be understood by the MOs-based Lieb lattice model. The TB band fitting of CB1–CB3 within Lieb lattice model is shown in Supplementary Fig. 5. Due to the small dE between the LUMOs of corner and edge ligands, the CB1–CB3 formed Dirac-flat bands (Fig. 3c, d) are close to the ideal ones (Fig. 1c).
Ferromagnetism in sp 2C-COF and sp 2N-COF
In the experiments, the FM phase, with an easy axis perpendicular to the COF plane, i.e., along the z direction in Fig. 3a, b, is found in sp2C-COF after iodine oxidization at a low temperature (8.1 K)27. This is unusual as the sp2C-COF not containing any transition-metal atoms can evoke FM ordering, which inspires us to further understand the magnetic properties of these COFs after iodine doping. Our extensive test calculations indicate that the most important role of iodine doping is to introduce extra holes in sp2C-COF and sp2N-COF. The nh depends on the iodine concentration, but it is insensitive to the adsorption sites of iodine atoms. Besides of the role of hole doping, ionized iodine dopants can introduce some localized defect levels around EF. In the following, we will discuss the mechanism of nh-dependent magnetic phases in the main text and leave the results of iodine doping in Supplementary Fig. 6. According to the Mermin–Wagner theorem32, for a 2D system there is no realistic long-range order at finite temperature. Hence our discussion should be constrained on finite-size 2D lattices.
Figure 4 shows the band structures of COFs under a typical doping concentration of nh = 0.5 holes per unitcell (u.c.). The VB1 in both sp2C-COF and sp2N-COF can become partially occupied. According to the Stoner model33, the FM of iterative carriers can arise from partially occupied (except for the case of half-filling34) band with a sufficient narrow bandwidth, i.e., the smaller the bandwidth is, the larger the effective electronic interactions will be. Since VB1 in sp2C-COF has a rather narrower bandwidth than that in sp2N-COF, the Stoner criterion for realizing a FM ordering can be more easily satisfied in sp2C-COF than in sp2N-COF. It is found that the spontaneous spin polarization exists in both COFs under nh = 0.5 holes per u.c., and the COFs becomes metallic. The calculated exchange energy splitting of VB1 in sp2C-COF (sp2N-COF) is about 85 (40) meV. The spin density distributions of these two systems are mainly localized around the corner sites (see Supplementary Fig. 7). For the sp2C-COF, the magnetic moment is 0.5 μB per u.c. at nh = 0.5 holes per u.c., i.e., 1.0 μB per hole. The FM state is energetically more favorable than the AFM one by \(\Delta E = E_{{\mathrm{FM}}} - E_{{\mathrm{AFM}}} = - 5.0\) meV [or than the nonmagnetic (NM) one by −6.9 meV] (see Supplementary Table 1). For sp2N-COF, the magnetic moment is ~0.378 μB per unitcell and the FM ordering is slightly more stable than the AFM one (\(\Delta E = - 1.5\;{\mathrm{meV}}\)). Therefore, the stronger FM stability in sp2C-COF is induced by its narrower VB1 bandwidth, which originates from the larger dE and torsion angles in sp2C-COF than in sp2N-COF.
Fig. 4: Ferromagnetism in sp2C-COF and sp2N-COF.
DFT-calculated spin-polarized band structures of monolayer (a) sp2C-COF and (b) sp2N-COF with nh = 0.5 holes per u.c. DFT-calculated exchange energy splitting (Es) is 85 meV and 40 meV for sp2C-COF and sp2N-COF, respectively. c Calculated average magnetization per site as a function of temperature for monolayer sp2C-COF and sp2N-COF by Monto Carlo simulations. The transition temperatures (Tc) are marked with dashed vertical lines.
Although the magnetism originates from the many-body effects, which are usually underestimated by the mean-field approximation (MFA) within the DFT methods, the spin-polarized ground states of similar weakly-correlated systems have been correctly evaluated35,36,37,38,39,40,41 and confirmed in the experiments42,43. The onsite electron-electron (e–e) interactions (U) in the sp2C-COF and sp2N-COF (sp2-hybridized C and N systems) are very weak, compared to the traditional magnetic (3d- or 4f-) systems. On the other hand, the widely acceptable remedy within the DFT formalism to include the electron correlations is DFT + U method. We have performed a GGA + U calculation with U = 1.5 eV on the 2p orbitals of C and N atoms, and our calculations show that the enhanced U can slightly increase the exchange splitting Es (see Supplementary Fig. 8) and hence further enhance the stability of FM ground-states. Therefore, the conclusion of FM ordering in hole-doped COFs obtained from DFT calculations can be enhanced under DFT + U calculations.
To estimate the FM transition temperature (Tc), we have performed Monte Carlo (MC) simulations within the 2D Heisenberg model with uniaxial anisotropy (See "Methods" section for details), which is successfully applied to simulate the Tc of many 2D systems44,45. The reduced exchange integral J and magnetic moment m are 2.5 (1.31) \({\mathrm{meV}} \cdot {\upmu}_{\mathrm{B}}^{ - 2}\) and 0.5 (0.378) μB for sp2C-COF (sp2N-COF), respectively. The magnetic anisotropic energy (MAE = Eout-of-plane − Ein-plane), determined by the energy difference between the out-of-plane and in-plane spin configurations, is ~−0.1 meV for both COFs, which indicates that the magnetic easy axis is perpendicular to the COF planes, agreeing with the experimental observations. Although the value of MAE is small, it is very important for these two COFs to overcome the spin fluctuation and maintain a long-range magnetic ordering at finite temperature. The calculated T-dependent magnetic moment (per TPPy) is shown in Fig. 4c. Here, both sp2C-COF and sp2N-COF are under nh = 0.5 holes per u.c. It can be seen that the in-plane components of magnetization are zero at arbitrary temperature for both COFs, but their out-of-plane components are nonzero at low temperature. The calculated Tc is 5.3 K (close to 8.1 K of experiments) and 1.7 K for sp2C-COF and sp2N-COF, respectively.
Magnetic phase transitions in sp 2C-COF
Besides of nh = 0.5 holes per u.c., we have further explored the possible magnetic phase transitions as a function of nh at the DFT-level calculations. As shown in Fig. 3, the VB1 of sp2C-COF is contributed by the MO of corner ligands, which has a small effective inter-site (ligand) hopping t1 = 0.1 eV. Meanwhile, the effective onsite (ligand) U is usually in the range of 0.5~1.5 eV for organic systems46,47, which is significantly larger than t1. The calculated \(\Delta E\) as a function of nh is shown in Fig. 5a. When nh < 0.3 holes per u.c. (yellow region), the magnetic moment on each ligand is too small to evoke a preferred magnetic ordering, giving rise to a paramagnetic (PM) phase. When 0.3 < nh < 0.75 holes per u.c., the FM state becomes to be more favorable (\(\Delta E < 0\), green region) according to the Stoner's criterion. Interestingly, when \(\Delta E\) reaches a maximum value of −5.7 meV at nh = 0.55 hole per u.c., it becomes to be reduced gradually and finally reaches \(\Delta E = 0\) at nh = 0.75 holes per u.c. When 0.75 < nh < 1.0 holes per u.c., the effect of Fermi surface nesting plays a dominated role in determining its magnetic ground-state, as shown in Supplementary Fig. 9, giving rise to a Néel AFM state (blue region).
Fig. 5: Magnetic phase transitions as a function of hole doping concentrations.
a DFT-calculated \(\Delta E\left( {\Delta E = E_{{\mathrm{FM}}} - E_{{\mathrm{AFM}}}} \right)\) as a function of nh for monolayer sp2C-COF. b nh-dependent magnetic phase diagram for monolayer sp2C-COF, calculated by the MC simulations with 2D Heisenberg models with the uniaxial anisotropy.
Employing the Monte Carlo simulations within the 2D Heisenberg model with the uniaxial anisotropy, we further show nh dependent magnetic phase diagram of sp2C-COF in Fig. 5b. It can be seen that the FM or AFM phase can survive at finite temperatures under different nh. Interestingly, our calculated magnetic phase diagram of sp2C-COF system agrees well with that by J. E. Hirsch based on a similar model system (at the MFA-level calculations) with 6t1 < U < 10t134. Remarkably, the Tc of FM ordering can reach a maximum value of ~5.7 K for the sp2C-COF at nh~0.55 holes per u.c. Finally, we predict that at an extremely heavy hole doping concentration (nh = 3.0 holes per u.c.), the system can reach a half-metallic Dirac semimetal state, as shown in Supplementary Fig. 10.
t is emphasized that a Heisenberg model is adopted to simulate the Tc of various magnetic phases. Due to finite-size effect, a series of phase transition temperatures can be obtained. Taking such temperature as a proxy to a real Tc, we can conclude that the critical temperature of sp2C-COF is higher than that of sp2N-COF. Moreover, the more accurate phase diagram may need more complex methods beyond MFA, e.g., quantum Monte Carlo48,49 or dynamical mean-field theory50,51,52, which is out of the scope of the current study.
In conclusion, based on TB modeling and first-principles calculations, we predict that the sp2C-COF and sp2N-COF, which have been synthesized in the recent experiments, are actually the first two material realizations of organic-ligand-based Lieb lattices. The lattice distortion and inequivalent corner and edge sites make the electronic structures of these two systems deviate from the ideal Dirac-flat bands in a Lieb lattice. Interestingly, these two COF systems can be converted to the FM and AFM states at certain nh due to the less dispersive top valence bands. Our discovery not only can extend our understanding on the electronic and topological properties of distorted Lieb lattices, but also offers a possible way to realize exotic magnetic states in frustrated organic lattices, which opens the possibility to design the novel organic spintronic devices.
DFT calculations
Our ab initio calculations for the electronic structures of periodic structures were carried out within the framework of the Perdew-Burke-Ernzerhof generalized gradient approximation (PBE-GGA)53, as embedded in the Vienna ab initio simulation package code54,55. All the calculations were performed with a plane-wave cutoff energy of 450 eV. For the sp2C-COF, we adopted the experimental lattice constants as starting structure to perform geometric optimization. sp2N-COF holds a very similar structure as sp2C-COF based on our calculation. Both COFs were fully relaxed without any constraint until the force on each atom was <0.01 eV Å−1. The final lattice constants are a = b = 24.634 Å, c = 3.569 Å, α = 79.407°, β = 100.593° and γ = 88.299° for the sp2C-COF and a = b = 24.282 Å, c = 3.465 Å, α = 78.492°, β = 101.508° and γ = 88.490° for sp2N-COF. To eliminate the interlayer interaction, we introduced a vacuum layer of 18 Å thickness for monolayer calculations. The Brillouin zone k-point sampling was set with a 3 × 3 × 13 mesh for bulk unitcell, and a 3 × 3 × 1 one for monolayer calculations, respectively.
Chern number and Berry curvature calculations
The Chern number for one energy band can be calculated by the integral of the Berry curvature over the whole BZ
$$C_n = \frac{1}{{2{\uppi}}}\mathop {\smallint }\limits_{{\mathrm{BZ}}} {\mathrm{\Omega }}_n\left( {\mathbf{k}} \right)d^2{\mathbf{k}}$$
The Berry curvature Ωn(k) of the nth band can be determined as29,56,57
$${\mathrm{\Omega }}_n\left( {\mathbf{k}} \right) = - 2{\mathrm{Im}}\mathop {\sum }\limits_{m \ne n} \frac{{\psi _n\left( {\mathbf{k}} \right){\mathrm{|}}v_x{\mathrm{|}}\psi _m\left( {\mathbf{k}} \right)\psi _m\left( {\mathbf{k}} \right){\mathrm{|}}v_y{\mathrm{|}}\psi _n\left( {\mathbf{k}} \right)}}{{\left( {\varepsilon _{n{\mathbf{k}}} - \varepsilon _{m{\mathbf{k}}}} \right)^2}}$$
where \(\psi _n\left( {\mathbf{k}} \right)\) and εnk are the spinor Bloch wave function and eigenvalue of the nth band at k point, and vx(y) is the velocity operator. The total Chern number29 of a system is the sum of Cn over all the occupied bands (n):
$$C = \mathop {\sum }\limits_n C_n = \frac{1}{{2{\uppi}}}\mathop {\sum }\limits_n \mathop {\smallint }\limits_{{\mathrm{BZ}}} {\mathrm{\Omega }}_n\left( {\mathbf{k}} \right)d^2{\mathbf{k}}$$
2D Monte Carlo simulations
We have performed the Monte Carlo simulations based on the 2D Heisenberg model with uniaxial anisotropy, which is written as
$$H_{\mathrm{M}} = - \mathop {\sum }\limits_{\langle{i,j}\rangle} J{\mathbf{m}}_i \cdot {\mathbf{m}}_j - \mathop {\sum }\limits_i D\left( {m_i^z} \right)^2$$
where mi and \(m_i^z\) represent the magnetic moment and its z component on the ith moment, and i, j confines the NN neighboring mi and mj. The J is the exchange integral (\(J = - \Delta E/8m^2,{\mathrm{where}}\Delta E = E_{{\mathrm{FM}}} - E_{{\mathrm{AFM}}}\)) with the assumption of \(\left| {{\mathbf{m}}_i} \right| = m\). D is the reduced onsite MAE, and D = −MAE/m2. A 100 × 100 supercell containing 10,000 local magnetic moments are adopted, and each simulation lasts 109 loops to relaxation and another 109 loops to collect the physical quantities. In each loop, one moment is rotated to a random direction.
The data that support the findings of this study are available from the corresponding author upon reasonable request.
The MC code used for simulations is available by written request to the corresponding authors.
An amendment to this paper has been published and can be accessed via a link at the top of the paper.
Weeks, C. & Franz, M. Topological insulators on the Lieb and perovskite lattices. Phys. Rev. B 82, 085310 (2010).
Shen, R., Shao, L. B., Wang, B. & Xing, D. Y. Single Dirac cone with a flat band touching on line-centered-square optical lattices. Phys. Rev. B 81, 041410 (2010).
Dauphin, A., Müller, M. & Martin-Delgado, M. A. Quantum simulation of a topological Mott insulator with Rydberg atoms in a Lieb lattice. Phys. Rev. A 93, 043611 (2016).
Julku, A., Peotta, S., Vanhala, T. I., Kim, D.-H. & Törmä, P. Geometric origin of superfluidity in the Lieb-lattice flat band. Phys. Rev. Lett. 117, 045303 (2016).
Article ADS PubMed CAS Google Scholar
Lieb, E. H. Two theorems on the Hubbard model. Phys. Rev. Lett. 62, 1201–1204 (1989).
Article ADS MathSciNet CAS PubMed Google Scholar
Tamura, H., Shiraishi, K., Kimura, T. & Takayanagi, H. Flat-band ferromagnetism in quantum dot superlattices. Phys. Rev. B 65, 085324 (2002).
Tang, E., Mei, J.-W. & Wen, X.-G. High-temperature fractional quantum Hall states. Phys. Rev. Lett. 106, 236802 (2011).
Neupert, T., Santos, L., Chamon, C. & Mudry, C. Fractional quantum Hall states at zero magnetic field. Phys. Rev. Lett. 106, 236804 (2011).
Wang, Y. F., Gu, Z. C., Gong, C., De & Sheng, D. N. Fractional quantum Hall effect of hard-core bosons in topological flat bands. Phys. Rev. Lett. 107, 146803 (2011).
Sheng, D. N., Gu, Z., Sun, K. & Sheng, L. Fractional quantum Hall effect in the absence of Landau levels. Nat. Commun. 2, 385–389 (2011).
Drost, R., Ojanen, T., Harju, A. & Liljeroth, P. Topological states in engineered atomic lattices. Nat. Phys. 13, 668–671 (2017).
Slot, M. R. et al. Experimental realization and characterization of an electronic Lieb lattice. Nat. Phys. 13, 672–676 (2017).
Guzmán-Silva, D. et al. Experimental observation of bulk and edge transport in photonic Lieb lattices. N. J. Phys. 16, 063061 (2014).
Mukherjee, S. et al. Observation of a localized flat-band state in a photonic Lieb lattice. Phys. Rev. Lett. 114, 245504 (2015).
Vicencio, R. A. et al. Observation of localized states in Lieb photonic lattices. Phys. Rev. Lett. 114, 245503 (2015).
Taie, S. et al. Coherent driving and freezing of bosonic matter wave in an optical Lieb lattice. Sci. Adv. 1, e1500854–e1500854 (2015).
Article ADS PubMed PubMed Central CAS Google Scholar
Cote, A. P. et al. Porous, crystalline, covalent organic frameworks. Science 310, 1166–1170 (2005).
Wan, S. et al. Covalent organic frameworks with high charge carrier mobility. Chem. Mater. 23, 4094–4097 (2011).
Ascherl, L. et al. Molecular docking sites designed for the generation of highly crystalline covalent organic frameworks. Nat. Chem. 8, 310–316 (2016).
Vyas, V. S. et al. A tunable azine covalent organic framework platform for visible light-induced hydrogen generation. Nat. Commun. 6, 8508 (2015).
Huang, N. et al. Multiple-component covalent organic frameworks. Nat. Commun. 7, 12325 (2016).
Mateo-Alonso, A. Pyrene-fused pyrazaacenes: from small molecules to nanoribbons. Chem. Soc. Rev. 43, 6311 (2014).
Granda, J. M., Grabowski, J. & Jurczak, J. Synthesis, structure, and complexation properties of a C3-symmetrical triptycene-based anion receptor: selectivity for dihydrogen phosphate. Org. Lett. 17, 5882–5885 (2015).
Zhou, T.-Y., Xu, S.-Q., Wen, Q., Pang, Z.-F. & Zhao, X. One-step construction of two different kinds of pores in a 2D covalent organic framework. J. Am. Chem. Soc. 136, 15885–15888 (2014).
Pachfule, P. et al. Diacetylene functionalized covalent organic framework (COF) for photocatalytic hydrogen generation. J. Am. Chem. Soc. 140, 1423–1427 (2018).
Diercks, C., Kalmutzki, M. & Yaghi, O. Covalent organic frameworks—organic chemistry beyond the molecule. Molecules 22, 1575 (2017).
Article PubMed Central CAS Google Scholar
Jin, E. et al. Two-dimensional sp2 carbon–conjugated covalent organic frameworks. Science 357, 673–676 (2017).
Kane, C. L. & Mele, E. J. Quantum spin Hall effect in graphene. Phys. Rev. Lett. 95, 226801 (2005).
Thouless, D. J., Kohmoto, M., Nightingale, M. P. & den Nijs, M. Quantized Hall conductance in a two-dimensional periodic potential. Phys. Rev. Lett. 49, 405–408 (1982).
Fukui, T., Hatsugai, Y. & Suzuki, H. Chern numbers in discretized Brillouin zone: efficient method of computing (spin) Hall conductances. J. Phys. Soc. Jpn. 74, 1674–1677 (2005).
Bredas, J. L., Calbert, J. P., da Silva Filho, D. A. & Cornil, J. Organic semiconductors: a theoretical characterization of the basic parameters governing charge transport. Proc. Natl Acad. Sci USA. 99, 5804–5809 (2002).
Mermin, N. D. & Wagner, H. Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models. Phys. Rev. Lett. 17, 1133–1136 (1966).
Stoner, E. C. Collective electron ferromagnetism. Proc. R. Soc. A Math. Phys. Eng. Sci. 165, 372–414 (1938).
Article ADS MATH Google Scholar
Hirsch, J. E. Two-dimensional Hubbard model: numerical simulation study. Phys. Rev. B 31, 4403–4419 (1985).
Son, Y.-W., Cohen, M. L. & Louie, S. G. Half-metallic graphene nanoribbons. Nature 444, 347–349 (2006).
Mounet, N. et al. Two-dimensional materials from high-throughput computational exfoliation of experimentally known compounds. Nat. Nanotechnol. 13, 246–252 (2018).
Fu, Y. & Singh, D. J. Applicability of the strongly constrained and appropriately normed density functional to transition-metal magnetism. Phys. Rev. Lett. 121, 207201 (2018).
Liu, Z., Wang, Z.-F., Mei, J.-W., Wu, Y.-S. & Liu, F. Flat Chern band in a two-dimensional organometallic framework. Phys. Rev. Lett. 110, 106804 (2013).
Peng, H. et al. Origin and Enhancement of Hole-Induced Ferromagnetism in First-Row d0 Semiconductors. Phys. Rev. Lett. 102, 017201 (2009).
Liu, L., Yu, P. Y., Ma, Z. & Mao, S. S. Ferromagnetism in GaN:Gd: A Density Functional Theory Study. Phys. Rev. Lett. 100, 127203 (2008).
Pan, H. et al. Room-temperature ferromagnetism in carbon-doped ZnO. Phys. Rev. Lett. 99, 127201 (2007).
Ruffieux, P. et al. On-surface synthesis of graphene nanoribbons with zigzag edge topology. Nature 531, 489–492 (2016).
Magda, G. Z. et al. Room-temperature magnetic order on zigzag edges of narrow graphene nanoribbons. Nature 514, 608–611 (2014).
Huang, C. et al. Toward intrinsic room-temperature ferromagnetism in two-dimensional semiconductors. J. Am. Chem. Soc. 140, 11519–11525 (2018).
Xu, C., Feng, J., Xiang, H. & Bellaiche, L. Interplay between Kitaev interaction and single ion anisotropy in ferromagnetic CrI3 and CrGeTe3 monolayers. npj Comput. Mater. 4, 57 (2018).
Yoshitake, M. et al. Reflectance spectra of the 1:1 salts of bis(methylenedithio)tetrathiafulvalene (BMDT-TTF): estimation of the on-site Coulomb energy. Bull. Chem. Soc. Jpn. 61, 1115–1119 (1988).
Tosatti, E., Fabrizio, M., Tóbik, J. & Santoro, G. E. Strong correlations in electron doped phthalocyanine conductors near half filling. Phys. Rev. Lett. 93, 117002 (2004).
Ma, N. et al. Anomalous quantum-critical scaling corrections in two-dimensional antiferromagnets. Phys. Rev. Lett. 121, 117202 (2018).
Manousakis, E. The spin-½ Heisenberg antiferromagnet on a square lattice and its application to the cuprous oxides. Rev. Mod. Phys. 63, 1–62 (1991).
Nomura, Y., Sakai, S., Capone, M. & Arita, R. Unified understanding of superconductivity and Mott transition in alkali-doped fullerides from first principles. Sci. Adv. 1, e1500568 (2015).
Lichtenstein, A. I., Katsnelson, M. I. & Kotliar, G. Finite-temperature magnetism of transition metals: an ab initio dynamical mean-field theory. Phys. Rev. Lett. 87, 067205 (2001).
Kotliar, G. et al. Electronic structure calculations with dynamical mean-field theory. Rev. Mod. Phys. 78, 865–951 (2006).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).
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).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996).
Fang, Z. The anomalous Hall effect and magnetic monopoles in momentum space. Science 302, 92–95 (2003).
Yao, Y. et al. First principles calculation of anomalous Hall conductivity in ferromagnetic bcc Fe. Phys. Rev. Lett. 92, 037204 (2004).
B.C., X.Z., and D.L. acknowledge the support from NSFC (Nos. 11574118). J.W. and B.H. acknowledge the support from Science Challenge Project (No. TZ2016003), NSFC (No. 11574024) and NSAF (No. U1930402). B.C. also acknowledges the support from Shandong Provincial Natural Science Foundation (No. ZR2019MA64) and the NSFC (No. 11404188). Computations were performed at Tianhe2-JK at CSRC. The authors thank Dr. Z. Liu (Tsinghua University), Dr. S. Yin and Dr. M. Zhao (Shandong University) for helpful discussions.
School of Physics, National Demonstration Center for Experimental Physics Education, Shandong University, Jinan, 250100, China
Bin Cui, Xingwen Zheng, Desheng Liu & Shijie Xie
Beijing Computational Science Research Center, Beijing, 100193, China
Jianfeng Wang & Bing Huang
Bin Cui
Xingwen Zheng
Jianfeng Wang
Desheng Liu
Shijie Xie
Bing Huang
B.C. and B.H. directed the project. B.C. and X.Z. calculated the results. B.C., J.W., D. L., S. X. and B.H. analyzed the results. B.C., J.W. and B.H. wrote the manuscript. All authors discussed the results and commented.
Correspondence to Bin Cui or Bing Huang.
Peer review information Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Peer Review File
Cui, B., Zheng, X., Wang, J. et al. Realization of Lieb lattice in covalent-organic frameworks with tunable topology and magnetism. Nat Commun 11, 66 (2020). https://doi.org/10.1038/s41467-019-13794-y
DOI: https://doi.org/10.1038/s41467-019-13794-y
Geometric characterization of anomalous Landau levels of isolated flat bands
Yoonseok Hwang
Jun-Won Rhim
Bohm-Jung Yang
Nature Communications (Nat Commun) ISSN 2041-1723 (online)