Wigner Crystals

I’d like to explain a conjecture about Wigner crystals, which we came up with in a discussion on Google+. It’s a purely mathematical conjecture that’s pretty simple to state, motivated by the picture above. But let me start at the beginning.

Electrons repel each other, so they don’t usually form crystals. But if you trap a bunch of electrons in a small space, and cool them down a lot, they will try to get as far away from each other as possible—and they can do this by forming a crystal!

This is sometimes called an electron crystal. It’s also called a Wigner crystal, because the great physicist Eugene Wigner predicted in 1934 that this would happen.

Only since the late 1980s have we been able to make electron crystals in the lab. Such a crystal can only form if the electron density is low enough. The reasons is that even at absolute zero, a gas of electrons has kinetic energy. At absolute zero the gas will minimize its energy. But it can’t do this by having all the electrons in a state with zero momentum, since you can’t put two electrons in the same state, thanks to the Pauli exclusion principle. So, higher momentum states need to be occupied, and this means there’s kinetic energy. And it has more if its density is high: if there’s less room in position space, the electrons are forced to occupy more room in momentum space.

When the density is high, this prevents the formation of a crystal: instead, we have lots of electrons whose wavefunctions are ‘sitting almost on top of each other’ in position space, but with different momenta. They’ll have lots of kinetic energy, so minimizing kinetic energy becomes more important than minimizing potential energy.

When the density is low, this effect becomes unimportant, and the electrons mainly try to minimize potential energy. So, they form a crystal with each electron avoiding the rest. It turns out they form a body-centered cubic: a crystal lattice formed of cubes, with an extra electron in the middle of each cube.

To know whether a uniform electron gas at zero temperature forms a crystal or not, you need to work out its so-called Wigner-Seitz radius. This is the average inter-particle spacing measured in units of the Bohr radius. The Bohr radius is the unit of length you can cook up from the electron mass, the electron charge and Planck’s constant:

\displaystyle{ a_0=\frac{\hbar^2}{m_e e^2} }

It’s mainly famous as the average distance between the electron and a proton in a hydrogen atom in its lowest energy state.

Simulations show that a 3-dimensional uniform electron gas crystallizes when the Wigner–Seitz radius is at least 106. The picture, however, shows an electron crystal in 2 dimensions, formed by electrons trapped on a thin film shaped like a disk. In 2 dimensions, Wigner crystals form when the Wigner–Seitz radius is at least 31. In the picture, the density is so low that we can visualize the electrons as points with well-defined positions.

So, the picture simply shows a bunch of points x_i trying to minimize the potential energy, which is proportional to

\displaystyle{ \sum_{i \ne j} \frac{1}{\|x_i - x_j\|} }

The lines between the dots are just to help you see what’s going on. They’re showing the Delauney triangulation, where we draw a graph that divides the plane into regions closer to one electron than all the rest, and then take the dual of that graph.

Thanks to energy minimization, this triangulation wants to be a lattice of equilateral triangles. But since such a triangular lattice doesn’t fit neatly into a disk, we also see some ‘defects’:

Most electrons have 6 neighbors. But there are also some red defects, which are electrons with 5 neighbors, and blue defects, which are electrons with 7 neighbors.

Note that there are 6 clusters of defects. In each cluster there is one more red defect than blue defect. I think this is not a coincidence.

Conjecture. When we choose a sufficiently large number of points x_i on a disk in such a way that

\displaystyle{ \sum_{i \ne j} \frac{1}{\|x_i - x_j\|} }

is minimized, and draw the Delauney triangulation, there will be 6 more vertices with 5 neighbors than vertices with 7 neighbors.

Here’s a bit of evidence for this, which is not at all conclusive. Take a sphere and triangulate it in such a way that each vertex has 5, 6 or 7 neighbors. Then here’s a cool fact: there must be 12 more vertices with 5 neighbors than vertices with 7 neighbors.

Puzzle. Prove this fact.

If we think of the picture above as the top half of a triangulated sphere, then each vertex in this triangulated sphere has 5, 6 or 7 neighbors. So, there must be 12 more vertices on the sphere with 5 neighbors than with 7 neighbors. So, it makes some sense that the top half of the sphere will contain 6 more vertices with 5 neighbors than with 7 neighbors. But this is not a proof.

I have a feeling this energy minimization problem has been studied with various numbers of points. So, there either be a lot of evidence for my conjecture, or some counterexamples that will force me to refine it. The picture shows what happens with 600 points on the disk. Maybe something dramatically different happens with 599! Maybe someone has even proved theorems about this. I just haven’t had time to look for such work.

The picture here was drawn by Arunas.rv and placed on Wikicommons on a Creative Commons Attribution-Share Alike 3.0 Unported license.

22 Responses to Wigner Crystals

  1. HaroldW says:

    Why is the potential energy \frac{1}{\|x_i - x_j\|^2}  and not \frac{1}{\|x_i - x_j\|^1} ?

    • jessemckeown says:

      I was going to ask why, in the 3-d case, the cloud doesn’t condense into a thin shell on the boundaries of the box? Because that’s what, in school, we were told happens to extra charge on a round conductor…

    • John Baez says:

      Sorry, I was rushing and made a mistake. I’ll fix it.

    • John Baez says:

      The main reason all the charge doesn’t go to the boundary of the disk is that we’re working in two dimensions yet using the potential \phi = 1/r that’s the solution of the Poisson equation

      \nabla^2 \phi = \frac{1}{4 \pi} \delta

      in three dimensions. If we minimized potential energy for the potential that solved Poisson’s equation in two dimensions, the charge particles would indeed cluster on the boundary of the disk…

      … at least in the limit where the number of particles approaches infinity.

      In the 3-dimensional case, for finite numbers of particles confined to a ball minimizing the potential energy

      \displaystyle{ \sum_{i \ne j} \frac{1}{\|x_i - x_j\|} ,}

      the particles don’t always lie on the surface of the ball! In the 2-dimensional case, with the potential energy appropriate to 2 dimensions

      \displaystyle{ - \sum_{i \ne j} \ln \|x_i - x_j\| , }

      I think the case of 7 particles would be a good one to try, to see this effect.

  2. witzar01 says:

    Puzzle solution: Let v, e, f be respectively the number of vertices, edges and faces of the triangulation, and let d be the difference between number of vertices with 5 edges and vertices with 7 edges. By counting edges using faces we get e = 3f/2, and by counting vertices using faces we get v = (3f + d)/6. Plugging it all into Euler characteristics (v – e + f = 2) we get the equation (3f + d)/6 – 3f/2 + f = 2, where all f’s magically cancel out and we are left with d = 12.

  3. jessemckeown says:

    Sub-conjectures, then: electrons on the boundary always have four neighbours each (so pasting in a second copy of the whole makes a sphere, and introduces no 5- or 7-neighbour electrons); and any configurations with 4 or 8 neighbours can relax further.

    • John Baez says:


      Theorem. Given a triangulated disk, if each vertex in the interior has 5, 6 or 7 neighbors, and each vertex on the boundary has 4, then there are six more vertices with 5 neighbors than with 7.

      Proof. Take two copies of the triangulated disk and attach them in the obvious way, attaching each point on the boundary of one to the corresponding point on the boundary of the other, to form a triangulated sphere. Now the vertices that were on the boundary all have 6 neighbors. Thus, every vertex in the triangulated sphere has 5, 6, or 7 neighbors. However, because the Euler characteristic of a sphere is 2, there must be exactly twelve more vertices with 5 neighbors than with 7. (See the proof earlier in these comments.) Thus, in each copy of the triangulated disk, there must be six more vertices with 5 neighbors than with 7.

      • Tamfang says:

        Of course you mean “because the Euler characteristic of a sphere is 2,” and each 5-way vertex accounts for 1/6 of this.

      • John Baez says:

        Right—I fixed that.

        A reasonable way to think about it is to imagine all the triangles as equilateral (even if they’re not really). Then a vertex with 6 neighbors is flat, a vertex with 5 neighbors has an angle deficit of 60° a vertex with 7 neighbors has an angle excess of 60°, and so on. This, combined with the fact that any triangulation of a 2-sphere needs to have a total angle deficit of 720°, gets the job done. It’s a discretized version of the Gauss—Bonnet theorem.

  4. Tamfang says:

    Related question: when N>>120 balls are packed in S3 (positively curved 3space), presumably they form regions of face-centred cubic and/or body-centred cubic, but what are the defects like? Are they concentrated around 120 points, or along arcs, or what?

  5. A says:

    I’d imagine you want to look into Burgers vectors if you haven’t already, since there are conservation laws related to them that can be applied to defects in crystals. Decent summary can be found in this book:

    Click to access UploadFile_6472.pdf

  6. John Baez says:

    Here is a paper on minimizing

    \displaystyle{ \sum_{i \ne j} \|x_i - x_j\|^{-1} }

    for a collection of points on the disk:

    • A. Worley, Minimal energy configurations for charged particles on a thin conducting disc: the perimeter particles.

    One interesting thing I hadn’t guessed is that the particles often lie on concentric shells.

    Abstract. The lowest energy configurations for N equal charged particles confined to a thin conducting disc have been investigated in detail up to N = 160 and in outline for further values up to N = 500. For all values of N up to 160 the particle configurations can be described in terms of concentric shells. The number of perimeter particles p appears to be simply related to N and to the mean radius of the outermost internal shell. Justification for these relations is obtained from a simple model based on the well-known distribution of continuous charge on a conducting disc.

    In the N → ∞ limit the density of particles in the unit disk is proportional to

    \displaystyle{  \frac{1}{\sqrt{ 1- r^2}}  }

    Other interesting things:

    In the range N ≤ 80, the configuration for N + 1 is, with one exception, the same as that for N with the addition of one particle to an existing shell or the creation of a new shell with a particle placed at the centre. The single exception is the pair N = 55 (with shell structure 5-13-37) and N = 56 (with shell structure 1-6-12-37). For the range 81 ≤ N ≤ 160 there are four exceptions to the usual pattern: N = 97–98, 117–118, 150–151 and 152–153.

    The first case where the particles don’t lie on concentric shells is N = 185.

    • HaroldW says:

      An interesting contrast is “Topological defect motifs in two-dimensional Coulomb clusters” by Arunas Radzvilavičius and Egidijus Anisimovas where the particles are confined by a parabolic energy well, rather than a square well. They write, “Distribution of defects in a parabolic confinement is by large determined by a conflict between the circular boundary and the bulk-like interior where a hexagonal lattice is preferred. Small clusters (N ≤ 70) have no bulk and exhibit a shell structure defined by the circular symmetry of the confinement.”

      Their Figure 5 (N=520) is similar to the figure at the top of your post.

      • John Baez says:

        There are very nice results in this paper – for example, a description of some common defect patterns. Such patterns should appear in a large 2d Wigner crystal for a large variety of potentials (i.e. parabolic potential, square well potential on a disc, etc.).

You can use Markdown or HTML in your comments. You can also use LaTeX, like this: $latex E = m c^2 $. The word 'latex' comes right after the first dollar sign, with a space after it.

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

This site uses Akismet to reduce spam. Learn how your comment data is processed.