# Eigenvalue equation

1. Sep 18, 2007

### jonny81

Hello,

I want to prove that the function $\mathcal{A}$ in the 1D case satisfy

$\mathcal{A}=\frac{48}{m}\sum_{j=1}^\infty \frac{\sin^2(qj/2)}{j^5}=\frac{12}{m}\left[2\zeta(5)-\text{Li}_5(e^{iq})-\text{Li}_5(e^{-iq})\right],$

with $\text{Li}_n(z)$ the polylogarithm function, and the matrix $\mathcal{A}$ in the 2D case satisfy

$\mathcal{A}=\frac{3}{m}\sum_{j\neq 0}\left[1-\frac{5}{|\vec{r}_j^0|^2}\begin{pmatrix}(x_j^0)^2&x_jy_j\\x_jy_j&(y_j^0)^2\end{pmatrix}\right]\frac{\left(1-e^{i\vec{q}\vec{r}_j^0}\right)}{|\vec{r}_j^0|^5}$

Can somebody help me, please, to do this? In 2D the equilibrium positions $\vec{r}_i^0$ form a triangular lattice with basic lattice vectors $a_1 = (1, 0)$ and $a_2 = (1,\sqrt{3})/2$.

Starting point:

$m\vec{\ddot x}_i=\frac{3}{2}\sum_{i\neq j}^{N}\left[\frac{5(\vec{r}_i^0 - \vec{r}_j^0)^2\left(\vec{x}_j-\vec{x}_i\right)}{|\vec{r}_i^0 - \vec{r}_j^0|^7}+\frac{\left(\vec{x}_i-\vec{x}_j\right)}{|\vec{r}_i^0 - \vec{r}_j^0|^5}\right]$

with the ansatz

$\vec{x}_i(t)=\epsilon_\lambda(\vec{q})e^{i\left(\vec{q}\vec{r}_i^0-\omega_\lambda(\vec{q})t\right)}$
$\vec{\ddot x}_i(t)=-\omega^2\epsilon_\lambda(\vec{q})e^{i\left(\vec{q}\vec{r}_i^0-\omega_\lambda(\vec{q})t\right)}$

follows

$-\omega^2_\lambda(\vec{q})\epsilon_\lambda(\vec{q})=\mathcal{A}\epsilon_\lambda(\vec{q})$

Thanks!

Last edited: Sep 18, 2007
2. Sep 18, 2007

### mjsd

new to PF?
what have you tried? is this a condense matter problem?

3. Sep 18, 2007

### jonny81

Yes!

with the ansatz in the equation of motion, i obtain

$\mathcal{A}=\frac{3}{2m}\sum_{i\neq j}\left[1-\frac{5}{|\vec{r}_i^0-\vec{r}_j^0|^2}\left(\vec{r}_i^0-\vec{r}_j^0\right)^2 \right]\frac{\left(1-e^{i\vec{q}(\vec{r}_i^0-\vec{r}_j^0)}\right)}{|\vec{r}_i^0-\vec{r}_j^0|^5}$

What I have to do next? I do not make any differentiation about 1D and 2D.

I consider a dipolar crystal.

4. Sep 18, 2007

### mjsd

your question is asked in some kind of a strange way... I don't even know what you are asking...or what are you trying to prove? just sub in the ansatz and check? what are you having problems with? it is difficult for anyone to help without knowing what you are stuck at. Besides in this forum you are meant to show us what you have tried