Spin-current pump in silicene*

Project supported by the National Natural Science Foundation of China (Grant Nos. 11274059, 11574045, and 11704165).

Jada Marcellino John Tombe1, Wang Mei-Juan1, Wang Sa-Ke2, Wang Jun1, †
Department of Physics, Southeast University, Nanjing 210096, China
Department of Fundamental Courses, Jinling Institute of Technology, Nanjing 211169, China

 

† Corresponding author. E-mail: jwang@seu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11274059, 11574045, and 11704165).

Abstract

We report a theoretical study of pumped spin currents in a silicene-based pump device, where two time-dependent staggered potentials are introduced through the perpendicular electric fields and a magnetic insulator is considered in between the two pumping potentials to magnetize the Dirac electrons. It is shown that giant spin currents can be generated in the pump device because the pumping can be optimal for each transport mode, the pumping current is quantized. By controlling the relevant parameters of the device, both pure spin currents and fully spin-polarized currents can be obtained. Our results may shed a new light on the generation of pumped spin currents in Dirac-electron systems.

1. Introduction

Spintronics is currently a very active research field in condensed matter physics, aiming to exploit both the charge and spin degrees of freedom of electrons in the new generation of spin-based devices.[1] Since graphene has many important properties with high application potentials, such as the gate voltage-controlled carrier conduction and high mobilities, spin currents generated or manipulated in graphene have attracted a great deal of research attention.[26] Especially, the weak spin–orbital interaction and the lack of nuclear spin in graphene will enable the maintenance of spin coherence over longer distances than that in conventional semiconductors. An experiment[7] on spin injection into single layer graphene found that the spin coherence length is about 1 μm in dirty samples and a much longer coherence length is expected in a clean system. This is also attributed to the properties of the massless Dirac electrons with linear energy dispersion.

Silicene, similar to graphene, is a monolayer of silicon atoms instead of carbon atoms assembled into a two-dimensional (2D) honeycomb lattice and also, electrons in silicene around the Fermi energy are linearly dispersed. Recently, this material has been experimentally synthesized and theoretically explored.[813] Different from graphene, silicene has a prominent buckled structure that the A and B silicon atoms are not coplanar, which allows one to control the mass of the Dirac electrons by a perpendicular electric field and leads to some peculiar transport properties of electrons in silicene. Therefore, the generation and manipulation of spin currents of Dirac electrons in the silicene system are similarly very attractive.

Quantum spin and charge pump is a striking topic in the context of quantum transport through nanostructures.[1416] The nature of the quantum parameter pump arises from the geometric Berry phases and quantum interference effects. Generally, a nonzero pumped current can be generated in a pump device with at least two pumping parameters out of phase. In the adiabaticity that the characteristic time of the parameter variations is much shorter than the dwell time of the carriers, the pumped current is proportional to the geometric area circled by the time-dependent parameters.[16] Quantum parameter pump is a significant and effective means to generate spin currents.[1722]

In this work, we propose to generate spin currents in a two-parameter pump device based on silicene, where the pumping potentials are simulated by the time-dependent perpendicular electric field applied to the system, while a magnetic insulator is assumed in between the two pumping potentials to induce possible magnetization. It is shown that the device can pump out spin currents, including both pure spin currents and fully spin-polarized charge current. Due to the quantum interference effect, the pumped spin currents can be modulated by controlling the magnetization parameters. When the Fermi energy is residing in the energy gap of the pumping potentials, the pumped currents are maximized.

2. Model

Let us start from a silicene-based pump device as schematically shown in Fig. 1, where the two staggered pumping potentials Δ1 and Δ2 are assumed from the time-dependent perpendicular electric field Ez with Δ1,2 = Ezl and l being the perpendicular distance between the A and B silicon atom planes, and a magnetic insulator[23] (say, EuO) is deposited on the silicene between the two pumping potentials with a length L0. Thus, a constant magnetization M is assumed to be induced in the silicene. The whole pump device is connected to the outside world through the left and right electrodes but without any bias. The sizes of the two pumping potentials are marked by their lengths L along the current direction (the x-axis) and the transverse direction (along the y-axis) is assumed homogeneous or to have a translational symmetry.

Fig. 1. (color online) Schematic of a two-terminal pump device based on a monolayer silicene. The two pumping sources, Δ1 = Δ0 cosωt and Δ1 = Δ0 cos(ωt + ϕ), are the staggered potentials with a phase lag ϕ between them. The pumped current direction is assumed along the x axis and the device is connected to electrodes through the left and right leads. L is the length of the pump region and L0 is the length of the magnetic insulator M inserted between Δ1 and Δ2.

A low-energy continuum Hamiltonian is adopted here to describe the pump device as[24]

where σx,y,z is the lattice spin operator, Θ1 = Θ(x)Θ(Lx), and Θ2 = Θ(xL0L)Θ(2L + L0x), with Θ(x) being a Heaviside step function. The first term describes the free Dirac electrons of silicene, the second term stands for the homogeneous magnetization assumed from the proximity effect of the magnetic insulator, sz is the realistic spin operator, and the magnetization M is set along the z axis. kx,y is the 2D momenta and η is the valley index. V0 denotes the local controllable potential and σ0 is a unit matrix. Δ1 = Δ0 cos(ωt) and Δ2 = Δ0 cos(ωt + ϕ) are homogeneous in the finite-length L region, where ω is the pumping frequency, Δ0 is the pumping strength, and ϕ is the pumping phase shift. The magnetized silicene length between Δ1 and Δ2 is assumed to be L0. A staggered potential breaking the inversion symmetry of silicene can open an energy gap for the Dirac electrons and the system becomes a quantum valley Hall insulator. Here, it comes from the assumed perpendicular electric field applied onto the silicene plane. The pumping potentials with opposite signs will generate a topological interface state between them, which in turn will lead to a quantized pumping current[25] when the Fermi energy lies in the energy gap due to Δ1,2.

In the adiabatic limit ω → 0, the pumped current can be calculated directly from the BPT formula as[26]

where T = 2π/ω is the pump cycle, is the instantaneous scattering matrix with α(β) = L, R being the left (right) lead index, which is often expressed as
Here and are left (right) reflection and transmission matrices, respectively. The trace is over the transverse modes and ± = ↑(↓) are the spin-up (↑) and spin-down (↓) indices. Since we consider the adiabatic limit, the transverse momentum ky is assumed to conserve in the instantaneous scattering process. It is also noted that in the following paragraph, we take only one valley into account (η = 1) because of the valley degeneracy.

In order to work out the scattering coefficients, one can construct a single scattering process to solve them directly or choose the following combining method to obtain them. We have the formulae of two relevant coefficients r± and through two individual scattering events at Δ1 or Δ2 potentials,

with
which is the dynamic phase of electrons passing through the magnetization region. ri and ti (i = 1, 2) are the reflection and transmission amplitudes of a single pumping potential Δ1 or Δ2, and V0 is assumed to be the overall potential profile of the magnetized region that can be modulated by an external voltage gate. For instance, we could utilize the following scattering process to solve r1 and t1:
Here three wavefunctions are respectively in the three regions of the first pumping potential Δ1. kx is the momentum of electrons in the normal silicene region \, and the counterpart in the potential region is given by (ℏvF = 1). Similarly, r2 and t2 can be obtained by using the same procedure by solving the scattering event at the second pumping potential Δ2. Then, and would be calculated by considering the injection from the right region conjugate to the process of Eq. (7).

3. Results

With these preparations above, we can directly calculate the spin-resolved pumped current IL↑(↓) (± = ↑(↓)) flowing through the left lead of the pump device in Fig. 1 according to Eq. (3). We set ℏvF = 1, so the Fermi energy E, the momentum kx(ky), and the magnetization strength M have the same energy unit for simplicity. Since the quantized pumped current is focused on, in that the Fermi energy resides in the energy gap opened by Δ1 and Δ2, and the transverse momentum ky does not bring any qualitative effect on the results, we can, for simplicity, consider a pure one-dimensional case ky = 0.

We first plot the pumped current (IL) versus the pumping phase ϕ and the results are shown in Fig. 2. The curves deviate from the usual sine function behavior of an adiabatic quantum parameter pump I ∼ sinϕ, since the pumped current obtained here is quantized as the Fermi energy is taken as EΔ0, i.e., IL↑,↓ = ±e/T. One can see that in Fig. 2(a), the spin-up and spin-down pumped currents are almost opposite, which implies that the pump device can pump out pure spin currents. Whereas in Fig. 2(b), the magnetization is changed from M = 0.003 eV to M = 0.01 eV, both spin species currents are almost the same so that pure spin current almost vanishes.

Fig. 2. (color online) Spin-resolved pumped current–phase relationship ILϕ. The magnetization strength is (a) M = 0.003 eV and (b) M = 0.01 eV. Other parameters are E = 0.001 eV, Δ0 = 0.01 eV, L = 1000a, and L0 = 1000a with a being the lattice constant.

The quantization of the pumped currents is ascribed to the possible topological interface state[25] forming between the two pumping potentials, which is evolving with time. As mentioned above, the staggered potential will result in a topological valley Hall phase appearing in the pumping region, but this topological phase state has no boundary or edge state. However, such an interface state with zero energy will appear in the bridging region when Δ1 and Δ2 have the opposite signs. It is this topological interface state that leads to a nonzero pumped current as well as its quantization when the Fermi energy locates in the energy gap. Otherwise, an insulator would not conduct any current. When ϕ, the pumped current is deviated from the quantized value e/T in Fig. 2. The reason is that the effective energy gap is quite small so that the Fermi energy locates outside of the energy gap.

Since the pumping effect is in nature a quantum interference effect, the pumping results can be modulated by the dynamic phase of the electrons. Therefore, it is shown that different magnetization M will lead to different spin-dependent pumped currents IL↑/IL↓ and this is attributed to the different phase φ± in Eq. (6). So we can obtain pure spin currents in the studied pump device. According to this reasoning, one can also utilize the length L0 of the magnetization region to control the pumping output and instead, we choose to modulate the overall potential profile V0 to change the pumped spin currents.

In Fig. 3, the spin-dependent currents are presented as a function of V0 with different lengths L0. It is seen that the spin degeneracy is broken by M and the spin-up and spin-down momenta in the magnetized silicene are split and changed by V0, so that one can obtain a nonzero spin current Is = IL↑IL↓, and even a fully spin-polarized current with IL↑ = 0 or IL↓ = 0. Since V0 is varied, the electron momentum in the magnetization region would change accordingly so as to the dynamic phase φ±. From this, one can also see that even the variation of the Fermi energy can modulate the pumped spin currents too, because E and V0 play the same roles in determining φ±. Similarly, it is not strange that when L0 is different, the periodicity of the pumped spin currents with V0 should also be different as shown in Figs. 3(a) and 3(b); the larger the L0, the shorter the V0 periodicity.

Fig. 3. (color online) Pumped current IL as a function of the gate voltage V0. The length of the magnetization region is (a) L0 = 1000a and (b) L0 = 500a. Other parameters are E = 0.001 eV, Δ0 = 0.01 eV, ϕ = π/2, L = 1000a, and M = 0.005 eV.

It is worth pointing out that it is not uncommon to obtain pumped spin currents in a device with magnetization involved, nevertheless, the studied model can provide an optimal pumping effect of spin currents, i.e., for each transport model, the pumped spin is maximized. Moreover, when the Fermi energy E locates in the energy gap ∼Δ0, the pumping is quantized[25] for all transverse modes ky. Although the magnetization introduced in the pumping device here is limited between the two pumping devices, this is not compulsory and it can extend to all devices. The key point to spin currents is the spin-splitted momentum in the normal silicene region without pumping potentials resulting in an energy gap of silicene. We have discussed the optimal spin currents pumped in a silicene-based device by using the time-dependent electric fields that induce the staggered potentials. Actually, such a conclusion is also valid in a graphene-based device by using the time-dependent strains, which can play the same role as the staggered potentials, however, the transport direction is quite important since the graphene cannot have a transport gap along all directions, which is necessary for the quantized pump.

4. Conclusion

We have investigated the pumped spin currents in a silicene-based pump device. The pumping potentials are from the perpendicular electric fields that result in the staggered potentials in the silicene and the magnetization is assumed to be induced by the magnetic proximity effect. According to the adiabatic pumping theory, we calculated the pumped currents flowing in the device and found that the currents are quantized when the Fermi energy resides in the energy gap opened by the staggered potentials. The pure spin currents and fully spin-polarized currents can be obtained by modulating the relevant parameters such as the Fermi energy and the magnetization strength. Our findings are helpful to generate spin currents in the Dirac-electron system.

Reference
[1] Žutić I Fabian J Das Sarma S 2004 Rev. Mod. Phys. 76 323
[2] Shakouri K Masir M R Jellal A Choubabi E B Peeters F M 2013 Phys. Rev. B 88 115408
[3] Szaszkó-Bogár V Peeters F M Földi P 2015 Phys. Rev. B 91 235311
[4] Akturk A Goldsman N 2008 J. Appl. Phys. 103 053702
[5] Han W Kawakami R K Gmitra M Fabian J 2014 Nat. Nanotechnol. 9 794
[6] Bolotin K I Sikes K J Jiang Z Klima M Fudenberg G Hone J Kim P Stormer H L 2008 Solid State Commun. 146 351
[7] Tombros N Tanabe S Veligura A Jozsa C Popinciuc M Jonkman H T van Wees B J 2008 Phys. Rev. Lett. 101 46601
[8] Lalmi B Oughaddou H Enriquez H Kara A Vizzini S Ealet B Aufray B 2010 Appl. Phys. Lett. 97 223109
[9] Vogt P de Padova P Quaresima C Avila J Frantzeskakis E Asensio M C Resta A Ealet B Le Lay G 2012 Phys. Rev. Lett. 108 155501
[10] Lin C-L Arafune R Kawahara K Tsukahara N Minamitani E Kim Y Takagi N Kawai M 2012 Appl. Phys. Exp. 5 45802
[11] Cahangirov S Topsakal M Aktürk E Sahin H Ciraci S 2009 Phys. Rev. Lett. 102 236804
[12] Liu C C Feng W Yao Y 2011 Phys. Rev. Lett. 107 76802
[13] Ezawa M 2012 Phys. Rev. Lett. 109 55502
[14] Makhlin Yu Mirlin A D 2001 Phys. Rev. Lett. 87 276803
[15] Zhou F Spivak B Altshuler B 1999 Phys. Rev. Lett. 82 608
[16] Brouwer P W 1998 Phys. Rev. B 58 R10135
[17] Zhang Q Chan K S Lin Z 2011 Appl. Phys. Lett. 98 032106
[18] Zhang Q Chan K S Lin Z Liu J F 2013 Phys. Lett. A 377 632
[19] Wang J Chan K S Lin Z 2014 Appl. Phys. Lett. 104 013105
[20] Abdollahipour B Mohammadkhani R 2014 J. Phys.: Condens. Matter 26 85304
[21] Chen M N Sheng L Shen R Sheng D N Xing D Y 2015 Phys. Rev. B 91 125117
[22] Deng W Y Luo W Geng H Chen M N Sheng L Xing D Y 2015 New J. Phys. 17 103018
[23] Haugen H Huertas-Hernando D Brataas A 2008 Phys. Rev. B 77 115406
[24] Luo W Sheng L Wang B G Xing D Y 2016 Sci. Rep. 6 31325
[25] Wang J Liu J F 2017 Phys. Rev. B 95 205433
[26] Büttiker M Thomas H Prêtre A 1994 Z. Phys. B: Condens. Matter 94 133