Bipolar Magnetic and Thermospin Transport Properties of Graphene Nanoribbons with Zigzag and Klein Edges

Magnetic nanoribbons based on one-dimensional materials are potential candidates for spin caloritronics devices. Here, we constructed ferromagnetic graphene nanoribbons with zigzag and Klein edges (N-ZKGNRs, N� 4–21) and found that the NZKGNRs are in the indirect-gap bipolar magnetic semiconducting state (BMS). Moreover, when a temperature difference is applied through the nanoribbons, spin-dependent currents with opposite flow directions and opposite spin directions are generated, indicating the occurrence of the spin-dependent Seebeck effect (SDSE). In addition, the spin-dependent Seebeck diode effect (SDSD) also appeared in these devices. More importantly, we found that the BMS with a larger bandgap is promising for generating the SDSD, while the BMS with a smaller bandgap is promising for generating the SDSE. *ese findings show that ZKGNRs are promising candidates for spin caloritronics devices.


Introduction
Spin caloritronics, which focuses on the interaction between the spin, charge, and heat in materials, has attracted intensive interest because it plays an important role in the development of fundamental science and novel low-powerconsumption technologies [1][2][3][4][5][6]. In this field, an important effect is the spin-dependent Seebeck effect (SDSE) [7][8][9], which features different Seebeck thermopowers for different spins in spin-polarized systems. In recent years, there have been many reports concerning the SDSE in graphene nanoribbons (GNRs) with armchair or zigzag edges [3,10,11]. e original work was done by Zeng et al. [3], and they demonstrated the SDSE in magnetized zigzag GNRs (ZGNRs). Ni et al. [10] then proposed a new spin caloritronics device based on the ZGNR heterojunction and observed the SDSE. In previous work, we also observed the SDSE in GNRs and silicene nanoribbons (SiNRs) with armchair or zigzag edges [12][13][14]. As we have known, armchair and zigzag/Klein edges can be formed by cutting the graphene honeycomb lattice along the <1 1 10> and <2 1 1 0> directions, respectively. However, studies of the SDSE in ZGNRs with reconstructed edges (i.e., the Klein edge) and the pentagon-heptagon edge are rare. Actually, researchers have already experimentally observed these two types of reconstructed edges for zigzag GNRs [15,16], and they have theoretically demonstrated that these edge reconstructions have drastic influences on the band structures and magnetic states of ZGNRs [17,18]. Hence, in this study, we investigated systematically the electronic structures and thermal spin-dependent transport characteristic of ZKGNRs by ab initio calculations combined with the nonequilibrium Green's function approach. We found that ZKGNRs can be used in the stable bipolar magnetic semiconductor [19][20][21], in which the valence bands and conduction bands approach the Fermi level through opposite spin channels, and the SDSE can also be obtained.
e results indicate that ZKGNRs are promising for application in spin caloritronics devices. zigzag edge and one Klein edge (see Figures 1(a) and 1(b)), and both edges are saturated by two H atoms. Here, N denotes the number of carbon rows across the GNR width, and N � 4-21. We then built two-probe spin caloritronics devices based on the N-ZKGNRs (N � 4 as an example), as shown in Figure 1(c). e left and right contacts are semiinfinite ZKGNRs, and the central scattering region contains five units of ZKGNRs. We focus on the spin currents driven by the temperature difference (ΔT) between the source temperature T L and drain temperature T R , that is, All of the calculations were performed with the Atomistix Toolkit (ATK) package [22,23], which uses spin density functional theory in combination with the nonequilibrium Green's function method. e geometry optimization and electronic structure calculations were performed with the double-zeta-polarized (DZP) basis set, and the exchange-correlation potential was treated by the generalized gradient approximation method [24,25]. e cutoff energy was 75 hartree, and a Monkhorst-Pack 1 × 1 × 100 k-mesh was chosen. In the Landauer-Büttiker formalism, the spin-dependent currents of the devices were obtained with the following equation [26]: where e is the electron charge, h is the Planck constant, and ) is the average Fermi-Dirac distribution of the left (right) electrode: where μ (R) and T L(R) are the chemical potential and temperature of the left (right) electrode, respectively, and k B is the Boltzmann constant. T ↑(↓) (E) is the spin-dependent transport coefficient: where G r(a) is Green's function that is retarded in the central region and Γ L(R) is the coupling matrix for the left (right) electrode. In addition, the calculation methods of the spinup thermopower (S up ), spin-down thermopower (S dn ), total charge thermopower (S ch ), and net spin thermopower (S sp ) can be found in our previous work [14]. It is obvious that all the E FM-AFM and E FM-NM are less than zero, confirming the FM ground state in the N-ZKGNRs. Next, we investigated the band structures of the N-ZKGNRs (N � 4-21) (Figure 3). e band structure of the nanoribbons noticeably changes as N increases. For N � 4 and 5, the conduction band minimum (CBM) states are related to the spin-down states and located at the Γ point, which are above the Fermi level (E F ), and the spin-up states below E F serve as the valence band maximum (VBM) states and are located at the Z point. e valence and conduction bands possess opposite spin polarization as they approach E F . Moreover, the spin-dependent bands of these ZKGNRs have finite gaps around E F . ese characteristics indicate that 4-ZKGNR and 5-ZKGNR are indirect-gap BMSs (IBMSs) [27][28][29]. As the nanoribbon width increases, the ferromagnetic configuration is maintained, whereas the band structure near the Fermi level remarkably changes. For N � 6, the CBM is related to the spindown states and located at the Γ point, and the spin-up states below E F serve as the VBM and are located at the Γ-Z line.

Results and Discussion
ese characteristics indicate that 6-ZKGNR is also an indirect-gap BMS. For N � 7-21, both the VBM of the spin-up channel and the CBM of the spin-down channel are located at the Γ-Z line with the same point. ese characteristics indicate that these ZKGNRs are direct-gap BMSs (DBMSs). e corresponding bandgaps are shown in Figure 4. In summary, as the width increases, the ZKGNR may have two different states, that is, the IBMS state and the DBMS state ( Figure 4).
To explore the SDSE, we turn our attention to the ZKGNR thermal spin transport properties. e spin-dependent currents through the N-ZKGNR (taking N � 4 and 21 as examples) devices versus T L and ΔT are shown in Figure 5. For the 4-ZKGNR device, there are no spin-up currents (I up ) and spin-down currents (I dn ) when T L < 300 K for the three values of ΔT. is suggests that no thermalinduced spin-dependent currents are generated in this T L range, irrespective of the temperature difference (ΔT). In other words, there is a threshold temperature (T th ) at approximately 300 K for I up and I dn . When T L > T th , both I up and I dn sharply increase with increasing T L . However, they flow in opposite directions, that is, I up is negative and I dn is positive. ere is no doubt that this is caused by the SDSE [30][31][32][33][34]. Furthermore, the spin-dependent current is larger for higher ΔT. e curves of spin-dependent current versus ΔT are shown in Figure 5(b), with T L set to 300, 350, and 400 K. For ΔT > 0, the curves clearly indicate that the spindependent currents are approximately symmetric about the zero-current axis and robust over a large range of temperature difference. erefore, the SDSE can also be confirmed by the curves of spin-dependent current versus ΔT. For ΔT < 0, both I up and I dn are approximately equal to zero.
ese characteristics indicate that the perfect SDSD occurs in the 4-ZKGNR device. e 21-ZKGNR device has similar characteristics to the 4-ZKGNR device, but it has smaller T th (T th � 80 K) and larger I up and I dn , as shown in Figure 5(c). In addition, as shown in Figure 5(d), both I up and I dn of the 21-ZKGNR are larger than zero, indicating a weak SDSD.

Advances in Materials Science and Engineering
To illustrate the underlying mechanism of these phenomena, we plot the corresponding spin-dependent transmissions of the N-ZKGNRs (N � 4, 21) in Figure 6. For the 4-ZKGNR (Figure 6(a)), the spin-up and spin-down transmissions are nearly symmetrically localized downside and upside of the Fermi level, which provides two independent spin-conducting channels. If a temperature difference is applied along the nanoribbons, the spin-down and spin-up current have opposite flow directions, unambiguously representing the occurrence of the SDSE, as shown in Figures 5(a) and 5(b). For the 21-ZKGNR (Figure 6(b)), the spin-dependent transmissions have similar characteristics to the 4-ZKGNR, except for a smaller bandgap, resulting in the smaller T th , as shown in Figure 5(c). erefore, the SDSE also occurs. It is important that we found that T th has a positive correlation with the bandgap, that is, T th is larger for larger bandgap. erefore, we can draw the following conclusion: the BMSs with larger bandgaps are promising candidates for generating the SDSD, while the BMSs with smaller bandgaps are promising candidates for generating the SDSE.
For further in-depth study into the spin thermal transport, we calculated the spin-up thermopower (S up ), spin-down thermopower (S dn ), total charge thermopower (S ch ), and net spin thermopower (S sp ) versus the device temperature at the Fermi level (E F � 0). As shown in     Figure 7, S up > 0, while S dn < 0, supporting the appearance of SDSE in these two devices. In particular, the S ch of 21-ZKGNR is nearly equal to zero, indicating a perfect SDSE. Furthermore, the S sp of 4-ZKGNR is larger than the S sp of 21-ZKGNR because 4-ZKGNR has a larger band gap.

Conclusion
We have investigated the electronic structures and thermal spin-dependent transport properties of a series of N-ZKGNRs by first-principle calculations combined with the nonequilibrium Green's function. First, we found that as the nanoribbon width parameter N is increased from 4 to 21, the N-ZKGNR transforms from an indirect-gap BMS to a directgap BMS. Second, the SDSE and SDSD appear if we produce a temperature difference through these ZKGNRs. Finally, the BMSs with larger bandgaps are promising candidates for generating the SDSD, while the BMSs with smaller bandgaps are promising candidates for generating the SDSE. ese findings strongly suggest that N-ZKGNRs are promising materials for spin caloritronics devices.

Data Availability
Data are available upon request.

Conflicts of Interest
e authors declare no conflicts of interest.