Electronic Structure and Minimal Models for Flat and Corrugated CuO Monolayers: An Ab Initio Study

Abstract

$\mathrm{CuO}$ atomic thin monolayer ($\mathrm{mlCuO}$) was synthesized recently. Interest in the $\mathrm{mlCuO}$ is based on its close relation to ${\mathrm{CuO}}_2$ layers in typical high temperature cuprate superconductors. Here, we present the calculation of the band structure, the density of states and the Fermi surface of the flat $\mathrm{mlCuO}$ as well as the corrugated $\mathrm{mlCuO}$ within the density functional theory (DFT) in the generalized gradient approximation (GGA). In the flat $\mathrm{mlCuO}$, the $\mathrm{Cu}$-$3d_{x^2-y^2}$ band crosses the Fermi level, while the $\mathrm{Cu}$-$3d_{xz,yz}$ hybridized band is located just below it. The corrugation leads to a significant shift of the $\mathrm{Cu}$-$3d_{xz,yz}$ hybridized band down in energy and a degeneracy lifting for the $\mathrm{Cu}$-$3d_{x^2-y^2}$ bands. Corrugated $\mathrm{mlCuO}$ is more energetically favorable than the flat one. In addition, we compared the electronic structure of the considered $\mathrm{CuO}$ monolayers with bulk $\mathrm{CuO}$ systems. We also investigated the influence of a crystal lattice strain (which might occur on some interfaces) on the electronic structure of both $\mathrm{mlCuO}$ and determined the critical strains of topological Lifshitz transitions. Finally, we proposed a number of different minimal models for the flat and the corrugated $\mathrm{mlCuO}$ using projections onto different Wannier functions basis sets and obtained the corresponding Hamiltonian matrix elements in a real space.

1. Introduction

Copper oxides stand apart from other transition metal compounds. First of all, they attract much attention because of their high temperature superconductivity (HTSC) and already existing applications [1,2,3,4,5] as catalysts [6], photocells [7] and thin-film transistors [8]. $\mathrm{CuO}$ is an exceptional member of the generally rocksalt family ($\mathrm{MnO}$ to $\mathrm{CuO}$), as it deviates both structurally and electronically from others. Unlike other members of the $3d$ transition oxides, which crystallize in the cubic rocksalt structure (with possible rhombohedral distortions), Tenorite (CuO) crystallizes in the lower symmetry monoclinic (C2/c) crystal structure [9], albeit the cubic crystal structure is also possible [10]. Thus far, the bulk compound $\mathrm{CuO}$ has been thoroughly studied using ab initio calculations: DFT+U [11,12,13,14,15], DFT with hybrid functional [4] and Charge Transition Level Approach [16].

The relatively recent interest in the $\mathrm{CuO}$ monolayer arose in part because one would expect a superconducting phase to occur here by analogy with the typical representative of HTSC cuprates ${\mathrm{La}}_2{\mathrm{CuO}}_4$. There, the superconductivity occurs in two-dimensional layers formed by ${\mathrm{CuO}}_2$ plaquettes. The $\mathrm{CuO}$ monolayer consists of the same plaquettes. However, the plaquettes in the cuprates are connected by vertices, whereas in the monolayer they are connected by faces. This fact leads to a difference in the chemical composition—the number of copper and oxygen atoms is identical in the monolayer ($\mathrm{CuO}$), while there are two oxygen atoms per copper atom in the cuprates (${\mathrm{CuO}}_2$).

In general, the electronic properties of the copper oxides are thoroughly studied. As everyone knows, ${\mathrm{La}}_2{\mathrm{CuO}}_4$ has the $\mathrm{Cu}$-$d_{x^2-y^2}$ orbital at the Fermi level [5]. In Ref. [17] the authors studied isolated ${\mathrm{CuO}}_2$ monolayer using DFT and showed that all $\mathrm{Cu}$ orbitals, except for $d_{z^2}$, have states near the Fermi level; the $\mathrm{Cu}$-$d_{xz}$ and $d_{yz}$ orbitals and the $\mathrm{O}$-$p_z$ orbital have most of the states near the Fermi level, leading to $\pi$ bonds in the entire monolayer. The bulk $\mathrm{CuO}$ with the monoclinic structure, known as a p-type semiconductor, has significantly three dimensional electronic structure with mainly $\mathrm{Cu}$-$3d$ ($\mathrm{Cu}$-$3d_{x^2-y^2}$) states near the Fermi level [11]. DFT study in Ref. [14] claims that the bulk $\mathrm{CuO}$ with the cubic structure is an indirect gap semiconductor; its valence band consists mainly of $\mathrm{O}$-$2p$ and $\mathrm{Cu}$-$3d$ orbitals. We can conclude that the electronic structure of the ${\mathrm{CuO}}_2$ monolayer and ${\mathrm{La}}_2{\mathrm{CuO}}_4$ is quite similar to each other, while the one of the flat $\mathrm{mlCuO}$ partially resembles them, but has some qualitatively differences, such as an extremum presence in the $\mathsf{\Gamma }-M$ direction (as the reader can observe later).

Moreover, a number of experiments and theoretical studies were carried out on various structural modifications of the $\mathrm{CuO}$ monolayer: $\mathrm{CuO}$ monolayer in a graphene pores and freestanding $\mathrm{CuO}$ monolayer [18], $\mathrm{CuO}$ monolayer on a graphene substrate [19], different combinations of $\mathrm{CuO}$ monolayers as an interface between bilayer graphene and finding thermodynamically stable freestanding $\mathrm{CuO}$ monolayer using the evolutionary algorithm [20]. In Ref. [18], the authors showed that freestanding the perfectly flat $\mathrm{CuO}$ monolayer can be corrugated in some cases (Figure 1b). This perfectly flat $\mathrm{mlCuO}$ can be easily constructed from the bulk $\mathrm{CuO}$ with the cubic structure. In Figure 1a, we show the $2\times 2\times 1$ supercell of the cubic $\mathrm{mlCuO}$. It matches with the perfectly flat $\mathrm{mlCuO}$ structure (Figure 1b, left) when rotated by ${45}^{\circ }$. The corrugated $\mathrm{mlCuO}$ can be constructed from the bulk $\mathrm{CuO}$ with the monoclinic structure, though in the monoclinic system atoms are much more displaced; see Figure 1c and compare it with Figure 1b, right. Thus, the corrugated state of $\mathrm{mlCuO}$ can be described as a transitional one relative to the $\mathrm{CuO}$ systems with the cubic and the monoclinic structures.

The description of the electronic structure of $\mathrm{mlCuO}$ systems is also far from being complete—the authors were able to find some data in Ref. [18,20], but no more than that. Besides, at the moment the authors are not aware of any works where a minimal model has been formulated neither for an isolated perfectly flat monolayer $\mathrm{CuO}$, nor for more complex crystal structures such as the corrugated monolayer, monolayer on a substrate or monolayer as an interface. Thus, it seems necessary to obtain on a more systematic basis the densities of states, the band structures and the Fermi surfaces for all listed monolayer $\mathrm{CuO}$ systems and to formulate a minimal model for them with the corresponding Hamiltonian parameter values.

In this work, we solve a task of proposing and comparing different minimal models for the $\mathrm{CuO}$ monolayer systems as a necessary first step of any further theoretical investigations.

2. Crystal Structure and Calculation Details

To calculate the band structure, the density of states (DOS) and the Fermi surface, we used the density functional theory with the full-potential linear augmented plane-wave framework, as implemented in WIEN2k [21] together with the generalized gradient approximation by Perdew, Burke and Ernzerhof [22], to the exchange-correlation functional.

Figure 2 shows the crystal structures of the systems discussed in this paper. The flat $\mathrm{CuO}$ monolayer space group is a 123 (P4/mmm). The lattice parameter is $a=2.69$ Å [18]. Atoms occupy the following positions: $\mathrm{Cu}$ 1a $(0,0,0)$ and $\mathrm{O}$ 1c $(0.5,0.5,0)$. We used a 20 Bohr vacuum gap. In order to construct the corrugated $\mathrm{CuO}$ monolayer, we doubled the unit cell and made these new additional $\mathrm{Cu}$ and $\mathrm{O}$ atoms unequivalent to the original ones by applying a small shift about 0.5 Bohr in z direction only for them. Resulting system has 59 (Pmmn) [origin choice 2] space group. Next, we did a set of structural relaxations with 10, 20 and 40 Bohr vacuum gaps. There was no difference between 20 and 40 Bohr vacuum gap cases, so in all further calculations we used the 20 Bohr gap. After relaxation neighboring $\mathrm{Cu}$ atoms shifted in z direction about $\pm 0.31$ Å relative to their original positions, while neighboring $\mathrm{O}$ atoms barely shifted at all. Since the final corrugated structure of $\mathrm{CuO}$ monolayer has a doubled unit cell and is rotated by ${45}^{\circ }$ relative to the flat $\mathrm{mlCuO}$, we cannot directly compare their calculated electronic structures. Thus, we use additional system—the flat $\mathrm{mlCuO}$ with a doubled (and rotated) unit cell—in order to make a proper comparison. Moreover, for the flat $\mathrm{mlCuO}$ we did a ${45}^{\circ }$ rotation of a local coordinate system in order to use a typical orbital convention, such as in cuprate compounds.

All calculations were nonmagnetic and converged self-consistently on a grid of 24 × 24 × 1 k-points in the irreducible Brillouin zone using the Monkhorst-Pack method [23]. We used energy convergence limit 0.1 mRy, force convergence limit 0.5 mRy/a.u. for optimization, $\mathrm{RKmax}=7$, $\mathrm{Gmax}=12$, energy separation $-6.0$ Ry. In Figure 2d, we show the Brillouin zone with the k-path used in the band structure analysis.

3. Results and Discussion

3.1. Electronic Structure

Figure 3 shows the DFT (GGA) band structures, the densities of states, bands with their orbital characters and the Fermi surfaces. The first row (a–c) of Figure 3 shows the results for the flat $\mathrm{mlCuO}$. The band structure in Figure 3a shows that there is an isolated set of bands in the range from $-8$ eV to $2.3$ eV resembling typical $\mathrm{Cu}$-based HTSC ${\mathrm{La}}_2{\mathrm{CuO}}_4$. It has only the $\mathrm{Cu}$-$3d$ and the $\mathrm{O}$-$2p$ states. The electronic bands of the flat $\mathrm{mlCuO}$ at the Fermi level are formed by the $\mathrm{Cu}$-$3d_{x^2-y^2}$ states (with small addition of the hybrid $\mathrm{O}$-$2p$ states) in consistence with the known results [11]. Note that there is a second band that almost crosses the Fermi level—it is only 0.02 eV lower. It includes the $\mathrm{Cu}$-$3d_{xz,yz}$ states hybridized with the $\mathrm{O}$-$2p_z$. The Fermi surface has a hole pocket around the X point.

The second row (d–f) of Figure 3 shows the results for the $\mathrm{mlCuO}$ with a doubled unit cell. Due to ${45}^{\circ }$ rotation of a unit cell, its band structure differs from that of the flat $\mathrm{mlCuO}$ in a much more complex way than simply by a large number of bands. The two bands crossing the Fermi level originate from two $\mathrm{Cu}$ atoms in the unit cell. They are formed by the $\mathrm{Cu}$-$3d_{x^2-y^2}$ states. Regarding them, we can note degeneracy lifting in the $\mathsf{\Gamma }-M$ direction. As in the flat $\mathrm{mlCuO}$, there are the $\mathrm{Cu}$-$3d_{xz,yz}$ states (four bands) just below the Fermi level. The Fermi surface has two hole pockets around the M point.

The third row (g–i) of Figure 3 shows the results for the corrugated $\mathrm{mlCuO}$. Its electronic structure is rather similar to that in the $\mathrm{mlCuO}$ with a doubled unit cell, but there are two notable differences. First, a significant shift of the $\mathrm{Cu}$-$3d_{xz,yz}$ bands to $-0.7$ eV. Second, a noticeably larger degeneracy lifting for the $\mathrm{Cu}$-$3d_{x^2-y^2}$ bands in the $\mathsf{\Gamma }-M$ direction. On top of that, the corrugated $\mathrm{mlCuO}$ total energy turns out to be lower than the flat $\mathrm{mlCuO}$ one by $0.07$ eV. In other words, the corrugated state appears to be more favorable and if we have the flat $\mathrm{mlCuO}$ as a topmost layer of some surface, it will most likely be corrugated.

Seeing such a significant shift of the $\mathrm{Cu}$-$3d_{xz,yz}$ bands of the corrugated $\mathrm{mlCuO}$, we wondered if it was possible to raise these states to the Fermi level only via the lattice strain. Significant lattice strain is observed at all sorts of interfaces where there is a mismatch between the lattice parameters. Besides, as we already mention, the corrugated state appears to be more energetically favorable.

To clarify this issue, we carried out a series of calculations where we varied the lattice parameter a from 0% to 10% for the flat $\mathrm{mlCuO}$ and from 0% to 35% for the corrugated $\mathrm{mlCuO}$. The corresponding results can be observed in Figure 4. For the flat $\mathrm{mlCuO}$, the lattice deformation $\Delta a=0.7\%$ leads to a topological Lifshitz transition with the appearance of a new hole pocket around the X point. Clearly, it is a very minor lattice parameter change of the order of experiment accuracy. For the corrugated $\mathrm{CuO}$ monolayer, such a transition requires a much larger lattice deformation; it appears only at $\Delta a=35\%$. Of course, such a strain is too large, and we bring it here only as an illustration. However, for the corrugated $\mathrm{mlCuO}$ case, we want to note a presence of what seems to be a flat band in the $\mathsf{\Gamma }-X$ direction near the Fermi level. It is likely that a flat band at the Fermi level can be obtained using a reasonable lattice strain and a hole doping.

Let us now return to the idea we proposed in the introduction: the corrugated of $\mathrm{mlCuO}$ can be described as a transitional system in between the $\mathrm{CuO}$ with the cubic and monoclinic structures (Figure 1). Table 1 shows $\mathrm{Cu}$-$\mathrm{Cu}$ and $\mathrm{Cu}$-$\mathrm{O}$ distances in the bulk $\mathrm{CuO}$ systems and the $\mathrm{CuO}$ monolayers. In the bulk $\mathrm{CuO}$ with the cubic structure, $\mathrm{CuO}$ layers are significantly stretched—by about $11.5\%$ relative to the flat $\mathrm{mlCuO}$, which is close to the considered lattice strain $\Delta a=10\%$ (Figure 4c). We studied the monolayer made from the bulk $\mathrm{CuO}$ by simply adding a vacuum between layers (getting slab+vacuum) in z direction and using no relaxation (Figure 1a). The final crystal and band structures were nearly identical to the flat $\mathrm{mlCuO}$ ones with lattice strain $\Delta a=10\%$ (as in Figure 4c,d); thus, we do not include them.

There are more surprising results for the bulk $\mathrm{CuO}$ with the monoclinic structure. Its $\mathrm{CuO}$ layers are also stretched, but to a noticeably lesser extent—by about $5.1\%$ relative to the corrugated $\mathrm{mlCuO}$. Its crystal structure has more pronounced corrugation pattern as compared to the corrugated $\mathrm{mlCuO}$ (Figure 1c). To calculate the monoclinic $\mathrm{CuO}$ monolayer, we again made a slab+vacuum in z direction and used no relaxation; moreover, we did a ${45}^{\circ }$ rotation of a local coordinate system to remain consistency with the orbital convention chosen for the corrugated $\mathrm{mlCuO}$. The final space group was 13 (P2/c). The results are given in Figure 5. The most surprising result is that the monoclinic $\mathrm{mlCuO}$ has the band gap at the Fermi level, which opens because of symmetry lowering without using DFT+U or hybrid potentials, as (but for the bulk compound) in Ref. [4,11]. There are mainly the $\mathrm{Cu}$-$3d_{x^2-y^2}$ (with the hybrid $\mathrm{O}$-$2p$) states near the Fermi level similar to the corrugated $\mathrm{mlCuO}$. However, it is rather difficult to compare these results in detail with the results for the corrugated $\mathrm{mlCuO}$ (Figure 3g,h) due to the different space group. Let us conclude here that a more complex corrugation pattern leads to a more complex bands structure.

3.2. Minimal Models

The next step in our study is to reveal a good minimal model for both flat and corrugated $\mathrm{mlCuO}$. In order to do that, we constructed a set of models using the maximally localised Wannier functions (MLWF) within wannier90 package [24]. We examined the following set of models: single-band model with Wannier projected $\mathrm{Cu}$-$d_{x^2-y^2}$, three-band model with $\mathrm{Cu}$-$d_{x^2-y^2},d_{xz},d_{yz}$, five-band model with $\mathrm{Cu}$-$d_{x^2-y^2},d_{xz},d_{yz}$, $\mathrm{O}$-$p_x,p_y$ and eight-band model with $\mathrm{Cu}$-d, $\mathrm{O}$-p. The corresponding band structures are depicted in Figure 6. Note that there are actually twice the number of bands for the corrugated $\mathrm{mlCuO}$ and the $\mathrm{mlCuO}$ with a doubled unit cell due to unit cell doubling in contrast to the flat $\mathrm{mlCuO}$. The eight-band model resulting band structure is in the excellent agreement with GGA calculated one; in other cases, the agreement is fine at the Fermi level. Thus, the single-band $\mathrm{Cu}$-$d_{x^2-y^2}$ model for the flat and the corrugated $\mathrm{mlCuO}$ can be used as a minimal model.

We also present on-site energies and hopping integrals for one- and three-band models for the flat (Table A1), the corrugated $\mathrm{mlCuO}$ and the $\mathrm{mlCuO}$ with a doubled unit cell (Table A2 and Table A3). The corresponding hopping schemes are illustrated in Figure 7. We used hoppings up to the second coordination sphere for the flat $\mathrm{mlCuO}$ and up to the fifth coordination sphere for the corrugated one. These numbers of neighbors are the minimum required to obtain good agreement between the model Hamiltonian band structure and the initial one. We also attach the values of the Hamiltonian matrix elements in a real space for all the models (Figure 6) as machine-readable data files to Supplementary Materials. So, depending on a task, an interested reader can use the Hamiltonian of the appropriate complexity.

4. Conclusions

We investigated the electronic properties of the flat and the corrugated $\mathrm{CuO}$ monolayers—DOS, the band structures and the Fermi surfaces. The flat $\mathrm{mlCuO}$ is similar to the bulk $\mathrm{CuO}$ with the monoclinic crystal structure and typical $\mathrm{Cu}$-based HTSC, e.g., ${\mathrm{La}}_2{\mathrm{CuO}}_4$ and has the $\mathrm{Cu}$-$3d_{x^2-y^2}$ states at the Fermi level (with a small addition of the hybrid $\mathrm{O}$-$2p$ states). There is a second band just below the Fermi level—that it is only 0.02 eV lower. It includes the $\mathrm{Cu}$-$3d_{xz,yz}$ states hybridized with $\mathrm{O}$-$2p_z$.

The corrugation effect leads to a significant shift of the $\mathrm{Cu}$-$3d_{xz,yz}$ bands to $-0.7$ eV and a degeneracy lifting for the $\mathrm{Cu}$-$3d_{x^2-y^2}$ bands. The corrugated $\mathrm{mlCuO}$ is energetically more favorable than the flat one by $0.07$ eV and is more likely formed as a topmost layer on some surfaces.

It is possible to create a topological Lifshitz transition via lattice strain: for the flat $\mathrm{mlCuO}$, a slight stretching of the lattice parameter $\Delta a=0.7\%$ already leads to the appearance of a hole pocket around the X point; this effect for the corrugated $\mathrm{mlCuO}$ occurs at $\Delta a=35\%$. It is interesting to note the presence of what seems to be a rather flat band in the $\mathsf{\Gamma }-X$ direction near the Fermi level. Probably, it will be possible to create it using a reasonable strain and a hole doping.

There is a significant mismatch in the lattice parameters of the considered $\mathrm{CuO}$ monolayers and the known bulk $\mathrm{CuO}$ systems (with the cubic and the monoclinic crystal structures). We conclude that $\mathrm{CuO}$ layers in the bulk $\mathrm{CuO}$ are stretched relative to the monolayer systems—by about $11.5$% (cubic bulk vs. flat monolayer) and by $5.1$% (monoclinic bulk vs. corrugated monolayer). The monolayer made from experimental bulk $\mathrm{CuO}$ with the monoclinic crystal structure turns out to have the band gap in our DFT calculation without using DFT+U or hybrid potentials. Clearly, the reason for this is the more complex corrugation patterns of its crystal structure.

We also suggested a set of minimal models for the flat and the corrugated $\mathrm{CuO}$ monolayers. The simplest model includes only the $\mathrm{Cu}$-$d_{x^2-y^2}$ states and agrees well with the GGA calculated band structure at the Fermi level (the latter is also correct for the other models). For the one- and three-band models, we obtained the values of the corresponding Hamiltonian matrix elements in a real space; therefore, depending on the needs, the Hamiltonian of the appropriate complexity can be used.

Before proceeding to more complex systems, such as a $\mathrm{CuO}$ monolayer on a substrate or a $\mathrm{CuO}$ monolayer as an interface, it was necessary to perform calculations for the original system and understand the features of its electronic structure. So, the next step might be the investigation of such complex systems or using more advanced methods such as DFT+DMFT.

Author Contributions

Conceptualization, I.A.N., M.M.K. and S.G.O.; calculations, A.A.S. and L.V.B.; writing, A.A.S. and I.A.N.; review and editing, I.A.M., M.M.K. and S.G.O.; funding acquisition, I.A.N. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

L.V.B. would like to thank the Irkutsk Supercomputer Center of SB RAS for providing the access to HPC-cluster «Akademik V.M. Matrosov» (Irkutsk Supercomputer Center of SB RAS, Irkutsk: ISDCT SB RAS; http://hpc.icc.ru, accessed 21 November 2022).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Table

Figure 1. The [Forumla omitted. See PDF.] with the cubic ([Forumla omitted. See PDF.] supercell) (a) and monoclinic (c) crystal structure, top and side view (following Ref. [14]). Freestanding perfectly flat and corrugated [Forumla omitted. See PDF.] monolayer crystal structure (b) (following Ref. [18]). Blue denote [Forumla omitted. See PDF.] atoms, red is [Forumla omitted. See PDF.] atoms. The figures are on the same spatial scale.

Figure 2. The flat [Forumla omitted. See PDF.] (a), the [Forumla omitted. See PDF.] with a doubled unit cell (b), the corrugated [Forumla omitted. See PDF.] ([Forumla omitted. See PDF.] supercell) (c) crystal structures considered in this work, top and side view; brown denote [Forumla omitted. See PDF.] atoms, red is [Forumla omitted. See PDF.] atoms. Brillouin zone for the doubled unit cell and corrugated [Forumla omitted. See PDF.] (d).

Figure 3. DFT (GGA) calculated DOS, the band structures, the band structures with the orbital characters and the Fermi surface of the flat [Forumla omitted. See PDF.] (a–c), the [Forumla omitted. See PDF.] with a doubled unit cell (d–f), the corrugated [Forumla omitted. See PDF.] (g–i). Zero corresponds to the Fermi level.

Figure 4. DFT (GGA) calculated Fermi surface with [Forumla omitted. See PDF.] (zero strain), [Forumla omitted. See PDF.], [Forumla omitted. See PDF.] and the band structure comparison for the flat [Forumla omitted. See PDF.] (a–d); the Fermi surface with [Forumla omitted. See PDF.] (zero strain), [Forumla omitted. See PDF.] and the band structure comparison for the corrugated [Forumla omitted. See PDF.] (e–g). Zero corresponds to the Fermi level.

Figure 5. DFT (GGA) calculated band structure, the band structure with the orbital characters and the corresponding Brillouin zone (a–c) of the monoclinic [Forumla omitted. See PDF.]. Zero corresponds to the Fermi level.

Figure 6. Comparison DFT (GGA) band structure with a Wannier projected one: [Forumla omitted. See PDF.]-[Forumla omitted. See PDF.] (a), [Forumla omitted. See PDF.]-[Forumla omitted. See PDF.] (b), [Forumla omitted. See PDF.]-[Forumla omitted. See PDF.], [Forumla omitted. See PDF.]-[Forumla omitted. See PDF.] (c), [Forumla omitted. See PDF.]-d, [Forumla omitted. See PDF.]-p (d) for the [Forumla omitted. See PDF.] with a doubled unit cell; same for the corrugated [Forumla omitted. See PDF.] (e–h). Zero corresponds to the Fermi level.

Figure 7. Hoppings schemes for the flat [Forumla omitted. See PDF.] (left) and the [Forumla omitted. See PDF.] with a doubled unit cell (right). A [Forumla omitted. See PDF.] supercell is displayed for both structures. Yellow transparent square denotes a single unit cell.

Table 1

$\mathrm{Cu}$-$\mathrm{Cu}$ and $\mathrm{Cu}$-$\mathrm{O}$ distances in the bulk and the monolayer $\mathrm{CuO}$ systems.

System	d($\mathbf{Cu}$-$\mathbf{Cu}$), Å	d($\mathbf{Cu}$-$\mathbf{O}$), Å	Ref.
Cubic bulk $\mathrm{CuO}$	$3.00$	$2.12$	[10]
Monoclinic bulk $\mathrm{CuO}$	$2.90$	$1.96$	[9]
Flat $\mathrm{mlCuO}$	$2.69$	$1.90$	This work, ref. [18]
Corrugated $\mathrm{mlCuO}$	$2.76$	$1.93$	This work

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ma16020658/s1.

References

1. Anisimov, V.I.; Zaanen, J.; Andersen, O.K. Band theory and Mott insulators: Hubbard U instead of Stoner I. Phys. Rev. B; 1991; 44, pp. 943-954. [DOI: https://dx.doi.org/10.1103/PhysRevB.44.943] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9999600]

2. Ruiz, E.; Alvarez, S.; Alemany, P.; Evarestov, R.A. Electronic structure and properties of Cu₂O. Phys. Rev. B; 1997; 56, pp. 7189-7196. [DOI: https://dx.doi.org/10.1103/PhysRevB.56.7189]

3. Ghijsen, J.; Tjeng, L.H.; van Elp, J.; Eskes, H.; Westerink, J.; Sawatzky, G.A.; Czyzyk, M.T. Electronic structure of Cu₂O and CuO. Phys. Rev. B; 1998; 38, pp. 11322-11330. [DOI: https://dx.doi.org/10.1103/PhysRevB.38.11322] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9946011]

4. Heinemann, M.; Eifert, B.; Heiliger, C. Band structure and phase stability of the copper oxides Cu₂O, CuO, and Cu₄O₃. Phys. Rev. B; 2013; 87, 115111. [DOI: https://dx.doi.org/10.1103/PhysRevB.87.115111]

5. Pickett, W.E. Electronic structure of the high-temperature oxide superconductors. Rev. Mod. Phys.; 1989; 61, pp. 433-512. [DOI: https://dx.doi.org/10.1103/RevModPhys.61.433]

6. Reitz, J.B.; Solomon, E.I. Propylene Oxidation on Copper Oxide Surfaces: Electronic and Geometric Contributions to Reactivity and Selectivity. J. Am. Chem. Soc.; 1998; 120, pp. 11467-11478. [DOI: https://dx.doi.org/10.1021/ja981579s]

7. Nakaoka, K.; Ueyama, J.; Ogura, K. Photoelectrochemical Behavior of Electrodeposited CuO and Cu₂O Thin Films on Conducting Substrates. J. Electrochem. Soc.; 2004; 151, C661. [DOI: https://dx.doi.org/10.1149/1.1789155]

8. Kim, S.Y.; Ahn, C.H.; Lee, J.H.; Kwon, Y.H.; Hwang, S.; Lee, J.Y.; Cho, H.K. p-Channel Oxide Thin Film Transistors Using Solution-Processed Copper Oxide. ACS Appl. Mater. Interfaces; 2013; 5, pp. 2417-2421. [DOI: https://dx.doi.org/10.1021/am302251s]

9. Åsbrink, S.; Norrby, L.J. A refinement of the crystal structure of copper(II) oxide with a discussion of some exceptional e.s.d.’s. Acta Crystallogr. Sect. B; 1970; 26, pp. 8-15. [DOI: https://dx.doi.org/10.1107/S0567740870001838]

10. Schmahl, N.G.; Eikerling, G.F. Über Kryptomodifikationen des Cu(II)-Oxids. Z. Für Phys. Chem.; 1968; 62, pp. 268-279. [DOI: https://dx.doi.org/10.1524/zpch.1968.62.5_6.268]

11. Ekuma, C.; Anisimov, V.; Moreno, J.; Jarrell, M. Electronic structure and spectra of CuO. Eur. Phys. J. B; 2014; 87, 23. [DOI: https://dx.doi.org/10.1140/epjb/e2013-40949-5]

12. Wu, D.; Zhang, Q.; Tao, M. LSDA+U study of cupric oxide: Electronic structure and native point defects. Phys. Rev. B; 2006; 73, 235206. [DOI: https://dx.doi.org/10.1103/PhysRevB.73.235206]

13. Nolan, M.; Elliott, S.D. The p-type conduction mechanism in Cu₂O: A first principles study. Phys. Chem. Chem. Phys.; 2006; 8, pp. 5350-5358. [DOI: https://dx.doi.org/10.1039/b611969g] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19810413]

14. Cao, H.; Zhou, Z.; Yu, J.; Zhou, X. DFT study on structural, electronic, and optical properties of cubic and monoclinic CuO. J. Comput. Electron.; 2018; 17, pp. 21-28. [DOI: https://dx.doi.org/10.1007/s10825-017-1057-9]

15. Grant, P.M. Electronic properties of rocksalt copper monoxide: A proxy structure for high temperature superconductivity. J. Phys. Conf. Ser.; 2008; 129, 012042. [DOI: https://dx.doi.org/10.1088/1742-6596/129/1/012042]

16. Cipriano, L.A.; Di Liberto, G.; Tosoni, S.; Pacchioni, G. Band Gap in Magnetic Insulators from a Charge Transition Level Approach. J. Chem. Theory Comput.; 2020; 16, pp. 3786-3798. [DOI: https://dx.doi.org/10.1021/acs.jctc.0c00134] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32427487]

17. Yazdani, A.; Barakati, B. A first-principles study on electronic structure and crystal field effect of layered La₂CuO₄ as composed of CuO₂ and La₂O₂ monolayers. Phys. E Low-Dimens. Syst. Nanostructures; 2021; 125, 114395. [DOI: https://dx.doi.org/10.1016/j.physe.2020.114395]

18. Yin, K.; Zhang, Y.Y.; Zhou, Y.; Sun, L.; Chisholm, M.F.; Pantelides, S.T.; Zhou, W. Unsupported single-atom-thick copper oxide monolayers. 2D Mater.; 2016; 4, 011001. [DOI: https://dx.doi.org/10.1088/2053-1583/4/1/011001]

19. Kano, E.; Kvashnin, D.G.; Sakai, S.; Chernozatonskii, L.A.; Sorokin, P.B.; Hashimoto, A.; Takeguchi, M. One-atom-thick 2D copper oxide clusters on graphene. Nanoscale; 2017; 9, pp. 3980-3985. [DOI: https://dx.doi.org/10.1039/C6NR06874J]

20. Kvashnin, D.G.; Kvashnin, A.G.; Kano, E.; Hashimoto, A.; Takeguchi, M.; Naramoto, H.; Sakai, S.; Sorokin, P.B. Two-Dimensional CuO Inside the Supportive Bilayer Graphene Matrix. J. Phys. Chem. C; 2019; 123, pp. 17459-17465. [DOI: https://dx.doi.org/10.1021/acs.jpcc.9b05353]

21. Blaha, P.; Schwarz, K.; Tran, F.; Laskowski, R.; Madsen, G.K.H.; Marks, L.D. WIEN2k: An APW+lo program for calculating the properties of solids. J. Chem. Phys.; 2020; 152, 074101. [DOI: https://dx.doi.org/10.1063/1.5143061] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32087668]

22. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett.; 1996; 77, pp. 3865-3868. [DOI: https://dx.doi.org/10.1103/PhysRevLett.77.3865] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10062328]

23. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B; 1976; 13, 5188. [DOI: https://dx.doi.org/10.1103/PhysRevB.13.5188]

24. Pizzi, G.; Vitale, V.; Arita, R.; Blügel, S.; Freimuth, F.; Géranton, G.; Gibertini, M.; Gresch, D.; Johnson, C.; Koretsune, T. et al. Wannier90 as a community code: New features and applications. J. Phys. Condens. Matter.; 2020; 32, 165902. [DOI: https://dx.doi.org/10.1088/1361-648X/ab51ff] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31658458]