Ferromagnetic topological crystalline insulating phase in the π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi$$\end{document}-stacked graphene nanobelts under a small pressure

This work presents a topological material based on the graphene nanobelts. It is intrinsically magnetic. The band inversion in this system is partial and anomalous, as it was recently found in the Au-doped graphene and 2D molecular crystal. The metallic edge states with the Chern number equal to 1 exist in the ferro- and antiferromagnetic phase. The edge states become gapless at a small pressure, around 0.87 GPa, applied along the stacking axis. Their dispersions are not symmetric for the spin up and down. The graphene nanobelts are π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi$$\end{document}-stacked with the AB-type, and the zigzag edges are terminated with the O and H atoms at the opposite sides, in the reversed order along the stacking axis. The propagation of the topological edge states is parallel to the O⋯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\cdots$$\end{document}H-interaction direction, across the stacked belts.


Introduction
Graphene has shown many appealing properties, since its discovery in 2004 and the Nobel Prize for Novoselov and Geim in 2010 [1,2]. Due to the existence of the Diraccone states, the topological insulating (TI) phase has been searched there. However, without success, because the spin-orbit coupling (SOC) in this material is negligible. On the other hand, a light doping of graphene with the hydrogen atoms strengthens the SOC by four orders of magnitude [3,4]. Recently, the status of the search for the TI phase in graphene has changed with a prediction of the TCI phase in the twisted graphene bilayers [5] with the negligible SOC effects, both the intrinsic and Rashba type. The SOC-free TCI phase was found also at the domain-wall of the AB-and BA-stacked graphene bilayers [6][7][8] and in the 2D molecular crystal [9].
Here, the TCI phase is found in the -stacked magnetic graphene nanoribbons at a small pressure of about 0.87 GPa, at which the metallic edge states become gapless. The belts are -stacked in the AB-type order, the same as in graphite. The zigzag edges are terminated with the oxygen and hydrogen atoms on the opposite sides, with an alternation of the O and H atoms in the -stacking direction-here called the AB-reversed order. It is plausible, that the vertical O ⋯ H-interactions play a role in the formation of the topological states, where the mechanism would be the same as in the stacked 2D molecular system [9].
The band inversion is essential for the topological phase, but not sufficient. In the conventional TIs, the time-reversal properties are responsible for the BI [10]. In contrast, the TCIs are ruled by the crystalline symmetry, the SOC should be in principle absent there [11,12]. These materials are supposed to be built of the light elements only [11]. Very recently, the new type of the BI-namely the anomalous BI-was reported in graphene with the substitutional N 3 -centers additionally doped with the Au atoms [13]. The anomalous character of the BI in the above system shows up as the inversion which occurs between the conduction band minimum (CBM) and the bands laying deeper below the valence band top (VBT). Our previous work, on the other hand, reports the partial BI in the 2D molecular crystal [9]. It will be demonstrated in this work, that a combination of both types of the BI, the partial and anomalous between the VBT and the bands laying above the CBM, occurs in the -stacked graphene nanobelts.
The magnetism in the graphene nanobelts originates from the edge termination. It is known for a few years that a saturation of the zigzag edges with the double and single bonds on the opposite sides leads to the formation of the magnetic moments at the edges [14,15]. We proposed to use this effect in the random-access memory devices [16]. Here, the saturation of the zigzag edges with O and H is responsible for two effects: magnetism and the topological phase. The second phenomenon occurs when the interlayer distance is slightly decreased. We think that these systems are possible to synthesize, since the chemical manipulation at graphene edge is nowadays an easy task [17].
To date, a very few TIs among native magnetic materials, such as undoped crystals, are known. One of them is the antiferromagnetic Heusler compound [18]. Usually, the conventional TIs are doped with magnetic elements. Some examples of the transition metals are V, Cr, Mn, Fe and Co in chalcogenides [19][20][21][22][23][24][25][26][27][28]. It is desired to know whether and how the doping might change the topological states. It was believed that the nonmagnetic doping does not influence the topological properties [29]. In contrast, the effect of magnetic impurities and magnetic field opens the gap at the Dirac cone of the surface states in some 3D TIs [21][22][23], challenging the robustness of the topological phase. While some other TIs with the nonmagnetic dopants pass the above test [24][25][26]. Recently, it has been shown that the nonmagnetic impurities such as indium in Bi 2 Se 3 also open the surface band gaps [28], in contrast to previous belief [29]. Thus, a discussion on a role of the magnetic anisotropy of the add atoms at the surface and in the bulk, and the position of the impurity band, for a protection of the topological surface states, is still actual [19,27]. Moreover, the topological properties are correlated with the mechanisms of the long-range magnetic order, and this fact influences the Curie temperature [20,27]. Since the impurities might harm the topological properties, and the spin-polarized states at the surfaces usually do not coexist with the bulk magnetism, a generation of the spin-polarized current through the bulk of the TI was proposed [30].
In view of the controversy on the topological protection and doping, and in contrast with the doped TIs, the intrinsically magnetic and SOC-free system with the asymmetric metallic surface states are presented. These states become gapless under a small pressure applied along the stacking axis. In the infinite 2D system, the band gap lies above the two bands above the Fermi level. In the terminated system, the magnetic edge states could be accessed with the electron excitation via the absorption of light. The proposed systems pave the way for a design of naturally magnetic TCIs based on graphene and the hydrogen interactions. In the antiferromagnetic (AF) phase, this system is similar to the 3D TCIs composed of the AF stacks of the 2D layers, predicted by the symmetry considerations [31,32]. However, with the difference that the dimensions of the systems presented here are lower. The ferromagnetic (FM) TCIs, on the other hand, were to date only predicted in the spinel structure materials belonging to the class of the filling-enforced semimetals [33].

The atomic, electronic, and magnetic structure
The atomic structure and the symmetry of the proposed system are presented in Fig. 1a-c. The side and top views at the AB-stacked graphene nanobelts with the alternated (named reversed) O and H terminations are drawn in Fig. 1a, b, respectively. This 2D nanomaterial contains two -stacked armchair ribbons in the elementary cell, and transforms with the inversion operation, as shown in Fig. 1c. The optimized separations between the belts do not differ much for the FM and AF phases and equal to 3.439 Åand 3.437 Å, respectively.
Let us further investigate the symmetry of these systems. The magnetic symmetry groups (MSG) are of the types I-IV, where the type I is equivalent to the nonmagnetic symmetry groups and can be described solely by the unitary symmetry operations. Other types, II-IV, have an equal number of the unitary and anti-unitary elements: M = G + A . The anti-unitary part A can be written as a product of the spacial and time-reversal operations Θ . If belongs to G, then our system is of the type II. In the opposite case ( does not belong to G), our system is of type IV (called sometimes IIIb) when is a pure translation. In case is not a pure translation, our system is of type III (called sometimes IIIa). For a reference see for example works [32,33]. The atoms of our elementary cell transform with the superposition of two operations: the rotation by 180 • around the stacking (z-axis) and the fractional translation along this axis (by 1 2 of the lattice constant c). Equivalently, one can say that there is a glide mirror symmetry (pmg2). The systems studied here can be classified as 3.4 MSG IV ( P a 112 ) [33].
The band structures, calculated with the density functional theory (DFT), with the SOC included, are presented for the FM and AF cases in Fig. 1d-g. The Fermi level, positioned at the zero energy, is around 0.7-0.9 eV below the VBT. Both magnetic phases are metallic, and they change the character to semiconducting after the n-type doping with one electron per each oxygen atom (two electrons per the elementary cell). The FM phase is energetically more favorable than the AF phase, by 31 meV per the elementary cell.
In Fig. 1d, e, the colors indicate a contribution of the oxygen states to the Bloch functions. The oxygen component predominates the valence band along the stacking axis, while the bands along the belt axis are rather built by carbons. The direct bandgap in the FM phase is 0.18 eV and located at the k-point on the line between the Γ and X symmetry points. In contrast, the band gap in the AF phase is indirect and equals to 0.2 eV. Its location is between the VBT at the Γ point and the CBM at the same k-point as in the FM phase.
The absolute magnetization, per one oxygen in the elementary cell, is 0.55 B for the FM phase and 0.28 B for the AF phase. The projections of the band structures onto the majority and minority spin-component are in Fig. 1f, g. The absence of the states which mix the large and small spinor components, in the whole Brillouin zone, is a signature of the SOC-free character of the band structures. The magnetic moments are mainly localized at the zigzag edges, as seen from the spin maps for the FM and AF phases in Figure S1a-j in the supporting information (SI). In the AF phase, the local spin moments sit mainly at the oxygen atoms, i.e., the magnetization is localized in the planes of the graphene belts. The local moments in the FM phase are similar to those in the AF phase when comparing the belts planes. Additionally, the FM phase possesses more rich spin structure, which is shifted away from the belts planes above and below the carbons adjacent to hydrogens, as presented in Figure S1c-d and S1e-g in SI.
Interestingly, the systems presented in this work have a very nice property which is desired for the photovoltaics. Namely, the carbon rings terminated with the atomic groups which contain the O ⋯ H or N ⋯ H-interactions possess the property of the space separation of the charge carriers [34]. In other words, the electrons and hole move along different paths through the material. Such separation reduces the unwanted charge recombination of carriers, which is one of the key factors lowering the conversion efficiency of the solar cells. The atomic structures (a-c) and the band structures (d-g) of the graphene nanobelts. The FM or AF case is denoted above each electronic structure. The projections at the oxygen-localized Wannier functions d, e are in color legends. If 100% of the band was built of the oxygen-localized functions then the band color would be red (the component factor was 1). The projection at spin-z component f, g is also in color legends. The absence of the green color in f, g indicates the negligible SOC Similar property has been suggested to appear in the organic lead halide perovskites [35,36] and covalent organic frameworks [37]. It is a nice coincidence that this effect occurs also in the O= and H-terminated graphene nanobelts. This fact is clear from Figures S2 and  S3 in SI, which present the maps of the Bloch's functions at the VBT and CBM for the FM and AF phase, respectively.

The topological properties
We assume the 2D system to be finite along the belt axis (the Y Cartesian axis). The 2D band structures and the 1D edge states of the systems truncated along the belts (the X Cartesian axis) are analyzed in Fig. 2. In particular, the dispersions are drawn in Fig. 2a, b for the FM phase and in Fig. 2d, e for the AF phase. The propagation direction of the edge states is along the stacking axis (the Z Cartesian axis), which is also the direction of the O ⋯ H-interactions. For the FM phase, two spin-polarized and not degenerate In the AF phase, the state within the gap which appears to be single is in fact doubly degenerate for the two spin components. Both magnetic orders are characterized with the parabolic dispersions of the edge states, falling within the band gap of the 1D structure. This property could be one of the prerequisites of the topological crystalline insulating phase, under the condition that the edge states become gapless under a small pressure. Such effect occurs in the studied systems, as will be clear later. Figure 2c, f displays the spacial localization of the wave functions of the two edge states, in the FM and AF phase. These states are well localized at the cut edges of the belts. In the AF case, both edge states localize at the two ends of the same belt. This is, in fact, every second belt of the stack, because two belts are in the elementary cell. In the FM case, the edge states on the left-and right-hand side of the belt localize at the alternating belts of the stack. Moreover, the wave function of the edge state which is more occupied by the electrons is less localized. It extends more to the interior of the belt, in comparison with the less occupied edge state (the one which is higher in the energy).
Let us switch to the band inversion. It is almost absent for both magnetic phases, if one looks only at Fig. 1d-g. However, an enlargement of the color scale reveals more subtle effects. The detailed characteristics is presented in Fig. 3. The contribution of the oxygen-localized states to the Bloch functions, which is printed up to 10% , shows a signature of the band inversion. The two bands of the CBM are composed of the O-states in much smaller part than the VBT states. The two states just above the CBM contain much more substantial portion of the oxygen states. This type of the BI is similar to the earlier published two cases: The first type of the BI, called anomalous, was reported for the Au-doped N 3 -centers in graphene, where the CBM exchanged its character with the bands laying deeper below the VBT [13]. The second type of the BI, called partial, was reported by us for the 2D molecular crystal, where the bands inversion was mixed with the non-inversed components [9]. It is worth to notice that these types of BI are numerically observed in the three reported cases which are the -stacked systems.
The maps of the square-root function of the Bloch densities-equivalent to the absolute values of the Bloch functions-are presented in Figures S2 and S3 in SI for the FM and AF phases, respectively. The VBT Blochs are mainly composed of the states which are localized within the planes of the belts, close to the O= terminated edge. In contrast, the CBM Bloch functions are composed of all carbon states which are distant from the belt edges, and the maxima of these states are shifted a bit away from the belt's plane. This shift is a signature of the p z -orbital character of the CBM.
Since our systems are magnetic, there is a question on the proper topological invariant, i.e., whether one could apply the Z 2 invariant [38]. In fact, the band gaps are positioned above the Fermi levels of the magnetic phases. One needs to occupy two more bands in order to move the Fermi level inside the band gap. However, if all states up to this gap are occupied, the magnetism will vanish. The only way to access the magnetic topological edge states would be the excitation of electrons via the absorption of a photon. Therefore, when analyzing the symmetry of the entangled valence band manifold, these unoccupied two states must be included. This way, we arrive to the same symmetry case as of the nonmagnetic system, which can be described by the Z 2 invariant.
The parities of the Bloch's functions, with respect to the inversion symmetry operation (or equivalently a superposition of C 2 and fractional translation), were examined for the manifold of the entangled valence states-in a number of 70 bands per spinor-printed in Figure S4 in SI. The product ( ) of the parities ( i ) within this manifold is given in Table 1 for the Γ-point and the three TRIM points ( TRIM = time-reversal invariant momentum ): X = (1∕2, 0, 0) , Z = (0, 0, 1∕2) and M = (1∕2, 0, 1∕2) . The parities of the two bands at the VBT and CBM are also included in Table 1. The list of the characters of the magnetic double point group C 2v (2) is in Table S1 in SI. The irreducible representations assigned to the considered states are enclosed in SI in Tables S2-S5 for the ferromagnetic phase and Tables S6-S9 for the antiferromagnetic phase. In conclusion, both systems presented in this work are characterized by the Chern number [38] equal to 1.
The pressure evolution of the 2D band structures and the corresponding 1D edge states is presented for the FM phase in Fig. 4. The larger overlap of the p z orbitals, with the applied pressure, causes three effects: (1) the bandgaps become smaller, (2) the local magnetic  Table 1 The product of the parities of the Bloch functions belonging to the manifold of the entangled valence bands at four high symmetry points, counted per large spinor component The parities of the two bands at VBT and CBM are denoted VBT (1,2) and CBM (1,2) . The Z 2 invariant indicates the topological phase in 2D. It is given as a product of all 's, which is The pressure evolution of the 2D band structure and the 1D truncated band structure, for the FM phase. The vertical distances d z between the nanobelts, as well as the pressures, are given above the corresponding panels. The valence band maximum of the 1D truncated structure is marked with the green dashed line. The edge states are marked with red. The number of cells, N, in the truncated finite direction (X-axis) is given above the 1D structures. The last row presents enlarged regions close to the Γ-point. The energy scale is common for all three choices of N states move not monotonically with both the pressure and the thickness of the truncated structure. The AF magnetic topological crystalline phase could be viewed as an example of the TCIs earlier predicted from the purely symmetrical considerations [31,32]. The Θ model [32] (with the glide mirror symmetry, as in our case) was discussed for a layered structure of the antiferromagnetically coupled planes, where the spins are perpendicular to the stacked planes. Such 3D structure was shown to possess the gapless edge states protected by the Π symmetry defined as a matrix of C 4 operations of two magnetic sublayers. In the AF phase studied here, this kind of stacking could be achieved on top of the 2D structures, if the belts are periodically repeated along the previously nonperiodic direction, i.e., the Y-axis. Building such 3D structure is plausible, due to the existence of the hydrogen bonds between the alternated O= and H-zigzag terminations.

Discussion
There is a quest for robust magnetic topological insulators for spintronic applications. The ways could be twofold: either to search over the intrinsically magnetic systems like Heusler compounds [18,39], or via doping of the TIs with magnetic elements [19][20][21][22][23][24][25][26][27][28]. Doping of the TIs with the transition metals might lead to an opening of the gap in the topological surface states. But this is not a rule. It has been postulated that the magnetic anisotropy plays a role in that effect [27]. Namely, in 3D TIs, the out-ofsurface magnetic polarization was believed to open the gap, as in the case of Mn-doped Bi 2 Te 3 . While the magnetic moments parallel to the surface were believed not to harm the surface states, as for the Co doping. With respect to the magnetic anisotropy, the doping at the surface or in the bulk usually makes a difference [27]. Recently, the above rule about the anisotropy and the gapless states has been contradicted, and the independence of the topological protection on magnetization was evidenced [28].
The 2D system presented in this work is intrinsically magnetic, with the magnetic moments localized at the plane spanned by the stacking axis and the zigzag edge of the belt. The belt has its thickness. In our case, it is sixteen carbons along the armchair chain of graphene, terminated with O and H. The belts are stacked with the AB-type reversed order. Such that each oxygen is located vertically between the carbons adjacent to the hydrogens of the belt below and belt above. This is important from the point of view of the magnetic interactions. The FM phase, with respect to the coupling between the belts in the stack, is energetically more favorable, by 31 meV, than the AF phase. The magnetic coupling occurs despite a quite large distance between the stacked layers, and the oxygen atoms are vertically separated by the two stacking distances, around 6.88 Å. However, the spin maps of the FM phase show the spindensity shifted from the belts plane close to the carbons adjacent to hydrogens. If we focus on the magnetic anisotropy, all O-centered moments are oriented out of the 2D surface. While the carbons along the armchair host much weaker magnetic moments than those of oxygens. The carbon moments are in the alternated spin order across the belt. This is a kind of a spin wave, schematically presented in Fig. 5a and in detail in Figure S5a in SI. After stacking the belts, the moments at the vertically overlapping carbons are parallel to each other. Therefore, weakly antiferromagnetic exchange across the belt Fig. 5 The schematic view at: a the magnetic structure of the stacked nanobelts (the red and green arrows represent the positions of the local spins up and down), b the edge states propagation between the belts (the O ⋯ H-interaction direction, which is along the stacking Z-axis), c the edge states propagation in the 2D molecular crystal from Ref. [9] (also the O ⋯ H-interaction direction, which is one of the planar axes in this case) drives the ferromagnetic exchange between the belts and also between the vertically coupled oxygens of every second belt. Moreover, the triangular symmetry of the system disfavors the AF order between the oxygens along the zigzag edge. If the AF order along the zigzag was existent, then the magnetic system would be frustrated, as shown in Figure S5b in SI. But it does not occur. The DFT calculations initialized at the AF-zigzag configuration converge to the spin unpolarized result, with the total energy higher than the FM phase, by 43 meV. The long-range ferromagnetic order in these belts might be robust even above the room temperature. This expectation is based on the similarities with the systems without the intentional O= and H-saturation [40][41][42].
The topological edge states form when one cuts through the belts, and not between the belts. Hence, the edge states propagate between the stacked layers. On the other hand, the local magnetic moments are parallel to the topological edge. Thus, the graphene belts would confirm an observation about the correlation between magnetic anisotropy and opening or not of the gap within the surface states [27]. This fact, however, is coincidental because our system should not be viewed as the TI material with the impurity states. This is a system where an intrinsic magnetism coexists with the crystal topological phase at a small pressure.
It is interesting to note that the propagation direction of the edge states along the stacking axis is parallel to the O ⋯ H-interaction direction. Whereas the bandgap location and the band inversion-although only partial-are between the states propagating along the belt axis. Similar situation was in the case of the 2D molecular crystal, namely in the -stacked phenalene derivatives planarly interacting via the hydrogen bonds [9]. The edge states in that system also propagated along the O ⋯ H-interaction direction, whereas the bandgap and band inversion were between the states propagating along the stacking direction. The aforementioned systems and their topological directions are schematically drawn in Fig. 5b, c, for these and the earlier reported TCIs, respectively. It is remarkable, that the carriers at the topological states move between the weakly interacting parts of the system, and not through the covalent bonds, for both systems. Additionally, not the -stacking axis, but the hydrogen bonds pin the propagation direction of the edge states in these numerically just-discovered TCIs.
From the symmetry-group point of view, the presented 2D systems are similar to the antiferromagnetic 3D TCIs constructed of the stacked 2D layers [31,32], while the ferromagnetic case might be similar to the filling-enforced semimetals [33]. However, with the difference that the access to the edge states in our case is via the excitation of electrons, because the band gap is almost 1 eV above the Fermi level. Our systems are metallic when they are charge undoped.
In summary, the intrinsically magnetic, the FM and AF coupled, 2D TCIs are predicted under a small pressure in the AB-type -stacked O= and H-terminated zigzag graphene nanobelts. The edge states dispersions are nondegenerate for the opposite spins. This is the second system just found [9] where the structures based on the benzene rings, -stacking and O ⋯ H-interactions are characterized with the metallic edge states. These states become gapless under a small pressure applied along the stacking axis. It would be desirable to take these and other hydrogeninteracting systems under a magnifying glass and conduct the experimental investigations.

Methods
The density functional theory [43,44] calculations were performed with the Quantum ESPRESSO (QE) package [45], which is based on the plane waves and the pseudopotentials for atomic cores. The norm-conserving pseudopotentials were used, and the energy cutoff of 60 Ry for the plane waves was set. The exchange-correlation functional was gradient corrected with the Perdew-Burke-Ernzerhof parametrization [46]. The uniform k-mesh by Monkhorst and Pack [47] with the 5-point grid along the stack and 10-point grid along the belt axis was sufficient to converge the band structure. The vacuum space of 20 Å was set between the 2D slabs in the nonperiodic direction. The geometry optimization was performed with the Broyden-Fletcher-Goldfarb-Shanno algorithm [48] until forces were smaller than 1 meV∕Å.
For an accurate interpolation of the band structures, the wannier90 package [49] was used. It enables finding the maximally localized Wannier functions for the composite bands [50,51]. The tight-binding model for the edge states analysis was constructed with the DFT Bloch functions. They were used for the generation of the Wannier functions [50]. Then, the Kohn-Sham Hamiltonian [43,44] between these well-localized functions was built. In the tight-binding simulations with the pythTB code [52], the cutoff for the nonvanishing matrix elements was set to 0.005 eV. The parities of the Bloch functions at the TRIM points were obtained with the bands.x utility from the QE package.