Hi all,(adsbygoogle = window.adsbygoogle || []).push({});

I'm having a few problems with crystal dynamics of a simple monatomic chain.

Taking the dispersion relation:

[itex]\omega^2 = \frac{4k}{m}\left(\sin^2 \left( \frac{\kappa a}{2}\right)\right)[/itex]

Where k=spring constant, m=mass, [itex]\kappa[/itex]=wavevector, a= lattice constant and [itex]\omega[/itex]= frequency.

Now assume we have a periodic chain made up of four ions and applying the boundary conditions we get:

[itex]\kappa = \frac{2 \pi n}{a L}[/itex]

where n is just an integer and L is the length of the chain (here L=4).

Substituting this into the dispersion relation and assuming mass and springs are unity we can calculate the four frequencies for the four admissible modes as [itex]\omega^2=[/itex]0,2,2 and 4.

Now solving numerically, we construct the dynamical matrix using the same values for masses and springs as:

(2, -1, 0, -1; -1, 2, -1, 0; 0, -1, 2, -1; -1, 0, -1, 2)

Diagonalizing this gives the eigenvalues ([itex]\omega^2[/itex]) identical to the ones above calculated analytically and the following eigenvectors (extensions from equilibrium):

(0.5, 0.5, 0.5, 0.5)

(0.7071, 0, -0.7071, 0)

(0, 0.7071, 0, -0.7071)

(-0.5, 0.5, -0.5, 0.5)

I've tried analytically to obtain these eigenvectors but I just can't see how. I recall the equation for extension from equilibrium as:

[itex] u_s = U exp[\frac{i 2 \pi n s}{L}] [/itex]

Again L=4, U is the amplitude and for normalization it is [itex]L^{-1/2}[/itex], [itex]u_s[/itex] is the extension of site s. If I just substitute the values site numbers 1 to 4 in along with mode numbers n= 1 to 4. I do not reproduce the eigenvectors found from diagonalization.

Thanks for any help in advance.

Sebastian

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Monatomic Linear Chain - Comparison of Numerical and Analytical Results

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

**Physics Forums | Science Articles, Homework Help, Discussion**