\section{Introduction}
One of the ways to search for New Physics is
precision measurement of the electric dipole moment $d_l$ of a charged lepton ($l=e,\mu,\tau$). The  value of  $d_l$ predicted by the Standard Model (SM) is too small for experimental measurement. Therefore, the observation of electric dipole moment or its manifestation would directly demonstrate the existence of New Physics.

The manifestation of the lepton electric dipole moment can be sought in the process of $l\bar l$ pair production in $e^+e^-$ annihilation. The general form of  $\gamma l \bar l$ vertex can be represented as
\begin{align}\label{genvert}
	&\Gamma^\mu=-ie\Bigg\{ F_1^l(k^2)\gamma^\mu +\dfrac{\sigma^{\mu\nu}k_\nu}{2m_l} \Big [iF_2^l(k^2) +F_3^l(k^2)\gamma_5 \Big ] +\Big( \gamma^\mu - \dfrac{2k^\mu \, m_l}{k^2}\Big)\gamma_5 F_4^l(k^2) \Bigg\}\,,
\end{align}
where $m_l$ is the lepton mass, $e<0$ is the electron charge, $k$ is the 4-momentum of the photon, $F_1^l(k^2)$ is the Dirac form factor, $F_2^l(k ^2)$ is the Pauli form factor, $F_3^l(k^2)$ is the electric dipole form factor, $F_4^l(k^2)$ is the anapole form factor, $\sigma^{\mu\nu} =i/2[\gamma^\mu,\gamma^\nu]$, $\gamma_5=-i\gamma^0\gamma^1\gamma^2\gamma^3$. In the limit $k^2 \to 0$ these form factors are
\begin{align*}
	& F_1^l(0)=1\, ,\quad F_2^l(0)=\mu'_l\,\dfrac{2m_l}{e}\, ,\quad F_3^l(0)=d_l\,\dfrac{2m_l}{e}\,,\quad
	F_4^l(0)=0\,,
\end{align*}
where $\mu'_l$ is the anomalous magnetic moment. It follows from \eqref{genvert} that a violation of $\mbox{P}$ and $\mbox{T}$ parities leads to the appearance of $F_3^l(k^2)$, while $F_4^l(k^2)$ is related to  violation of $\mbox{P}$-parity alone. Assuming that the $\mbox{CPT}$ theorem holds, violation of $\mbox{T}$-parity is equivalent to violation of $\mbox{CP}$-parity. Thus, $d_l$ occurs due to $\mbox{CP}$ violations.

For all charged leptons, predictions of  $\mu'_l$  in  SM \cite{Aoyama:2019ryr,Volkov:2019phy,Aoyama:2020ynm,Eidelman:2007sb}  can be experimentally verified \cite{Fan:2022eto,Muong-2:2021ojo,DELPHI:2003nah}. For $d_l$ the situation is essentially different. An estimate of $d_l$ in  SM  gives $d_l=(m_l/\mbox{MeV})\times 1.6\times 10^{-40}\,e\,\mbox{cm} $ \cite{Booth:1993af,Mahanta:1996er}. Hence $F_3^e(0) < F_3^\mu(0) < F_3^\tau(0)\approx 5 \times10^{-23}\ll 1$. The sensitivity of modern experiments does not allow one to  measure  $F_3^\tau(0)$ with an accuracy of $10^{-23}$. Therefore, extracting a non-zero value of $F_3^\tau(0)$ from the experiment would be a discovery of New Physics.  In Refs.~ \cite{Barr:1988mc,Grifols:1990ha,Escribano:1993pq,Escribano:1996wp,Taylor:1997tn,OPAL:1998dsa,L3:1998gov,DELPHI:2003nah,L3:2004ful,Blinov:2008mu,Grozin:2008nw} upper limits were set to $|F_3^\tau(k^2)|$, and in Refs.~\cite{ARGUS:2000riz,Belle:2002nla,Belle:2021ybo} upper limits were set to $\mbox{Re} \, F_3^\tau(k^2)$ and $\mbox{Im}\, F_3^\tau(k^2)$ separately.

In our work,  the processes $e^+e^-\rightarrow \tau^+\pi^-\nu_\tau$,
 $e^+e^-\rightarrow \pi^+\tau^-\bar\nu_\tau$, $e^+e^-\rightarrow \tau^+\rho^-\nu_\tau$,  and $e^+e^-\rightarrow \rho^+\tau^-\bar\nu_\tau$
are studied for longitudinally polarized electron beams, and  the   process  $e^+e^-\rightarrow \pi^+\pi^-\nu_\tau\bar\nu_\tau$ is discussed for unpolarized electrons. The asymmetric with respect to $\mbox{CP} $ transformation parts of the corresponding cross sections  are obtained for $e^+e^-$ invariant masses  $\sqrt{s}\ll m_Z$.  To derive these results, it is sufficiently to consider   $\gamma\tau^+\tau^-$ vertex  in the form
\begin{align}\label{tauvert}
	&\Gamma^\mu=-ie\Bigg [ \gamma^\mu +\dfrac{\sigma^{\mu\nu}k_\nu}{2M} F_3^\tau(k^2)\gamma_5\Bigg ]\,,
\end{align}
where $M$ is the $\tau$ lepton mass and $k^2=s$.

Measurement of $\mbox{CP}$ odd parts of the cross sections  can diminish the upper limits of both $\mbox{Re}\,F_3^\tau(s)$ and $\mbox{Im}\,F_3^\tau(s)$. Moreover,  polarization of beams simplifies the measurement of $\mbox{Re}\,F_3^\tau(s)$. At present, there are no corresponding experiments with    longitudinally polarized electrons. However, such  experiments are planned to be carried out at the Super-Charm-Tau factory (SCTF) \cite{TauFactory}. This  collider, having high luminosity $\sim 10^{35}\,\mbox{cm}^{-2}\,\mbox{s}^{-1}$ and $\sqrt{s}$ from
$3$ to $5-7\,\mbox{GeV}$,  will become an intense source of  $\tau$ leptons.

\section{$e^+e^-\to \tau^+\tau^-$}
To study the effect of polarization, let us consider a longitudinally polarized electron beam and an unpolarized positron beam. Since $\sqrt{s}\ll m_Z$, we neglect the contribution of the Z boson. Then the cross section $d\sigma_0$ of the process $e^+e^-\to \tau^+\tau^-$ in the center-of-mass frame is
\begin{align}
	&d\sigma_0=\dfrac{\beta\alpha^2}{4s}\,\left |\phi_\tau^\dagger\,H\,\chi_\tau\right|^2\,d\Omega_{\bm q}\,,\nonumber\\
	&H=\bm\sigma\cdot\bm e_\lambda-\dfrac{(\bm\sigma\cdot\bm q)(\bm e_\lambda\cdot\bm q)}{E(E+M)}+i\dfrac{b}{M}\,(\bm e_\lambda\cdot\bm q)\,,\nonumber\\
	&b=F_3^\tau(s)\,,\quad s=4E^2\,,\quad \beta=q/E\,,\quad \bm e_\lambda=\frac{1}{\sqrt{2}}(\bm e_x+i\lambda\bm e_y)\,.
\end{align}
Here $\alpha$ is the fine structure constant charge, $E$ is the electron energy, $\bm q$ is the momentum of the $\tau^-$ lepton, the vector $\bm e_z$ is directed along the electron momentum, $\lambda$ is the electron helicity, $\phi_\tau$ and $\chi_\tau$ are two-component spinors entering, respectively, into the positive-frequency and negative-frequency solutions of the Dirac equation for the $\tau$ lepton,
\begin{align}\label{wfin}
	&U_{\bf q}=\sqrt{\dfrac{E_q+M}{2E_q}}
	\begin{pmatrix}
	\phi_\tau\\
	\dfrac{\bm\sigma\cdot\bm q}{E_q+M}\,\phi_\tau\end{pmatrix}\, ,\quad
V_{\bf -q}=\sqrt{\dfrac{E_q+M}{2E_q}}
\begin{pmatrix}
\dfrac{-\bm\sigma\cdot\bm q}{E_q+M}\,\chi_\tau\\
\chi_\tau\end{pmatrix}\,,
	\end{align}
where $E_q=\sqrt{\bm q^2+M^2}$.
The cross section $d\sigma_0$ summed over polarizations of $\tau^+$ is
\begin{align}
	&d\sigma_0=\dfrac{\beta\alpha^2}{4s}\,\left[ 1-\dfrac{ q_\perp^2}{2E^2}+\bm\zeta\cdot\bm Z\right]\,d\Omega_{\bm q}\,,\nonumber\\
	&\bm Z=\mbox{Im}\,b\,\dfrac{q_\perp^2\,\bm q}{ME(E+M)}-\mbox{Im}\,b\,\dfrac{\bm q_\perp}{M}-\,\mbox{Re}\,b\,\dfrac{[\bm q_\perp\times\bm \Lambda]}{M}\nonumber\\
	&+\dfrac{M}{E}\bm \Lambda+\dfrac{(\bm q\cdot\bm \Lambda)\,\bm q}{E(E+M)}\,,
\end{align}
where $\bm q_\perp=\bm q- \bm \Lambda (\bm q \cdot \bm \Lambda)$, $\bm \zeta$ is the spin of $\tau^-$, $\bm\Lambda =\lambda\bm e_z$,  terms quadratic in $b$ are omitted. It is seen that the  linear in $b$ terms contribute only to the $\bm \zeta$-dependent part of the cross section. Measurement of $\tau$ decay products is a way to measure polarization, which in turn makes it possible to measure $b$. The total cross section summed over the $\tau^-$ polarization is
 \begin{align}\label{sig0}
 	&\sigma_0=\dfrac{\pi\beta\alpha^2}{3E^2}\,\left( 1+\dfrac{M^2}{2E^2}\right)\,.
 \end{align}

\section{$e^+e^-\to \tau^+\pi^-\nu_\tau$ , $e^+e^-\to \tau^-\pi^+{\bar\nu}_\tau$}
Consider the cross section of the process $e^+e^-\to \tau^+\tau^-$ followed by the decay  $\tau^-\to \pi^-\nu_\tau$. Taking into account the smallness of the $\tau$  lepton  width,  $\Gamma_\tau\approx 2.27\,\mbox{meV}$ \cite{Workman:2022ynf},  we can make the substitutions
\begin{align}
&\dfrac{1}{{\hat q}-M}\rightarrow  \,\dfrac{2E_q}{{\cal E}^2-E_q^2+i\Gamma_\tau M}\,\sum_{\mu}U_{\bm q,\mu}\bar U_{\bm q,\mu} \,,\nonumber\\
&\dfrac{4E_q^2}{({\cal E}^2-E_q^2)^2+\Gamma_\tau^2 M^2}\rightarrow \dfrac{2\pi E_q}{M\Gamma_\tau}\,\delta({\cal E}-E_q)\,,
\end{align}
where ${\cal E}=q^0$. After that, the cross section of the process $e^+e^-\to \tau^+\pi^-\nu_\tau$ can be represented as
\begin{align}
&d\sigma_\pi^{(-)}(\bm k)=B_\pi\dfrac{\beta\alpha^2(E+M)\,d\Omega_{\bm q}\,d\bm k}{4\pi sM^2\,\omega_k}\,R^{(-)}\,\delta(E-\omega_k-|\bm q-\bm k|)\,,\nonumber\\
&R^{(-)}=\left |\phi^+_\nu\left[1+\dfrac{\bm\sigma\cdot\bm q}{E+M}\right]\,H\,\chi_\tau\right|^2\,,
\end{align}
where $B_\pi\approx 10.8\%$ \cite{Workman:2022ynf} is the branching ratio of $\tau\to \pi\nu$  decay, $\omega_k=\sqrt{\bm k^2+m_\pi^2 }$ and $\bm k$ are the pion energy and momentum, respectively, $\phi_\nu$ is a two-component spinor in the Dirac  spinor $U_{\bm Q}$ for neutrino, $\bm Q=\bm q-\bm k$ and $(\bm \sigma\cdot\bm Q)\phi_\nu=-Q\phi_\nu$.
We also neglected the pion mass $m_\pi$ compared to $M$ and took into account that the  matrix element of the decay $\tau\to \pi\nu$ is proportional to $\bar U_{\bm Q}\,U_{\bm q }$ \cite{Okun:1982ap}. After summing over polarizations $\chi_\tau$, we find
\begin{align}
&R^{(-)}=A_0-\bm n_{\bm Q}\bm A\,,\quad \bm n_{\bm Q}=\bm Q/Q\,,\nonumber\\
&A_0=\dfrac{2E}{E+M}\left[ 1-\dfrac{ q_\perp^2}{2E^2}+\dfrac{(\bm q\cdot\bm \Lambda)}{E}- \mbox{Im}\,b\,\dfrac{q_\perp^2}{E^2}\right]\,,\nonumber\\
&\bm A=\dfrac{2}{(E+M)}\Bigg\{\dfrac{M^2}{E}\bm \Lambda+\left[1+\dfrac{(\bm q\cdot\bm \Lambda)}{E}-\dfrac{ q_\perp^2}{2E^2}\right]\,\bm q\nonumber\\
&-\mbox{Im}\,b\,\bm q_\perp-\mbox{Re}\,b\,[\bm q\times\bm \Lambda]\Bigg\}\,.
\end{align}
Integrating over the angles of the vector $\bm q$, we obtain
\begin{align}\label{sm}
	&d\sigma_\pi^{(-)}(\bm k)=B_\pi\dfrac{\alpha^2}{4E^3}\,\Big[\xi_0+\xi_1\cos\theta+\xi_2(3\cos^2\theta-1)\Big]\,d\omega_k\,d\Omega_{\bm k}\,,\nonumber\\
	&\xi_0=\dfrac{1}{3}\left(1+\dfrac{M^2}{2E^2}\right)+\dfrac{\mbox{Im}\,b}{3}\left(1-\dfrac{2\omega_k}{E}\right)\,,\nonumber\\
	&\xi_1=\dfrac{M^2}{4E\omega_k}-\dfrac{1}{2}\left(1-\dfrac{2\omega_k}{E}\right)\,,\nonumber\\
	&\xi_2=\dfrac{M^2(M^2-4E\omega_k)}{32E^2\omega_k^2}+\dfrac{1}{12}\left(1+\dfrac{M^2}{2E^2}\right)+\dfrac{\mbox{Im}\,b}{3}\left(1-\dfrac{\omega_k}{2E}-\dfrac{3M^2}{8E\omega_k}\right)\,.
\end{align}
Here $\cos\theta=\bm k\cdot\bm\Lambda/k=\lambda k_z/k$, the available pion energy range is determined by the relation $|2\omega_k-E|\leq q$. It is seen that the only term, which contains the helicity $\lambda$, is the term proportional to $\xi_1$. However, this term is independent of $b$. Moreover, the cross section $d\sigma_\pi^{(-)}(\bm k)$ in \eqref{sm} contains only the imaginary part of $b$, which does not require electron polarization to be measured.
For the cross section $d\sigma_\pi^{(+)}(\bm k)$ of the process $e^+e^-\to \tau^-\pi^+{\bar\nu}_\tau$, we similarly find
\begin{align}\label{sp}
	&d\sigma_\pi^{(+)}(\bm k)=B_\pi\dfrac{\alpha^2}{4E^3}\,\Big[\xi'_0+\xi'_1\cos\theta+\xi'_2(3\cos^2\theta-1)\Big]\,d\omega_k\,d\Omega_{\bm k}\,,\nonumber\\
	&\xi'_0=\dfrac{1}{3}\left(1+\dfrac{M^2}{2E^2}\right)-\dfrac{\mbox{Im}\,b}{3}\left(1-\dfrac{2\omega_k}{E}\right)\,,\nonumber\\
	&\xi'_1=-\xi_1=-\dfrac{M^2}{4E\omega_k}+\dfrac{1}{2}\left(1-\dfrac{2\omega_k}{E}\right)\,,\nonumber\\
	&\xi'_2=\dfrac{M^2(M^2-4E\omega_k)}{32E^2\omega_k^2}+\dfrac{1}{12}\left(1+\dfrac{M^2}{2E^2}\right)-\dfrac{\mbox{Im}\,b}{3}\left(1-\dfrac{\omega_k}{2E}-\dfrac{3M^2}{8E\omega_k}\right)\,.
\end{align}
We define the differential asymmetry $dA_\pi$ as follows
\begin{align}\label{api}
	&dA_\pi=\dfrac{d\sigma_\pi^{(-)}(\bm k)-d\sigma_\pi^{(+)}(-\bm k)}{2\sigma_0}\,,
\end{align}
where $\sigma_0$ is defined in \eqref{sig0}.
The asymmetry is proportional to $\mbox{Im}\,b$, so its measurement is important from the experimental point of view. We have
\begin{align}\label{da}
	&dA_\pi=\dfrac{B_\pi{\mbox{Im}\,b}\,d\omega_k\,d\Omega_{\bm k}}{4\pi q(1+M^2/2E^2)}\,\left[1-\dfrac{2\omega_k}{E}+(3\cos^2\theta-1)\left(1-\dfrac{\omega_k}{2E}-\dfrac{3M^2}{8E\omega_k}\right)\right]\,.
\end{align}
After integration over angles, the asymmetry is
\begin{align}\label{da1}
	&dA_\pi=\dfrac{{B_\pi\mbox{Im}\,b}\,d\omega_k}{ q(1+M^2/2E^2)}\,\left(1-\dfrac{2\omega_k}{E}\right)\,.
\end{align}
As it should be, after integration over the pion energy, the asymmetry vanishes. Therefore, we define the total asymmetry $A_\pi$  as $dA_\pi$ \eqref{da1} integrated over $\omega_k$ from $(E-q)/2$ to $E/2$ (half of the allowed energy range),
\begin{align}\label{apiE}
	&A_\pi=\dfrac{B_\pi\,\mbox{Im}\,b}{ 2(2+M^2/E^2)}\,\sqrt{1-\dfrac{M^2}{E^2}}\,.
\end{align}

\section{$e^+e^-\to \tau^+\rho^-\nu_\tau$ , $e^+e^-\to \tau^-\rho^+{\bar\nu}_\tau$}
To measure $\mbox{Re}\,b$, consider the decay of one $\tau$ lepton into $\rho$ meson with the momentum
$\bm p$, energy $\varepsilon_p=\sqrt{\bm p^2+m_\rho^2}$ and 4-polarization vector $f=(f_0,\,\bm f)$, where $m_\rho $ is the mass of $\rho$ meson. We define the differential asymmetry $dA_\rho$ as
\begin{align}\label{arho}
	&dA_\rho=\dfrac{d\sigma_\rho^{(-)}(\bm p,\,\bm f)-d\sigma_\rho^{(+)}(-\bm p,\,-\bm f)}{2\sigma_0}\,,
\end{align}
where $d\sigma_\rho^{(-)}(\bm p,\,\bm f)$ is the cross section of the process $e^+e^-\to \tau^+\rho^-\nu_\tau $ and $d\sigma_\rho^{(+)}(-\bm p,\,-\bm f)$ is the cross section of the process $e^+e^-\to \tau^-\rho^+{ \bar\nu}_\tau$. Using the  matrix element of decay $\tau\to \rho\nu$ \cite{Okun:1982ap},  we obtain as a result of straightforward calculations
\begin{align}\label{darhopol}
	&dA_\rho=\dfrac{3B_\rho\,d\varepsilon_p\,d\Omega_{\bm p}\,\,[C_1\,\mbox{Re}\,b+C_2\,\mbox{Im}\,b]}{2\pi pM^2(1+M^2/2E^2)(2+M^2/m_\rho^2)(1-m^2_\rho/M^2)^2}\,,\nonumber\\
	&C_1=\left[\dfrac{q}{p}\varepsilon_p\,P_2(x_0)-E P_1(x_0)\right]\,([\bm\Lambda\times\bm f]\cdot\bm p)\,f_0\,,\nonumber\\ &C_2=\dfrac{qp}{3}\left[\left(2+\dfrac{\varepsilon_p}{E}\right)f_0^2+\left(1-\dfrac{\varepsilon_p}{E}\right)\bm f^2-(\bm\Lambda\cdot\bm f)^2\right]\nonumber\\
	&+P_1(x_0)\,\Big[\dfrac{1}{2}[p^2-(\bm\Lambda \cdot\bm p)^2](\bm f^2-f_0^2) \nonumber\\
	&+Ef_0[(\bm\Lambda\cdot\bm p)(\bm\Lambda\cdot\bm f)-(\bm p\cdot\bm f)]
	+\dfrac{2q^2}{5E}f_0[(\bm\Lambda\cdot\bm p)(\bm\Lambda\cdot\bm f)-2(\bm p\cdot\bm f)]\Big]\nonumber\\
	&+\dfrac{q}{p}\,P_2(x_0)\,\Bigg\{\left(\dfrac{\varepsilon_p}{2E}(\bm f^2-f_0^2)-f_0^2\right)\left[(\bm\Lambda\cdot\bm p)^2-\dfrac{p^2}{3}\right]+\left[(\bm f\cdot\bm p)^2-\dfrac{p^2\bm f^2}{3}\right]\nonumber\\
	&-(\bm f\cdot\bm\Lambda)\left[(\bm\Lambda\cdot\bm p)(\bm f\cdot\bm p)-\dfrac{p^2}{3}(\bm f\cdot\bm\Lambda)\right]\Bigg\}\nonumber\\
	&+P_3(x_0)\,\dfrac{q^2}{Ep^2}f_0\,\Big\{(\bm f\cdot\bm p)(\bm\Lambda\cdot\bm p)^2-\dfrac{p^2}{5}[(\bm f\cdot\bm p)+2(\bm f\cdot\bm\Lambda)(\bm p\cdot\bm\Lambda)]\Big\}\,,\nonumber\\
	&x_0=\dfrac{2E\varepsilon_p-M^2-m_\rho^2}{2qp}\,,
\end{align}
where $B_\rho\approx25.5\%$ \cite{Workman:2022ynf} is the branching ratio of $\tau\to\rho\nu$ decay and $P_n(x)$ are Legendre polynomials. It is seen that for a polarized electron beam and  a polarized $\rho$ meson the contribution of  $\mbox{Re}\,b$ does not vanish. Note that
\begin{align}
	& dA_\rho\Big|_{\lambda=+1}-dA_\rho\Big|_{\lambda=-1}=\dfrac{3B_\rho\,\mbox{Re}\,b\,C_1\,d\varepsilon_p\,d\Omega_{\bm p}}{\pi pM^2(1+M^2/2E^2)(2+M^2/m_\rho^2)(1-m^2_\rho/M^2)^2}\,.
\end{align}
Then we perform summation over polarizations of $\rho$ meson, using the formula
$$\sum_{pol} f^\mu f^\nu= -g^{\mu\nu}+\dfrac{p^\mu p^\nu}{m_\rho^2}\,,$$
and obtain
\begin{align}\label{darho1}
		&\sum_{pol}C_1=0\,,\quad \sum_{pol}C_2=C_{21}\,F+C_{22}\,,\nonumber\\
		&C_{21}=-\dfrac{qp}{3}+P_1(x_0)\,\left(-\dfrac{3}{2}m_\rho^2+\dfrac{7}{5}E\varepsilon_p
		-\dfrac{2\varepsilon_p}{5E}M^2\right)\nonumber\\
	&+\dfrac{q}{p}P_2(x_0)\left(\dfrac{3\varepsilon_p}{2E}m_\rho^2-\dfrac{5}{3}\varepsilon_p^2
	+\dfrac{2}{3}m_\rho^2\right)+P_3(x_0)\,\dfrac{3\varepsilon_pq^2}{5E}\,,\nonumber\\
&C_{22}=\dfrac{p(M^2-2m_\rho^2)[m_\rho^2 E+M^2(E-2\varepsilon_p)]}{6m_\rho^2E\,q}\,,\nonumber\\
&F=\dfrac{1}{3m_\rho^2}\left[3(\bm\Lambda\cdot\bm p)^2-	p^2\right]\,.
\end{align}
In this formula, the contribution $\propto\mbox{Re}\,b$ is absent even for the case of polarized electrons. After integration over the angles of the vector $\bm p$, the term  $\propto F$ vanishes, and  the asymmetry reads
\begin{align}\label{darhoeps1}
	&dA_\rho=\dfrac{6B_\rho\,d\varepsilon_p\,C_{22}\,\mbox{Im}\,b}{ pM^2(1+M^2/2E^2)(2+M^2/m_\rho^2)(1-m^2_\rho/M^2)^2}\,.
	\end{align}
The allowed region of the energy $\varepsilon_p$ is given by the relation
$$\left|2\varepsilon_p-\left(1+\dfrac{m_\rho^2}{M^2}\right)E\right|\leq   q\left(1-\dfrac{m_\rho^2}{M^2}\right)\,.      $$
After integration over $\varepsilon_p$ in Eq.~\eqref{darhoeps1}, the asymmetry vanishes. Therefore, we define the total asymmetry  $A_\rho$ as a result of integration of Eq.~ \eqref{darhoeps1}  over $\varepsilon_p$ in the range
$$(1+m_\rho^2/M^2)E-(1-m_\rho^2/M^2)q<2\varepsilon_p< (1+m_\rho^2/M^2)E\,.$$
One has
\begin{align}\label{arhoE}
	&A_\rho=\dfrac{B_\rho\,\mbox{Im}\,b}{ 2(2+M^2/E^2)}\,\sqrt{1-\dfrac{M^2}{E^2}}\Bigg(\dfrac{M^2-2m_\rho^2}{M^2+2m_\rho^2}\Bigg)\,.
\end{align}
Note that
 \begin{align}\label{arhoE1}
 	&A_\rho=\dfrac{B_\rho}{ B_\pi}\Bigg(\dfrac{M^2-2m_\rho^2}{M^2+2m_\rho^2}\Bigg)\,A_\pi\,.
 \end{align}

To determine  the polarization of  $\rho$ meson, it is possible to measure the main decay channel of  $\rho$ meson with momentum $\bm p$ into two pions with momenta $\bm k_1$ and $\bm k_2$. The corresponding asymmetry can be obtained from Eq.~
\eqref{darhopol} with obvious substitution
\begin{align}\label{darhopipi}
	&\frac{d\bm p}{2\varepsilon_p(2\pi)^3}\,f^\mu\,f^\nu\rightarrow  \,\dfrac{f^2_{\rho\pi\pi}\,d\bm k_1\,d\bm k_2\,(k_1-k_2)^\mu(k_1-k_2)^\nu}{ 4\omega_{1}\omega_{2}(2\pi)^6 [(s_1-m_\rho^2)^2+\Gamma_\rho^2m_\rho^2]}\,.
\end{align}
Here $s_1=(k_1+k_2)^2$, $\Gamma_\rho=149.1\,\mbox{MeV}$ \cite{Workman:2022ynf} is the $\rho$ meson width, $\omega_1=|\bm k_1|$, $\omega_2=|\bm k_2|$, the constant $f^2_{\rho\pi\pi}$ is
$$ f^2_{\rho\pi\pi}=\dfrac{48\pi\Gamma_\rho}{m_\rho\,(1-4m_\pi^2/m_\rho^2)^{3/ 2}}\,.$$

\section{$e^+e^-\to \tau^+\tau^-\to \pi^+\pi^-\nu_\tau{\bar\nu}_\tau$}
If each $\tau$ lepton decays into a pion and a neutrino (antineutrino), both the imaginary and real parts of $b$ can be measured even in the case of an unpolarized electron beam. It is this case that we consider in this section.  We define asymmetry as
\begin{align}\label{apipi}
	&dA_{\pi\pi}=\dfrac{d\sigma_{\pi\pi}(\bm k_1,\,\bm k_2)-d\sigma_{\pi\pi}(-\bm k_2,\,-\bm k_1)}{2\sigma_0}\,,
\end{align}
where $d\sigma_{\pi\pi}(\bm k_1,\,\bm k_2)$ is the cross section of the process $e^+e^-\to \pi^+\pi^-\nu_\tau{\bar\nu}_\tau$, $\bm k_1$ is the momentum of $\pi^-$ and $\bm k_2$ is the momentum  of $\pi^+$. For $dA_{\pi\pi}$ we obtain
\begin{align}\label{dapipi}
	&dA_{\pi\pi}=\dfrac{3B^2_{\pi}\,d\bm k_1\,d\bm k_2}{ (2\pi)^3M^4(1+M^2/2E^2)\omega_{1}\omega_{2}(E-\omega_{1})(E-\omega_{2})}\nonumber\\
	&\times\,\int\,d\Omega_{\bm q}\,\delta(E-\omega_{1}-|\bm q-\bm k_1|)\,\,\delta(E-\omega_{2}-|\bm q+\bm k_2|)\nonumber\\
	&\times\Bigg<\mbox{Re}\,b\,\Bigg\{(\bm q\cdot\bm\Lambda)\,\Big([\bm q\times\bm\Lambda]\cdot(\omega_{2}\bm k_1+\omega_{1}\bm k_2)\Big)\nonumber\\
	&-\dfrac{\Big(M^2+(\bm q\cdot\bm\Lambda)^2\Big)}{E}\,\Big(\bm q\cdot[\bm k_1\times\bm k_2]\Big)
	+E(\bm q\cdot\bm\Lambda)\big(\bm \Lambda\cdot[\bm k_1\times\bm k_2]\Big)\Bigg\}\nonumber\\
	&+
	\mbox{Im}\,b\,\dfrac{M^2}{2}\Bigg\{\dfrac{(\omega_{2}-\omega_{1})}{E}[M^2+(\bm q\cdot\bm\Lambda)^2]+(\bm q\cdot\bm\Lambda)\big(\bm \Lambda,\bm k_1+\bm k_2\Big)\Bigg\}\Bigg>\,.
\end{align}
Integrating over the angles of  vector $\bm q$, we find
\begin{align}\label{dapipi1}
	&dA_{\pi\pi}=\dfrac{3B^2_{\pi}\,d\bm k_1\,d\bm k_2}{ (2\pi)^3M^4(1+M^2/2E^2)\omega^2_{1}\omega^2_{2}\,\sqrt{(1-x^2)\,{ q}^2(q^2-{ P}^2)}}\nonumber\\
	&\times\Bigg<\mbox{Re}\,b\,\Big(\bm\Lambda\cdot[\bm k_1\times\bm k_2]\Big)\Bigg[2(q^2-P^2)(\bm N_2\cdot\bm\Lambda)\nonumber\\
	&+\dfrac{(\bm P\cdot\bm\Lambda)}{E}\Big(M^2+2P^2-q^2-\dfrac{Ea_2}{1-x}\Big)\Bigg]+
	\mbox{Im}\,b\,\dfrac{M^2}{2}\Bigg[\dfrac{(\omega_{2}-\omega_{1})}{E}
	\Big(M^2\nonumber\\
	&+(\bm N_3\cdot\bm\Lambda)^2(q^2-{ P}^2)+(\bm P\cdot\bm\Lambda)^2\Big)+(\bm P\cdot\bm\Lambda)\big(\bm \Lambda,\bm k_1+\bm k_2\Big)\Bigg]\Bigg>\,,
\end{align}
 Here the following notation is introduced:
 \begin{align}\label{not}
 	&\bm N_1=\frac{\bm n_1+\bm n_2}{2(1+x)}\,,\quad  \bm N_2=\frac{\bm n_1-\bm n_2}{2(1-x)}\,,\quad \bm N_3=\frac{[\bm n_2\times\bm n_1]}{\sqrt{1-x^2}}\,,\nonumber\\
 	&\bm P=a_1\,\bm N_1+a_2\,\bm N_2\,,\quad a_1=\dfrac{M^2(\omega_{1}-\omega_{2})}{2\omega_{1}\omega_{2}}\,,\quad a_2=2E-\dfrac{M^2(\omega_{1}+\omega_{2})}{2\omega_{1}\omega_{2}}\,,\nonumber\\
 	&x=(\bm n_1\cdot\bm n_2)\,,\quad \bm n_1=\dfrac{\bm k_1}{\omega_{1}}\,,\quad \bm n_2=\dfrac{\bm k_2}{\omega_{2}}\,.
 \end{align}
 It is also interesting to consider the asymmetry that remains after the integration over the angles of  vectors $\bm k_1$ and $\bm k_2$.
The corresponding result reads
 \begin{align}\label{dapipiEE}
 	&dA_{\pi\pi}=\dfrac{2B^2_{\pi}\,\mbox{Im}\,b\,(\omega_2-\omega_1)d\omega_1\,d\omega_2}{ Eq^2(1+M^2/2E^2)}\,.
 \end{align}
Here only the imaginary part of $b$ contributes. Integrating over $\omega_2$ in the region $|2\omega_2-E|<q$, we obtain the result
\begin{align}\label{dapipiE}
	&dA_{\pi\pi}=\dfrac{B^2_{\pi}\,\mbox{Im}\,b\,(1-2\omega_1/E)d\omega_1}{ q(1+M^2/2E^2)}\,,
\end{align}
which is consistent with \eqref{da1}. Integrating \eqref{dapipiE} over $\omega_2$ in the range $(E-q)/2<\omega_2 <  E/2$, we get
\begin{align}\label{apipiE}
	&A_{\pi\pi}=\dfrac{B_\pi^2\,\mbox{Im}\,b}{ 2(2+M^2/E^2)}\,\sqrt{1-\dfrac{M^2}{E^2}}\,.
\end{align}
 Naturally, $A_{\pi\pi}={B_\pi}A_\pi$.

\section{Discussion of the results}
In our work we have obtained the asymmetries which contain both $\mbox{Re}b$ and $\mbox{Im}\,b$, as well as asymmetries proportional only to $\mbox{Re}\,b$ or only to $\mbox{Im}\,b$. To measure the asymmetry, which is proportional only to $\mbox{Re}\,b$, it is necessary to have longitudinal polarization of electrons.

\begin{figure}[h]
	\centering
	\includegraphics[width=0.7\linewidth]{"dfig"}
	\caption[AsimFig]{Total asymmetry $A$ in units of $\mbox{Im}\,b$. Solid curve: $A_\pi$; dotted curve: $A_\rho$; dashed curve: $A_{\pi\pi}$.}
	\label{fig:plot-for-article}
\end{figure}
Asymmetries \eqref{api},\eqref{arho}, and \eqref{apipi} are the odd quantities  with respect to $\mbox{CP}$ transformations. Indeed, as a result of this transformation
\begin{align*}
	&d\sigma_\pi^{(-)}(\bm k)\rightarrow \,d\sigma_\pi^{(+)}(-\bm k)\,,\quad d\sigma_\pi^{(+)}(-\bm k)\rightarrow\, d\sigma_\pi^{(-)}(\bm k)\,;\\
	&d\sigma_\rho^{(-)}(\bm p,\,\bm f)\rightarrow\, d\sigma_\rho^{(+)}(-\bm p,\,-\bm f)\,,\quad d\sigma_\rho^{(+)}(-\bm p,\,-\bm f)\rightarrow\, d\sigma_\rho^{(-)}(\bm p,\,\bm f)\,;\\
	&d\sigma_{\pi\pi}(\bm k_1,\,\bm k_2)\rightarrow\, d\sigma_{\pi\pi}(-\bm k_2,\,-\bm k_1)\,,\quad d\sigma_{\pi\pi}(-\bm k_2,\,-\bm k_1)\rightarrow\, d\sigma_{\pi\pi}(\bm k_1,\,\bm k_2)\,.
\end{align*}
The term $\propto\gamma^\mu$ in \eqref{tauvert} is $\mbox{CP}$-even while the term
 $\propto\,b$  is  $\mbox{CP}$-odd. Therefore, asymmetries appear due to interference between $\mbox{CP}$-odd and  $\mbox{CP}$-even terms and are linear in $b$.
The energy dependence of total asymmetries is shown in Fig. \ref{fig:plot-for-article}. The solid curve corresponds to $A_\pi$, the dotted curve to $A_\rho$, and the dashed curve to $A_{\pi\pi}$. All curves have the same energy dependence and differ only in scale. Indeed, it follows from Eqs.~ \eqref{apiE}, \eqref{arhoE}, \eqref{apipiE}  that
$$A_\rho=1.06\,A_{\pi}\,,\quad A_{\pi\pi}=0.108\,A_{\pi}\,.$$

\section{Conclusion}
In conclusion, we have considered the processes $e^+e^-\rightarrow \tau^+\pi^-\nu_\tau$, $e^+e^-\rightarrow \pi^+\tau^-\bar\nu_ \tau$, $e^+e^-\rightarrow \tau^+\rho^-\nu_\tau$, and $e^+e^-\rightarrow \rho^+\tau^-\bar\nu_\tau $ with longitudinally polarized electrons, as well as the process $e^+e^-\rightarrow \pi^+\pi^-\nu_\tau\bar\nu_\tau$ with unpolarized electrons for the invariant masses of the initial electron and positron $\sqrt {s}\ll m_Z$.  We have  calculated the $\mbox{CP}$-odd  asymmetries $\propto\mbox{Re}\,b$ and $\propto\mbox{Im}\,b$, as well as the asymmetries proportional only to $\mbox{ Re}\,b$ or only to $\mbox{Im}\,b$. Measuring these quantities can improve the upper limits for $\mbox{Re}\, b$ and $\mbox{Im}\, b$. It is shown that to measure $\mbox{Im}\,b$, polarization is not needed, and to measure $\mbox{Re}\,b$, the polarization is not necessary, but simplifies the measurement.

\section*{Acknowledgements}
We are grateful to D.A. Epifanov and V.S. Vorobiev for valuable discussions.


\section{\textbf{Introduction}}
The spectral analysis of operators defined over sequence spaces has been treated by many
researchers worldwide. The literature contains the spectrum and fine spectrum of several classes of Toeplitz operators \cite{t3,bir,duren}, Ces{\`a}ro operators \cite{cesaro3,cesaro2,cesaro1}, Rhaly operators \cite{rhaly2,rhaly1}, operators generated by various difference equations \cite{diff1,diff2,diff3,patra2019some,patra2020spectrum}, etc. For a detailed review, one may refer to the survey articles \cite{baliarsingh2021survey,yecsilkayagil2019survey} and the references therein. In particular, numerous mathematicians have focused their research on an important class of tridiagonal matrix, known as Jacobi matrix  \cite{bairamov2001non,janas2003similarity,janas1998point,}. The Jacobi matrix $J$ is generated by the difference equation
\[(Jy)_n=a_{n-1}y_{n-1}+b_ny_n+a_ny_{n+1},~~ y_n\in \mathbb{C^N}\]
  with certain initial conditions where $\{a_n\}$ and $\{b_n\}$ are complex sequences. In this paper, we focus on the spectral properties of a class of Jacobi-like penta-diagonal band matrices defined over the sequence space $\ell_p (1 < p < \infty)$ where the entries in the non-zero bands form sequences with two limit points.

Let $\ell_{p}$ represents the Banach space of  $p$-absolutely summable sequences of real or complex numbers with the norm
\[ \norm{x}_{p}=\left( \sum_{n=1}^{\infty}|x_n|^p\right)^{\frac{1}{p}}.\]
Also let $\mathcal{D}_p$ denotes the set of all diagonal operators on $\ell_p.$ For any operator $T \in \mathcal{D}_p,$ $\mbox{diag}(T)$ represents the sequence in the diagonal of $T.$ In this work we investigate the spectral properties of a class of operators $T$ defined over $\ell_p$ represented by the following form:
\[T = S_r^2D_1 + D_2S_{\ell}^2 + D_3\]
where $S_r, S_{\ell} : \ell_p \to \ell_p$, denotes the right shift operator, left shift operator respectively and $D_1, D_2, D_3 \in \mathcal{D}_p$ with $\mbox{diag}(D_1) = \{c_n\},$ $\mbox{diag}(D_2) = \{b_n\}$ and $\mbox{diag}(D_3)=\{a_n\}$. We further assume that the subsequences $\{a_{2n-1}\},$ $\{b_{2n-1}\}$, $\{c_{2n-1}\}$ converges to the real numbers $r_1,$ $s_1,$ $s_1$ respectively and $\{a_{2n}\},$ $\{b_{2n}\}$, $\{c_{2n}\}$ converges to the real numbers $r_2,$ $s_2,$ $s_2$ respectively where $s_1 \neq 0$ and $s_2 \neq 0$.

Our focus is to investigate the spectral properties of the operator $T$ using compact perturbation technique. Let us consider another operator $T_0$ over $\ell_p$ defined by
\[T_0 = S_r^2D_1^\prime + D_2^\prime S_{\ell}^2 + D_3^\prime\]
where $D_1^\prime, D_2^\prime, D_3^\prime \in \mathcal{D}_p$ with
\[\mbox{diag}(D_1^\prime) =\mbox{diag}(D_2^\prime) =  \{s_1, s_2, s_1, s_2, \cdots\}, \ \mbox{diag}(D_3^\prime) = \{r_1, r_2, r_1, r_2, \cdots\}.\]
Using the properties of compact operators, it can be proved that $T - T_0$ is a compact operator over $\ell_{p}$. Both the operators $T$ and $T_0$ can be represented by the following penta-diagonal matrices
\[T=\begin{pmatrix}
	a_1 &0& b_1 &0&0&\cdots &  \\
	0 & a_2&0 & b_2& 0&\cdots&\\
	c_1&0&a_3&0&b_3&\cdots&\\
	0&c_2&0&a_4&0&\cdots\\
	0&0&c_3&0&a_5&\cdots\\
	\vdots  &\vdots & \vdots  & \vdots &\vdots&\ddots

	\end{pmatrix},
	T_0=
		\begin{pmatrix}
			r_1 & 0 &s_1&0&0&\cdots &  \\
			0 & r_2 & 0& s_2&0& \cdots&\\
			s_1&0&r_1&0&s_1&\cdots&\\
			0&s_2&0&r_2&0&\cdots&\\
			0&0&s_1&0&r_1&\cdots\\
			\vdots  & \vdots  & \vdots &\vdots&\vdots& \ddots
		\end{pmatrix}.
	\]

  We obtain the spectrum, fine spectrum and the sets of various spectral subdivisions of the operator $T_0.$ It is interesting to note that the spectrum of $T_0$ is given by \[[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2],\] which is also its essential spectrum with no eigenvalues. Later we investigate how the spectrum of $T$ and $T_0$ are related. Few results on the essential spectrum of $T$ (which is identical to the essential spectrum of $T_0$) being devoid of its eigenvalues are also derived. This helps us to characterize the point spectrum of $T.$ Theory of difference equations plays an important role in our study. We use various results on the asymptotic behaviour of solutions of difference equations to demonstrate the findings of our paper. For more details on difference equations one can refer \cite{difference}.

The remainder of paper is organized as follows: section 2 is devoted to introduce some terminologies and results which are relevant to our work. Section 3 contains the results on the spectrum and fine spectrum of $T_0$ over $\ell_p.$ The spectral properties of $T$ are discussed in section 4.

 \section{\bf Preliminaries}
  Let $X$ and $Y$ are Banach spaces and for any operator $T:X \to Y$, $N(T)$ and $R(T)$ denote the null space and range space of $T$ respectively. The operator $T^*:Y^* \to X^*$ is called the adjoint operator and defined by
 \begin{equation*}
 	(T^{*}f)(x)=f(Tx) \quad \text{for all}~ f \in Y^{*} \text{and}~ x\in X
 \end{equation*}
 where $X^*$, $Y^*$ are the dual spaces of $X$ and $Y$ respectively. $B(X)$ denotes the set of all bounded linear operators from $X$ to itself. For any $T\in B(X)$, the resolvent set $ \rho(T, X)$ of $T$ is the set of all $\lambda $ in the complex plane such that $(T-\lambda I)$ has a bounded inverse in $X$ where $I$ is the identity operator defined over $X$. The complement of resolvent set in the complex plane $\mathbb{C} $ is called the spectrum of $T$ and it is denoted by $ \sigma(T, X)$. The spectrum $\sigma(T,X)$ can be partitioned into three disjoint sets which are
 \begin{enumerate}
 	\item[(i)] the point spectrum, denoted by $\sigma_p(T, X)$, is the set of all such $\lambda \in \mathbb{C}$ for which $(T-\lambda I)^{-1}$ does not exist. An element $\lambda\in\sigma_p(T , X)$ is called an eigenvalue of $T$,
 	\item[(ii)] the continuous spectrum, denoted by $\sigma_c(T , X)$, is the set of all such $\lambda \in\mathbb{C}$ for which $(T-\lambda I)^{-1}$ is exists, unbounded and the domain of $(T-\lambda I)^{-1}$ is dense in $X$ but $R(T-\lambda I)\neq X$,
 	\item[(iii)] the residual spectrum, denoted by $\sigma_r(T, X)$, is the set of all such $\lambda \in\mathbb{C}$ for which $(T-\lambda I)^{-1}$ exists (and may be bounded or not) but the domain of $(T-\lambda I)^{-1}$ is not dense in $X$.
 \end{enumerate}
These three disjoint sets are together known as fine spectrum and their union becomes the whole spectrum. There are some other important subdivisions of the spectrum such as approximate point spectrum $\sigma_{app}(T, X)$, defect spectrum $\sigma_{\delta}(T, X)$ and compression spectrum $\sigma_{co}(T, X)$, defined by

\begin{eqnarray*}
	\sigma_{app}(T, X)&=&\{\lambda \in \mathbb{C}:~~ (T - \lambda I)~~ \text{is not bounded below}\},\\
	\sigma_{\delta}(T,X)&=&\left\lbrace \lambda \in \mathbb{C}: (T-\lambda I ) \text{ is not surjective} \right\rbrace,\\
	\sigma_{co}(T,X)&=&\left\lbrace  \lambda \in \mathbb{C}:\overline{R(T-\lambda I)} \neq X\right\rbrace.
\end{eqnarray*}
 The sets which are defined above also forms subdivisions of spectrum of $T$  which are not necessarily disjoint as follows
 \begin{eqnarray*}
 	\sigma(T,X)&=&\sigma_{app}(T,X) \cup \sigma_{co}(T,X),\\
 	\sigma(T,X)&=&\sigma_{app}(T,X) \cup \sigma_{\delta}(T,X).
 \end{eqnarray*}

   An operator $T \in B(X)$ is said to be Fredholm operator if $R(T)$ is closed and  $\dim(N(T))$, $\dim(X/R(T))$ are finite. In this case the number
   \begin{equation*}
   	 \dim(N(T))- \dim(X/R(T))
   \end{equation*} is called the index of the Fredholm operator $T$. The essential spectrum of $T$ is defined by the set \[ \sigma_{ess}(T,\ell_{p})=\left\lbrace \lambda \in \mathbb{C}:(T-\lambda I)~ \text{is not a Fredholm operator}\right\rbrace.\]
   If $T$ is a Fredholm operator and $K \in B(X)$ is a compact operator then $T+K$ is also a Fredholm operator with same indices. Since compact perturbation does not effect the Fredholmness and index of a Fredholm operator, we have \[\sigma_{ess}(T, X)=\sigma_{ess}(T+K, X).\]

For any isolated eigenvalue $\lambda$ of $T,$ the operator $P_T$ which is defined by
\[P_T(\lambda) = \frac{1}{2 \pi i} \int_\gamma  (\mu I - T)^{-1} d \mu,\]
is called the Riesz projection of $T$ with respect to $\lambda$ where $\gamma$ is positively orientated circle centred at $\lambda$ with sufficiently small radius such that it excludes other spectral values of $T.$ An eigenvalue $\lambda$ of $T$ is said to be a discrete eigenvalue if it is isolated and the rank of the associated Riesz projection is finite. The rank of the Riesz projection is called the algebraic multiplicity of $\lambda$. The set of all such eigenvalues with finite multiplicities is called the discrete spectrum of $T$ and it is denoted by $\sigma_d(T,X).$ This type of eigenvalues sometimes referred as eigenvalues with finite type.

In the following proposition, we mention some inclusion relation of spectrum of a bounded linear operator and its adjoint operator.
\begin{proposition} \label{p3.1} \cite[p.195]{basar2012summability}
	If $X$ is a Banach space and $T\in B(X)$, $T^*\in B(X^*)$ then the spectrum and subspectrum of $T$ and $T^*$ are related by the following relations:
	\begin{itemize}
	\item[(a)] $\sigma(T^*,X^*)=\sigma(T,X).$
	\item[(b)] $\sigma_c(T^*,X^*)\subseteq\sigma_{app}(T,X).$
	\item[(c)] $\sigma_{app}(T^*,X^*)=\sigma_{\delta}(T,X).$
	\item[(d)] $\sigma_{\delta}(T^*,X^*)=\sigma_{app}(T,X).$
	\item[(e)] $\sigma_p(T^*,X^*)=\sigma_{co}(T,X).$
	\item[(f)] $\sigma_{co}(T^*,X^*)\supseteq\sigma_{p}(T,X).$
	\item[(g)] $ \sigma(T,X)=\sigma_{app}(T,X)\cup\sigma_{p}(T^*,X^*)$=$\sigma_p(T,X)\cup\sigma_{app}(T^*,X^*).$
	\end{itemize}
\end{proposition}
Here we record few lemmas related to the boundness of an infinite matrix defined over sequence spaces, which are useful to our research.
\begin{lemma}\cite[p. 253]{choudhary1989functional} \label{l1}
The matrix $A=(a_{nk})$ gives rise to a bounded linear operator $T \in B(\ell_{1})$ from $\ell_1$ to itself if and only if the supremum of $\ell_1$ norms of  the columns of $A$ is bounded.
\end{lemma}
\begin{lemma}\cite[p. 245]{choudhary1989functional} \label{l2}
The matrix $A=(a_{nk})$ gives rise to a bounded linear operator $T \in B(\ell_\infty)$  from $\ell_\infty$ to itself if and only if the supremum of $\ell_1$  norms of the rows of $A$ is bounded.
\end{lemma}
\begin{lemma}\cite[p. 254]{choudhary1989functional}\label{L3.5}
The matrix $A=(a_{nk})$ gives rise to a bounded linear operator $T\in B(\ell_p)(1<p<\infty)$ if $T \in B(\ell_1)\cap B(\ell_{\infty}).$
\end{lemma}

\section{\bf Spectra of $T_0$}

  It is already mentioned that we study the spectral properties of $T$ by using the spectral properties of $T_0$ and compact perturbation technique. In this section we derive the spectrum and fine spectrum of $T_0$. The notation $\|T\|_p$ denotes the operator norm of an operator $T \in B(\ell_{p})$ where $1 \leq p \leq \infty$.
\begin{theorem} The operator
	$T_0:\ell_p \rightarrow \ell_p$ is a bounded linear operator which satisfies the following inequality
  $$	\left(\frac{|r_1|^p+|r_2|^p+|s_1|^p+|s_2|^p}{2}\right)^{\frac{1}{p}}\leq \norm{T_0}_{p}\leq \left( 3^{p-1}\left(|r_1|^p+2|s_1|^p+|r_2|^p+2|s_2|^p\right)\right)^{\frac{1}{p}}.$$
\end{theorem}
\begin{proof}
As linearity of $T_0$ is trivial, we omit it. Let $e=(1,1,0,0,...) \in\ell_{p}$. Then $T_{0}(e)=(r_1,r_2,s_1,s_2,0,...)$ and one can observe that
\begin{equation*}
\frac{\norm{T_0(e)}_{p}}{\norm{e}_{p}}=\left( \frac{|r_1|^p+|r_2|^p+|s_1|^p+|s_2|^p}{2}\right)^{\frac{1}{p}}.
\end{equation*}
This proves \[ \left( \frac{|r_1|^p+|r_2|^p+|s_1|^p+|s_2|^p}{2}\right)^{\frac{1}{p}}\leq \norm{T_0}_{p}.\]

\noindent Also, let $x=\{x_n\}\in \ell_p$ and $x_n=0$ if $n \leq0$.
Then,
\begin{align*}
\norm{T_0(x)}_{p}^p=&\sum_{n=1}^{\infty}|s_{1}x_{2n-3}+r_1x_{2n-1}+s_1x_{2n+1}|^p+	\sum_{n=1}^{\infty}|s_2x_{2n-2}+r_2x_{2n}+s_2x_{2n+2}|^p\\
\leq& \sum_{n=1}^{\infty}(|s_{1}x_{2n-3}|+|r_1 x_{2n-1}|+|s_1 x_{2n+1}|)^p+	\sum_{n=1}^{\infty}(|s_2x_{2n-2}|+|r_2x_{2n}|+|s_2x_{2n+2}|)^p\\
\end{align*}
 By Jensen’s inequality we get,
\begin{align*}
\norm{T_0(x)}_{p}^p \leq & 3^{p-1}\sum_{n=1}^{\infty}\left( |s_1x_{2n-3}|^p+|r_1x_{2n-1}|^p+|s_1x_{2n+1}|^p\right)\\ +& 3^{p-1}\sum_{n=1}^{\infty}\left( |s_2x_{2n-2}|^p+|r_2x_{2n}|^p+|s_2x_{2n+2}|^p\right)\\
\leq& 3^{p-1}\left(|r_1|^p+2|s_1|^p+|r_2|^p+2|s_2|^p\right){\norm{x}_p^p}.
\end{align*}
This implies,
$$\norm{T_0} \leq \left( 3^{p-1}\left(|r_1|^p+2|s_1|^p+|r_2|^p+2|s_2|^p\right)\right)^{\frac{1}{p}}.$$ This completes the proof.
\end{proof}
The following theorem proves the non-existence of eigenvalues of the operator $T_0$ in $\ell_p$.
\begin{theorem}\label{t1.2}
The point spectrum of $T_0$ is given by $\sigma_p(T_0,\ell_p)=\emptyset.$
\end{theorem}
\begin{proof}
Consider $(T_0 - \lambda I)x=0$ for $\lambda\in\mathbb{C}$ and $x=\{x_n\}\in \mathbb{C^N}$. This gives the following system of equations

	\begin{eqnarray*}
		\begin{rcases*}
		(r_{1} - \lambda) x_{1}+s_{1} x_{3}&=0\\
		(r_{2} - \lambda) x_{2}+s_{2} x_{4}&=0 \\
		s_1 x_1 + (r_1 - \lambda) x_3 + s_1 x_5 &= 0\\
		s_2x_2+(r_2 - \lambda) x_4+s_2x_6&=0 \\
		                       &\vdots \\
    	s_1 x_{2n-1}+(r_1 - \lambda) x_{2n+1}+s_1 x_{2n+3}&=0\\
		s_2 x_{2n}+(r_2-\lambda)x_{2n+2}+s_2 x_{2n+4}&=0\\
		                             &\vdots
	\end{rcases*}
\end{eqnarray*}
We assume that $x_1 \neq 0$ and $x_2 \neq 0$ otherwise $x_n = 0$ for all $n \in \mathbb{N}.$ Let us consider two sequences $\{y_n\}$ and $\{z_n\}$ where  $y_n=x_{2n-1}$ and $z_n=x_{2n}$, $n \in \mathbb{N}$ respectively. Then the system of equations of $(T_0-\lambda I)x=0$ reduces to
 \begin{equation}\label{3.1}
  y_{n}+p_1y_{n+1}+y_{n+2}=0,

 \end{equation}
 \begin{equation} \label{3.2}
 z_{n}+p_2z_{n+1}+z_{n+2}=0,\\

 \end{equation}
where $p_1=\frac{r_1-\lambda}{s_1}$, $p_2=\frac{r_2-\lambda}{s_2}$, $n \in \mathbb{N}\cup \{0\}$ and $y_0=z_0=0$. The general solution of the difference equation (\ref{3.1}) is given by
\begin{equation} \label{eqq1}
y_n = \begin{cases}
& (c_1 + nc_2) (-1)^n, \ \ \mbox{if} \ p_1 = 2\\
& c_1+ n c_2, \ \ \mbox{if} \ p_1 = -2\\
& c_1 \alpha_1^n + c_2 \alpha_2^n, \ \ \mbox{if} \ p_1 \notin \{-2,2\}
\end{cases}
\end{equation}
where $c_1,$ $c_2$ are arbitrary constants and $\alpha_1,$ $\alpha_2$ are the roots of the polynomial
\begin{equation} \label{ch1}
y^2+p_1y+1=0
\end{equation}
which is called the characteristic polynomial of (\ref{3.1}). The following two equalities
\begin{equation*}
\alpha_1\alpha_2=1  \text{ and }  \alpha_1+\alpha_2= -p_1
\end{equation*}
are useful.
There are three cases to be considered.\\

\noindent Case 1: If $p_1=2$ (i.e., $\lambda =r_1-2s_1 $). In this case the general solution of (\ref{3.1}) is
 \begin{equation*}
 	y_{n}=(c_1+c_2n)(-1)^{n}, \ n \in \mathbb{N}\cup\{0\}
 \end{equation*}
with the initial condition $y_{0}=0$ which gives $c_1 = 0.$ This reduces the solution as $y_n = nc_2 (-1)^n.$ This also implies $c_2 = -y_1$ and the solution in this case is
\[y_n = ny_1 (-1)^{n+1}, ~~ n \in \mathbb{N}.\]

\noindent Case 2: If $p_1=-2$ (i.e., $\lambda =r_1+2s_1 $). Similar as Case 1, the solution reduces to
\[y_n = ny_1,~~ n \in \mathbb{N}.\]

\noindent Case 3: If $p_1 \notin \{-2, 2\}.$ The general solution of (\ref{3.1}) is given by
 \begin{equation*}
 	y_{n}=c_1 \alpha_1^n+c_2 \alpha_2^n.
 \end{equation*}
With the help of initial condition $y_0=0$ and by using the equalities $\alpha_1\alpha_2=1$, $\alpha_1+\alpha_2= -p_1$, one can obtain that $c_2 = -c_1$ and the solution reduces to
\[y_n = \frac{\alpha_1^n - \alpha_2^n}{\alpha_1 - \alpha_2} y_1,~~ n \in \mathbb{N}.\]
If $y_1 \neq 0$ then  $\{y_n\} \notin \ell_p$ in Case 1 and Case 2. In Case 3 $\{y_n\} \in \ell_p$ if and only if $|\alpha_1| <1$ and $|\alpha_2|<1$ which can not be the case since $\alpha_1 \alpha_2  = 1.$ Hence in all the three cases $\{y_n\} \in \ell_p$ if and only if $y_1 = 0$ and this leads to the trivial solution of the difference equation (\ref{3.1}). Hence there is no non-trivial solution of ({\ref{3.1}}).

Similarly for the difference equation (\ref{3.2}), the general solution $\{z_n\}$ is of the form
\begin{equation} \label{eqq2}
z_n = \begin{cases}
& (d_1 + nd_2) (-1)^n, \ \ \mbox{if} \ p_2 = 2\\
& d_1+ n d_2, \ \ \mbox{if} \ p_2 = -2\\
& d_1 \beta_1^n + d_2 \beta_2^n, \ \ \mbox{if} \ p_2 \notin \{-2,2\}
\end{cases}
\end{equation}
where $d_1,$ $d_2$ are arbitrary constants and $\beta_1,$ $\beta_2$ are the roots of the polynomial
\begin{equation} \label{ch2}
z^2+p_2z+1=0
\end{equation}
which is called the characteristic polynomial of difference equation (\ref{3.2}). In a similar way, it can be proved that $\{z_n\} \in \ell_p$ if and only if $z_1 = 0$ and this leads to the trivial solution of (\ref{3.2}). Hence, there does not exist any non-trivial solution of the system $(T_0 - \lambda I)x = 0$ such that $x \in \ell_p.$ This proves the required result.
\end{proof}

\begin{remark}
The solution $x = \{x_n\}$ of the system $Tx = \lambda x,$ which are obtained in terms of the sequences $\{y_n\}$ and $\{z_n\}$ in the equations (\ref{eqq1}) and (\ref{eqq2}) respectively, actually depends on the unknown $\lambda.$ Therefore, instead of writing $x_n(\lambda),$ we write $x_n$ for the sake of brevity throughout this paper except in Theorem \ref{t4.12} where the dependency of the solutions on $\lambda$ is vital.
\end{remark}

It is well known that the adjoint operator of $T_0$ is $T_0^*$ which is defined over sequence space $\ell_p^*$ where $\ell_p^*$ denotes the dual space of $\ell_p$ which is isomorphic to $\ell_q$ where $\frac{1}{p}+ \frac{1}{q}=1$.
\begin{corollary}\label{c1.3}
	The point spectrum of adjoint operator $T_0^*$ over the sequence space $\ell_p^*$ is given by
	 $\sigma_p(T_0^*,\ell_p^*)=\emptyset.$
\end{corollary}
\begin{proof}
It is well known that the adjoint operator $ {T_0}^*:\ell_p^* \to \ell_p^* $, is represented by transpose of the matrix $T_0$. Since $T_0$ is represented by a symmetric matrix, using the same argument as Theorem \ref{t1.2}, it is easy to prove that $\sigma_p(T_0^*,\ell_p^*)=\emptyset.$
\end{proof}
\begin{corollary}
	The residual spectrum of $T_0$ over the sequence space $\ell_{p}$ is given by
$\sigma_r(T_0,\ell_p)=\emptyset$.
\end{corollary}
\begin{proof}
	We know that the operator $T$ has a dense range if and only if $T^*$ is one to one \cite[p.197]{basar2012summability}. Using this we have the following relation
\[\sigma_r(T_0,\ell_p)=\sigma_p(T_0^*,\ell_p^*)\setminus\sigma_p(T_0,\ell_p).\]
Hence, $\sigma_r(T_0,\ell_p)=\emptyset.$
\end{proof}

Following that, we obtain the spectrum of $T_0$.
\begin{theorem}
	The spectrum of $T_0$ is given by
	\begin{equation*}
	\sigma(T_0,\ell_p) =\left[r_{1}-2s_1, r_1+2s_1\right] \cup [r_2-2s_2, r_2+2s_2].
	\end{equation*}
\end{theorem}

\begin{corollary}
The continuous spectrum of $T_0$ is given by
	$$\sigma_c(T_0,\ell_p)=\left[r_{1}-2s_1, r_1+2s_1\right] \cup [r_2-2s_2, r_2+2s_2].$$
	\end{corollary}
	\begin{proof}
		It is evident that $\sigma(T_0,\ell_{p})$ is the disjoint union of $\sigma_p(T_0,\ell_{p})$,  $\sigma_r(T_0,\ell_{p})$ and  $\sigma_c(T_0,\ell_{p})$, we have
		\begin{eqnarray*}
			\sigma(T_0,\ell_p)=\sigma_c(T_0,\ell_p).
		\end{eqnarray*}
		Hence, $ \sigma_c(T_0,\ell_p)=\left[r_{1}-2s_1, r_1+2s_1\right] \cup [r_2-2s_2, r_2+2s_2]$.
	\end{proof}

\begin{corollary}\label{ess}
Essential spectrum of $T_0$ is given by $$\sigma_{ess}(T_0,\ell_{p})=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2].$$
\end{corollary}
\begin{proof}
It is well-known that $\sigma_c(T_0, \ell_p) \subseteq \sigma_{ess}(T_0,\ell_p)$ and we have $$\sigma_{ess}(T_0,\ell_p)\subseteq \sigma(T_0,\ell_{p}) = \sigma_c(T_0,\ell_{p})\subseteq \sigma_{ess}(T_0, \ell_{p}).$$ The desired result is obvious.
\end{proof}
Using the relations which are mentioned in Proposition \ref{p3.1} we can easily obtain the following results.
\begin{corollary}
The compression spectrum, approximate point spectrum and defect spectrum of $T_0$ are as follows
\item[(i)] $\sigma_{co}(T_0,\ell_p)=\emptyset$.
\item[(ii)]$\sigma_{app}(T_0,\ell_p)=\left[r_{1}-2s_1, r_1+2s_1\right] \cup [r_2-2s_2, r_2+2s_2]$.
\item[(iii)] $\sigma_{\delta}(T_0,\ell_p)=\left[r_{1}-2s_1, r_1+2s_1\right]\cup[r_2-2s_2, r_2+2s_2]$.
\end{corollary}

\section{\bf Spectra of $T = T_0 + K$}

This section contains the spectral properties of $T$ which can be expressed as $T=T_0+K$  where $K=T-T_0$ is represented by the following matrix
\begin{equation*}
	K=
	\begin{pmatrix}
		a_1-r_1 & 0 &b_1-s_1&0&0&\cdots &  \\
		0 & a_2-r_2 & 0& b_2-s_2&0& \cdots&\\
		c_1-s_1&0&a_3-r_1&0&b_3-s_1&\cdots&\\
		0&c_2-s_2&0&a_4-r_2&0&\cdots&\\
		0&0&c_3-s_1&0&a_5-r_1&\cdots\\
		\vdots  & \vdots  & \vdots &\vdots&\vdots& \ddots
	\end{pmatrix}.
\end{equation*}
The following result proves the compactness of $K$ on $\ell_{p}$.

\begin{theorem}\label{theorem1.6}
  The operator $K$ is a compact operator on $\ell_{p}$.

\end{theorem}
\begin{proof}
	The operator $K$ on $\ell_{p}$ can be represented by the following infinite matrix
	 \begin{equation*}
		K=
		\begin{pmatrix}
			u_1 & 0 &v_1&0&0&\cdots &  \\
			0 & u_2 & 0& v_2&0& \cdots&\\
			w_1&0&u_3&0&v_3&\cdots&\\
			0&w_2&0&u_4&0&\cdots&\\
			\vdots  & \vdots  & \vdots &\vdots&\vdots& \ddots
		\end{pmatrix},
	\end{equation*}
where, $\{u_n\}$, $\{v_n\}$ and $\{w_n\}$ are null sequences, which are defined as follows,

\begin{equation*}
	u_n=\begin{cases*}
		a_n-r_1,~~ \mbox{$n$ is odd}\\
		a_n-r_2, ~~ \mbox{$n$ is even},
	\end{cases*}~~~~\qquad
	v_n=\begin{cases*}
		b_n-s_1,~~ \mbox{ $n$ is odd}\\
		b_n-s_2,~~ \mbox{$n$ is even}
	\end{cases*}
\end{equation*}
and
\begin{equation*}
w_n=\begin{cases*}
	c_n-s_1,~~ \mbox{ $n$ is odd}\\
	c_n-s_2,~~ \mbox{$n$ is even}.
\end{cases*}
\end{equation*}
Let $x=\{x_1,x_2,x_3...\}\in \ell_p$.
 We construct a sequence of compact operators $\left\lbrace K_n\right\rbrace $ such that for $i \in \mathbb{N}$, \\
\begin{equation*}
(K_n(x))_i=\begin{cases}
	(Kx)_i,~~ i=1,2,...n\\
	0,~~ \text{otherwise.}
\end{cases}
\end{equation*}

For $n \geq2$,
\begin{align*}
\norm{(K-K_n)x}_{p}
= & \left( \sum_{k=n-1}^{\infty}|w_kx_k+u_{k+2}x_{k+2}+v_{k+2}x_{k+4}|^p\right) ^{\frac{1}{p}}\\
\leq & \left( \sup_{k\geq n-1}|w_k|\right)\norm{x}_p + \left( \sup_{k \geq n-1}|u_k|\right)\norm{x}_p+ \left( \sup_{k\geq n-1}|v_k|\right) \norm{x}_p
\end{align*}
This implies,
$$ \|{K-K_n}\|_p ~\leq  \sup_{k \geq n-1}|w_k|+\sup_{k \geq n-1 }|u_k|+\sup_{k \geq n-1}|v_k|.$$
Thus, $\{K_n\}$ converges to $K$ as $n\to\infty$ in operator norm and hence $K$ is a compact operator over $\ell_p$.
\end{proof}

Next we derive an inclusion relation between $\sigma(T_0, \ell_{p})$ and $\sigma(T,\ell_{p})$.
\begin{theorem}\label{t4.11}
The spectrum of $T$ satisfies the following inclusion relation $$\sigma(T_0,\ell_p) \subseteq \sigma(T,\ell_p)$$ and $\sigma(T,\ell_p)\setminus \sigma(T_0,\ell_p)$ contains finite or countable number of eigenvalues of $T$ of finite type with no accumulation point in $\sigma(T,\ell_p)\setminus \sigma(T_0,\ell_p)$.
\end{theorem}

\begin{corollary}\label{ess1}
	$\sigma_{ess}(T, \ell_{p})=\sigma(T_0, \ell_{p})=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]$.
\end{corollary}
\begin{proof}
	As we are aware of that compact perturbation does not effect the Fredholmness and index of a Fredholm operator. Therefore, $\sigma_{ess}(T_0,\ell_{p})=\sigma_{ess}(T,\ell_{p})$. Hence, by using Corollary \ref{ess}, we have
	\begin{equation*}
		\sigma_{ess}(T, \ell_{p})=\sigma(T_0, \ell_{p})=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2].
	\end{equation*}
\end{proof}
We now focus on the point spectrum of $T$. First we analyze the eigenvalues of $T$ lying in $\sigma_{ess}(T, \ell_{p})=\sigma(T_0, \ell_p)$, in particular we derive sufficient conditions for the absence of point spectrum on $\sigma_{ess}(T, \ell_p)$. In Theorem \ref{poT}, sufficient conditions are provided in terms of the rate of convergence of the sequences $\{a_{2n-1}\}, \ \{a_{2n}\},  \ \{b_{2n-1}\}, \ \{b_{2n}\}, \ \{c_{2n-1}\}$ and $\{c_{2n}\}$. Sufficient conditions of absence of point spectrum on $\sigma_{ess}(T, \ell_{p})$ are also provided in Theorem \ref{sc} in terms of the entries of the matrix $T$.

\begin{theorem}\label{poT}
If the convergence of the sequences $\{a_{2n-1}\}, \ \{a_{2n}\}, \ \{b_{2n-1}\}, \ \{b_{2n}\}, \\ \{c_{2n-1}\}$ and $\{c_{2n}\}$ are exponentially fast then $$\sigma_{ess}(T, \ell_{p}) \cap \sigma_{p}(T,\ell_{p})= \emptyset.$$
\end{theorem}

In the next theorem, we apply transfer matrix approach as discussed in \cite{janas2003similarity, janas1998point} . This enables us to examine the sufficient condition for the absence of point spectrum in essential spectrum of $T$ in terms of the entries of matrix $T$.
\begin{theorem}\label{sc}
	If \(\lambda \in \sigma_{ess}(T,\ell_{p})\) satisfies either of the following conditions\\
	(i) $ \sum_{n=1}^{\infty} \prod_{j=1}^n \left[\frac{1}{2} \left(P_j(\lambda)-\sqrt{P_j(\lambda)^2-\left| \frac{2c_{2j+1}}{b_{2j+1}}\right|^2}\right) \right]^{\frac{1}{2}}=+\infty$ \\
	\\ or\\
	(ii) $ \sum_{n=1}^{\infty} \prod_{j=1}^n \left[\frac{1}{2} \left( Q_j(\lambda)-\sqrt{Q_j(\lambda)^2-\left| \frac{2c_{2j+2}}{b_{2j+2}}\right|^2}\right) \right]^{\frac{1}{2}}=+\infty$\\
	where,
	$$  P_j(\lambda)= \left| \frac{c_{2j-1}}{b_{2j+1}}\right| ^2 +\left| \frac{a_{2j+1-\lambda}}{b_{2j+1}}\right|^2+1,~~Q_j(\lambda)= \left| \frac{c_{2j}}{b_{2j+2}}\right| ^2 +\left| \frac{a_{2j+2-\lambda}}{b_{2j+2}}\right|^2+1, $$
	then \(\lambda \notin \sigma_{p}(T,\ell_{p})\).
\end{theorem}
Now we focus our study on the point spectrum of $T$.
Under the sufficient conditions as mentioned in previous two results, we have
$\sigma_p(T, \ell_{p})\cap \sigma(T_0, \ell_{p})=\emptyset$. In this case, all the eigenvalues of $T$ are lying outside the set $\sigma(T_0, \ell_p)$.
To characterize the eigenvalues, let $Tx=\lambda x$, $x \in \mathbb{C^N}$ and $\lambda \in \sigma(T_0,\ell_{p})^c$ where $\sigma(T_0,\ell_{p})^c$ denotes the complement of $\sigma(T_0,\ell_{p})$. From equations (\ref{point_diff}) and (\ref{point_diff1}) in Theorem \ref{poT}, we have the following system
 \begin{equation} \label{point_diff4}
 	c_{2n-1}y_n+(a_{2n+1}-\lambda) y_{n+1}+b_{2n+1}y_{n+2}=0,
 \end{equation}
 \begin{equation}\label{point_diff5}
 	c_{2n}z_n+(a_{2n+2}- \lambda)z_{n+1}+b_{2n+2}z_{n+2}=0,
 \end{equation}
where $n \in \mathbb{N}\cup\{0\}$ with $y_0=z_0=0$ and $y_n=x_{2n-1}$, $z_n=x_{2n}$.
Clearly each of the difference equations (\ref{point_diff4}) and (\ref{point_diff5}) have two fundamental solutions. Let $\{y_n^{(1)}(\lambda),y_n^{(2)}(\lambda)\}$ and $\{z_n^{(1)}(\lambda),z_n^{(2)}(\lambda)\}$ are the sets of fundamental solutions of the equations  (\ref{point_diff4}) and (\ref{point_diff5}) respectively. Under this setting we have the following result:
\begin{theorem}\label{t4.12}
If either of the sufficient conditions mentioned in Theorem \ref{poT} and Theorem \ref{sc} hold true then the point spectrum of $T$ is given by
\[\sigma_p(T,\ell_p)=\left\lbrace \lambda \in \mathbb{C}: y_{0}^{(1)}(\lambda) =0\right\rbrace \cup \left\lbrace \lambda \in \mathbb{C}:{z_0}^{(1)}(\lambda)=0\right\rbrace .\]
\end{theorem}

\begin{remark}\label{remark1}
	The adjoint operator $ {T}^*:\ell_p^*\to \ell_p^* $, is represented by transpose of the matrix $T$ and dual of $\ell_p$ is isomorphic to $\ell_q$ where $\frac{1}{p}+ \frac{1}{q}=1$ and $1<q<\infty$. Similar as $T$, the operator $T^*$ can also be written as
	\begin{equation*}
		T^*=T_0+K^t,
	\end{equation*}
where $K^t$ denotes the transpose of $K$ and $K^t$ is also a compact operator. Since $\sigma(T, \ell_{p})=\sigma(T^*, \ell_{p}^*)$, Theorem \ref{t4.11} implies $$\sigma(T_0, \ell_{p})\subseteq \sigma(T^*, \ell_{p}^*),$$ and using similar argument of the proof of Theorem \ref{t4.11} it can be obtain that $\sigma(T^*,\ell_{p}^*)\setminus \sigma(T_0,\ell_p)$ contains finite or countable number of eigenvalues of $T^*$ of finite type with no accumulation point in $\sigma(T^*,\ell_{p}^*)\setminus \sigma(T_0,\ell_p)$. Assuming similar hypothesis on the rate of convergence of sequences in Theorem \ref{poT},
 we can prove that $$\sigma_{ess}(T^*, \ell_{p}^*)\cap \sigma_p(T^*,\ell_{p}^*)=\emptyset$$
 and this implies, the point spectrum of $T^*$ is lying outside of the region \[[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2].\]
Now similar as Theorem \ref{t4.12},
let $\{g_n^{(1)}(\lambda), g_n^{(2)}(\lambda)\}$ and $\{h_n^{(1)}(\lambda), h_n^{(2)}(\lambda)\}$ are the sets of fundamental solutions of the following difference equations respectively
\begin{eqnarray*}
	b_{2n-1}g_n+(a_{2n+1}-\lambda)g_{n+1}+c_{2n+1}g_{n+2}=0,\\
	b_{2n}h_n+(a_{2n+2}-\lambda)h_{n+1}+c_{2n+2}h_{n+2}=0,
\end{eqnarray*}
 which are obtained from $T^*f= \lambda f$,
$f \in\ell_{p}^*$ and $g_n(\lambda)=f_{2n-1}(\lambda)$, $h_n(\lambda)=f_{2n}(\lambda)$. Also, $g_0(\lambda)=h_0(\lambda)=0$.
This leads us to the following result
\begin{equation*}
\sigma_p(T^*, \ell_{p}^*)=\left\lbrace \lambda \in \mathbb{C}:g_0^{(1)}(\lambda)=0\right\rbrace \cup \left\lbrace \lambda \in \mathbb{C}:h_0^{(1)}(\lambda)=0\right\rbrace.
\end{equation*}
Eventually, we obtain that
$$\sigma(T^*,{\ell_p^*})=[r_1-2s_1, r_1+2s_1]\cup [r_2-2s_2, r_2+2s_2] \cup S_2. $$
where,
\begin{equation*}
	S_2=\left\lbrace \lambda \in \mathbb{C}:g_0^{(1)}(\lambda)=0\right\rbrace \cup \left\lbrace \lambda \in \mathbb{C}:h_0^{(1)}(\lambda)=0\right\rbrace.
\end{equation*}
Since, $\sigma(T,\ell_{p})=\sigma(T^*,\ell_{p}^*)$ and $S_1,S_2$ both sets are disjoint from $[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]$ we have, $S_1=S_2$.

\end{remark}
 Using the observations in Remark \ref{remark1} and Proposition \ref{p3.1}, we can summarize all the results of spectrum and various spectral subdivisions of the operator $T$ in the following theorem,
	\begin{theorem}
		If the convergence of the sequences $\{a_{2n-1}\}, \ \{a_{2n}\}, \ \{b_{2n-1}\}, \ \{b_{2n}\},\\ \ \{c_{2n-1}\}$ and $\{c_{2n}\}$ are exponentially fast then the following hold,
	\begin{itemize}
		\item[(i)] The spectrum of $T$ on $\ell_p$ is
		\[\sigma(T,\ell_p)=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.\]
		\item[(ii)] The point spectrum of $T$ on $\ell_{p}$ is
		 \[\sigma_p(T,\ell_p)=\left\lbrace \lambda \in \mathbb{C}: y_{0}^{(1)}(\lambda) =0\right\rbrace \cup \left\lbrace \lambda \in \mathbb{C}:{z_0}^{(1)}(\lambda)=0\right\rbrace .\]
		 \item[(iii)] The residual spectrum of $T$ on $\ell_{p}$ is
		 	\[\sigma_r(T, \ell_p)=\emptyset.\]
		 \item[(iv)] The continuous spectrum of $T$ on $\ell_{p}$ is
		 \[\sigma_c(T,\ell_{p})=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2].\]
		 \item[(v)] The essential spectrum of $T$ on $\ell_{p}$ is
		 \[\sigma_{ess}(T,\ell_{p})=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2].\]
		 \item[(vi)] The discrete spectrum of $T$ on $\ell_{p}$ is
		 \[\sigma_d(T, \ell_{p})=\left\lbrace \lambda \in \mathbb{C}: y_{0}^{(1)}(\lambda) =0\right\rbrace \cup \left\lbrace \lambda \in \mathbb{C}:{z_0}^{(1)}(\lambda)=0\right\rbrace.\]
		\item [(vii)] The compression spectrum of $T$ on $\ell_{p}$ is \[\sigma_{co}(T,\ell_p)=\left\lbrace \lambda \in \mathbb{C}: y_{0}^{(1)}(\lambda) =0\right\rbrace \cup \left\lbrace \lambda \in \mathbb{C}:{z_0}^{(1)}(\lambda)=0\right\rbrace.\]
		\item[(viii)] The approximate spectrum of $T$ on $\ell_{p}$ is \[\sigma_{app}(T,\ell_p)=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.\]
		\item[(ix)] The defect spectrum of $T$ on $\ell_{p}$ is \[\sigma_{\delta}(T,\ell_p)=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.\]
	\end{itemize}
\end{theorem}
\begin{proof} The proofs of the above statements are given below.
	\begin{enumerate}
			\item[(i)] It is well known that $\sigma_p(T,\ell_p)\subseteq \sigma(T,\ell_p)$ and $\sigma(T_0,\ell_p)\subseteq \sigma(T,\ell_p)$. This implies, $$[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1\subseteq \sigma(T,\ell_p).$$ Also by using Theorem (\ref{t4.11}) we get, $$\sigma(T,\ell_p)\subseteq[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.$$ Hence,
		$$\sigma(T,\ell_p)=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.$$
		\item[(ii)]  Result has been proved in Theorem \ref{t4.12}.
		\item [(iii)] We already aware of that $\sigma_r(T,\ell_p)=\sigma_p(T^*,\ell_p^*)\setminus\sigma_p(T,\ell_p)$.
		Hence, $$\sigma_r(T, \ell_p)=\emptyset.$$
		\item [(iv)] Spectrum of an operator is the disjoint union of point spectrum, residual spectrum and continuous spectrum. By using this result we can obtain the desired result.
		\item [(v)]
		The required result has been proved in Corollary \ref{ess1}.

		\item [(vi)] We already proved that the point spectrum of $T$ is disjoint from $[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]$, and by Theorem \ref{t4.11} we have, every element of $\sigma_{p}(T, \ell_{p})$ is of finite type. Hence,
		\[\sigma_d(T, \ell_{p})=\left\lbrace \lambda \in \mathbb{C}: y_{0}^{(1)}(\lambda) =0\right\rbrace \cup \left\lbrace \lambda \in \mathbb{C}:{z_0}^{(1)}(\lambda)=0\right\rbrace.\]
		\item [(vii)] By part (e) of Proposition \ref{p3.1}, the desired result is obvious.
		\item [(viii)] Clearly, $$\sigma_{app}(T, \ell_{p})\subseteq \sigma(T,\ell_p)=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.$$ Also, we know that point spectrum is always a subset of approximate point spectrum.  By using this fact and with the help of part (g) of Proposition \ref{p3.1}, we have
		\[[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1 \subseteq \sigma_{app}(T, \ell_{p}).\]
		Hence,\[\sigma_{app}(T, \ell_{p})=\sigma(T, \ell_{p})=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.\]
		\item [(ix)] By part (c) of Proposition \ref{p3.1}, we have
		\[\sigma_{app}(T^*, \ell_{p}^*)=\sigma_{\delta}(T,\ell_{p}).\]
		Clearly, $$\sigma_{app}(T^*, \ell_{p}^*)\subseteq \sigma(T^*,\ell_{p}^*)=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.$$
		As $\sigma_{p}(T^*, \ell_{p}^*)\subseteq\sigma_{app}(T^*, \ell_{p}^*)$ thus
		$S_1\subseteq\sigma_{app}(T^*, \ell_{p}^*)$.  By using this fact and with the help of part (g) of Proposition \ref{p3.1}, we have
		\[[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1 \subseteq \sigma_{app}(T^*, \ell_{p}^*).\] Therefore, $\sigma_{app}(T^*, \ell_{p}^*)=\sigma(T,\ell_{p})$. Hence,
		\[\sigma_{\delta}(T,\ell_p)=[r_1-2s_1,r_1+2s_1]\cup[r_2-2s_2,r_2+2s_2]\cup S_1.\]
	\end{enumerate}
\end{proof}

\section*{\bf Disclosure statement}
No potential conflict of interest was reported by the authors.

\bibliographystyle{plain}