0
$\begingroup$

Introduction

I am trying to obtain the eigenvectors of a unitary matrix $M(k)$ which depends on a parameter k.

This matrix $M(k)$ has dimension 6, and while for general matrices of dimension 6 it's not possible to write down algebraic expressions for their eigenvalues, because the characteristic polynomial will also be of order 6, for the case of $M(k)$ it's characteristic polynomial is such that it is possible to write its eigenvalues in algebraic form.

When I use the Eigenvectors[] function in Mathematica, it gives me eigenvectors which

  1. Cannot be written algebraically, i.e. it involves expressions with roots and # when it should be possible to avoid this since the eigenvalues are algebraic.
  2. Discontinuous when I don't think they should be for this matrix.

So basically I don't trust how Eigenvectors[] is working for this matrix. I would like to use another way of calculating the eigenvectors.

Question

I have tried computing the eigenvectors $v$ of a matrix $M(k)$ of dimension 6 in Mathematica by using Solve[] on this equation

$(M(k) - aI)v = 0$

where a is an eigenvalue of $M(k)$.

Why does Mathematica only give me the trivial solution $v=0$?

I checked that the determinant of $(M(k) - aI)$ is zero, and so if I denote an eigenvector $v$ as $v = (v_1,v_2,v_3,v_4,v_5,v_6)$ then there should be a solution where $v_2, v_3, v_4, v_5,$ and $v_6$ are written solely in terms of $v_1$, but this solution doesn't appear.

I tried doing the exact same procedure as above for a different matrix of dimension 2 and IT DID GIVE ME the non-trivial solutions, so I'm not sure why it can't do it for $M(k)$.

This is the matrix

$ \frac{1}{\sqrt{2}}\begin{bmatrix} 0 & 0 & 1 & i & 0 & 0\\ 0 & 0 & 0 & 0 & ie^{-ik} & e^{-ik} \\ 0 & 0 & 0 & 0 & e^{i\frac{2\pi}{3}} & ie^{i\frac{2\pi}{3}} \\ ie^{i\frac{2\pi}{3}} & e^{i\frac{2\pi}{3}} & 0 & 0 & 0 & 0\\ e^{ik}e^{-i\frac{2\pi}{3}}& ie^{ik}e^{-i\frac{2\pi}{3}} & 0 & 0 & 0 & 0 \\ 0 & 0 & ie^{-i\frac{2\pi}{3}} & e^{-i\frac{2\pi}{3}} & 0 & 0 \\ \end{bmatrix}$

{{0,0,1/Sqrt[2],I/Sqrt[2],0,0},
 {0,0,0,0,(I E^(-I k))/Sqrt[2],E^(-I k)/Sqrt[2]},
 {0,0,0,0,E^((2 I π)/3)/Sqrt[2],(I E^((2 I π)/3))/Sqrt[2]},
 {(I E^((2 I π)/3))/Sqrt[2],E^((2 I π)/3)/Sqrt[2],0,0,0,0},
 {E^(I k-(2 I π)/3)/Sqrt[2],(I E^(I k-(2 I π)/3))/Sqrt[2],0,0,0,0},
 {0,0,(I E^(-((2 I π)/3)))/Sqrt[2],E^(-((2 I π)/3))/Sqrt[2],0,0}} 

I need to find the eigenvectors of this matrix

$\endgroup$
8
  • 2
    $\begingroup$ NullSpace[ M - a I ] ? $\endgroup$
    – LouisB
    Apr 26, 2019 at 12:01
  • $\begingroup$ This hasn't worked but thanks. It returns an empty bracket $\endgroup$
    – Ben Pepper
    Apr 26, 2019 at 12:11
  • 2
    $\begingroup$ Root objects are algebraic numbers. $\endgroup$
    – John Doty
    Apr 26, 2019 at 12:41
  • 3
    $\begingroup$ Why do you want a solution in radical terms? Root objects are much easier to work with. One thing you might find useful about them is their well-defined ordering. $\endgroup$
    – John Doty
    Apr 26, 2019 at 14:20
  • 2
    $\begingroup$ Eigenvectors[M[k]] // ToRadicals gives explicit expressions without Root objects. $\endgroup$
    – Roman
    Jan 13, 2021 at 12:17

1 Answer 1

3
$\begingroup$

The 6×6 matrix $m$ looks like three 2×2 matrices, so re-ordering the matrix may be helpful. We can find the first eigenvector, in terms of the Root expressions, this way

μ = Eigenvalues[m];

rowOrder = {4, 5, 1, 6, 2, 3};
xfrm = IdentityMatrix[6][[rowOrder]];
s = xfrm.m;
v = First@NullSpace[s - μ[[1]] xfrm];
m.v - μ[[1]] v // Simplify

(*  {0, 0, 0, 0, 0, 0}  *)

So, $v$ is an eigenvector for the first eigenvalue, $\mu_1$. The transformation that re-ordered the rows preserved the null space. If someone can find a transformation that preserves the eigenvalues and re-orders the matrix to have a banded structure, Mathematica may (or may not) be able to find the eigenvalues that use radicals instead of Root.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.