*Article* **On Symmetry Properties of The Corrugated Graphene System**

#### **Mihal Pudlak <sup>1</sup> , Jan Smotlacha <sup>2</sup> and Rashid Nazmitdinov 2,3,\***


Received: 20 February 2020; Accepted: 24 March 2020; Published: 3 April 2020

**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.

**Keywords:** graphene; ripple; transport; symmetry

#### **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<sup>z</sup>* orbital orientation, and is leading to shift of the on-site energies of *p<sup>z</sup>* 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–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).

**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 *θ*<sup>0</sup> < *x* < ∞. A ripple is modelled by an arc of a circle (Region II) of radius *R*, defined as −*R* cos *θ*<sup>0</sup> < *x* < *R* cos *θ*0. At *θ*<sup>0</sup> = 0, the ripple is a half of the nanotube, while at *θ*<sup>0</sup> = *π*/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.

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* 0 point.

The Hamiltonian of the nanotube has the following form in the effective mass approximation [12]

$$H\_{\mathbf{r}} = \gamma \left(\pi\_{\mathbf{x}} \mathfrak{R}\_{\mathbf{x}} + \pi\_{\mathbf{y}} \mathfrak{R}\_{\mathbf{y}}\right) \otimes I - \lambda\_{\mathbf{y}} \pi\_{\mathbf{y}} \otimes \sigma\_{\mathbf{x}}(\vec{r}) - \mathfrak{J}\_{\mathbf{x}} \pi\_{\mathbf{x}} \otimes \sigma\_{\mathbf{y}} \,. \tag{1}$$

Here, the operators *π*ˆ *<sup>x</sup>* = −*i ∂ <sup>R</sup>∂θ* , *π*ˆ *<sup>y</sup>* = −*i ∂ ∂y* , *σx*(~*r*) = *σ<sup>x</sup>* cos *θ* − *σ<sup>z</sup>* sin *θ*, and *ξ<sup>x</sup>* = 2*δγp*/*R*, *λ<sup>y</sup>* = *δγ*0/4*R*. The Pauli matrices *σx*,*<sup>y</sup>* act in the spin space. The matrices *τx*,*<sup>y</sup>* act on the sublattice degree of freedom. The Pauli matrix *τ<sup>i</sup>* is called the "pseudospin", to distinguish it from the real electron spin.

The following notations are used: *γ* = − √ 3*V π ppa*/2 = √ 3*γ*0*a*/2, *γ* 0 = √ 3(*V σ pp* − *V π pp*)*a*/2 = *γ*1*a*, *p* = 1 − 3*γ* 0/8*γ*; where *V σ pp* and *V π pp* 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 *γ*<sup>0</sup> = −*V π pp* ≈ 3 eV and *γ*<sup>1</sup> ≈ 8 eV (see, e.g., [12]). The intrinsic source of the spin–orbit coupling *δ* = ∆/3*eπσ* is defined by means of the quantity

$$
\Delta = i \frac{\Im \hbar}{4m^2 c^2} \left< \mathbf{x} \middle| \frac{\partial V}{\partial \mathbf{x}} p\_y - \frac{\partial V}{\partial y} p\_\mathbf{x} \middle| y \right> , \tag{2}
$$

where *V* is the atomic potential and *eπσ* = *e π* <sup>2</sup>*<sup>p</sup>* − *e σ* 2*p* . Here, the energies *e π* 2*p* and *e σ* 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 = \kappa E\_{\pm \nu} \quad \kappa = \pm 1,\tag{3}$$

where *κ* = +1(−1) is associated with the conductance (valence) band, and the energies *E*<sup>±</sup> are defined as

$$E\_{\pm} = \sqrt{t\_m^2 + t\_y^2 + \lambda\_y^2 + \lambda\_x^2 \pm 2\sqrt{\lambda\_x^2 \left(t\_m^2 + \lambda\_y^2\right) + t\_y^2 \lambda\_y^2}}.\tag{4}$$

Here, *t<sup>m</sup>* = *mγ*/*R*, *t<sup>y</sup>* = *γky*, *λ<sup>x</sup>* = *γ*(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

$$\boldsymbol{\Phi}\_{m\boldsymbol{k}y}^{s=\pm 1}(\boldsymbol{\theta}, \boldsymbol{y}) = e^{im\boldsymbol{\theta}} e^{i\boldsymbol{k}\_y y} \boldsymbol{N}\_{\pm} \begin{bmatrix} \kappa (\cos \theta / 2A\_{\pm} - \sin \theta / 2B\_{\pm}) \\ \kappa (\sin \theta / 2A\_{\pm} + \cos \theta / 2B\_{\pm}) \\ \cos \theta / 2\mathbb{C}\_{\pm} - \sin \theta / 2D\_{\pm} \\ \sin \theta / 2\mathbb{C}\_{\pm} + \cos \theta / 2D\_{\pm} \end{bmatrix} \tag{5}$$

where

$$D\_{\pm} = \frac{\lambda\_y \lambda\_x \pm \sqrt{\lambda\_x^2 (t\_m^2 + \lambda\_y^2) + t\_y^2 \lambda\_y^2}}{i t\_m \lambda\_x - t\_y \lambda\_y},\tag{6}$$

$$A\_{\pm} = \frac{1}{E\_{\pm}} \left( t\_m - it\_y + i(\lambda\_y + \lambda\_x)D\_{\pm} \right), \tag{7}$$

$$B\_{\pm} = \frac{1}{E\_{\pm}} \left[ (t\_{\mathcal{W}} - it\_{\mathcal{Y}}) D\_{\pm} + i \left( \lambda\_{\mathcal{Y}} - \lambda\_{\mathcal{X}} \right) \right],\tag{8}$$

$$\mathbb{C}\_{\pm} = 1,\tag{9}$$

and *N*<sup>±</sup> is a normalisation constant

$$\mathcal{N}\_{\pm}^{2} = \frac{t\_{y}^{2}\lambda\_{y}^{2} + t\_{m}^{2}\lambda\_{x}^{2}}{2\left[\left(\lambda\_{y}\lambda\_{x} \pm \sqrt{\lambda\_{x}^{2}\left(t\_{m}^{2} + \lambda\_{y}^{2}\right) + t\_{y}^{2}\lambda\_{y}^{2}}\right)^{2} + t\_{y}^{2}\lambda\_{y}^{2} + t\_{m}^{2}\lambda\_{x}^{2}\right]}. \tag{10}$$

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* 0 ) electron dispersion is determined by the Dirac-type Hamiltonian

$$H\_f = \gamma(\tau\_\mathbf{x}\hat{k}\_\mathbf{x} + \tau\_\mathbf{y}\hat{k}\_\mathbf{y}) \otimes I,\tag{11}$$

where again the Pauli matrices *τx*,*<sup>y</sup>* act on the sublattice degrees of freedom, *I* is 2 × 2 unity matrix acting in the spin space, with *<sup>k</sup>***<sup>ˆ</sup>** <sup>=</sup> <sup>−</sup>*i*(*∂*/*∂x*, *<sup>∂</sup>*/*∂y*). The eigenvalues and eigenstates of the flat graphene Hamiltonian are

$$E = \kappa \gamma \sqrt{k\_x^2 + k\_{y'}^2}, \kappa = \pm 1,\tag{12}$$

$$\Psi\_{\mathbf{k}}^{\sigma}(\mathbf{x}, \mathbf{y}) = \frac{1}{2} \begin{pmatrix} \kappa e^{-i\varphi} \\ 1 \end{pmatrix} \otimes \begin{pmatrix} 1 \\ \sigma i \end{pmatrix} e^{i\mathbf{k} \cdot \mathbf{r}}, \qquad \sigma = \pm \,\,\,\tag{13}$$

where *e* <sup>−</sup>*i<sup>ϕ</sup>* = (*k<sup>x</sup>* <sup>−</sup> *iky*)/ q *k* 2 *<sup>x</sup>* + *k* 2 *y* , *k* = (*kx*, *ky*), *r* = (*x*, *y*), and *k* = q *k* 2 *<sup>x</sup>* + *k* 2 *y* . 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*ˆ *<sup>t</sup>* = *τ<sup>y</sup>* ⊗ *iσyC*

The spin–orbit coupling implies that one of the symmetries should be related to the time-reverse symmetry operator T = *iσyC* 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*ˆ *<sup>t</sup>* = *τ<sup>y</sup>* ⊗ T commutes with the both Hamiltonians

$$\left[\hat{\mathbb{S}}\_{t\prime}H\_f\right] = \mathbf{0} \,, \quad \left[\hat{\mathbb{S}}\_{t\prime}H\_r\right] = \mathbf{0} \,. \tag{14}$$

Let us investigate the properties of this operator with respect to the eigenfunctions of the Hamiltonian *H<sup>f</sup>* ( *Hr*) described in Section 2. For the sake of convenience, we introduce the following equivalent definitions: Ψ*<sup>σ</sup> k* (*x*, *y*) ≡ h*r*|*σ*, *k*i, *σ* = (+/−) ⇔ *σ* = (↑ / ↓). As a result, for the wave function, associated with the flat graphene sheet, we have

$$
\begin{aligned}
\langle \hat{S}\_t | \uparrow \pm k \rangle &=& \pm e^{i\varphi} | \downarrow \mp k \rangle \\
\langle \hat{S}\_t | \downarrow \pm k \rangle &=& \mp e^{i\varphi} | \uparrow \mp k \rangle
\end{aligned}
\quad \Rightarrow S\_t^2 = 1,\tag{15}
$$

i.e., the operator *S*ˆ *<sup>t</sup>* has two eigenvalues +1 and −1. Since the phase *e <sup>i</sup><sup>ϕ</sup>* does not affect our results, hereafter, we omit it in our calculations.

Any ket |*ψ*i can be expressed as

$$\left|\psi\right> = \frac{1}{2}\left[\left(1 + \hat{S}\_{l}\right)\left|\psi\right> + \left(1 - \hat{S}\_{l}\right)\left|\psi\right>\right] = \left|\psi\_{+}\right> + \left|\psi\_{-}\right>\,\tag{16}$$

where

$$|\psi\_{+}\rangle = \frac{1 + \mathcal{S}\_{l}}{2}|\psi\rangle , \quad |\psi\_{-}\rangle = \frac{1 - \mathcal{S}\_{l}}{2}|\psi\rangle , \tag{17}$$

with the property *S*ˆ *t* |*ψ*±i = ±|*ψ*±i. In our particular case, we can form four types of the wave functions:

$$|\psi\_{\pm}\rangle = \frac{1 \pm \hat{S}\_{l}}{2} |\uparrow + k\rangle = \frac{1}{2} (|\uparrow + k\rangle \pm |\downarrow - k\rangle) ,\tag{18}$$

$$\left|\phi\_{\pm}\right\rangle = \frac{1 \pm \hat{S}\_{\text{f}}}{2} \left|\downarrow + k\right\rangle = \frac{1}{2} (\left|\downarrow + k\right\rangle \mp \left|\uparrow - k\right\rangle) \,. \tag{19}$$

Thus, for the plane graphene sheet the full set of the operator *S*ˆ *<sup>t</sup>* 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 *ky*. In this case, we introduce the equivalent notations Φ*<sup>s</sup> m*,*ky* (*θ*) ≡ h*θ*|*s*, *<sup>m</sup>*, *<sup>k</sup>y*i. The action of the operator *<sup>S</sup>*<sup>ˆ</sup> *<sup>t</sup>* on the wave function in Equation (5), associated with the nanotube, yields

$$\begin{aligned} \left| \hat{S}\_{t} \right| \left| s=+, \pm m, \pm k\_{y} \right\rangle &= \left| iB\_{+}^{\*} \right| \left| s=+, \mp m, \mp k\_{y} \right\rangle \\ \left| \hat{S}\_{t} \right| \left| s=-, \pm m, \pm k\_{y} \right\rangle &= \left| iB\_{-}^{\*} \right| \left| s=-, \mp m, \mp k\_{y} \right\rangle \end{aligned} \right\} \Rightarrow S\_{t}^{2} = 1 \,,\tag{20}$$

since |*B*±| <sup>2</sup> = 1 [see Equation (8)]. Applying the same arguments [see Equations (16) and (17)], we obtain

$$|\chi\_{\pm}\rangle = \frac{1 \pm \hat{S}\_{l}}{2}|\mathbf{s} = + \,, m\_{\prime}k\_{\underline{y}}\rangle = \frac{1}{2}(|\mathbf{s} = + \,, m\_{\prime}k\_{\underline{y}}\rangle \pm i\mathcal{B}\_{+}^{\*}|\mathbf{s} = + \,, -m\_{\prime} - k\_{\underline{y}}\rangle) , \tag{21}$$

$$|\langle \boldsymbol{\varrho} \pm \rangle = \frac{1 \pm \hat{S}\_l}{2} |\mathbf{s} = -\,, m, k\_y \rangle = \frac{1}{2} (|\mathbf{s} = -\,, m, k\_y \rangle \pm i \mathcal{B}\_-^\* |\mathbf{s} = -\,, -m, -k\_y \rangle) \,. \tag{22}$$

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 *θ*<sup>0</sup> (the boundary between Regions I and II), and *x* = *R* cos *θ*<sup>0</sup> (the boundary between Regions II and III), we define the unknown reflection and transmission amplitudes *r β <sup>α</sup>*, *t β <sup>α</sup>* (*α*, *β* = ↑, ↓). 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:

$$\mathbf{Y}\_{k\_x k\_y}^+(\mathbf{x}) + r\_L(\boldsymbol{\phi})\_\uparrow^\uparrow \mathbf{Y}\_{-k\_x k\_y}^+(\mathbf{x}) + r\_L(\boldsymbol{\phi})\_\downarrow^\uparrow \mathbf{Y}\_{-k\_x k\_y}^-(\mathbf{x}) = \tag{23}$$

$$a\_+\Phi^+\_{m+k\_\mathcal{Y}}(\theta) + b\_+\Phi^+\_{-m+k\_\mathcal{Y}}(\theta) + a\_-\Phi^-\_{m-k\_\mathcal{Y}}(\theta) + b\_-\Phi^-\_{-m-k\_\mathcal{Y}}(\theta); \quad \mathbf{x} = -\mathcal{R}\cos\theta\_0, \theta = -\phi/2.$$

The unknown coefficients *<sup>a</sup>*(+/−) , *<sup>b</sup>*(+/−) characterise transport properties of the electron transfer across the rippled region. The wave numbers *m*<sup>+</sup> and *m*<sup>−</sup> are determined by the equation

$$m\_{\pm} = \frac{R}{\gamma} \sqrt{E^2 - t\_y^2 - \lambda\_y^2 + \lambda\_x^2 \mp 2\sqrt{\lambda\_x^2(E^2 - t\_y^2) + \lambda\_y^2 t\_y^2}}\tag{24}$$

where *E* = *γk* (*E* > 0) is the electron energy. At the boundary between Regions II and III, we have the following conditions

$$a\_+\Phi^+\_{m\_+,k\_y}(\theta) + b\_+\Phi^+\_{-m\_+,k\_y}(\theta) + a\_-\Phi^-\_{m\_-,k\_y}(\theta) + b\_-\Phi^-\_{-m\_-,k\_y}(\theta) = \tag{25}$$

$$t\_L(\phi)^\uparrow\_\uparrow \Psi^+\_k(\mathbf{x}) + t\_L(\phi)^\uparrow\_\downarrow \Psi^-\_k(\mathbf{x});\ \mathbf{x} = \mathcal{R}\cos\theta\_{0\prime}, \ \theta = \phi/2.$$

The system of Equations (23) and (25) determine the coefficients *rL*(*φ*) *β <sup>α</sup>* and *tL*(*φ*) *β <sup>α</sup>*. Acting by the operator *S*ˆ *<sup>t</sup>* on these equations, we obtain:

• at the boundary between Region I and II

$$\begin{split} \Psi^{-}\_{-k\_{x}-k\_{y}}(\mathbf{x}) - r\_{L}(\boldsymbol{\phi})^{\uparrow\ast}\_{\uparrow} \Psi^{-}\_{k\_{x}-k\_{y}}(\mathbf{x}) + r\_{L}(\boldsymbol{\phi})^{\uparrow\ast}\_{\downarrow} \Psi^{+}\_{k\_{x}-k\_{y}}(\mathbf{x}) &= \\ \mathrm{i}a^{\ast}\_{+} \boldsymbol{B}^{\ast}\_{+}(m\_{+},k\_{y}) \boldsymbol{\Phi}^{+}\_{-m\_{+},-k\_{y}}(\boldsymbol{\theta}) + \boldsymbol{i}b^{\ast}\_{+} \boldsymbol{B}^{\ast}\_{+}(-m\_{+},k\_{y}) \boldsymbol{\Phi}^{+}\_{m\_{+},-k\_{y}}(\boldsymbol{\theta}) + \\ + \mathrm{i}a^{\ast}\_{-} \boldsymbol{B}^{\ast}\_{-}(m\_{-},k\_{y}) \boldsymbol{\Phi}^{-}\_{-m\_{-},-k\_{y}}(\boldsymbol{\theta}) + \boldsymbol{i}b^{\ast}\_{-} \boldsymbol{B}^{\ast}\_{-}(-m\_{-},k\_{y}) \boldsymbol{\Phi}^{-}\_{m\_{-},-k\_{y}}(\boldsymbol{\theta}); \quad \mathbf{x} = \mathrm{R}\cos\theta\_{0\prime}, \boldsymbol{\theta} = -\boldsymbol{\phi}/2; \end{split} \tag{26}$$

• at the boundary between Region II and III

$$\begin{array}{ll} ia\_{+}^{\*}B\_{+}^{\*}(m\_{+},k\_{y})\Phi\_{-m\_{+},-k\_{y}}^{+}(\theta) + ib\_{+}^{\*}B\_{+}^{\*}(-m\_{+},k\_{y})\Phi\_{m\_{+},-k\_{y}}^{+}(\theta) + \\ + ia\_{-}^{\*}B\_{-}^{\*}(m\_{-},k\_{y})\Phi\_{-m\_{-},-k\_{y}}^{-}(\theta) + ib\_{-}^{\*}B\_{-}^{\*}(-m\_{-},k\_{y})\Phi\_{m\_{-},-k\_{y}}^{-}(\theta) = \\ t\_{L}(\boldsymbol{\phi})\_{\uparrow}^{\uparrow\star}\Psi\_{-k}^{-}(\boldsymbol{x}) - t\_{L}(\boldsymbol{\phi})\_{\downarrow}^{\uparrow\star}\Psi\_{-k}^{+}(\boldsymbol{x}); \quad \boldsymbol{x} = -\mathcal{R}\cos\theta\_{0\prime}\ \theta = \boldsymbol{\phi}/2. \end{array} \tag{27}$$

Next, we consider the equations that determine the coefficients *rR*(*φ*) *β <sup>α</sup>* and *tR*(*φ*) *β <sup>α</sup>* 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

$$\begin{array}{c} \mathbf{Y}^{-}\_{-k\_x,-k\_y}(\mathbf{x}) + r\_R(\boldsymbol{\phi})^\perp\_\downarrow \mathbf{Y}^{-}\_{k\_x,-k\_y}(\mathbf{x}) + r\_R(\boldsymbol{\phi})^\perp\_\uparrow \mathbf{Y}^{+}\_{k\_x,-k\_y}(\mathbf{x}) = \\ \mathbb{1}\_+ \boldsymbol{\Phi}^+\_{m,-k\_y}(\boldsymbol{\theta}) + \tilde{b}\_+ \boldsymbol{\Phi}^+\_{-m,+-k\_y}(\boldsymbol{\theta}) + \mathbb{1}\_- \boldsymbol{\Phi}^-\_{m,-k\_y}(\boldsymbol{\theta}) + \tilde{b}\_- \boldsymbol{\Phi}^-\_{-m,-k\_y}(\boldsymbol{\theta}); \quad \mathbf{x} = \mathrm{R} \cos \theta\_0, \boldsymbol{\theta} = \boldsymbol{\phi}/2; \end{array} \tag{28}$$

• at the boundary between Regions I and II

$$\begin{array}{ll} \mathfrak{a}\_{+}\Phi^{+}\_{m\_{+}-k\_{y}}(\theta) + \tilde{b}\_{+}\Phi^{+}\_{-m\_{+}-k\_{y}}(\theta) + \mathfrak{a}\_{-}\Phi^{-}\_{m\_{-}-k\_{y}}(\theta) + \tilde{b}\_{-}\Phi^{-}\_{-m\_{-}-k\_{y}}(\theta) =\\ \mathfrak{t}\_{R}(\phi)^{\downarrow}\_{\downarrow}\Psi^{-}\_{-k\_{x}-k\_{y}}(\mathbf{x}) + \mathfrak{t}\_{R}(\phi)^{\downarrow}\_{\uparrow}\Psi^{+}\_{-k\_{x}-k\_{y}}(\mathbf{x}); \quad \mathbf{x} = -R\cos\theta\_{0}, \quad \theta = -\phi/2. \end{array} \tag{29}$$

From the comparison of Equations (26) and (27) with Equations (28) and (29), it follows that the coefficients *rR*(*φ*) *β <sup>α</sup>* and *tR*(*φ*) *β <sup>α</sup>* of Equations (28) and (29) can be expressed in the following form

$$r\_{\mathbb{R}}(\boldsymbol{\phi})\_{\downarrow}^{\downarrow} = -r\_{\mathbb{L}}(-\boldsymbol{\phi})\_{\uparrow}^{\uparrow\*} \, ; \quad r\_{\mathbb{R}}(\boldsymbol{\phi})\_{\uparrow}^{\downarrow} = r\_{\mathbb{L}}(-\boldsymbol{\phi})\_{\downarrow}^{\uparrow\*} \, ; \tag{30}$$

and

$$t\_R(\boldsymbol{\phi})\_\downarrow^\downarrow = -t\_L(-\boldsymbol{\phi})\_\uparrow^{\uparrow\*}; \quad t\_R(\boldsymbol{\phi})\_\uparrow^\downarrow = -t\_L(-\boldsymbol{\phi})\_\downarrow^{\uparrow\*}.\tag{31}$$

As a result, we obtain for the reflection probabilities

$$|r\_{\mathbb{R}}(\boldsymbol{\phi})\_{\downarrow}^{\downarrow}|^{2} = |r\_{L}(-\boldsymbol{\phi})\_{\uparrow}^{\uparrow}|^{2}; \quad |r\_{\mathbb{R}}(\boldsymbol{\phi})\_{\uparrow}^{\downarrow}|^{2} = |r\_{L}(-\boldsymbol{\phi})\_{\downarrow}^{\uparrow}|^{2} \,. \tag{32}$$

while for the transmission probabilities we have

$$|\left|t\_{\mathbb{R}}(\boldsymbol{\phi})\_{\downarrow}^{\downarrow}\right|^{2} = |t\_{L}(-\boldsymbol{\phi})\_{\uparrow}^{\uparrow}|^{2} \; ; \quad |t\_{\mathbb{R}}(\boldsymbol{\phi})\_{\uparrow}^{\downarrow}|^{2} = |t\_{L}(-\boldsymbol{\phi})\_{\downarrow}^{\uparrow}|^{2} \; . \tag{33}$$

We have similar probabilities for the electron motion with the spin up polarisation from the right to the left side

$$|r\_R(\boldsymbol{\phi})\_\uparrow^\uparrow|^2 = |r\_L(-\boldsymbol{\phi})\_\downarrow^\downarrow|^2; \quad |r\_R(\boldsymbol{\phi})\_\downarrow^\uparrow|^2 = |r\_L(-\boldsymbol{\phi})\_\uparrow^\downarrow|^2. \tag{34}$$

and, consequently, for the transmission probabilities

$$|t\_R(\boldsymbol{\phi})\_\uparrow^\uparrow|^2 = |t\_L(-\boldsymbol{\phi})\_\downarrow^\downarrow|^2; \quad |t\_R(\boldsymbol{\phi})\_\downarrow^\uparrow|^2 = |t\_L(-\boldsymbol{\phi})\_\uparrow^\downarrow|^2. \tag{35}$$

From these results, it follows that the operator *S*ˆ *<sup>t</sup>* 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*<sup>±</sup> → −*m*<sup>±</sup> and the sign of the *k<sup>y</sup>* component: *k<sup>y</sup>* → −*ky*. We conclude that the operator *S<sup>t</sup>* acts like a time-reversal operator in the single valley.

#### *3.2. The Operator S*ˆ *ch* = *τ<sup>x</sup>* ⊗ *σ<sup>y</sup>*

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<sup>y</sup>* = 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 = \pm (\mathcal{D} \pm \lambda\_{\mathbf{x}}) \, , \tag{36}$$

where D = q *λ*2 *<sup>y</sup>* + *t* 2 *<sup>m</sup>* and *t<sup>m</sup>* = *γ <sup>R</sup>m*. For the rippled (arc) piece (see Figure 1), there are four eigenenergies

$$E = \begin{cases} \mathcal{E}\_1 = \lambda\_\mathcal{X} + \mathcal{D} \\ \mathcal{E}\_2 = \lambda\_\mathcal{X} - \mathcal{D} \\ \mathcal{E}\_3 = -\lambda\_\mathcal{X} + \mathcal{D} \\ \mathcal{E}\_4 = -\lambda\_\mathcal{X} - \mathcal{D} \end{cases} \tag{37}$$

The connection between the energy and the quantum number *m* can be formulated in the form

$$m \Rightarrow m\_{\rm s} = \pm \frac{R}{\gamma} \sqrt{(sE - \lambda\_{\rm x})^2 - \lambda\_{\rm y}^2} \quad \text{s} = \pm 1 \,, \tag{38}$$

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 : \Longleftrightarrow \begin{cases} s = +1, & \mathcal{E}\_1 > 0 \\ s = +1, & \text{if } \mathcal{E}\_2 > 0 \\ s = -1, & \text{if } \mathcal{E}\_3 > 0 \end{cases} \tag{39}$$

On the other hand, at *E* < 0, the branches E<sup>3</sup> and E<sup>4</sup> have the quantum number *s* = +1, while the branch E<sup>2</sup> 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*

$$
\mathcal{E}\_1 \, \mathcal{E}\_2 \to \pm t\_{m\_{\pm 1}} \Longleftrightarrow s = +1 \, \text{s} \tag{40}
$$

$$
\mathcal{E}\_3 \to \pm t\_{m\_{-1}} \Longleftrightarrow s = -1 \,\, . \tag{41}
$$

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 Figures 2 and 3). The anticrossings yield the energy gaps = 2*λ<sup>y</sup>* indicated by the arrows (see the insets in Figures 2 and 3). As a result, the energy gaps give rise to evanescent modes at energies *λ<sup>x</sup>* − *λ<sup>y</sup>* < |*E*| < *λ<sup>x</sup>* + *λ<sup>y</sup>* in our system.

**Figure 2.** The spectrum in Equation (4) (*k<sup>y</sup>* = 0) versus the magnetic quantum number *m*. The non-quantised values ±*ms*=±<sup>1</sup> at the energy *<sup>E</sup>* = 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*<sup>1</sup> , *E*2, *E*3, *E*<sup>4</sup> 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<sup>1</sup> , E2) at *E* > 0. Similar anticrossing occurs at *E* < 0, when the pair (E3, E<sup>4</sup> ) has the same quantum number *s* = +1. These anticrossings are caused by the term *λ<sup>y</sup>* in the Hamiltonian in Equation (1), which creates the energy gaps 2*λ<sup>y</sup>* near the energy *E* = ±*λ<sup>x</sup>* (see [13,17]). The following parameters are used: *R* = 10 Å, *δ* = 0.01, *p* = 0.1, *γ* = (4.5 · 1.42) eV· Å, *γ* <sup>0</sup> = <sup>8</sup> 3 *γ*. These parameters define the values of the spin–orbit coupling strengths: *λ<sup>y</sup>* = *δγ*0/4*R* = 0.0043 eV, *λ<sup>x</sup>* = *γ*(1/2 + 2*δp*)/*R* = 0.32 eV (see Section 2).

**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<sup>1</sup> , E2) at *E* > 0 and for the pair (E3, E<sup>4</sup> ) at *E* < 0.

At *k<sup>y</sup>* = 0, the eigenvectors in Equation (5) of the Hamiltonian in Equation (1) transform to the following form (see details in [13,17]):

$$\Phi\_{j,m}(\theta) = \left(I \otimes e^{-i\frac{\theta y}{2}\theta}\right)\vec{v}\_{\dot{j}}(m) \ \textit{e}^{im\theta} \ , \dot{j} = 1 \ . \tag{42}$$

Here

$$
\vec{\upsilon}\_1 = \begin{pmatrix} -it\_m \\ -A\_- \\ -iA\_- \\ -t\_m \end{pmatrix}, \vec{\upsilon}\_2 = \begin{pmatrix} -it\_m \\ A\_+ \\ iA\_+ \\ -t\_m \end{pmatrix}, \vec{\upsilon}\_3 = \begin{pmatrix} -it\_m \\ A\_+ \\ -iA\_+ \\ t\_m \end{pmatrix}, \vec{\upsilon}\_4 = \begin{pmatrix} -it\_m \\ -A\_- \\ iA\_- \\ t\_m \end{pmatrix}, \tag{43}
$$

where

$$A\_{\pm} = \pm \lambda\_y + \sqrt{t\_m^2 + \lambda\_y^2}. \tag{44}$$

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* > *λ<sup>x</sup>* ⇒ *E* = E<sup>1</sup> ∨ *E* = E3, the base eigenfunctions are: ∼ exp(±i*m*+1*θ*)~*ν*1(±*tm*+<sup>1</sup> ); ∼ exp(±i*m*−1*θ*)~*ν*3(±*tm*−<sup>1</sup> ), 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<sup>y</sup>* = 0 is convenient to present in the form:

$$\Psi\_k^{\sigma}(\mathbf{x}) = \frac{1}{2} \begin{pmatrix} \pi \frac{k}{|k|} \\ \mathbf{1} \end{pmatrix} \otimes \begin{pmatrix} \mathbf{1} \\ \sigma i \end{pmatrix} e^{i\mathbf{k}\mathbf{x}} \Big/ \qquad \sigma = \pm \,. \tag{45}$$

As above, we use the following definitions: Ψ*<sup>σ</sup> k* (*x*) ≡ h*x*|*σ*, *k*i, *σ* = (+/−) ⇔ *σ* = (↑ / ↓). 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<sup>x</sup>* corresponds to the direction of the electron motion from the left to the right side of the considered system.

In the case *k<sup>y</sup>* = 0, there is one more operator *S*ˆ *ch* = *τ<sup>x</sup>* ⊗ *σ<sup>y</sup>* that commutes with both Hamiltonians in Equations (1) and (11), describing the flat and the rippled pieces of our system, respectively:

$$\left[\left[\hat{S}\_{\text{ch}\prime}H\_f\right]=0\,,\quad\left[\hat{S}\_{\text{ch}\prime}H\_r\right]=0\,,\tag{46}$$

For the wave function, associated with the flat graphene sheet, we have

$$
\begin{aligned}
\langle \mathcal{S}\_{\rm cl} | \uparrow \pm k \rangle &=& \pm | \uparrow \pm k \rangle \\
\langle \mathcal{S}\_{\rm cl} | \downarrow \pm k \rangle &=& \mp | \downarrow \pm k \rangle
\end{aligned}
\quad \Rightarrow \mathfrak{a} = \pm 1.
\tag{47}
$$

Thus, the operator *S*ˆ *ch* 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*ˆ *ch* 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*ˆ *ch* on the electron wave functions, associated with the rippled piece of our system, is as follows:

$$
\hat{S}\_{ch}\Phi\_{1,m}(\theta) = -\Phi\_{1,m}(\theta) \, , \tag{48}
$$

*S*ˆ *ch*Φ2,*m*(*θ*) = −Φ2,*m*(*θ*), (49)

$$
\hat{S}\_{\text{ch}} \Phi\_{\text{3,m}}(\theta) = + \Phi\_{\text{3,m}}(\theta) \, , \tag{50}
$$

$$
\hat{S}\_{\text{ch}} \Phi\_{4,m}(\theta) = +\Phi\_{4,m}(\theta) \,. \tag{51}
$$

In this case, the operator *S*ˆ *ch* 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.

	- 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*i]. 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*i 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*i]. 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*i 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<sup>y</sup>* = 0.

	- (a) The electron is moving from the left side (flat graphene sheet) with the spin down polarisation [Ψ − *k* (*x*) ≡ | ↓ +*k*i]. 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*i 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*i]. 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*i 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<sup>y</sup>* = 0.

	- 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*i]. In the rippled graphene region, there is only the transmission channel, defined by the wave function Φ2,*m*(*θ*). The wave function Ψ + −*k* (*x*) ≡ | ↑ −*k*i 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*i]. In the rippled graphene region, there is only the transmission channel, defined by the wave function Φ2,*m*(*θ*). The wave function Ψ − *k* (*x*) ≡ | ↓ +*k*i 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<sup>y</sup>* = 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<sup>y</sup>* = 0, the symmetry, associated with the operator *S*ˆ *ch*, 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 *<sup>k</sup><sup>y</sup>* <sup>6</sup><sup>=</sup> 0, the operator *<sup>S</sup>*<sup>ˆ</sup> *ch* 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*ˆ *<sup>t</sup> and S*ˆ *ch*

Note that these two symmetry operator commute if *k<sup>y</sup>* = 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:

$$
\hat{S}\_{l}|\psi\_{\pm}\rangle = \pm|\psi\_{\pm}\rangle \,, \quad \hat{S}\_{l}|\phi\_{\pm}\rangle = \pm|\phi\_{\pm}\rangle \,, \tag{52}
$$

$$
\mathcal{S}\_{\rm cl}|\psi\_{\pm}\rangle = +|\psi\_{\pm}\rangle \,, \quad \mathcal{S}\_{\rm cl}|\phi\_{\pm}\rangle = -|\phi\_{\pm}\rangle \,. \tag{53}
$$

For the rippled graphene piece, we have to consider only the case *k<sup>y</sup>* = 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

$$\Phi\_{\dot{j}} = \frac{1}{2} [\Phi\_{\dot{j},m}(\theta) + \Phi\_{\dot{j},-m}(\theta)], \quad \dot{j} = 1, \ldots, 4. \tag{54}$$

As a result, we obtain

$$
\pounds\_l \Phi\_{1,2} = \Phi\_{1,2} \prime \quad \pounds\_l \Phi\_{3,4} = -\Phi\_{3,4} \,, \tag{55}
$$

$$
\hat{S}\_{\text{cl}}\Phi\_{1,2} = -\Phi\_{1,2} \quad \hat{S}\_{\text{cl}}\Phi\_{3,4} = \Phi\_{3,4} \,. \tag{56}
$$

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*ˆ *<sup>t</sup>* (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*ˆ *ch* (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<sup>y</sup>* = 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**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
