1
$\begingroup$

Let $A$ be an invertible, symmetric and tridiagonal matrix of size $n \times n$. Assume that $A_{i,i}=a \neq 0$ for $i=1\dotsc n$ and all the elements in the sub- and super-diagonal of $A$ are $b \neq 0$. I would like to simplify the following Kronecker product: $e^{-A} \otimes e^{A}$.

I know that, given the Kronecker sum property of matrix exponential ($e^{A\oplus B}= e^{A}\otimes e^{B}$), the following holds:

\begin{equation} e^{-A} \otimes e^{A} = e^{-A \otimes I_n +I_n \otimes A}. \end{equation}

Since $A \otimes I_n$ and $I_n \otimes A$ commutes, using Zassenhaus formula,

\begin{equation} e^{-A} \otimes e^{A} = e^{-A \otimes I_n} e^{I_n \otimes A}=(e^{-A}\otimes I_n)(I_n \otimes e^{A}) \end{equation}

Given the above mentioned properties of matrix $A$, I was wondering whether it would be possible to further simplify this expression.

$\endgroup$
7
  • $\begingroup$ What does $e^A$ mean? Is there a missing hypothesis about where the entries of $A$ come from? $\endgroup$ Commented Oct 18, 2022 at 20:01
  • $\begingroup$ It is the matrix exponential as defined in the provided link. I am not sure what you mean be the second part of your question. The enteries can be any non zero real numbers. $\endgroup$
    – Mirar
    Commented Oct 18, 2022 at 20:10
  • 1
    $\begingroup$ The formula you obtained can also be proved more quickly using $(A\otimes B)(C\otimes D) = (AC)\otimes (BD)$. $\endgroup$ Commented Oct 18, 2022 at 21:51
  • $\begingroup$ Agree on the first comment. I just showed a different variation of it. $\endgroup$
    – Mirar
    Commented Oct 18, 2022 at 21:55
  • 3
    $\begingroup$ There is a sign error in your exponents, right? $\endgroup$ Commented Oct 19, 2022 at 1:27

1 Answer 1

6
$\begingroup$

One can diagonalize your $A = VDV^{-1}$ explicitly; the closed formulas are here for instance.

Once you have those matrices, you can write the orthogonal eigendecomposition $$ \exp(-A) \otimes \exp(A) = (V\otimes V) \,\, (\exp(-D)\otimes \exp(D)) \, \,(V\otimes V)^{-1}. $$ The diagonal matrix $\exp(-D)\otimes \exp(D)$ has elements $\exp(-\lambda_i + \lambda_j)$, $i,j=1,\dots,n$; when you replace it with the formulas for the eigenvalues, the $a$'s simplify out and you get a difference of two cosines that you can further manipulate using the sum-to-product formulas. I don't think it gets any simpler than that.

Alternatively, you can write $A = aI + bZ$, where $Z$ is the matrix with ones on the super- and subdiagonal, and write $\exp(A) = e^a \exp(bZ)$ (since the two summands commute) and $$ \exp(-A) \otimes \exp(A) = \exp(-bZ) \otimes \exp(bZ). $$ This shows in a simpler way that $a$ simplifies out from your expression.

$\endgroup$
3
  • $\begingroup$ +1 Fair enough. Thank you for the answer. $\endgroup$
    – Mirar
    Commented Oct 19, 2022 at 6:46
  • $\begingroup$ @Federico: Is it possible to find an explicit form for $exp(Z)$? where Z is the sum of two nilpotent matrices one upper and one lower triangular matrices. $\endgroup$
    – Mirar
    Commented Oct 19, 2022 at 8:33
  • $\begingroup$ @Mirar Yes, with the techniques described in the first half of the answer; just set $a=0$ (which is not restrictive since the answer does not depend on $A$ anyway). $\endgroup$ Commented Oct 19, 2022 at 8:51

Not the answer you're looking for? Browse other questions tagged .