# Spin Current Generation and Magnetic Response in Carbon Nanotubes by the Twisting Phonon Mode

## Abstract

We theoretically investigate spin current and magnetic response induced by the twisting phonon mode in carbon nanotubes via the spin-rotation coupling. An effective magnetic field due to the twisting mode induces both spin and orbital magnetizations. The induced spin and orbital magnetizations have both radial and axial components. We show that AC pure spin current is generated by the twisting phonon mode. The magnitude of the spin current and orbital magnetization for a (10,10) armchair nanotube is estimated as an example. We find that the AC pure spin current is detectable in magnitude when the frequency of the twisting mode is of the order of GHz, and that the orbital magnetization is found to be larger than the spin magnetization.

###### pacs:

61.48.De, 72.25.-b, 62.25.-g, 85.75.-dIntroduction. Carbon nanotubes (CNTs) have been studied extensively, and today various interesting physical characters of CNTs have been reported. CNTs can be regarded as rolls of graphene sheets, and therefore their electronic bands are closely related with those of graphene. In graphene, the electronic bands form linear dispersion called a Dirac cone. By rolling the graphene into a CNT, the band structure becomes a metal or semiconductor depending on its chirality. From the viewpoint of spin transport, the spin-orbit interaction of CNT is so weak that spin current cannot be generated for CNT alone by usual methods which require strong spin-orbit interaction.

In this paper, we propose a new possibility of the spin current generation by using a twisting mode, one of the characteristic phonon modes in
the CNT Kane *et al.* (1998); Kane and Mele (1997); Sánchez-Portal *et al.* (1999).
Our proposal of spin current generation by the twisting mode is based on a spin-rotation coupling.
The spin-rotation coupling can be derived in relativistic quantum mechanics Matsuo *et al.* (2011); de Oliveira and Tiomno (1962); Mashhoon (1988), and
it couples mechanical rotation with spins Matsuo *et al.* (2013a, 2011).
Physically
it can be regarded as the Barnett effect Barnett (1915); mechanical rotation serves as an
effective magnetic field, and couples with spin via Zeeman interaction.
In fact, this effective magnetic field couples with electrons not only via spins but also via orbitals, as seen from Ref. Matsuo *et al.* (2011). Thus,
generation of both spin and orbital magnetizations by the twisting mode is expected.
We show that in CNTs, AC pure spin current is generated by the twisting mode. We also show that the twisting mode induces an orbital magnetization.
We find that the orbital magnetization is larger than the spin magnetization for realistic parameters. We estimate the magnitude of the spin current and orbital magnetization for a (10,10) armchair nanotube as an example,
and we find that the AC pure spin current is detectable in magnitude when the frequency of the twisting mode is of the order of GHz.

Since it is shown that in CNTs and graphene, spins have a relatively long lifetime, they are promising as leads for spintronics devices.
Spin current injection into CNT were reported in Refs. Tombros *et al.* (2006); Yang *et al.* (2012), and
the methods of spin current injection are spin valves and non-local spin injection, both of which need a ferromagnet to inject spin. On the contrary, the generation of spin current by the twisting mode does not require a ferromagnet to generate spin current, and would pave the way for spintronics application of CNTs.

Twisting mode. Among the three types of phonon modes in CNTs, i.e. stretching, breathing, and twisting modes Suzuura and Ando (2002), we focus on the twisting mode, which represents a deformation along the circumference of the tube.
In the Cartesian coordinate shown in Fig. 1, the displacement of the twisting mode is given by
with the angular amplitude , wavenumber of the twisting mode , and frequency of the twisting mode . Here is the velocity of the mode given by with shear modulus and the mass of a carbon atom .
We note that the dispersion of the
twisting mode is almost independent of the chirality of CNT, since the phonon dispersion in a CNT is almost the same as that in graphene Sánchez-Portal *et al.* (1999).
The angular velocity is given by

(1) |

Then the angular velocity vector , which characterizes the local rotation of the CNT in Fig. 1, is

(2) |

with , and the tube radius .
This angular velocity works as an effective magnetic field Matsuo *et al.* (2013b)
due to spin-rotation coupling, where is the electron charge and is the electron mass.
Hence in the present case, the effective magnetic field and the corresponding vector potential are given as

(3) | |||||

(4) |

The component is along the axial direction, while the other components are along the radial direction. This effective magnetic field affects both the electronic spins and orbitals. First, it enters the Zeeman coupling:

(5) |

where is the Bohr magneton. Second, the orbital magnetic field enters the Hamiltonian via minimal coupling: , which affects electronic orbital motions.

Spin current. We will calculate the spin current in CNTs induced by the effective
Zeeman magnetic field due to the twisting mode. We consider the twisting phonon mode in the semi-infinite system of CNTs.
We use the spin diffusion equation and calculate spin current polarized along the -axis, using the formalism proposed in Ref. Matsuo *et al.* (2013a).
Here we retain only the -component of the Zeeman magnetic field and neglect the other components, i.e.
, because the twisting mode is assumed to have a long wavelength. We will discuss justification of this assumption later.
The spin-rotation coupling generates a difference of chemical potentials between up- and down-spins. The spin diffusion equation with spin-rotation coupling Matsuo *et al.* (2013a) is written as

(6) |

with the diffusion constant , the spin lifetime . The spin current polarized along the -axis is then given by

(7) |

where is the electrical conductivity. Note that the generated spin current is not accompanied by charge current. By substituting the angular velocity Eq. (1) into Eq. (6) and solving the differential equation, we obtain

(8) |

By substituting Eq. (8) into Eq. (7), we can numerically calculate the spin current. We set the realistic parameters of the CNT obtained in Ref. Yang *et al.* (2012): and at .
As an example, we consider the armchair CNT with radius , circumference , the group velocity of the twisting mode Suzuura and Ando (2002), and conductivity at Fischer *et al.* (1997). The twisting mode is assumed to be described by , and .
The result is shown in Fig. 2. The behavior can be classified into two regimes, (Fig. 2 (a)) and (Fig. 2 (b)) with spin diffusion length . When , the behavior of spin current is nearly the same as that of the twisting mode, since an influence of the reflection at the edge of CNT is negligible. This is because for , the diffusion and relaxation are much faster than the twisting mode. When , the behavior changes around
the edge of CNT: the signal is strongly enhanced. Sufficiently away from the edge, the behavior is nearly the same as that of the twisting mode and the magnitude of the spin current is proportional to the frequency of the mode.
For , the spin current diffuses slowly and hence can be reflected at the edge.

This theory requires presence of carriers, and therefore our theory applies to metallic CNT, doped semiconducting CNT, or semiconducting CNT at finite temperature. The resulting spin current depends on the chiralities and temperature through the physical quantities , and . When the frequency and wavenumber are and , the maximum of spin current is evaluated as , as shown in Fig. 2 (b), which is measurable in experiments. By exciting the twisting mode in CNT with the piezoelectric device, the spin current can be measured via inverse spin Hall effect. For example, by making a contact between Pt and CNT, one can convert the spin current into a voltage signal. When a ferromagnet is attached to CNT, the generated spin current can be injected into the ferromagnet. This spin current can then exert spin torque on the magnetzation, leading to magnetization reversal.

Orbital magnetization. Next we investigate the orbital magnetization of CNTs by the orbital magnetic field. Here we focus on its radial component perpendicular to the tube axis. We follow the method to derive the orbital susceptibility for other magnetic field directions in Refs. Ajiki and Ando (1993, 1992). For convenience in calculation, we introduce a flux penetrating along the tube axis. We also assume that the frequency is sufficiently low and the response to the orbital magnetic field is regarded to be static. We use the cylindrical coordinate in Fig. 1.

To calculate the orbital magnetization, we first calculate the shift of the energy levels up to the second order in the orbital magnetic field . The vector potential can be written as . The Hamiltonian near the point is given by

(9) |

with the envelope function near the point , the band parameter , and being a wave vector operator defined by . In the absence of the orbital magnetic field along radial direction, , the corresponding eigenvector is

(10) |

and the energy eigenvalue is where is the wave number in the circumferential direction with circumference length , is the band index, is a dimensionless parameter with the flux quantum , correspond to metal and semiconductor CNT, respectively, is the wave number along the -direction, and denotes the conduction and valence bands, respectively.

The envelope function near the point and the vector potential are then Fourier expanded in terms of in the range of . Then, the perturbation due to the magnetic field is written as

(11) |

with the magnetic length . The matrix elements of the perturbation between the eigenstates for are calculated as,

(12) |

Therefore, the magnetic field causes an energy shift for the eigenvalues to the second order in :

(13) |

At zero temperature, the shift of the total energy becomes

(14) |

where is a cutoff function, and is the length of CNT. The susceptibility per unit area is calculated as . Therefore, we obtain the susceptibility as

(15) | |||||

with , dimensionless parameter , and lattice constant . Adding the contribution from the point and multiplying the result by a spin-degeneracy factor of two, we obtain the susceptibilities for metalic and semiconducting CNT ( and ) as follows:

(16) | |||||

(17) |

The result of the susceptibility for the radial magnetic field is shown in Fig. 3(a). For comparison, we show the susceptibility for a uniform magnetic field perpendicular to the tube axis in Fig. 3(b), using the method in Ref. Ajiki and Ando (1993). It is seen that the magnitude of the diamagnetic susceptibility for the radial magnetic field is about ten times larger than that for the magnetic field perpendicular to the tube axis. This is related with the enhanced diamagnetic response of the graphene sheet for the out-of-plane orbital magnetic field, when the Fermi energy is close to the Dirac point McClure (1956). We note that the susceptibility in a radial magnetic field survives above room temperature, as is similar to that in a magnetic field perpendicular to the tube axis in Ref. Lu (1995).

Discussion. We discuss the magnitude of the effective magnetic field. For this purpose, we decompose the effective magnetic field Eq. (3) into the component along the tube axis and that along the radial direction :

(18) | ||||

(19) |

We evaluate them for the armchair nanotube as an example. The twisting mode is characterized by the group velocity , and the radius is Suzuura and Ando (2002). For the frequency of the twisting mode and amplitude of the twisting mode , we have . The magnitude of the effective magnetic field is and , and the ratio of and is

(20) |

The axial magnetic flux is given by where denotes the cross section of CNT. Hence, we have . Similarly, for the frequency , we have , , , and . Therefore, when the frequency of the twisting mode is between and , the ratio of is between and , and the axial magnetic flux is of the order from to . Hence, the radially-polarized AC spin current due to the Zeeman magnetic field is from to times smaller than that along the axial direction. Therefore, we can neglect the radial direction component , and the axial magnetic flux is very small.

Next, we roughly estimate four types of magnetizations: axial spin magnetization , radial spin magnetization , axial orbital magnetization , and radial orbital magnetization . The spin magnetization is given by the Pauli paramagnetic susceptibility with density of state . When the Fermi energy is , the density of state is . Then the magnetization per unit area is .

The orbital susceptibility by the radial magnetic field is given by Eqs. (16) and (17), and is shown in Figs. 3(c) and (d) for various values of the radii. It can be seen that for thicker CNT, the susceptibility becomes larger and dependence on becomes weaker. These behaviors are similar to those in response to uniform magnetic field perpendicular to the tube Ajiki and Ando (1993). For a rough estimate of the magnetization, we use the averaged value of the susceptibility from to , , which is common for metallic and semiconducting CNTs. The magnetization per unit area is . For the axial orbital magnetization , the susceptibility diverges at Ajiki and Ando (1993), and therefore the response is nonlinear in . Therefore, we use the results in Ref. Ajiki and Ando (1993) to calculate for the given value of

In the frequency range from to , the magnitudes of the four magnetizations are evaluated as follows: to , to , to , and to . We find that the orbital magnetizations are larger than the spin magnetizations in both directions, and axial magnetizations are larger than radial magnetizations . The axial orbital magnetization is the largest among the four in this setup, and is the main contribution. The magnetization can be measured by the magnet-optic Kerr effect or using SQUID.

Summary. In this paper, we have theoretically investigated spin current and magnetic response induced by the twisting phonon mode in carbon nanotubes via the spin-rotation coupling. We found that the AC spin current and orbital magnetizations are generated by the twisting phonon mode. The spin current is detectable in magnitude when the frequency is of the order of GHz. It was shown that the orbital magnetizations are larger than the spin magnetizations, and the axial orbital magnetization is the main contribution for realistic values of the parameters.

###### Acknowledgements.

We thank M. Matsuo, J. Ieda, K. Harii, and S. Maekawa for useful discussions. This work was supported by Grant-in-Aid for Young Scientists (B) (No. 23740236) , the “Topological Quantum Phenomena” (No. 25103709) Grant-in-Aid for Scientific Research on Innovative Areas (No.26103006), and MEXT Elements Strategy Initiative to Form Core Research Center (TIES) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan.### References

- C. L. Kane, E. J. Mele, R. S. Lee, J. E. Fischer, P. Petit, H. Dai, A. Thess, R. E. Smalley, A. R. M. Verschueren, S. J. Tans, and C. Dekker, EPL (Europhysics Letters) 41, 683 (1998).
- C. L. Kane and E. J. Mele, Phys. Rev. Lett. 78, 1932 (1997).
- D. Sánchez-Portal, E. Artacho, J. M. Soler, A. Rubio, and P. Ordejón, Phys. Rev. B 59, 12678 (1999).
- M. Matsuo, J. Ieda, E. Saitoh, and S. Maekawa, Phys. Rev. Lett. 106, 076601 (2011).
- C. de Oliveira and J. Tiomno, Il Nuovo Cimento 24, 672 (1962).
- B. Mashhoon, Phys. Rev. Lett. 61, 2639 (1988).
- M. Matsuo, J. Ieda, K. Harii, E. Saitoh, and S. Maekawa, Phys. Rev. B 87, 180402 (2013a).
- S. J. Barnett, Phys. Rev. 6, 239 (1915).
- N. Tombros, S. J. van der Molen, and B. J. van Wees, Phys. Rev. B 73, 233403 (2006).
- H. Yang, M. E. Itkis, R. Moriya, C. Rettner, J.-S. Jeong, D. S. Pickard, R. C. Haddon, and S. S. P. Parkin, Phys. Rev. B 85, 052401 (2012).
- H. Suzuura and T. Ando, Phys. Rev. B 65, 235412 (2002).
- M. Matsuo, J. Ieda, and S. Maekawa, Phys. Rev. B 87, 115301 (2013b).
- J. E. Fischer, H. Dai, A. Thess, R. Lee, N. M. Hanjani, D. L. Dehaas, and R. E. Smalley, Phys. Rev. B 55, R4921 (1997).
- H. Ajiki and T. Ando, J. Phys. Soc. Jpn. 62, 2470 (1993).
- H. Ajiki and T. Ando, J. Phys. Soc. Jpn. 62, 1255 (1992).
- J. W. McClure, Phys. Rev. 104, 666 (1956).
- J. P. Lu, Phys. Rev. Lett. 74, 1123 (1995).