Next Article in Journal
Taxonomy of Three-Qubit Mermin Pentagrams
Next Article in Special Issue
Kramers Degeneracy and Spin Inversion in a Lateral Quantum Dot
Previous Article in Journal
Hydrodynamics Interactions of Metachronal Waves on Particulate-Liquid Motion through a Ciliated Annulus: Application of Bio-Engineering in Blood Clotting and Endoscopy
Previous Article in Special Issue
Modulational Instability, Inter-Component Asymmetry, and Formation of Quantum Droplets in One-Dimensional Binary Bose Gases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Symmetry Properties of The Corrugated Graphene System

1
Institute of Experimental Physics, Watsonova 47, 04001 Košice, Slovakia
2
BLTP, Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia
3
BLTP, Dubna University, 141982 Dubna, Moscow Region, Russia
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(4), 533; https://doi.org/10.3390/sym12040533
Submission received: 20 February 2020 / Revised: 19 March 2020 / Accepted: 24 March 2020 / Published: 3 April 2020
(This article belongs to the Special Issue Symmetry and Mesoscopic Physics)

Abstract

:
The properties of the ballistic electron transport through a corrugated graphene system are analysed from the symmetry point of view. The corrugated system is modelled by a curved surface (an arc of a circle) connected from both sides to flat sheets. The spin–orbit couplings, induced by the curvature, give rise to equivalence between the transmission (reflection) probabilities of the transmitted (reflected) electrons with the opposite spin polarisation, incoming from opposite system sides. We find two integrals of motion that explain the chiral electron transport in the considered system.

1. Introduction

It appears that graphene possesses a remarkable stretchability. For example, the DFT and molecular dynamics simulations predict that it can be stretched up to about 20–30%, without being damaged [1]. The experimental measurements demonstrate a good agreement with the theoretical estimations, while the robust engineering results indicate on sample-wide elastic strain ∼6% [2]. Evidently, transforming the flat surface to the curved one, one creates the strain that affects the graphene properties. This fact suggests that, by altering the stretchability, one might tune electronic and transport properties of the graphene sheet.
Recent experimental techniques enable demonstrating evidently a spatial variation of graphene and its direct consequences. For example, ripples can be formed by means of the electrostatic manipulation without any change of doping [3]. Periodically rippled graphene can be fabricated by the epitaxial technique (e.g., [4]). In this case, in contrast to free-standing graphene, a strong modification of the electronic structure of graphene is observed, which gives rise to localised phonon [5] and plasmon [6] modes. Periodic nanoripples can be created as well by means of the chemical vapour deposition [7]. It is found that ripples, acting as potential barriers, yield the localisation of charged carriers [8]. The potential surface variations could reach the figure of 20–30 meV. Similar independent prediction has been done in the study of Klein collimation by the rippled graphene superlattice [9]. In this model, the hybridisation between the π - and σ -orbitals creates the potential barrier between the flat and curved graphene pieces. The barrier value could reach Δ ε 24 meV at the ripple radius R = 12 Å. This fact provides the confidence in the vitality and the validity of our model (outlined in Section 2) and following analyses its symmetry properties, presented in our paper.
Indeed, the lattice deformation changes the distance between ions, p z orbital orientation, and is leading to shift of the on-site energies of p z orbitals. This affects the effective Dirac equation that could simulate the low energy electron states as a result of a deformation-induced gauge field [10]. The surface curvature modulates also the hopping parameter in the tight-binding approach [10,11]. Moreover, it enhances as well the effect of the spin–orbit coupling [12], usually neglected in flat carbon-based systems.
Recently, it is predicted that rippled graphene could lead to the spin selectivity effect for the ballistic electrons [13,14] in virtue of the curvature-induced spin–orbit couplings (see details in [12,13,15,16,17]). As a result, it is shown that at the particular energy values the ballistic electrons with the one spin polarisation can travel through periodically repeated ripples without any reflection. At the same time, electrons with the opposite spin polarisation are fully reflected. Once we change the flow direction through the considered system, the situation becomes inverse. It is noteworthy that different experiments of a spin selective electron transmission through biomolecules has been discussed recently in Ref. [18]. In this review, the authors claimed that this phenomenon implies that chirality and spin may play an important role in biology.
In mesoscopic systems, symmetries are key points that allow to illuminate essential features of finite quantum systems (e.g., [19,20]). The basic goal of this paper is to elucidate the above discussed phenomenon from point of view of the symmetry properties of the considered system.

2. Basic Physics of The Corrugated Graphene

In our consideration, the corrugated graphene structure consists of a rippled graphene connected to two flat graphene sheets (see Figure 1).
For analysis of the curved graphene surface, we recapitulate the major results [17] obtained for armchair CNTs. In this case, only the interaction between nearest neighbour atoms is considered. The analysis is done in an effective mass approximation for the point K, in the vicinity of the Fermi level E = 0 . A similar approach can be applied for K point.
The Hamiltonian of the nanotube has the following form in the effective mass approximation [12]
H r = γ ( τ x π ^ x + τ y π ^ y ) I λ y τ y σ x ( r ) ξ x τ x σ y .
Here, the operators π ^ x = i R θ , π ^ y = i y , σ x ( r ) = σ x cos θ σ z sin θ , and ξ x = 2 δ γ p / R , λ y = δ γ / 4 R . The Pauli matrices σ x , y act in the spin space. The matrices τ x , y act on the sublattice degree of freedom. The Pauli matrix τ i is called the “pseudospin”, to distinguish it from the real electron spin.
The following notations are used: γ = 3 V p p π a / 2 = 3 γ 0 a / 2 , γ = 3 ( V p p σ V p p π ) a / 2 = γ 1 a , p = 1 3 γ / 8 γ ; where V p p σ and V p p π are the transfer integrals for σ - and π -orbitals, respectively, in a flat graphene. The distance d between atoms in the unit cell determines the length of the primitive translation vector a = 3 d 2.46 Å. For numerical illustration, we assume that γ 0 = V p p π 3 eV and γ 1 8 eV (see, e.g., [12]). The intrinsic source of the spin–orbit coupling δ = Δ / 3 ϵ π σ is defined by means of the quantity
Δ = i 3 4 m 2 c 2 x | V x p y V y p x | y ,
where V is the atomic potential and ϵ π σ = ϵ 2 p π ϵ 2 p σ . Here, the energies ϵ 2 p π and ϵ 2 p σ denote the energies of π - and σ -orbitals orbitals, respectively. We recall that σ -orbitals are localised between carbon atoms, while π -orbitals are directed perpendicular to the nanotube surface.
With the aid of the method discussed in [13,17], one obtains the eigenvalues of the Hamiltonian in Equation (1)
E = κ E ± , κ = ± 1 ,
where κ = + 1 ( 1 ) is associated with the conductance (valence) band, and the energies E ± are defined as
E ± = t m 2 + t y 2 + λ y 2 + λ x 2 ± 2 λ x 2 t m 2 + λ y 2 + t y 2 λ y 2 .
Here, t m = m γ / R , t y = γ k y , λ x = γ ( 1 / 2 + 2 δ p ) / R , and m is a magnetic quantum number (see details in [13,17]). Due to the curvature-induced spin–orbit coupling spin is no anymore a good quantum number. The eigenstates of Equation (1) are characterised by a quantum number s = ± 1 , and have the following form
Φ m , k y s = ± 1 ( θ , y ) = e i m θ e i k y y N ± κ ( cos θ / 2 A ± sin θ / 2 B ± ) κ ( sin θ / 2 A ± + cos θ / 2 B ± ) cos θ / 2 C ± sin θ / 2 D ± sin θ / 2 C ± + cos θ / 2 D ± ,
where
D ± = λ y λ x ± λ x 2 ( t m 2 + λ y 2 ) + t y 2 λ y 2 i t m λ x t y λ y ,
A ± = 1 E ± t m i t y + i ( λ y + λ x ) D ± ,
B ± = 1 E ± ( t m i t y ) D ± + i λ y λ x ,
C ± = 1 ,
and N ± is a normalisation constant
N ± 2 = t y 2 λ y 2 + t m 2 λ x 2 2 λ y λ x ± λ x 2 t m 2 + λ y 2 + t y 2 λ y 2 2 + t y 2 λ y 2 + t m 2 λ x 2 .
Generally, the relations | A ± | = | D ± | and | B ± | = | C ± | are fulfilled.
The solution for a flat graphene is well known (e.g., [21,22]). Near the center of each valley (the point K or K ) electron dispersion is determined by the Dirac-type Hamiltonian
H f = γ ( τ x k ^ x + τ y k ^ y ) I ,
where again the Pauli matrices τ x , y act on the sublattice degrees of freedom, I is 2 × 2 unity matrix acting in the spin space, with k ^ = i ( / x , / y ). The eigenvalues and eigenstates of the flat graphene Hamiltonian are
E = κ γ k x 2 + k y 2 , κ = ± 1 ,
Ψ k σ ( x , y ) = 1 2 κ e i φ 1 1 σ i e i k · r , σ = ± ,
where e i φ = ( k x i k y ) / k x 2 + k y 2 , k = ( k x , k y ) , r = ( x , y ) , and k = k x 2 + k y 2 . The spin degeneracy is taken into account. In our consideration, the states with the spin up and down are the eigenstates of the operator σ y . The above describe solutions are used to calculate the electron transmission through the corrugated graphene. We recall that, in rippled graphene, the symmetries related to the spin degree of freedom and to the angular momentum are not conserved [13,17].

3. Symmetries

To illuminate specific symmetries of our system, we have to identify the corresponding operators that act within one valley only. Evidently, these operators should act on the A an B sublattices of the honeycomb lattice.

3.1. The Operator S ^ t = τ y i σ y C

The spin–orbit coupling implies that one of the symmetries should be related to the time-reverse symmetry operator T = i σ y C with C the operator of complex conjugation (see, e.g., [23,24]). However, the operator T does not commute neither with the Hamiltonian of the flat graphene in Equation (11) or with the Hamiltonian of the ripple in Equation (1). Taking into account the “pseudospin” degree of freedom, we observe consequently that the operator S ^ t = τ y T commutes with the both Hamiltonians
[ S ^ t , H f ] = 0 , [ S ^ t , H r ] = 0 .
Let us investigate the properties of this operator with respect to the eigenfunctions of the Hamiltonian H f ( H r ) described in Section 2. For the sake of convenience, we introduce the following equivalent definitions: Ψ k σ ( x , y ) r | σ , k , σ = ( + / ) σ = ( / ) . As a result, for the wave function, associated with the flat graphene sheet, we have
S ^ t | ± k = ± e i φ | k S ^ t | ± k = e i φ | k } S t 2 = 1 ,
i.e., the operator S ^ t has two eigenvalues + 1 and 1 . Since the phase e i φ does not affect our results, hereafter, we omit it in our calculations.
Any ket | ψ can be expressed as
| ψ = 1 2 1 + S ^ t | ψ + 1 S ^ t | ψ = | ψ + + | ψ ,
where
| ψ + = 1 + S ^ t 2 | ψ , | ψ = 1 S ^ t 2 | ψ ,
with the property S ^ t | ψ ± = ± | ψ ± . In our particular case, we can form four types of the wave functions:
| ψ ± = 1 ± S ^ t 2 | + k = 1 2 ( | + k ± | k ) ,
| ϕ ± = 1 ± S ^ t 2 | + k = 1 2 ( | + k | k ) .
Thus, for the plane graphene sheet the full set of the operator S ^ t consists of the wave functions in Equations (18) and (19). These wave functions contain the equal mixture of the spin up and down states, associated with electrons that move in opposite directions of our structure.
Since for the curved graphene we use the eigenstates of the CNT, a complete set of quantum numbers consists of the magnetic quantum number m and the wave number k y . In this case, we introduce the equivalent notations Φ m , k y s ( θ ) θ | s , m , k y . The action of the operator S ^ t on the wave function in Equation (5), associated with the nanotube, yields
S ^ t | s = + , ± m , ± k y = i B + * | s = + , m , k y S ^ t | s = , ± m , ± k y = i B * | s = , m , k y } S t 2 = 1 ,
since | B ± | 2 = 1 [see Equation (8)]. Applying the same arguments [see Equations (16) and (17)], we obtain
| χ ± = 1 ± S ^ t 2 | s = + , m , k y = 1 2 ( | s = + , m , k y ± i B + * | s = + , m , k y ) ,
| φ ± = 1 ± S ^ t 2 | s = , m , k y = 1 2 ( | s = , m , k y ± i B * | s = , m , k y ) .
Having the symmetry properties of the wave functions of the different elements of our structure, we are ready to elucidate the symmetry properties of the transmission and reflection probabilities. Equating the wave functions Ψ , Φ at points x = R cos θ 0 (the boundary between Regions I and II), and x = R cos θ 0 (the boundary between Regions II and III), we define the unknown reflection and transmission amplitudes r α β , t α β ( α , β = , ). In these amplitudes, the upper (bottom) index denotes the spin polarisation of the incoming (outgoing) (reflected and transmitted) electron. For the sake of discussion, in our model, electrons move from the left to the right direction. As a result, at the boundary between Regions I and II for the electron, which moves from the left side with the spin up polarisation, we have:
Ψ k x , k y + ( x ) + r L ( ϕ ) Ψ k x , k y + ( x ) + r L ( ϕ ) Ψ k x , k y ( x ) =
a + Φ m + , k y + ( θ ) + b + Φ m + , k y + ( θ ) + a Φ m , k y ( θ ) + b Φ m , k y ( θ ) ; x = R cos θ 0 , θ = ϕ / 2 .
The unknown coefficients a ( + / ) , b ( + / ) characterise transport properties of the electron transfer across the rippled region. The wave numbers m + and m are determined by the equation
m ± = R γ E 2 t y 2 λ y 2 + λ x 2 2 λ x 2 ( E 2 t y 2 ) + λ y 2 t y 2 ,
where E = γ k ( E > 0 ) is the electron energy. At the boundary between Regions II and III, we have the following conditions
a + Φ m + , k y + ( θ ) + b + Φ m + , k y + ( θ ) + a Φ m , k y ( θ ) + b Φ m , k y ( θ ) =
t L ( ϕ ) Ψ k + ( x ) + t L ( ϕ ) Ψ k ( x ) ; x = R cos θ 0 , θ = ϕ / 2 .
The system of Equations (23) and (25) determine the coefficients r L ( ϕ ) α β and t L ( ϕ ) α β . Acting by the operator S ^ t on these equations, we obtain:
  • at the boundary between Region I and II
    Ψ k x , k y ( x ) r L ( ϕ ) * Ψ k x , k y ( x ) + r L ( ϕ ) * Ψ k x , k y + ( x ) = i a + * B + * ( m + , k y ) Φ m + , k y + ( θ ) + i b + * B + * ( m + , k y ) Φ m + , k y + ( θ ) + + i a * B * ( m , k y ) Φ m , k y ( θ ) + i b * B * ( m , k y ) Φ m , k y ( θ ) ; x = R cos θ 0 , θ = ϕ / 2 ;
  • at the boundary between Region II and III
    i a + * B + * ( m + , k y ) Φ m + , k y + ( θ ) + i b + * B + * ( m + , k y ) Φ m + , k y + ( θ ) + + i a * B * ( m , k y ) Φ m , k y ( θ ) + i b * B * ( m , k y ) Φ m , k y ( θ ) = t L ( ϕ ) * Ψ k ( x ) t L ( ϕ ) * Ψ k + ( x ) ; x = R cos θ 0 , θ = ϕ / 2 .
Next, we consider the equations that determine the coefficients r R ( ϕ ) α β and t R ( ϕ ) α β for the electron that moves from the right to the left side with the down spin polarisation. Namely, we have:
  • at the boundary between Regions II and III
    Ψ k x , k y ( x ) + r R ( ϕ ) Ψ k x , k y ( x ) + r R ( ϕ ) Ψ k x , k y + ( x ) = a ˜ + Φ m + , k y + ( θ ) + b ˜ + Φ m + , k y + ( θ ) + a ˜ Φ m , k y ( θ ) + b ˜ Φ m , k y ( θ ) ; x = R cos θ 0 , θ = ϕ / 2 ;
  • at the boundary between Regions I and II
    a ˜ + Φ m + , k y + ( θ ) + b ˜ + Φ m + , k y + ( θ ) + a ˜ Φ m , k y ( θ ) + b ˜ Φ m , k y ( θ ) = t R ( ϕ ) Ψ k x , k y ( x ) + t R ( ϕ ) Ψ k x , k y + ( x ) ; x = R cos θ 0 , θ = ϕ / 2 .
From the comparison of Equations (26) and (27) with Equations (28) and (29), it follows that the coefficients r R ( ϕ ) α β and t R ( ϕ ) α β of Equations (28) and (29) can be expressed in the following form
r R ( ϕ ) = r L ( ϕ ) * ; r R ( ϕ ) = r L ( ϕ ) * ;
and
t R ( ϕ ) = t L ( ϕ ) * ; t R ( ϕ ) = t L ( ϕ ) * .
As a result, we obtain for the reflection probabilities
| r R ( ϕ ) | 2 = | r L ( ϕ ) | 2 ; | r R ( ϕ ) | 2 = | r L ( ϕ ) | 2 ,
while for the transmission probabilities we have
| t R ( ϕ ) | 2 = | t L ( ϕ ) | 2 ; | t R ( ϕ ) | 2 = | t L ( ϕ ) | 2 .
We have similar probabilities for the electron motion with the spin up polarisation from the right to the left side
| r R ( ϕ ) | 2 = | r L ( ϕ ) | 2 ; | r R ( ϕ ) | 2 = | r L ( ϕ ) | 2 ,
and, consequently, for the transmission probabilities
| t R ( ϕ ) | 2 = | t L ( ϕ ) | 2 ; | t R ( ϕ ) | 2 = | t L ( ϕ ) | 2 .
From these results, it follows that the operator S ^ t does not involve the other valley, i.e., QED. However, it interchanges the sign of the vector k and the electron spin polarisation in the flat part of the considered graphene structure. In the rippled graphene region, it changes m ± m ± and the sign of the k y component: k y k y . We conclude that the operator S t acts like a time-reversal operator in the single valley.

3.2. The Operator S ^ c h = τ x σ y

We recall that it was found in Refs. [13,14] that the electron scattering in the superlattice, created by periodically repeated elements, has a curious behaviour. Note that in this way we mimic periodically rippled graphene as a set curvatures between flat graphene areas (e.g., see Figure 1 and the corresponding discussion in Ref. [14]). One element of the superlattice gives rise to the dominance of the electron transmission with a certain spin polarisation. While this effect is small for a few ripples, it defines the perfect transmission for electrons with the one spin polarisation and the perfect reflection for electrons with the opposite spin polarisation in the case of hundreds of ripples. As a result, we find the optimal angle values for the ripple that ensures the perfect transmission at relatively large number of elements N 1 in the superlattice. Note that the transmission depends on the ripple radius and the spin–orbit coupling strengths (see details in [14]). It is noteworthy that transmitted electrons with different spin orientation choose different channels characterised by the quantum number s = ± 1 [13]. This fact implies the existence of the additional symmetry that is fundamental for this feature.
To give an insight into this symmetry, we consider the case k y = 0 , when the spectrum in Equation (4) and the eigenspinors are particular simple. In this case, Equations (3) and (4) are reduced to the form
E = ± ( D ± λ x ) ,
where D = λ y 2 + t m 2 and t m = γ R m . For the rippled (arc) piece (see Figure 1), there are four eigenenergies
E = E 1 = λ x + D E 2 = λ x D E 3 = λ x + D E 4 = λ x D
The connection between the energy and the quantum number m can be formulated in the form
m m s = ± R γ ( s E λ x ) 2 λ y 2 , s = ± 1 ,
that determines four possible values of the quantum number m. Here, we introduce the additional quantum number s that characterises our eigenstates. Note that in Equation (38) the sign of the quantum number s depends on the sign of the energy E. In particular, the following relations take place
E > 0 : s = + 1 , E 1 > 0 s = + 1 , if E 2 > 0 s = 1 , if E 3 > 0
On the other hand, at E < 0 , the branches E 3 and E 4 have the quantum number s = + 1 , while the branch E 2 is characterised by the quantum number s = 1 .
At a fixed electron energy E ( E > 0 ) , Equations (37), (38), and (39), establish the connection between the energies E j ( j = 1 , 2 , 3 ) and the magnetic quantum number m with the quantum number s
E 1 , E 2 ± t m + 1 s = + 1 ,
E 3 ± t m 1 s = 1 .
We recall that the angular momentum is not conserved. As a result, the eigenfunction is the mixture of the eigenfunctions in Equation (5) at a given energy (see Figure 2). The energy branches in Equation (37) with the same quantum number s repel each other, while there is a crossing of the branches with the different s (see Figure 2 and Figure 3). The anticrossings yield the energy gaps = 2 λ y indicated by the arrows (see the insets in Figure 2 and Figure 3). As a result, the energy gaps give rise to evanescent modes at energies λ x λ y < | E | < λ x + λ y in our system.
At k y = 0 , the eigenvectors in Equation (5) of the Hamiltonian in Equation (1) transform to the following form (see details in [13,17]):
Φ j , m ( θ ) = I e i σ y 2 θ ν j ( m ) e i m θ , j = 1 , . . 4 .
Here
ν 1 = i t m A i A t m , ν 2 = i t m A + i A + t m , ν 3 = i t m A + i A + t m , ν 4 = i t m A i A t m ,
where
A ± = ± λ y + t m 2 + λ y 2 .
The choice of the components of the wave function in Equation (42) depends on the energy interval in Equation (37), available for electrons. For example (see Figure 2), for E > λ x E = E 1 E = E 3 , the base eigenfunctions are: exp ( ± i m + 1 θ ) ν 1 ( ± t m + 1 ) ; exp ( ± i m 1 θ ) ν 3 ( ± t m 1 ) , respectively. As discussed above, the quantum number s, characterising the components of the wave function in Equation (42), is associated with the specific energy (see Equations (40) and (41)). The wave function in Equation (13) of the flat graphene at k y = 0 is convenient to present in the form:
Ψ k σ ( x ) = 1 2 τ k | k | 1 1 σ i e i k x , σ = ± .
As above, we use the following definitions: Ψ k σ ( x ) x | σ , k , σ = ( + / ) σ = ( / ) . The symbol τ = ( + ) is ascribed to the valence (conduction) band. We study the case E 0 ; the opposite case can be analysed with the same method. As above, the positive value of the wave number k k x corresponds to the direction of the electron motion from the left to the right side of the considered system.
In the case k y = 0 , there is one more operator S ^ c h = τ x σ y that commutes with both Hamiltonians in Equations (1) and (11), describing the flat and the rippled pieces of our system, respectively:
[ S ^ c h , H f ] = 0 , [ S ^ c h , H r ] = 0 .
For the wave function, associated with the flat graphene sheet, we have
S ^ c h | ± k = ± | ± k S ^ c h | ± k = | ± k } α = ± 1 .
Thus, the operator S ^ c h has the eigenvalue α = + 1 ( 1 ) , acting on the wave function of the electron, traveling from the left side of our structure with spin up (down) polarisation. The eigenvalues are equal in value but opposite in sign if the operator S ^ c h acts on the wave function of the electron, traveling from the right side of our structure with spin up (down) polarisation. On the other hand, the action of the operator S ^ c h on the electron wave functions, associated with the rippled piece of our system, is as follows:
S ^ c h Φ 1 , m ( θ ) = Φ 1 , m ( θ ) ,
S ^ c h Φ 2 , m ( θ ) = Φ 2 , m ( θ ) ,
S ^ c h Φ 3 , m ( θ ) = + Φ 3 , m ( θ ) ,
S ^ c h Φ 4 , m ( θ ) = + Φ 4 , m ( θ ) .
In this case, the operator S ^ c h has the eigenvalue α = + 1 ( 1 ) as well, providing us the connection between the eigenstates of the flat and rippled pieces of the considered system.
Now, we are ready to discuss the list of eigenfunctions, responsible for the transport properties in different energy range.
  • E > λ x .
    1.
    Quantum numbers: α = + 1 , s = 1 . These quantum numbers determine the available set of the wave functions: Ψ k + ( x ) , Ψ k ( x ) , Φ 3 , m ( θ ) . If the corresponding symmetries are responsible for the transport properties, there are only the following options.
    (a)
    The electron is moving from the left side of our structure (flat graphene sheet) with the spin up polarisation [ Ψ k + ( x ) | + k ]. In this case, in the rippled graphene region (see Figure 2), there is one open channel, defined by the wave function Φ 3 , m . The wave function Ψ k ( x ) | k describes the reflection with the electron spin–flip.
    (b)
    The electron is moving from the right side (flat graphene sheet) with the spin down polarisation [ Ψ k ( x ) | k ]. In this case, in the rippled graphene region, there is only the transmission channel, defined by the wave function Φ 3 , m . The wave function Ψ k + ( x ) | + k describes the reflection with the electron spin–flip.
    As a result, we expect the equivalence between the left/right transmission probabilities with the opposite spin polarisations. Indeed, this expectation is consistent with Equation (33), obtained from the different arguments at k y = 0 .
    2.
    Quantum numbers: α = 1 , s = + 1 . The available set of the wave functions: Ψ k ( x ) , Ψ k + ( x ) , Φ 1 , m ( θ ) . In this case, the symmetries dictate the following options.
    (a)
    The electron is moving from the left side (flat graphene sheet) with the spin down polarisation [ Ψ k ( x ) | + k ]. In this case, in the rippled graphene region (see Figure 2), there is one open channel, defined by the wave function Φ 1 , m ( θ ) . The wave function Ψ k + ( x ) | k describes the reflections with the electron spin–flip.
    (b)
    The electron is moving from the right side (flat graphene sheet) with the spin up polarisation [ Ψ k + ( x ) | k ]. In this case, in the rippled graphene region for this electron there is only the transmission channel, defined by the wave function Φ 1 , m ( θ ) . The wave function Ψ k ( x ) | + k describes the reflections with the electron spin–flip.
    Again, we expect the equivalence between the left/right transmission probabilities with the opposite spin polarisations. Indeed, this expectation is consistent with Equation (35), obtained from different arguments at k y = 0 .
  • E < λ x .
    1.
    Quantum numbers: α = 1 , s = + 1 . In this case, the available set includes the following wave functions: Ψ k + ( x ) , Ψ k ( x ) , Φ 2 , m ( θ ) . The symmetries dictate the following options.
    (a)
    The electron is moving from the left side (flat graphene sheet) with the spin down polarisation [ Ψ k ( x ) | + k ]. In the rippled graphene region, there is only the transmission channel, defined by the wave function Φ 2 , m ( θ ) . The wave function Ψ k + ( x ) | k describes the reflections with the electron spin–flip.
    (b)
    The electron is moving from the right side (flat graphene sheet) with the spin up polarisation [ Ψ k + ( x ) | k ]. In the rippled graphene region, there is only the transmission channel, defined by the wave function Φ 2 , m ( θ ) . The wave function Ψ k ( x ) | + k describes the reflections with the electron spin–flip.
    Again, the expected equivalence between the left/right transmission probabilities with the opposite spin polarisations is consistent with Equation (35), obtained from different arguments at k y = 0 .
    2.
    Quantum numbers: α = + 1 , s = 1 . In this case, the available set includes the following wave functions: Ψ k + ( x ) , Ψ k ( x ) , Φ 3 , m ( θ ) . This situation is completely equivalent to the case discussed at E > λ x , Point 1.
    Thus, at k y = 0 , the symmetry, associated with the operator S ^ c h , determines the following transport properties through the rippled graphene piece: (i) at the transmission, it preserves the electron spin polarisation, while forbids the spin–flip; and (ii) the reflection occurs only with the spin–flip.
For k y 0 , the operator S ^ c h does not commute with the Hamiltonians. In this case, the discussed symmetry is broken. It results in constraint release on the reflections and transmissions mechanisms in our system (i.e., [25]).

3.3. The Relation Between the Operators S ^ t and S ^ c h

Note that these two symmetry operator commute if k y = 0 . Evidently, they have the common basis, while having different eigenvalues. Let us analyse this situation in details. One can readily see that the eigenfunctions in Equations (18) and (19), associated with the flat graphene piece, are common eigenstates for the both operators:
S ^ t | ψ ± = ± | ψ ± , S ^ t | ϕ ± = ± | ϕ ± ,
S ^ c h | ψ ± = + | ψ ± , S ^ c h | ϕ ± = | ϕ ± .
For the rippled graphene piece, we have to consider only the case k y = 0 . In this case, it is convenient to construct the common basis from the set in Equation (42). Taking into account the properties in Equations (48)–(51), we introduce the following superpositions
Φ j = 1 2 [ Φ j , m ( θ ) + Φ j , m ( θ ) ] , j = 1 , . . , 4 .
As a result, we obtain
S ^ t Φ 1 , 2 = Φ 1 , 2 , S ^ t Φ 3 , 4 = Φ 3 , 4 ,
S ^ c h Φ 1 , 2 = Φ 1 , 2 , S ^ c h Φ 3 , 4 = Φ 3 , 4 .
Thus, the eigenfuctions in Equation (54) form the complete set for the both symmetry operators in the case of the rippled graphene piece.

4. Summary

Evidently, symmetries play an essential role in our understanding different phenomena in mesoscopic physics. In graphene physics, they become especially apparent in the transport properties of the corrugated systems. We find two symmetry operators that explain the chiral behaviour of the ballistic electron transport through the rippled graphene. This unusual behaviour has emerged due to the curvature-induced spin–orbit coupling. In particular, the symmetry operator S ^ t (see Section 3.1) elucidates the equivalence between the transport characteristics of the ballistic electrons travelling from opposite sides of our system that have different type of polarisations. This operator acts as a time-reversal operator in the single valley system, considered in our paper. The other symmetry operator S ^ c h (see Section 3.2) enables us to explain the selection of open energy channels for the ballistic electrons travelling through the rippled graphene subsystem at the direct incident of the electron flow ( k y = 0 ). This symmetry explains the dominance of different electron spin polarisations that depend on the direction of the electron flow. This selection becomes increasingly important at multiple periodic repetition of the corrugated graphene structure, considered in our paper (see also [14]). From our preliminary analysis, it follows that similar symmetry preserves if we consider the down side of the CNT. This problem is, however, beyond the scope of the present studies and will be discussed in detail in forthcoming paper.

Author Contributions

All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partially supported by the Slovak Academy of Sciences in the framework of VEGA Grant No. 2/0009/19, and by Votruba-Blokhintsev program BLTP, JINR.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumar, S.; Parks, D.M. Strain Shielding from Mechanically Activated Covalent Bond Formation during Nanoindentation of Graphene Delays the Onset of Failure. Nano Lett. 2015, 15, 1503–1510. [Google Scholar] [CrossRef] [Green Version]
  2. Cao, K.; Feng, S.; Han, Y.; Gao, L.; Hue, L.T.; Xu, Z.; Lu, Y. Elastic straining of free-standing monolayer graphene. Nat. Commun. 2020, 11, 284. [Google Scholar] [CrossRef] [Green Version]
  3. Alyobi, M.M.; Barnett, C.J.; Rees, P.; Cobley, R.J. Modifying the electrical properties of graphene by reversible point-ripple formation. Carbon 2019, 143, 762–768. [Google Scholar] [CrossRef]
  4. Vázquez de Parga, A.L.; Calleja, F.; Borca, B.; Passeggi, M.C.G.; Hinarejos, J.J.; Guinea, F.; Miranda, R. Periodically Rippled Graphene: Growth and Spatially Resolved Electronic Structure. Phys. Rev. Lett. 2008, 100, 056807. [Google Scholar] [CrossRef] [Green Version]
  5. Maccariello, D.; Al Taleb, A.; Calleja, F.; Vázquez de Parga, A.L.; Perna, P.; Camarero, J.; Gnecco, E.; Farías, D.; Miranda, R. Observation of Localized Vibrational Modes of Graphene Nanodomes by Inelastic Atom Scattering. Nano Lett. 2016, 16, 2–7. [Google Scholar] [CrossRef]
  6. Politano, A.; Chiarello, G. Plasmon modes in graphene: Status and prospect. Nanoscale 2014, 6, 10927–10940. [Google Scholar] [CrossRef]
  7. Ni, G.X.; Zheng, Y.; Bae, S.; Kim, H.R.; Pachoud, A.; Kim, Y.S.; Tan, C.L.; Im, D.; Ahn, J.H.; Hong, B.H.; et al. Quasi-Periodic Nanoripples in Graphene Grown by Chemical Vapor Deposition and Its Impact on Charge Transport. ACS Nano 2012, 6, 1158–1164. [Google Scholar] [CrossRef] [Green Version]
  8. Vasić, B.; Zurutuza, A.; Gajić, R. Spatial variation of wear and electrical properties across wrinkles in chemical vapour deposition graphene. Carbon 2016, 102, 304–310. [Google Scholar] [CrossRef]
  9. Pudlak, M.; Nazmitdinov, R.G. Klein collimation by rippled graphene superlattice. J. Phys. Cond. Matter 2019, 31, 495301. [Google Scholar] [CrossRef] [Green Version]
  10. Guinea, F.; Katsnelson, M.I.; Vozmediano, M.A.H. Midgap states and charge inhomogeneities in corrugated graphene. Phys. Rev. B 2008, 77, 075422. [Google Scholar] [CrossRef] [Green Version]
  11. Katsnelson, M.; Geim, A. Electron scattering on microscopic corrugations in graphene. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2008, 366, 195–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Ando, T. Spin-Orbit Interaction in Carbon Nanotubes. J. Phys. Soc. Jpn. 2000, 69, 1757–1763. [Google Scholar] [CrossRef]
  13. Pudlak, M.; Pichugin, K.N.; Nazmitdinov, R.G. Cooperative phenomenon in a rippled graphene: Chiral spin guide. Phys. Rev. B 2015, 92, 205432. [Google Scholar] [CrossRef] [Green Version]
  14. Smotlacha, J.; Pudlak, M.; Nazmitdinov, R.G. Spin transport in a rippled graphene periodic chain. J. Phys. Conf. Ser. 2019, 1416, 012035. [Google Scholar] [CrossRef] [Green Version]
  15. Izumida, W.; Sato, K.; Saito, R. Spin–Orbit Interaction in Single Wall Carbon Nanotubes: Symmetry Adapted Tight-Binding Calculation and Effective Model Analysis. J. Phys. Soc. Jpn. 2009, 78, 074707. [Google Scholar] [CrossRef]
  16. Del Valle, M.; Margańska, M.; Grifoni, M. Signatures of spin-orbit interaction in transport properties of finite carbon nanotubes in a parallel magnetic field. Phys. Rev. B 2011, 84, 165427. [Google Scholar] [CrossRef] [Green Version]
  17. Pichugin, K.N.; Pudlak, M.; Nazmitdinov, R.G. Spin-orbit effects in carbon nanotubes—Analytical results. Eur. Phys. J. B 2014, 87, 124. [Google Scholar] [CrossRef] [Green Version]
  18. Michaeli, K.; Kantor-Uriel, N.; Naaman, R.; Waldeck, D.H. The electron’s spin and molecular chirality—How are they related and how do they affect life processes? Chem. Soc. Rev. 2016, 45, 6478–6487. [Google Scholar] [CrossRef]
  19. Nazmitdinov, R.G. From Chaos to Order in Mesoscopic Systems. Phys. Part. Nucl. Letts. 2019, 16, 159–169. [Google Scholar] [CrossRef]
  20. Birman, J.; Nazmitdinov, R.; Yukalov, V. Effects of symmetry breaking in finite quantum systems. Phys. Rep. 2013, 526, 1–91. [Google Scholar] [CrossRef] [Green Version]
  21. Foa Torres, L.E.F.; Roche, S.; Charlier, J.C. Introduction to Graphene-Based Nanomaterials: From Electronic Structure to Quantum Transport; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar] [CrossRef]
  22. Katsnelson, M.I. Graphene: Carbon in Two Dimensions; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar] [CrossRef]
  23. Sakurai, J.J. Modern Quantum Mechanics, Revised Edition; Addison-Wesley: Reading, MA, USA, 1994. [Google Scholar]
  24. Bohr, A.; Mottelson, B. Nuclear Structure; Number v. 1 in Nuclear Structure; World Scientific: Singapore, 1998. [Google Scholar]
  25. Busa, J.; Pudlak, M.; Nazmitdinov, R.G. On Electron Scattering through a Single Corrugated Graphene Sructure. Phys. Part. Nucl. Lett. 2019, 16, 729–733. [Google Scholar] [CrossRef]
Figure 1. The corrugated graphene system. There are two flat surfaces: Region I, defined in the intervals < x < R cos θ 0 ; and Region III, defined in the intervals R cos θ 0 < x < . A ripple is modelled by an arc of a circle (Region II) of radius R, defined as R cos θ 0 < x < R cos θ 0 . At θ 0 = 0 , the ripple is a half of the nanotube, while at θ 0 = π / 2 the ripple does not exist. The angle ϕ = π 2 θ 0 . Here, we have < y < . We keep the translational invariance along the y-axis, which is chosen as the symmetry and the quantisation axis.
Figure 1. The corrugated graphene system. There are two flat surfaces: Region I, defined in the intervals < x < R cos θ 0 ; and Region III, defined in the intervals R cos θ 0 < x < . A ripple is modelled by an arc of a circle (Region II) of radius R, defined as R cos θ 0 < x < R cos θ 0 . At θ 0 = 0 , the ripple is a half of the nanotube, while at θ 0 = π / 2 the ripple does not exist. The angle ϕ = π 2 θ 0 . Here, we have < y < . We keep the translational invariance along the y-axis, which is chosen as the symmetry and the quantisation axis.
Symmetry 12 00533 g001
Figure 2. The spectrum in Equation (4) ( k y = 0 ) versus the magnetic quantum number m. The non-quantised values ± m s = ± 1 at the energy E = 0.2 eV (thin horizontal line that mimics the Fermi energy) are indicated at the crossing of the energy branches with different s. Symbols E 1 , E 2 , E 3 , E 4 are used to guide the eyes on the formal solutions (straight lines) defined by Equation (36), irrespective of the sign of the quantum number s. In contrast, there is anticrossing of the energy branches with the same quantum number s = + 1 for the pair ( E 1 , E 2 ) at E > 0 . Similar anticrossing occurs at E < 0 , when the pair ( E 3 , E 4 ) has the same quantum number s = + 1 . These anticrossings are caused by the term λ y in the Hamiltonian in Equation (1), which creates the energy gaps 2 λ y near the energy E = ± λ x (see [13,17]). The following parameters are used: R = 10 Å, δ = 0.01 , p = 0.1 , γ = ( 4.5 · 1.42 ) eV · Å, γ = 8 3 γ . These parameters define the values of the spin–orbit coupling strengths: λ y = δ γ / 4 R = 0.0043 eV, λ x = γ ( 1 / 2 + 2 δ p ) / R = 0.32 eV (see Section 2).
Figure 2. The spectrum in Equation (4) ( k y = 0 ) versus the magnetic quantum number m. The non-quantised values ± m s = ± 1 at the energy E = 0.2 eV (thin horizontal line that mimics the Fermi energy) are indicated at the crossing of the energy branches with different s. Symbols E 1 , E 2 , E 3 , E 4 are used to guide the eyes on the formal solutions (straight lines) defined by Equation (36), irrespective of the sign of the quantum number s. In contrast, there is anticrossing of the energy branches with the same quantum number s = + 1 for the pair ( E 1 , E 2 ) at E > 0 . Similar anticrossing occurs at E < 0 , when the pair ( E 3 , E 4 ) has the same quantum number s = + 1 . These anticrossings are caused by the term λ y in the Hamiltonian in Equation (1), which creates the energy gaps 2 λ y near the energy E = ± λ x (see [13,17]). The following parameters are used: R = 10 Å, δ = 0.01 , p = 0.1 , γ = ( 4.5 · 1.42 ) eV · Å, γ = 8 3 γ . These parameters define the values of the spin–orbit coupling strengths: λ y = δ γ / 4 R = 0.0043 eV, λ x = γ ( 1 / 2 + 2 δ p ) / R = 0.32 eV (see Section 2).
Symmetry 12 00533 g002
Figure 3. The same as in Figure 2. Solid lines are associated with the positive energy states ( E > 0 ), while the negative energies are denoted by dashed lines (see text). Once the energy changes the sign, it affects the sign of the corresponding quantum number s. There are anticrossing of the energy branches with the same quantum number s = + 1 for the pair ( E 1 , E 2 ) at E > 0 and for the pair ( E 3 , E 4 ) at E < 0 .
Figure 3. The same as in Figure 2. Solid lines are associated with the positive energy states ( E > 0 ), while the negative energies are denoted by dashed lines (see text). Once the energy changes the sign, it affects the sign of the corresponding quantum number s. There are anticrossing of the energy branches with the same quantum number s = + 1 for the pair ( E 1 , E 2 ) at E > 0 and for the pair ( E 3 , E 4 ) at E < 0 .
Symmetry 12 00533 g003

Share and Cite

MDPI and ACS Style

Pudlak, M.; Smotlacha, J.; Nazmitdinov, R. On Symmetry Properties of The Corrugated Graphene System. Symmetry 2020, 12, 533. https://doi.org/10.3390/sym12040533

AMA Style

Pudlak M, Smotlacha J, Nazmitdinov R. On Symmetry Properties of The Corrugated Graphene System. Symmetry. 2020; 12(4):533. https://doi.org/10.3390/sym12040533

Chicago/Turabian Style

Pudlak, Mihal, Jan Smotlacha, and Rashid Nazmitdinov. 2020. "On Symmetry Properties of The Corrugated Graphene System" Symmetry 12, no. 4: 533. https://doi.org/10.3390/sym12040533

APA Style

Pudlak, M., Smotlacha, J., & Nazmitdinov, R. (2020). On Symmetry Properties of The Corrugated Graphene System. Symmetry, 12(4), 533. https://doi.org/10.3390/sym12040533

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop