# Zigzag nanoribbons of two-dimensional hexagonal crystals:

magnetic, topological and thermoelectric properties

###### Abstract

Effects of electron-electron and spin-orbit interactions on the ground-state magnetic configuration and on the corresponding thermoelectric and spin thermoelectric properties in zigzag nanoribbons of two-dimensional hexagonal crystals are analyzed theoretically. Thermoelectric properties of quasi-stable magnetic states are also considered. Of particular interest is the influence of Coulomb and spin-orbit interactions on the topological edge states and on the transition between the topological insulator and conventional gap insulator states. It is shown that the interplay of both interactions has also a significant impact on transport and thermoelectric characteristics of the nanoribbons. The spin-orbit interaction additionally determines an in-plane magnetic easy axis. Thermoelectric properties of the nanoribbons with in-plane magnetic moments are compared with those of the nanoribbons with edge magnetic moments oriented perpendicularly to their plane. Nanoribbons with ferromagnetic alignment of the edge moments are shown to reveal spin thermoelectricity, in addition to the conventional one.

###### pacs:

73.63.-b,75.75.-c,73.22.Pr,71.30.+h## I Introduction

There is currently a great interest in strictly two-dimensional (2D) crystalline materials, like graphene,Novoselov2004 () silicene, Aufray2010 (); Kara2012 () germanene, Bishnoi2013 () and others. A common feature of these materials is a 2D hexagonal lattice of atoms – carbon, silicon, and germanium, respectively. However, both silicene and germanene have buckled atomic structure, where the two triangular sublattices are slightly displaced along the normal to the atomic plane, whereas all carbon atoms in graphene are localized in a common plane. This buckling of atomic structure leads to new properties of silicene and germanene, Ezawa2012 (); Ezawa2013a (); Ezawa2013b (); Lado2014 () which make them different from graphene.

Electronic states around the charge neutrality points of these materials are usually described by the relativistic Dirac model.Neto2009 (); Katsnelson2012 () This model, however, neglects the electron-electron interaction. Moreover, some deviations of the exact electronic structure from the conical dispersion relation characteristic of the Dirac model appear for states distant from the Dirac points. In other words, the exact electronic structure is more complex, and either tight binding models or ab-inito techniques are usually employed to calculate the corresponding band structure, especially in the case of nanoribbons.Lopez2013 (); Zhang2010 (); Lou2009 (); Lou2011 (); Guan2012 (); Lou2013 (); Chen2014 (); Ding2013 () It is already well known that the Coulomb interaction between electrons leads to edge magnetism in nanoribbons with zigzag edges.Yazyev2010 () In the ground state of the system, magnetic moments localized at one edge are opposite to those at the other edge.Yazyev2010 (); Son2006 (); Pisani2007 () In some situations, however, one can stabilize a quasi-stable ferromagnetic state, where all edge moments are oriented in parallel.Haugen2008 (); Kim2008 (); Zheng2012 (); Li2014 (); Wang2014 () Spin-orbit interaction gives rise to a magnetic anisotropy, and when this coupling is sufficiently strong, an in-plane orientation of the edge moments is energetically more favorable.Lado2014 () Thus, the interplay of spin-orbit and Coulomb interactions is of fundamental importance for the edge magnetism.

The spin-orbit interaction also leads to topologically protected edge states in the spin-orbit gap,Kane2005 () and these states are responsible for topological-insulator properties of these materials.Hasan2010 (); Qi2011 () However, Coulomb interaction has a significant impact on these properties. Accordingly, the interplay of spin-orbit and Coulomb interactions is important not only for the edge magnetism, but also for the topological state of these materials. Therefore, in this paper we consider this issue in the case of zigzag nanoribbons of two-dimensional buckled hexagonal crystals, and especially the role of Coulomb interaction and spin-orbit coupling. We focus on the thermoelectric and spin thermoelectric effects of these materials in the ballistic transport regime, and on their dependence on the Coulomb and spin-orbit interactions. First, we calculate energy of various magnetic configurations in zigzag nanoribbons, and find that the antiparallel configuration with in-plane magnetic moments is the ground state structure for large enough spin-orbit coupling. For small spin-orbit coupling, the in-plane and perpendicular-to-plane configurations are degenerate. Then, we calculate the electronic states, with a special emphasis on the edge states responsible for the topological properties. Having determined the magnetic state and also the electronic edge states, we calculate thermoelectric properties with the main focus on the role of topological edge states and transition from the topological to conventional insulator state.

Generally, the thermoelectric and especially the spin thermoelectric effects in systems of reduced dimensionality are currently of great interest,Hochbaum2008 (); Harman2002 (); Duarte2009 (); Walter2011 (); Liebing2011 (); Uchida2008 (); Swirkowicz2009 (); Misiorny2014 () mainly due to a hope to find an efficient way to convert dissipated heat into electrical energy. In fact, thermoelectric properties of nanoribbons have been studied theoretically in a couple of papers.Pan2012 (); Zberecki2013 (); Yan2013 (); Zberecki2014 (); Zberecki2014a () However, the role of Coulomb interaction and topological edge states in the gap, has not been studied thoroughly enough yet. In a recent paperWierzbicki2015 () we have analyzed the influence of topological states on the thermoelectric properties, with the main focus on the role of a staggered exchange field Ezawa2012 (); Ezawa2013a (); Ezawa2013b () and of an electric field perpendicular to the atomic plane. Here we consider in more details the role of Coulomb interaction, which has a significant impact on the topological properties.Matthes2014 ()

In section 2 we present the Hamiltonian used to describe the 2D materials under consideration, and also introduce briefly the method used to calculate transmission function and then all the thermoelectric and transport coefficients. In section 3 we present numerical results on the edge magnetism, topologically protected states, and thermoelectric properties of the corresponding zigzag nanoribbons. Summary and final conclusions are in section 4.

## Ii Theoretical description

### ii.1 Tight-binding Hamiltonian

The zigzag nanoribbons of hexagonal 2D crystals, like silicene or germanene, which are analyzed in this paper consist of zigzag atomic chains. The system can be described by Hamiltonian of the following general form:

(1) |

where the first term, , is the tight-binding Hamiltonian, the term represents the spin-orbit coupling, while the last term, , stands for the electron-electron Coulomb interaction.

The tight-binding Hamiltonian is assumed in the form

(2) |

where () is the creation (annihilation) operator for an electron with spin at the lattice point (). The spin quantization axis is taken along the -axis, which is normal to the atomic plane. The tight-binding term is restricted here to electron hopping between nearest-neighbor sites, with the corresponding hopping parameter . For instance, in silicene.Ezawa2012 ()

The intrinsic spin-orbit interaction in Eq. (1) follows from electron hopping between next-nearest-neighbors in the hexagonal lattice, and can be written in the form,Ezawa2013b ()

(3) |

where is the spin-orbit interaction parameter, while () when the hopping path from second-neighbor sites to in the hexagonal lattice is clockwise (anticlockwise) with respect to the positive -axis.

The electron-electron interaction will be taken into account in the Hubbard form within the mean-field approximation (MFA). The form of the MFA Hubbard Hamiltonian depends on whether the spin quantization axis is parallel to the local magnetization orientation or not. We assume, similarly as in Ref.[Lado2014, ], that in a hypothetical state the magnetic moments are tilted from the normal orientation by an angle , and will calculate the corresponding total energy. The ground state configuration will be thus the one with the lowest energy. Assuming the quantization axis parallel to the local magnetic moments, the Hubbard term in the mean-field approximation can be written in the form

(4) |

where and are the particle number operators for a site and for electrons with spin and with respect to the spin quantization axis tilted from the normal to the plane towards a nanoribbon edge by an angle . In turn, is the on-site Coulomb repulsion parameter, . We also note, that the Coulomb interaction of electrons at different lattice sites is omitted.

The MFA Hubbard Hamiltonian should be now written in the quantization axis common for all terms of the Hamiltonian. This can be done by rotating the quantization axis back to the orientation normal to the plane, which is achieved by the transformationSwirkowicz1997 ()

(5) |

where () for ().

The spin-dependent mean number of electrons, , at the site , required to evaluate the Hubbard interaction term (4), is determined self-consistently from the band-structure of the corresponding nanoribbon,

(6) |

where is the Fermi-Dirac distribution function, with denoting the chemical potential, while is the spin-dependent electron density at site for a Bloch-wave of energy . The parameter denotes the length of the elementary cell of a zigzag nanoribbon. The chemical potential of a pristine nanoribbon is equal to the corresponding Fermi energy . By doping with donor or acceptor impurities, or by applying an external gate voltage, one can shift up or down the Fermi level, and its position with respect to will be described by . In the following and will denote spin projection on the quantization axis parallel to the orientation of the magnetic moments.

To identify the edge-states, we evaluate the mean value of , where the -axis is in the nanoribbon plane and normal to the nanoribbon axis, is the -coordinate of the nanoribbon center, while is the width of the nanoribbon. For (), the states are localized at the left (right) edge of the nanoribbon.

### ii.2 Electronic transmission and general formulas for thermoelectric coefficients

Ballistic electronic transmission in quasi-one-dimensional structures is usually determined in terms of the Green’s functions formalism and an expression derived by Caroli et al.Caroli71 () An alternative approach is based on matching wave functions in the scattering region to the Bloch modes of ideal bulk leads.Khomyakov2005 () The total transmission, , as a function of electron energy is then defined as the sum of transmissions for all incident Bloch waves for a given wave-number , and is a piece-wise constant function equal to the number of energy bands for a given energy . In turn, the spin transmission, , is equal to the sum of the mean values of the Pauli operator . More detailed description of the transmission calculation is presented in Ref. [Wierzbicki2015, ]. The total transmission, , and spin transmission, , are related to the spin-dependent transmission, (), as and , respectively.

Knowledge of the transmission function allows to calculate all transport characteristics in the linear response regime, including also the thermoelectric properties. The spin resolved transport and thermoelectric coefficients can be then determined from the moments of the electronic current series expansion around the Fermi level,

(7) |

For a piece-wise constant transmission function, calculation of can be reduced to calculation of the so-called incomplete Fermi-Dirac integral,Wierzbicki2015 ()

(8) |

for which there exists an efficient numerical procedure.goano ()

The spin-dependent conductance is proportional to the moment , where is the electron charge. The total conductance is then equal to , and the spin conductance is . The heat conductance is defined at vanishing charge current and can be written as , where stands for temperature. In turn, the spin-dependent thermoelectric coefficient (thermopower) can be expressed in the form

(9) |

for . Then, one can define the conventional (charge) Seebeck coefficient (thermopower) as and also the corresponding spin Seebeck coefficient (spin thermopower) as .Swirkowicz2009 () When the spin thermopower is absent, one can write the conventional Seebeck coefficient () in the standard form , with for . In turn, the electronic contribution to the thermal conductance is given by .

Finally, one can define the dimensionless figures of merit for the conventional and spin thermoelectric effects,

(10a) | |||

(10b) |

where is the phonon contribution to the heat conductance. In the absence of spin thermoelectricity, only the dimensionless figure of merit for the conventional thermoelectricity is relevant (), and ,

## Iii Numerical results

In this section we present some numerical results on the basic thermoelectric properties. Before this, however, we consider various stable and quasi-stable magnetic states of the nanoribbons. Then, we consider separately the antiferromagnetic (AFM) states, with magnetic moments at one edge being antiparallel to the magnetic moments at the other edge), and ferromagnetic (FM) states, with parallel alignment of all the edge moments.

### iii.1 Stability diagram

We have performed numerical calculations of the total energy of nanoribbons with the edge magnetic moments tilted by an angle from the orientation normal to the nanoribbon plane, for ranging from to . The calculations have been carried out for both FM and AFM magnetic configurations. The numerical results clearly show that the states for are neither stable nor quasi-stable. The only stable or quasi-stable states are those for (moments normal to the plane) and (moments in the nanoribbon plane). This applies to both FM and AFM configurations. From this analysis, we have constructed the stability diagram of the FM and AFM states. Due to the spin-orbit interaction, both AFM and FM configurations have lower energy when the edge moments are oriented in the nanoribbon plane than when they are along the normal to the plane, and the energy difference depends on the Hubbard parameter and spin-orbit parameter . Thus, for each magnetic configuration one can distinguish between the perpendicular () and parallel () states.Lado2014 () In the former case the edge moments are normal to the nanoribbon plane, while in the latter one they are in the nanoribbon plane. As a result, one finds the following four states, FM, AFM, FM and AFM. Numerical results show that the AFM state is the ground state of the system, i.e. the state of the lowest energy for .

In Fig. 1 we show the energy of the states FM, AFM, and FM, measured from the energy of the AFM ground state. The excess energy, , of a particular state over the AFM ground state is presented there as a function of for indicated values of the spin-orbit parameter . One should note, that only nonmagnetic solutions have been found for low values of , while magnetic phases appear for larger values of . ¿From Fig. 1a follows, that the perpendicular configuration of magnetic moments is not favorable in the presence of spin-orbit coupling, which is in agreement with the results obtained by Lado et al.Lado2014 () Both AFM and FM states correspond to energies higher than that of the AFM state, and the energy difference strongly increases with the spin-orbit parameter. As one might expect, the AMF and AMF states are practically degenerate in the limit of small values of . In turn, energy of the FM state is considerably higher than the energy of the AMF state, and the difference increases roughly linearly with increasing . Interestingly, the energy difference between the perpendicular phases AFM and FM for a given decreases with increasing , so these two configurations become almost degenerate for strong spin-orbit coupling and moderate values of the Hubbard parameter . Remarkably different behavior can be observed for the FM state, as shown in Fig. 1b. For high values of and moderate values of , this state is practically degenerate with the AFM one, but the energy of FM state notably increases when the spin-orbit interaction becomes small. Finally, in the limit of , the states AFM and AFM are fully degenerate for all relevant values of (Fig. 1a), whereas the FM states are quasi-stable and correspond to higher values of energy.

Below we will analyze transport and thermoelectric properties of the AFM ground state as well as of the quasi-stable states AFM, FM, and FM. Though the latter states correspond to energy larger than the ground-state one, they can be stabilized by some external forces, like a weak magnetic field, exchange coupling to an adjacent magnetic subsystems, etc.

### iii.2 Phonon heat conductance

To calculate some of the thermoelectric parameters, like thermoelectric efficiency (figure of merit), one needs to know not only the electronic contribution to the heat conductance, but also the corresponding phonon term . In a recent paper Yang et al Yang2014 () have determined the phonon contribution to the heat conductance of narrow germanene and silicene nanoribbons, utilizing the Landauer formula (ballistic transport approximation) and the phonon transmission functions calculated by an ab-initio method. We checked that presented in Ref.[Yang2014, ] scales linearly with the nanoribbons width, and then extrapolated up to zigzag nanoribbons containing zigzag chains (16 atoms in the zigzag chain), which are considered in this paper. We obtained the following values for and : for germanene, and for silicene. Thus, for the model calculations we take in Eqs (10a) and (10b).

### iii.3 Antiferromagnetic phase

Now, we analyze ballistic transport and thermoelectric properties in the AFM configuration in more details. The corresponding electronic states for the AMF (left panel) and AMF (right panel) configurations are presented in Fig. 2 for three different values of the spin-orbit parameter. The corresponding transmission functions are presented in Fig. 3. The lowest value of assumed in these figures is comparable to the spin-orbit parameter typical for siliceneEzawa2012b () (note the factor in the definition of spin-orbit coupling in this reference). However, larger values of the parameter can appear in other related two-dimensional crystals, like germaneneLiu2011 () () or staneneLiu2011 () (). Moreover, results for larger values of the spin-orbit parameter are very interesting from the point of view of the topological properties. Therefore, in numerical calculations of the band structure, transmission, and thermoelectric coefficients we assume the following three values of : meV (strong spin-orbit coupling), meV, and meV (weak spin-orbit coupling). In turn, for the electron-electron interaction we assume eV.

To emphasize behavior of the edge states, which are responsible for topological properties, we displayed in Figs. 2 and 3 only a small energy region near the Fermi level of the corresponding charge-neutral nanoribbon. Moreover, only the region of wavevectors , in which the states are strongly localized near the left (L) or right (R ) edges of the nanoribbons is presented there. As follows from these figures, the electronic spectra reveal an insulating gap for small values of , and this gap appears in both AFM and AFM phases. Moreover, width of the gap is practically the same in both cases. Note, such a gap is absent in the limit of , when zero-energy topologically protected states appear.Wierzbicki2015 () Furthermore, the valence and conduction states localized at the left and right edges are degenerate in both phases. These states are additionally spin degenerate.

Qualitatively different behavior can be observed for larger values of . The edge states are then still degenerate, but in the AFM phase they cross the Fermi level, so the system acquires properties of a topological insulator, whereas the AFM phase still exhibits properties of a conventional insulator. Thus, at a certain value of , there is a transition from the conventional to topological insulator in the AFM phase, while no such a transition appears in the AFM phase. Moreover, one may conceive that by applying an appropriate external magnetic field it would be possible to switch between the AFM and AFM states, and thus also to change the character of the system from conventional to topological insulator, and vice versa. Opening of an energy gap near the Fermi level at the transition from the AFM state to the AFM can be clearly seen in the corresponding transmission function presented in Fig. 3.

Numerical results for the whole range of the angle (not presented), i.e. including also the hypothetical unstable magnetic states, show that for a fixed value of the width of energy gap strongly depends on the angle as well as on the Hubbard parameter . Generally, the gap increases with and , independently of the nanoribbon width. Interestingly, the system is non-magnetic for and , but for higher values of it exhibits the antiferromagnetic arrangement with quite considerable edge moments. Transition to the AFM state opens an energy gap, and the gap width monotonically increases with .

Thermoelectric properties in the AFM phase were already analyzed in our earlier paper for a small value of .Wierzbicki2015 () Therefore, we focus here mainly on the AFM ground state, though some new features in the thermelectric properties in the AFM state appear for large spin-orbit coupling, especially those associated with the transition from the conventional to topological insulator state. The energy gap in the AFM phase has a significant influence on the transport and thermoelectric properties of the system. The electrical conductance and the thermal conductance due to electrons are practically equal to zero in a wide energy region (corresponding to the gap) near the Fermi level , and then they strongly increase at the gap edges, see Fig. 4. This figure shows both and as a function of the chemical potential, and for both AFM and AFM (insets) phases. Interestingly, both and in the AFM state only weakly depend on . This is rather obvious as the corresponding gap in the spectrum (see Figs. 2 and 3) only weakly depends on the spin-orbit coupling parameter. Contrary, owing to the transition to the topological insulator state in the AFM phase, and closure of the gap for large enough spin-orbit coupling, both and are constant and nonzero for chemical potentials in the region around (see the insets in Fig. 4. It is also worth to note, that in some narrow regions of the chemical potential, the electrical and thermal conductances in the AFM and AFM phases are smaller for higher values of .

Due to the difference in the edge states, both AFM and AFM phases display also remarkably different thermoelectric properties, especially the Seebeck coefficient (thermopower), see Fig. 5. Due to the gap in the AFM phase, see Fig. 3, the thermopower considerably increases for chemical potentials in the vicinity of the Fermi level , achieving very high values at low temperatures. Since the gap weakly depends on the spin-orbit coupling, the thermopower also weakly depends on . The situation in the AFM state is different. Now, the gap disappears for large values of , so the thermopower remains zero in a broad region of chemical potentials around , see the inset in Fig. 5. Then, smaller peaks appear due to the narrow peaks in transmission around meV, as shown in Fig. 3.

The presented results show that nanoribbons with relatively strong spin-orbit coupling and exhibiting AFM magnetic configuration may be attractive for applications for conversion of thermal into electrical energy. The corresponding thermoelectric efficiency, described by the figure of merit , is shown in Fig. 6. Unfortunately, , so the thermoelectric efficiency is below that anticipated for practical applications. However, one may conceive some ways of reducing the phonon heat conductance in order to increase .

### iii.4 Ferromagnetic phase

Now, we analyze ferromagnetic arrangement of the edge moments and compare the results with those obtained for the corresponding antiferromagnetic phase, presented and discussed in the preceding section. Both FM and FM configurations are quasi-stable states with the corresponding energies above the energy of the AFM ground state. The ferromagnetic configurations, however, can be stabilized externally by a magnetic field or due to proximity effects.

The calculated band structure near the Fermi level is presented in Fig. 7 for both FM and FM phases, whereas the corresponding spin-dependent transmission functions are shown in Fig. 8. As in the AFM configurations, electronic spectra of both FM and FM phases are similar for small values of , and both phases are in the topological insulator state. The situation changes for larger values of . According to Fig. 7, the FM phase remains in the topological insulator state, independently of the spin-orbit coupling strength assumed in this figure. However, for higher values of the degeneracy of states localized at right and left edges of the nanoribbon is lifted. Significantly different qualitative behavior is obtained for the FM phase. Though the system is in the topological insulator phase for very small values of , its character changes for stronger spin-orbit interaction, and the system transforms into a conventional insulator with a relatively wide energy gap. It is interesting to note, that the range of -wavevectors, where the edge states appear in the FM phase, is much narrower than in the corresponding FM state as well as in the previously discussed AFM phases (see the dotted vertical lines in Fig. 7). Moreover, the states localized at the left and right edges are roughly degenerate.

The opening of an energy gap in the FM phase with increasing and absence of such a gap in the FM state are well seen in Fig. 8, where the spin-dependent transmission function is presented as a function of energy. It is worth to note that the transmission function corresponding to a particular spin orientation, up or down, is asymmetrical near the Fermi level, i.e. transmission in the spin-up channel exhibits a peak above , while the appropriate peak in the spin-down channel appears below the Fermi energy. This practically applies to both FM and FM configurations. The width of the energy gap in the FM state increases with increasing . However, no gap appears in the FM phase. Thus, the situation is remarkably different from that in the AFM configurations discussed above.

In the ferromagnetic phases, the two spin channels contribute differently to electronic transport. If the spin relaxation time is long and therefore spin mixing is negligible, the spin-up and spin-down channels can be treated as independent. One can then distinguish contributions of individual spin channels to the electrical conductance (and also to the Seebeck coefficient). The spin resolved electrical conductance is presented in Fig. 9a as a function of chemical potential. The wide central region with the conductance close to zero can be clearly seen for the FM phase. This region corresponds to the energy gap. The conductance strongly increases near the gap edges, where it also remarkably depends on the spin orientation. In turn, the conductance in the FM phase is finite and constant in the whole region of chemical potential shown in Fig. 9a (see the inset), except the two peaks corresponding to the enhanced transmission in narrow regions of chemical potential for the spin-up () and spin-down () channels, which appear near meV.

The thermal conductance due to electrons is shown in Fig. 9b for two values of the spin-orbit parameter. Interestingly, in some regions of chemical potential the thermal conductance is smaller for strong spin-orbit coupling, see Fig. 9b. This appears mainly because the energy gap in the FM state is wider for stronger spin-orbit coupling, which effectively suppresses in a broader region of the chemical potential. This difference is especially pronounced in the FM state (see the inset in Fig. 9b), where both peaks, i.e. the one below and another one above the Fermi level, are notably reduced for the larger value of the spin-obit parameter .

In addition to the conventional thermopower, one can now also define the spin thermopower. Both the charge and spin Seebeck coefficients are presented in Fig. 10 for two values of the parameter . Due to the energy gap in the FM state, the conventional Seebeck coefficient is considerably enhanced in the vicinity of the Fermi level. The enhancement strongly depends on and increases for higher values of the spin-orbit coupling. Similarly, the spin thermopower also exhibits a very pronounced peak for a narrow region of chemical potential near , and then remains rather small for chemical potentials outside this region. The peak in is very narrow, much narrower than the relevant energy gap. Furthermore, height of the peak strongly depends on the spin-orbit parameter. Remarkably different behavior has been obtained for the FM phase which exhibits typical topological insulator features. Now, the and Seebeck coefficients are rather independent of the chemical potential in the vicinity of the Fermi level, as shown in the insets to Fig. 10). The small peaks which appear for higher values of follow from the narrow peaks in the corresponding transmission functions, see Fig. 8.

The corresponding figures of merit for the conventional and spin thermoelectricity are shown in Figs. 11a and 11b, respectively. Significant values of the figures of merit occur only in certain regions of the chemical potential, where the corresponding thermopowers reach local maxima. Both of them, however, are relatively small, similarly as it was in the case of AFM configurations. However, especially small is the spin figure of merit, even if the corresponding spin thermopower is relatively high. This is because the spin thermoelectric efficiency in the narrow region of chemical potential around , where the relatively high peak appears in the spin thermopower, is negligible (compare Fig. 10b and Fig. 11b). This suppression of appears due to vanishing spin conductance in this region of chemical potential (note, the transmission functions for both spin orientations are equal this region, see Fig. 8).

## Iv Summary and conclusions

We have analyzed the influence of Coulomb and spin-orbit interactions in nanoribbons of 2D buckled hexagonal crystals with zigzag edges on the edge magnetic moments and on the topological and thermoelectric properties. We have shown that the Coulomb interaction, taken in the form of Hubbard term in the mean field approximation, leads to formation of the edge magnetic moments, which for a sufficiently strong spin-orbit coupling are oriented in the nanoribbon plane in the ground state, and the moments at one edge are opposite to those at the other edge (AFM phase) – in agreement with Lado et al.Lado2014 () For small spin-orbit coupling, this configuration is almost degenerate with the state, in which magnetic moments are oriented perpendicularly to the nanoribbon plane (AFM phase). States with ferromagnetic arrangements of the moments (FM and FM phases) are quasi-stable and correspond to higher energies. Such states can be stabilized externally, e.g. by a magnetic field and/or due to proximity to a ferromagnetic contact.

As the AFM phase remains in the topological insulator state when the pin-orbit coupling increases, the AFM phase undergoes then a transition from the topological to conventional insulator, i.e. an energy gap opens in the corresponding edge states. This gap has a significant influence on transport and thermoelectric properties. First of all, the opening of the gap is associated with a nonzero thermopower in the vicinity of the gap edges. The ferromagnetic quasi-stable states, FM and FM, behave differently with increasing , as well. This behavior qualitatively resembles that of the antiferromagnetic phases, but now the gap in the FM state, which opens with increasing , is spin dependent. As a result, a nonzero spin thermopower may be observed, in addition to the conventional one.

###### Acknowledgements.

This work was supported by the National Science Center in Poland as Project No. DEC-2012/04/A/ST3/00372.## References

- (1) K.S Novoselov, A. K. Geim, S. V. Morozov, D. Jang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306, 666 (2004).
- (2) B. Aufray, A. Kara, S. Vizzini, H. Oughaddou, C. Léandri, B. Ealet, and G. Le Lay, Appl. Phys. Lett. 96, 183102 (2010).
- (3) A. Kara, H. Enriquez, A. P. Seitsonen, L. C. Lew Yan Voon, S. Vizzini, B. Aufray and H. Oughaddou, Surf. Sci. Rep. 67, 1 2012.
- (4) B. Bishnoi and B. Ghosh, RSC Adv. 3, 26153 (2013).
- (5) M. Ezawa, New J. Phys. 14, 033003 (2012).
- (6) M. Ezawa, Phys. Rev. B 87, 155415 (2013).
- (7) M. Ezawa and N. Nagaosa, Phys. Rev. B 88, 121401 (2013).
- (8) J.L. Lado and J. Fernandez-Rossier, Phys. Rev. Letters 113, 027203 (2014).
- (9) A.H.C. Neto, F. Guinea, N.M.R. Peres, K.S. Novoselov, and A.K. Geim, Rev. Mod. Phys. 81, 109 (2009).
- (10) M. I. Katsnelson, Graphene: Carbon in Two Dimensions (Cambridge Univ. Press, 2012).
- (11) A. Lopez-Bezanilla, J. Huang, P. R. C. Kent, and B. G. Sumpter, J. Phys. Chem. C 117, 15447 (2013).
- (12) J. M. Zhang, F. L. Zheng, Y. Zhang, and V. Ji, J. Mater. Sci. 45, 3259 (2010).
- (13) P. Lou and J. Y. Lee, J. Phys. Chem. C 113, 21213 (2009).
- (14) P. Lou, Phys. Chem. Chem. Phys. 13, 17194 (2011).
- (15) J. Guan, W. Chen, X. Zhao, G. Yu, X. H. and Ch. Sun, J. Mater. Chem. 22, 24166 (2012).
- (16) P. Lou, Phys. Stat. Sol. B 250, 1265 (2013).
- (17) W. Chen, H. Zhang, X. Ding, G. Yu, D. Liu and X. Huang, J. Mater. Chem. C 2, 7836 (2014).
- (18) X. Ding, G. Yu, X. Huang and W. Chen, Phys. Chem. Chem. Phys. 15, 18039 (2013).
- (19) O. V. Yazyev, Rep. Prog. Phys. 73, 056501 (2010).
- (20) Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature (London) 444, 347 (2006).
- (21) L. Pisani, J. A. Chan, B. Montanari, and N. M. Harrison, Phys. Rev. B 75, 064418 (2007).
- (22) H. Haugen, D. Huertas-Hernando, and A. Brataas, Phys. Rev. B 77, 115406 (2008).
- (23) W. Y. Kim and K. S. Kim, Nature Nanotechnology 3, 408 (2008).
- (24) X. H. Zheng, J. Lan, X. L. Wang, L. F. Huang, H. Hao and Z. Zeng, Appl. Phys. Lett. 101, 053101 (2012).
- (25) J. Li, Z. H. Zhang, D. Wang, Z. Zhu, Z. Q. Fan, G. P. Tang, and X. Q. Deng, Carbon 69, 142 (2014).
- (26) D. Wang, Z. H. Zhang, Z. Zhu, Z. L. Yu, B. K. Tao, S. Y. Hou, and C. Ye, Organic Electronics 15, 3406 (2014).
- (27) C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005).
- (28) M. Z. Hasan and C. L. Kane, Rev. Modern Phys. 82, 1057 (2010).
- (29) X.-L. Qi and S.-C. Zhang, Rev. Modern Phys. 83, 1057 (2011).
- (30) A. I. Hochbaum, R. Chen, R.D. Delgado, W. Liang, E.C. Garnett, M. Najarian, A. Majumdar, and P. Yang, Nature 451, 163 (2008).
- (31) T. C. Harman, P. J. Taylor, M. P. Walsh, and B. E. LaForge, Science 297, 2229 (2002).
- (32) N. B. Duarte, G. D. Mahan, and S. Tadigadapa, Nano Lett. 9, 617 (2009).
- (33) M. Walter, J. Walowski, V. Zbarsky, M. Münzenberg, M. Schäfers, D. Ebke, G. Reiss, A. Thomas, P. Peretzki, M. Seibt, J. Moodera, M. Czerner, M. Bachmann, and C. Heiliger, Nature Materials 10, 742 (2011).
- (34) N. Liebing, S. Serrano-Guisan, K. Rott, G. Reiss, J. Langer, B. Ocker, and H. W. Schumacher, Phys. Rev. Lett. 107, 177201 (2011).
- (35) K. Uchida, S. Takahashi, K. Harii, J. Ieda, W. Koshibae, K. Ando, S. Maekawa, and E. Saitoh, Nature 455, 778 (2008).
- (36) R. Swirkowicz, M. Wierzbicki, and J. Barnaś, Phys. Rev. B 80, 195409 (2009).
- (37) M. Misiorny and J. Barnaś, Phys. Rev. B 89, 235438 (2014).
- (38) L. Pan, H. J. Liu, X. J. Tan, H. Y. Lv, J. Shi, X. F. Tang, and G. Zheng, Phys. Chem. Chem. Phys. 14, 13588 (2012).
- (39) K. Zberecki, M. Wierzbicki, J. Barnaś, and R. Swirkowicz, Phys. Rev B 88, 115404 (2013).
- (40) Y. Yan, H. Wu, F. Jiang, and H. Zhao, Eur. Phys. Journ. B 86, 457 (2013).
- (41) K. Zberecki, R. Swirkowicz, and J. Barnaś, Phys. Rev B 89, 165419 (2014).
- (42) K. Zberecki, R. Swirkowicz, M. Wierzbicki, and J. Barnaś, Phys. Chem. Chem. Phys. 16, 12900 (2014).
- (43) M. Wierzbicki, J. Barnaś, and R. Swirkowicz, Phys. Rev. B 91, 165417 (2015).
- (44) L. Matthes and F. Bechstedt, Phys. Rev B 90,165431 (2014).
- (45) R. Swirkowicz, Acta Phys. Polon. A 90, 1215 (1997).
- (46) C. Caroli, R. Combescot, P. Nozieres and D. Saint-James, J. Phys. C 4, 916 (1971).
- (47) P.A. Khomyakov, G. Brocks, V. Karpan, M. Zwierzycki, P.J. Kelly, Phys. Rev. B 72, 035450 (2005).
- (48) M. Goano, ACM Transactions on Mathematical Software 21, 221 (1995).
- (49) K. Yang, S. Cahangirov, A. Cantarero, A. Rubio, R. D’Agosta, Phys. Rev. B 89, 125403 (2014).
- (50) M. Ezawa, Phys. Rev. Lett. 109, 055502 (2012).
- (51) C.-C. Liu, H. Jiang, Y. Yao, Phys. Rev. B 84, 195430 (2011).