Demystifier said:
Even if that was true (which was not), you missed that I mentioned also a much more general argument (actually a theorem) that the energy is conserved whenever the Hamiltonian does not have an explicit dependence on time, which in our case is true because the omegas in the Hamiltonian are time independent.
I'm still not understanding all these arguments. I guess one should move this discussion to the mechanics forum, but isn't LL discussing the case of small oscillations around a minimum of a potential. Then the eigenfrequencies are of course all real.
That's easy to see. Without loss of generality we can shift the coordinates such that the minimum of the potential is at ##\vec{x}=0##. Then we can approximate the potential in a neibourhood of the origin by
$$V(\vec{x})=\frac{1}{2} A_{jk} x_j x_k,$$
where ##A_{jk}## is a symmetric positive definite matrix, i.e., there's a orthogonal transformation, ##\vec{x}'=\hat{O} \vec{x}## such that ##A_{jk}' x_j' x_k'=m \sum_{j} \omega_j^2 x_j^{\prime 2}## with all ##\omega_j^2>0##, and then the solutions for the normal coordinates ##x_k'## are just
$$x_k'(t)=A_k \exp(\mathrm{i} \omega_k t) + B_k \exp(-\mathrm{i} \omega_k t).$$
So for this case LL seems right to me.
You can of course solve everything in the original coordinates. The EoMs read
$$m \ddot{x}_j=-\partial_{j} V=A_{jk} x_k \; \Rightarrow \; m\ddot{\vec{x}}=-\hat{A} \vec{x}.$$
This you can solve by the ansatz
$$\vec{x}(t)=\vec{a} \exp(\mathrm{i} \omega t).$$
This gives
$$-m \omega^2 \vec{a} = -\hat{A} \vec{a},$$
which means that ##\omega## must be chosen such that ##m \omega^2## is an eigenvalue of ##\hat{A}## and ##\vec{a}## the corresponding eigenvector.
Since ##\hat{A}=\hat{A}^{\text{T}}## is positive definite there are ##3## positive eigenvalues (which may all be different from each other or not) you have ##\omega \in \mathbb{R}##, and one can choose the eigenvectors ##\vec{a}_k## real and orthonormal. Then the general solution reads
$$\vec{x}(t) = \sum_{j=1}^3 \vec{a}_k [A_k \exp(\mathrm{i} \omega_k t)+B_k \exp(-\mathrm{i} \omega k t)].$$
This can be written in real form as
$$\vec{x}(t)=\sum_{j=1}^3 \vec{a}_k [C_k \cos(\omega_k t) + D_k \sin(\omega_k) t],$$
where ##C_k,D_k \in \mathbb{R}##.
If ##\hat{A}## is not positive definite, the approximation doesn't make sense, because then there are unbound solutions as discussed above.