Question about Angle-Preserving Operators
Solution 1:
If necessary, rescale each $x_i$ so they all have unit length and relabel the indices so that $|\lambda_1|\neq |\lambda_2|$.
Now, consider the vectors $v_1 = x_1 + x_2$ and $v_2 = x_1-x_2$.
First, I claim that $v_1$ and $v_2$ are orthogonal, for \begin{align*} \langle v_1, v_2\rangle &= \langle x_1+x_2,\; x_1-x_2\rangle \\\ &= \langle x_1, x_1\rangle + \langle x_1, x_2\rangle - \langle x_2,x_1\rangle - \langle x_2, x_2\rangle \\\ &= |x_1|^2-|x_2|^2 \\\ &=0\end{align*} where the last equality follows since both $x_1$ and $x_2$ are unit length.
Next, I claim that $Tv_1$ and $Tv_2$ are not orthogonal. For, \begin{align*} \langle Tv_1, Tv_2\rangle &= \langle \lambda_1 x_1 + \lambda_2 x_2,\; \lambda_1 x_1 - \lambda_2 x_2\rangle \\\ &= \lambda_1^2 \langle x_1, x_1\rangle +\lambda_1 \lambda_2 \langle x_1, x_2\rangle - \lambda_2\lambda_1 \langle x_2,x_1\rangle - \lambda_2^2 \langle x_2, x_2\rangle \\\ &= \lambda_1^2 |x_1|^2 - \lambda_2^2 |x_2|^2 \\\ &= \lambda_1^2 - \lambda_2^2\end{align*}
and this last line is $0$ iff $|\lambda_1| = |\lambda_2|$.
Solution 2:
This is not true.
Let $T=\begin{bmatrix} 1 & -2 \\ 0 & -1 \end{bmatrix}$. Then $x_1=\begin{bmatrix} 1 \\ 0 \end{bmatrix}$, $x_2=\begin{bmatrix} 1 \\ 1 \end{bmatrix}$ form a basis, and $T x_1 = x_1$, $T x_2 = - x_2$, hence the eigenvector requirement is satisfied.
However, take $x=\begin{bmatrix} 1 \\ 0 \end{bmatrix}$, $y=\begin{bmatrix} 0 \\ 1 \end{bmatrix}$, then $\langle x , y \rangle = 0$, but $\langle T x , T y \rangle = -2$. Since $\arccos$ is bijective on the domain $[-1,1]$, it follows that $\arccos \frac{\langle x , y \rangle}{\|x\| \|y\|} \neq \arccos \frac{\langle T x , T y \rangle}{\|T x\| \|T y\|}$.
Here is a slightly more satisfactory answer:
$T$ is angle preserving iff $T^T T = \sigma^2 I$, for some $\sigma >0$.
($\Rightarrow$) Suppose $T^T T = \sigma^2 I$. Then $T' = \frac{1}{\sigma} T $ is orthogonal. It follows that $\langle T'x, T'y \rangle = \langle x, y \rangle$, and $\|T'v \| = \|v\|$, hence $\frac{\langle T'x, T'y \rangle}{\|T'x \| \|T'y \|} = \frac{\langle Tx, Ty \rangle}{\|Tx \| \|Ty \|} = \frac{\langle x, y \rangle}{\|x \| \|y \|}$.
($\Leftarrow$) Let $e_k$ be the standard basis on $\mathbb{R}^n$, and suppose $T$ is angle preserving. It follows immediately that $\frac{\langle T e_i, T e_j \rangle}{\|Te_i \| \|Te_j \|} = \delta_{i,j}$, hence $\frac{T e_i}{\|Te_i \|}$ form an orthonormal basis. Let $ Q = \begin{bmatrix} \frac{T e_1}{\|Te_1 \|} & \cdots & \frac{T e_n}{\|Te_n \|} \end{bmatrix}$, and let $\Lambda = \mathbb{diag}(\|Te_1 \|, \cdots, \|Te_n \|)$. It should be clear that $Q$ is orthogonal, and since $Tx = \sum x_i \|Te_i \| \frac{T e_i}{\|Te_i \|} = Q \Lambda x$, we have $T = Q \Lambda$. Now, adapting @Jason DeVito's trick, we notice that $\langle e_i+e_j, e_i-e_j \rangle = 0$, and so since $T$ is angle preserving, we have $\langle T(e_i+e_j), T(e_i-e_j) \rangle = \langle \Lambda(e_i+e_j), \Lambda(e_i-e_j) \rangle = \langle \Lambda e_i, \Lambda e_i \rangle - \langle \Lambda e_j, \Lambda e_j \rangle= 0$. It follows that $ \|T e_i \|^2 = \|T e_j\|^2$, and so we can write $\Lambda = \sigma M$, with $M$ a diagonal matrix filled with $\pm1$ (note $M^2=I$) and $\sigma = \|T e_1 \|$. Hence $T = \sigma QM$, from which is follows that $T^T T = \sigma^2 I$.