Abstract

We performed the first-principles calculations to investigate the spin-dependent electronic transport properties of zigzag-edged germanium nanoribbons (ZGeNRs). We choose of ZGeNRs with odd and even widths of 5 and 6, and the symmetry-dependent transport properties have been found, although the mirror plane is absent in ZGeNRs. Furthermore, even- and odd- ZGeNRs have very different current-voltage relationships. We find that the even 6-ZGeNR shows a dual spin-filter effect in antiparallel (AP) magnetism configuration, but the odd 5-ZGeNR behaves as conventional conductors with linear current-voltage dependence. It is found that when the two electrodes are in parallel configuration, the 6-ZGeNR system is in a low resistance state, while it can switch to a much higher resistance state when the electrodes are in AP configuration, and the magnetoresistance of 270% can be observed.

1. Introduction

Graphene and graphene nanoribbons (GNRs) have attracted broad academic and industrial interests owing to their amazing properties [1, 2] and potential applications in nanodevices [35]. In particular, their applications in spintronics field [68] offer the most promising solutions for future high operating speed and energy-saving electronic devices. For instance, the zigzag graphene nanoribbon (ZGNR) is an anti-ferromagnetic insulator with considerable magnetic moment located at edge sites, which give rise to many unusual properties, such as half metallic [9], magnetoelectric effect [10, 11], and spin-filter effect [11, 12].

These interesting behaviors of graphene motivate the further exploration of honeycomb lattice with higher group IV elements. Walking down the periodic table from carbon to silicon (Si) to germanium (Ge) enhances the -orbital participation in their electronic behavior and thus develops strong electronic correlation. And there are many interesting magnetic and charge excitations behaviors that would be introduced. Theory predicts free standing silicene or germanene sheets and ribbons to be stable in a low puckered configuration [1316]. Different from graphene, in silicene or germanene, two adjacent atoms belonging to the same sublattice are not in the same plane owing to the buckled structure of silicene or germanene. Thus, the buckled structural feature would give rise to tunable spin-valley coupled band structures, which accounts for many exotic transport [17, 18] and superconducting phenomena [19] and makes silicene and germanene good candidates for the spintronics [2023]. The first-principle calculations show that the intrinsic carrier mobility of germanene sheet can even reach ~6 × 105 cm2 V−1 s−1, which is larger than that of graphene [24]. In particular, compared with Si, Ge has longer spin-diffusion length and the relative strong spin-orbit coupling would be helpful to the spin injection and the manipulation of spin [25, 26]. Thus, the developments of multifunctional germanene-based spintronics components and the corresponding nanodevices design are more worth studying thoroughly. A feasible approach to realize germanene-based spintronics is to construct the device junctions of germanium nanoribbons (ZGeNRs) and then systematically study the transport properties as fundamental and crucial components of germanene-based circuits. Recently, intriguing transport phenomena have been found in zigzag graphene nanoribbons (ZGNRs) [27], zigzag -graphyne nanoribbons [28], and silicene nanoribbons (ZSiNRs) with buckled structure [21], where symmetric and asymmetric zigzag ribbons are found to have completely different - characteristics despite the similarity of their band structure. More interestingly, bipolar spin diode behavior was observed in symmetry devices. Spin polarization and direction of the current through the device can be controlled through either the voltage or magnetic configuration of the electrodes. Such freedom in controlling the spin-polarized current is useful in design of spin devices such as rectification and amplification for spintronics. However, until now, no one is working on the spin transport property of symmetric and asymmetric ZGeNRs. Developments of multifunctional germanium-based spintronics devices that offer effective manipulation of spin-polarized current are an important supplement of exploration of spintronics.

So, in this work, we investigate the spin-dependent transport properties of ZGeNRs by first-principles calculations. Our results show that although the plane is absent in ZGeNRs due to the low puckered configuration, the ZGeNRs still exhibit symmetry-dependent transport properties. And the magnetoresistance effect is also observed in ZGeNRs. Our works suggest a route to manipulate spin-polarized current and to design the germanene-based spintronic devices.

2. Models and Theoretical Method

ZGeNRs can be classified by the number of Ge atoms along the width of the ribbon. We choose ZGeNRs with odd and even widths of 5 and 6, which are denoted here as 5-ZGeNR and 6-ZGeNR, respectively. The geometric structures from top view and cross section view of 6-ZGeNR are shown in Figures 1(a) and 1(b), respectively. The system is divided into three regions: left electrode, right electrode, and the central scattering region. Each electrode is described by a supercell with two repeated ZGeNR unit cells along transport direction, and the scattering region is ZGeNR with length of 6 unit cells. An external magnetic field controlling the magnetization of the left (right) ZGeNR electrode can be set to be 1 or −1, corresponding to magnetization along the or direction. In this paper, we use and to represent the magnetization of the left and right electrodes; that is, shows that the left and right electrodes are parallel (P) in spin orientation, while shows that the left and right electrodes are antiparallel (AP) in spin orientation. The geometric optimization and spin-resolved electronic transport properties are calculated by a developed ab initio software package Atomistix ToolKit [29, 30], which is based on the spin-polarized density-functional theory combined with the nonequilibrium Green functions. In the calculation, the Troullier-Martins norm-conserving pseudopotential, the local spin density approximation (LSDA) for exchange-correlation potential, and -point grid 1 × 1 × 100 are used. The real space grid techniques are used with the energy cutoff of 150 Ry as a required cutoff energy in numerical integrations and the solution of Poisson equation using fast Fourier transform (FFT). The geometry optimization is performed using quasi-Newton method until the absolute value of force acting on each atom is <0.05 eV/Å. The wave functions of Ge and H atoms are expanded by single-zeta polarized (SZP) basis set. The spin-resolved current is calculated using the Landauer formula:where (spin up) and (spin down), is the electrochemical potential of the left/right electrode, and the difference of them is , and is the spin-resolved transmission defined aswhere is the retarded (advanced) Green functions of the central region and is the contact broadening functions.

3. Results and Discussions

The spin configuration is crucial for the spin-filter effect, and many researchers have found that ZGNR, ZGNRs, and ZSiNR can show perfect spin-filter effect with AP configuration [12, 21, 27, 28]. So, in here, as shown in Figure 2, we consider firstly that the current-voltage (-) curves come with the spin-dependent transmission spectrum as a function of the electron energy and bias of 5-ZGeNR and 6-ZGeNR under the AP magnetic configuration. From the Landauer formula (1), we know that the spin-resolved current through central scattering region is decided by the transmission within the considered bias window. Therefore, the properties of - curves can be understood based on the transmission spectrum. In Figure 2(b), for spin-up state of 5-ZGeNR, the area within the bias window increases linearly under the whole bias window; thus a linear - curve can be observed in Figure 2(a). The spin-down state of 5-ZGeNR presents a similar property, as shown in Figures 2(c) and 2(d). Interestingly, as shown in Figures 2(e)2(h), the - curves of 6-ZGeNR can be modulated by the spin orientation of the electrodes and the direction of external . It is clearly seen in Figure 2(f) that the spin-up transmission within the bias window is always zero in the positive bias range, which leads to full suppression of spin-up current () under the positive bias as shown in Figure 2(e). Nevertheless, in the negative bias range, the variation of the current with the bias voltage goes through two different stages. When bias is less than 0.20 V, as shown in Figure 2(e), there are more and more transmission coefficients entering into the bias window with the increase of bias voltage, leading to a quick increase of current. When the bias voltage is above 0.2 V, the rate of the increase of transmission in the bias window is slower owing to the presence of regions with smaller transmission coefficient, which lower down the increase of the current as shown in Figure 2(e). The spin-down state of the AF configuration shows a completely opposite property. That is to say, in the positive bias range, only the spin-down current () is allowed to pass through the central scattering region, while in the negative bias range, only the spin-up current is possible. A device with such property is referred to as a bipolar spin diode [6] or a dual spin filter [9, 12] due to the unidirectional nature of the spin-dependent current, which is different from the traditional bipolar devices based on the polarization of the electron or hole carriers. In contrast, the current for 5-ZGeNR with the AP configuration is almost independent of spin and shows a linear metallic behavior. This interesting phenomenon is very similar to those of [21, 27, 28], which is possible because of a symmetry selection rule that was first proposed by Li et al. [27]. And, for symmetry one, the flexible control over spin current makes it possible to use the two-terminal device as a basic component for building multifunctional spintronic devices. Compared with other 2D materials nanoribbons such as ZGNR, ZGNRs, and ZSiNR spin devices, ZGeNR-based spin devices possess several unique features besides the common point that the spin polarization of the current can be simply tuned by switching either the magnetization of the electrodes or the direction. First, the intrinsic carrier mobility of germanene is much larger than that of graphene, making it a good candidate of ultrafast devices. Second, germanene is a much better spin transport medium because the spin-diffusion length in germanene is much longer than that in graphene and silicene.

The above results also illuminate high spin polarization and rectifying ratio () for 6-ZGeNR in the AP configuration. As shown in Figures 3(a) and 3(b), we calculate the corresponding spin polarization (SP) and rectifying ratio (), defined by the formulas SP = and , respectively. And we can clearly see that the maximum values of SP can reach 92%, and the maximum of rectifying ratio is 22 for the spin-up current of 6-ZGeNR.

These interesting spin transport properties can be understood from the transmission spectra and band structures of the left and right electrodes for the 5-ZGeNR and 6-ZGeNR at zero bias, as shown in Figure 4. From the band structures panels, we can find that the bands of spin-up and spin-down states around FL are split in each system. And also, owing to the different spin orientation, we can find that the subbands of spin-up and spin-down states are opposite on the left and right electrodes. In the central panel of Figure 4(a), we can see the broad transmission platforms for both spin-up and spin-down states in the 5-ZGeNR at zero bias. As the bias increases, the transmission probability remains to be one ( is the quantum conductance) around the Fermi level () for both spin states. Thus, the - curves are linear and display metallic behaviors in Figures 2(a) and 2(c). For 6-ZGeNR, the transmission spectra of spin-up and spin-down states are degenerate at zero bias, and a transmission gap (TG) appears around in central panel of Figure 4(b). And there is a large continuous transmission platform when energy is beyond the TG.

As we know, the transmission spectrum can be driven by spin matching of the left and right electrodes in a spin transport device [31]. Transmission is allowed only when both the spin and orbital symmetry of the energy bands on the two electrodes are matched; otherwise, it would be forbidden. To understand the origin of the different transmission behaviors around FL in 5-ZGeNR and 6-ZGeNR devices, as shown in Figures 4(c) and 4(d), we give the isosurface plots of the wave functions of and subbands for spin-up states of 5-ZGeNR and 6-ZGeNR at the -point, respectively. It is known that ZGeNRs have no definite parity with respect to the operation because of the buckled structure. From Figure 4(c), we cannot find any definite parity for the wave functions of and subbands in 5-ZGeNR. However, it is noted that the center axis of 6-ZGeNR (dashed line in Figure 4(d)) is a twofold axis, and 6-ZGeNR has symmetry with respect to this axis. From Figure 4(d), we can find that the orbitals of the and subbands have odd and even parity, respectively. Therefore, the TG is equal to the energy area with the parity mismatching of the orbitals with respect to the operation. Thus, although the plane is absent in ZGeNR, symmetry still plays an important role in the transport properties.

To understand the nature of the observed dual spin-filter effect, we plot the band structures of left and right electrodes and transmission spectra for 6-ZGeNR at bias of ±0.5 V, as shown in Figures 5(a) and 5(b). When the positive bias is applied, the energy bands are shifted downward (upward) for the left (right) electrode. Within the bias window, for 6-ZGeNR, as shown in Figure 5(a), the subband of left electrode overlaps with the subband of right electrode for spin-down state, since the and subbands have opposite parities (as in Figure 4(d)) with respect to the midplane mirror operation, the electron transmission from the subband of left electrode to the subband of right electrode is forbidden, and no electron transmission contributes to the current in Figure 2(g). Nevertheless, for the spin-up state, there is a wide spin-up transmission peak entering into the bias windows, which mainly originates from the coupling between the () subband of left electrode and the () subband of right electrode; thus a large spin-up current can be observed in Figure 2(e). When bias is 0.5 V, the shifting of energy bands of the left and right electrodes and transmission matching of different spin states are just opposite to the case of −0.5 V, which leads to zero current for spin-up and larger current for spin-down states. As a result, the dual spin-filter effect can be observed on the - curves of 6-ZGeNR device.

Moreover, in order to further understand the transport properties, we also show the transmission pathway at 0.24 eV for the spin-up and spin-down currents of 6-ZGeNR under the bias of ±0.5 V. The transmission pathways can show us the local bond contributions to the transmission coefficient. The pathways have the property that if the system is divided into 2 parts (), then the pathways across the boundary between and sum up to the total transmission coefficient . From Figures 5(c) and 5(d), it is shown that the spin-up electrons can transfer through the whole scattering region at −0.5 V and reach the other electrode, corresponding to a high transmission peak in middle panel of Figure 5(b), while for the spin-down electrons, the transmission channel is forbidden, and electrons cannot reach the right electrode, corresponding to the transmission coefficient which is zero at 0.24 eV in middle panel of Figure 5(b). At 0.5 V, the transmission pathways demonstrate that the current channels for the spin-up state are blocked. In contrast, the current channels of the spin-down states are unblocked, and electrons can transfer through the whole scattering region. Interestingly, we can clearly find in the GeNR-based devices that there are also two different current channels like GNR-based devices: bond current (current via a chemical bond) and hopping current (current through hopping between atoms of the same symmetry sublattice) [12, 32]. From Figures 5(c)5(f), we can clearly find that the hopping current is obviously weaker than the bond one under the same bias voltage.

By modulating the direction of spin polarization orientations of the left and right electrodes, giant magnetoresistance has been found in the ZGNRs and ZSiNRs [8, 21]. Thus, we wonder whether the same magnetoresistance behavior would also appear in ZGeNRs. Figure 6(a) shows the - curves for 6-ZGeNR in P and AP spin configurations, and Figure 6(b) presents the corresponding magnetoresistance ratio (MR) as a function of applied bias. We can see that the current of 6-ZGeNR in P configuration is much bigger than that in AP configuration at low bias voltages, and the magnetoresistance behavior appears. The magnetoresistance ratio at a certain bias is defined by MR = , where and are the currents in the P and AP configurations related to the same bias voltage, respectively. As shown in Figure 6(b), MR is gained at the bias voltage range from 0.01 to 0.1 V and its maximum reaches about 270%.

Moreover, since the width effect is always important in nanoribbons [3335], we also considered the dependence of the transport property on the width of the ZGNRs. In Figure 7(a), we show the - curves for -ZGeNRs ( = 5, 7, and 9). As can be seen, the current for odd ZGeNRs with AP spin configuration is almost independent of spin and linearly increases under low , leading to of nearly one. But for even -ZGeNRs ( = 6, 8, and 10), shown in Figure 7(b), the current-voltage (-) characteristics depend on interestingly not only spin, but also the direction of , and the up-spin electron substantially flows only in the negative regime of , while the down-spin electron flows only in the positive regime. Thus, it can be concluded that the width of ZGeNRs does not affect the qualitative transport properties and the dual spin-filter effects can also be observed in all widths of even ZGeNRs. The major reason is that the electronic transport properties for ZGeNRs are mainly determined by the parity of the and subbands of left and right electrodes; the parity of ZGeNRs is always maintained even when the width of the nanoribbon is changed. Meanwhile, it is interesting to note that both the spin polarization and rectifying ratio are increased with the increasing of the width as shown in Figures 7(c) and 7(d). For 10-ZGeNR, the spin splitting effect can reach 93%, and rectifying ratio is 23. This means that the performance of the GeNRs-based spin filter and spin rectifier would be modulated by controlling the ribbon width.

4. Summary

In conclusion, we have investigated the electronic transport properties of 5-ZGeNR and 6-ZGeNRs devices. Our results show that the symmetry-dependent transport characteristics are observed in ZGeNRs, although the mirror plane is absent, and the different symmetry of the wave functions of and subbands at the -point would determine the spin-dependent transport behavior of ZGeNRs. We also find that the even 6-ZGeNR shows a dual spin-filter effect in AP magnetism configuration, but the odd 5-ZGeNR behaves as conventional conductors with linear current-voltage dependence. Moreover, it is found that when the two electrodes are in P configuration, the 6-ZGeNR system is in a low resistance state, while it can switch to a much higher resistance when the electrodes are in AP configuration. And the order of the corresponding magnetoresistance ratio can reach 270%. We also find that the spin polarization and rectifying ratio are increased with the increasing of width in even ZGeNRs. These results may be useful in designing of high performance spin filter and spin rectifier based on GeNRs.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grants nos. 21103232 and 61306149), the Natural Science Foundation of Hunan Province (no. 14JJ3026), China Postdoctoral Science Foundation (nos. 2013M542130 and 2014M552145), Hong Kong Scholars Program (no. XJ2013003), and Shenghua Lieying Scholarship by the Central South University.