19. Dynamics of a One-Dimensional Crystal
Notation! In this lecture, I use for the spring constant ( is a wave number) and for frequency ( is a root of unity).
A good classical model for a crystal is to represent the atoms by balls held in place by light springs, representing valence bonds, between nearest neighbors. The simplest such crystal that has some realistic features is a single chain of connected identical atoms. To make the math easy, we’ll connect the ends of the chain to make it a circle. This is called “imposing periodic boundary conditions”. It is common practice in condensed matter theory, and makes little difference to the physics for a large system.
We’ll take the rest positions of the atoms to be uniformly spaced, apart, with the first atom at the atom at the final atom at the origin.
Away from the lowest energy state, we denote the position of the atom , so, as in our earlier discussion of oscillating systems, is the displacement from equilibrium (which we take to be along the linewe are not considering transverse modes of vibration at this time).
The Lagrangian of this circular chain system is:
We’re going to call the spring constant , we’ll need for something else. We’ll also call the frequency .
Looking for eigenstates with frequency , we find the set of equations
Taking a solution , with the understanding that may be complex, and at the end is just the real part of the formal solution, we find the eigenvalue equation for a chain of four atoms (the biggest Mathtype can handle!)
Actually we’d have a much bigger matrix, with lots of zeroes, but hopefully the pattern is already clear: for each diagonal element and ’s in two diagonal-slanting lines flanking the main diagonal (corresponding to the links between nearest neighbors) and finally ’s in the two far corners, these coming from the spring joining to 1 to complete the circle.
Notice first that if there is an eigenvector , since the sum of the elements in one row is zero (the means transpose, that is, it’s really a column vector, but row vectors are a lot easier to fit into the text here).
This eigenvector is just uniform displacement of the whole system, which costs zero energy since the system isn’t anchored to a particular place on the ring. We’ll assume, though, that the system as a whole is at rest, meaning the center of mass is stationary, and the atoms have well-defined rest positions as in the picture at
The Circulant Matrix: Nature of its Eigenstates
The matrix we’ve constructed above has a very special property: each row is identical to the preceding row with the elements moved over one place, that is, it has the form
Such matrices are called circulants, and their properties are well known. In particular, we’ll show that the eigenvectors have the form where
Recall the roots of the equation are points equally spaced around the unit circle,
The standard mathematical notation is to label these points as shown in the figure, but notice that .
Comparison with Raising Operators
Actually these matrices are related to the raising and lowering operator for angular momentum (and simple harmonic oscillators) in quantum mechanics. For example, the matrices would be for a spin 3/2, with four eigenstates.
The quantum mechanical raising and lowering matrices look like
They move the spin component up (and down) by one notch, except that on applying the raising operator, the top state is annihilated, similarly the lowering operator on the bottom state.
Our circular generalizations have one extra element:
This makes the matrices circulants, and gives them a “recycling” property: the top element isn’t thrown away, it just goes to the bottom of the pile.
(And bear in mind that the standard notation for a vector has the lowest index (0 or 1) for the top element, so when we bend the ladder into a circle, the “raising” operator actually moves to the next lower number, in other words, it’s a shift to the left.)
We’ll take this shift operator as our basic matrix:
It should be evident from this that the circulant matrix having top row is just the matrix
This generalizes trivially to matrices.
Finding the Eigenvectors
Now let’s look at the eigenvectors, we’ll start with those of Let’s call the eigenvalue .
Then for an eigenstate of the shift operator, the shifted vector must be just a multiple of the original vector:
Reading off the element by element equivalence of the two vectors,
The first three equalities tell us the eigenvector has the form , the last tells us that .
From our earlier discussion of circulant matrices, writing the smallest phase nontrivial root of unity as the roots of the equation are just this basic root raised to different powers: the roots are .
This establishes that the eigenvectors of have the form
where with corresponding eigenvalue the basic root raised to the power,
Try it out for : the eigenvalues are given by The corresponding eigenvectors are found to be For the case, there are different, linearly independent, vectors of this form, so this is a complete set of eigenvectors of .
They are also, of course, eigenvectors of all powers of , and therefore of all the circulant matrices! This means that all circulant matrices commute.
Eigenvectors of the Linear Chain
Let’s get back to our chain, with eigenfunction equation of motion the dimensional equivalent of
. We see the matrix is a circulant, so we know the eigenvectors are of the form , which we’ll now write
What does this mean for our chain system? Remember that the element of the eigenvector represents the displacement of the atom of the chain from its equilibrium position, that would be proportional to
The steady phase progression on going around the chain makes clear that this is essentially a (longitudinal) wave. (The actual particle displacement is the real part of the element, but there could be an overall complex factor fixing the phase.)
Allowed Wavenumbers from Boundary Conditions
The usual way of representing a wave on a line in physics is to have displacement proportional to , and is called the wavenumber. For our discretized system, the displacement parameter for the atom, at position , would therefore be proportional to
But we know this is an eigenvector of a circulant, so we must have and the allowed values of are
with an integer.
The circulant structure of the matrix has determined the eigenvectors, but not the eigenvalues .
Finding the Eigenvalues
The eigenvalues are found by operating on the eigenvector we just found with the matrix, meaning the dimensional generalization of
Applying the matrix to the column vector , and cancelling out the common factor, we have
(Of course, this same result comes from every row.)
The complete set of eigenvalues is given by inserting in the above expression
so is displacement of the system as a whole, as is .
Wavenumber values beyond repeats the eigenstates we already have, since
The eigenvalue equation is
To see the dynamics of this eigenstate , we need to multiply by the time dependence , then finally take the real part of the solution:
Notice that in the continuum limit, meaning large and small the atom displacement as a function of position has the form , in other words we’re looking at a sinusoidal wave disturbance with wavenumber here.
Now, is also a solution, but that is the same as so one must be careful not to overcount. The two frequencies correspond to waves going in opposite directions.
The Discrete Fourier Transform
It’s worth looking over this one more time from a slightly different perspective. In finding the energy of an oscillating continuous string, a standard approach is to analyze the motion of the string in terms of an infinite Fourier series of shorter and shorter wavelength oscillations, find the energy in each of these modes, and add to find the total energy. We’ll apply the same approach herebut with a difference. Since the waves only have meaning in our chain at a discrete set of uniformly spaced points, the set of waves needed to fully account for all possible motions is finite. In fact, it’s the same as the number of points. As we’ve discussed above, a wave with a higher wavenumber gives an identical set of displacements of the atoms as some lower one. So a complete Fourier analysis of the displacements at these equally spaces points only needs linear combinations of waves. This is the Discrete Fourier Transform (DFT).
Writing the complex (amplitude and phase) coefficient of the frequency eigenstate , the position of the atom in a superposition of such waves (with the standard normalization convention)
Given the positions of the atoms, the amplitude coefficients can be found using the inverse mapping:
gives establishing that we have the correct form for the inverse transformation.
The instantaneous configuration of the system is completely defined by the set and equally by the set All possible particle displacements at the equally spaced sites can be mapped into amplitudes of the distinct waves (eigenvectors).
(This DFT mapping is widely used in the time domain in signal processing: the signal amplitude is sampled, say every millisecond, then the data can be DFT’d to give the wave components down to a minimum frequency around one millisecond. A good quality voice signal would need a shorter time interval, maybe 0.2 milliseconds.)
Now, from ,
and again using
Back to our chain: for a physical configuration of the atoms, all the must be real, so from we see that (This reduces the number of apparent degrees of freedom in the representation to the correct . is real, , etc., and if there is a middle , it must be real.)
The kinetic energy of the chain particles,
We can find potential energy similarly:
and using the same routine as before,
Finally, Putting all this together, the Lagrangian can be written in terms of the transformed variables:
The equation of motion is then
This is of course the same result we found earlier, but it is perhaps worth seeing how it comes from the (mathematically equivalent) DFT analysis.
A Note on the Physics of These Waves
For wavelength long compared to the interparticle spacing, , these are like sound waves (and indeed they are what are called acoustic phonons in a crystal). As the wavelength shortens, the wave eigenstates are moving more slowly, remember the group velocity of a wavepacket goes as This is because there is some Bragg reflection of the waves by the lattice. At we have a standing wave. This is the highest energy mode, with even numbered sites all in sync with each other, and the odd numbered sites all half a cycle behind, so the restoring force experienced by an atom as a function of displacement is the maximum possible.