Eigenvectors and eigenvalues of nonsymmetric Tridiagonal matrix

Question

Hi, the question is following: We have one matrix $$\begin{pmatrix} -\beta & \Delta & 0 & 0 &\cdots & 0 & 0 & 0 \newline \beta & -(\beta+\Delta) & \Delta & 0 &\cdots & 0 & 0 & 0 \newline 0 & \beta & -(\beta + \Delta) & \Delta &\cdots & 0 & 0 & 0 \newline \vdots & \vdots & \vdots & \vdots & \ddots & \vdots & \vdots& \vdots \newline 0 & 0 & 0 & 0 &\cdots & \beta & -(\beta + \Delta) & \Delta \newline 0 & 0 & 0 & 0 & \cdots & 0 & \beta & -(\beta + \Delta) \end{pmatrix}. $$

Is it possible to find analytically all eigenvalues and eigenvectors of this matrix? What I found here is quite similar, but not exactly the same. The element in the last row and last column $-(\beta +\Delta)$ can be replaced by $-\Delta$ if it simplifies the solution. $\beta>0$, $\Delta>0$.

Thanks

Answer 1

As Denis remarked, the characteristic polynomial $P_n(\lambda)$ of the $n \times n$ matrix does satisfy a three-term linear recurrence $$P_{n+2}(\lambda) = (\lambda + \beta + \Delta) P_{n+1}(\lambda) - \beta \Delta P_n(\lambda)$$ with initial conditions $P_1(\lambda) = \lambda+\beta$, $P_2(\lambda) = \lambda^2+(2 \beta+\Delta)\lambda +\beta^2$, and then $$P_n(\lambda) = \frac{(P_2 - r_2 P_1)}{r_1 (r_1 - r_2)} r_1^n + \frac{(P_2 - r_1 P_1)}{r_2 (r_2 - r_1)} r_2^n $$ where $r_1$ and $r_2$ are the roots of $r^2 - (\lambda + \beta + \Delta) r + \beta \Delta$. However, I don't see how this implies a "closed form" for the eigenvalues.

EDIT: By scaling, we may as well assume $\Delta=1$. I don't know if this leads to a "closed form", but the eigenvalue $\lambda_0$ that is analytic in a neighbourhood of $\beta=0$ has some interesting regularities in its Maclaurin series:

$$\matrix{n=2 & -\beta^{2}+2 \beta^{3}-5 \beta^{4}+14 \beta^{5}-42 \beta^{6}\cr&+132 \beta^{7}-429 \beta^{8}+1430 \beta^{9}-4862 \beta^{10}\cr n=3 & -\beta^{3}+2 \beta^{4}-\beta^{5}-6 \beta^{6}+20 \beta^{7}\cr&-22 \beta^{8}-49 \beta^{9}+260 \beta^{10}-441 \beta^{11}\cr&-320 \beta^{12}+3652 \beta^{13}\cr n=4 & -\beta^{4}+2 \beta^{5}-\beta^{6}-8 \beta^{8}+26 \beta^{9}\cr&-28 \beta^{10}+10 \beta^{11}-100 \beta^{12}+442 \beta^{13}\cr&-729 \beta^{14}+532 \beta^{15}-1641 \beta^{16}\cr n=5 & -\beta^{5}+2 \beta^{6}-\beta^{7}-10 \beta^{10}+32 \beta^{11}\cr&-34 \beta^{12}+12 \beta^{13}-155 \beta^{15}+672 \beta^{16}\cr&-1089 \beta^{17}+782 \beta^{18}-210 \beta^{19}\cr n=6 & -\beta^{6}+2 \beta^{7}-\beta^{8}-12 \beta^{12}+38 \beta^{13}\cr&-40 \beta^{14}+14 \beta^{15}-222 \beta^{18}+950 \beta^{19}\cr&-1521 \beta^{20}+1080 \beta^{21}-287 \beta^{22}\cr n=7 & -\beta^{7}+2 \beta^{8}-\beta^{9}-14 \beta^{14}+44 \beta^{15}\cr&-46 \beta^{16}+16 \beta^{17}-301 \beta^{21}+1276 \beta^{22}\cr&-2025 \beta^{23}+1426 \beta^{24}-376 \beta^{25}\cr n=8 & -\beta^{8}+2 \beta^{9}-\beta^{10}-16 \beta^{16}+50 \beta^{17}\cr&-52 \beta^{18}+18 \beta^{19}-392 \beta^{24}+1650 \beta^{25}\cr&-2601 \beta^{26}+1820 \beta^{27}-477 \beta^{28}\cr n=9 & -\beta^{9}+2 \beta^{10}-\beta^{11}-18 \beta^{18}+56 \beta^{19}\cr&-58 \beta^{20}+20 \beta^{21}-495 \beta^{27}+2072 \beta^{28}\cr&-3249 \beta^{29}+2262 \beta^{30}-590 \beta^{31}\cr n=10 & -\beta^{10}+2 \beta^{11}-\beta^{12}-20 \beta^{20}+62 \beta^{21}\cr&-64 \beta^{22}+22 \beta^{23}-610 \beta^{30}+2542 \beta^{31}\cr&-3969 \beta^{32}+2752 \beta^{33}-715 \beta^{34}\cr}$$

All coefficients are integers, and it starts with $-\beta^n + 2 \beta^{n+1} - \beta^{n+2} - 2 n \beta^{2n} + (6n+2) \beta^{2n+1} - (6n+4) \beta^{2n+2} + (2n+2) \beta^{2n+3} - (6n+1)n \beta^{3n} \ldots$

Answer 2

Various commenters have pointed out that one can easily get asymptotic estimates for the eigenvalues as $n \rightarrow \infty$ (specifically, with $\beta, \Delta > 0$ they converge to $-2 (\beta + \Delta) \cos(\pi k/2 n)^2$ for $k = 1$ to $n$). OTOH, you specify in the comments that you require the eigenvalues in exact form. This is not possible. Even for $n = 5$, $\beta = 2$, and $\Delta = 1$, the characteristic polynomial of the matrix above is $$x^5 + 14 x^4 + 70 x^3 + 150 x^2 + 129 x + 32$$ whose splitting field has Galois group $S_5$, and hence there cannot be any exact formula (in radicals) for the roots by a theorem of Abel. This suggests that it is highly unlikely that there exists any "exact" formula for the roots, although "exact formula" is a somewhat nebulous notion, since of course one can define a function $f_{n,k}(\beta,\Delta)$ to be the $k$th smallest root of the corresponding characteristic polynomial.

Answer 3

There are such formulas, but beware: they are rather intimidating.

Look at this paper:

"EXPLICIT EIGENVALUES AND INVERSES OF TRIDIAGONAL TOEPLITZ MATRICES WITH FOUR PERTURBED CORNERS" by WEN-CHYUAN YUEH and SUI SUN CHENG2, ANZIAM J. 49(2008), 361–387.