# Electrons in One Dimension: the Peierls Transition

*Michael Fowler, UVa *

### Introduction: Looking for Superconductors and Finding Insulators

In 1964, Little suggested (Phys. Rev 134, A1416) that it might be possible to synthesize a room temperature superconductor using organic materials in which the electrons traveled along certain kinds of chains, effectively confined to one dimension.

The first satisfactory theory of “ordinary” superconductivity, that of Bardeen, Cooper and Schrieffer (BCS) had appeared a few years earlier, in 1957. The key point was that electrons became bound together in opposite spin pairs, and at sufficiently low temperatures these bound pairs, being boson like, formed a coherent condensate$\u2014$all the pairs had the same total momentum, so all traveled together, a supercurrent. The locking of the electrons into this condensate effectively eliminated the usual single-electron scattering by impurities that degrades ordinary currents in conductors.

But what could bind the electrostatically repelling electrons? The answer turned out to be lattice distortions, as first suggested by Fröhlich in 1950. An electron traveling through the crystal attracts the positive ions, the consequent excess of local positive charge attracts another electron. The strength of this binding, and hence the temperature at which the superconducting transition takes place, depends on the rapidity of the lattice response. This was confirmed by the isotope effect: lattice response time obviously depends on the inertia of the lattice, the BCS theory predicted that for a superconducting element with different isotopic varieties, the ratio of the superconducting transition temperatures for pure isotopes ${T}_{2}/{T}_{1}$ was equal to $\sqrt{{M}_{1}/{M}_{2}}$, ${M}_{1},{M}_{2}$ being the ion masses, the lighter isotope having the higher transition temperature. This was indeed the case.

Little’s idea was that the build up of positive charge by a passing electron could be speeded up dramatically if instead of having to move ions, it need only rearrange other electrons. Unfortunately, there were no obvious three-dimensional candidate materials. However, if the conduction electrons moved along a one-dimensional chain, polarizable side chains might be attached, and rearrangement of the electronic charge distribution in these side chains would respond very rapidly to a passing conduction electron, building up a local positive charge. If this worked, order of magnitude arguments suggested possible enhancement of the transition temperature by a factor $\sqrt{M/m}$ over ordinary superconductors, $m$ being the electron mass.

In the 1970’s, various organic materials were synthesized
and tested, beginning with one called TTF-TCNQ, in which a set of polymer-like long
molecules donated electrons to another set, leaving one-dimensional conductors with
partially filled bands (see later), seemingly good candidates for
superconductivity. Unfortunately, on
cooling, these materials surprisingly became *insulators* rather than superconductors! This was the first example of a *Peierls transition*, a widespread
phenomenon in quasi one-dimensional systems.

The basic mechanism of the Peierls transition can be understood with a simple model. It is a nice example of applied second-order perturbation theory, including the degenerate case. We examine the model and the result below.

It should be added that in some newer materials the Peierls transition is (unexpectedly) suppressed under high pressure, and superconductivity has in fact been observed in organic salts, but so far only at transition temperatures around one Kelvin: Little’s dream is not yet realized.

### Second-Order Perturbation Theory: a Periodic Potential in One Dimension

To understand how a one-dimensional conductor might turn into an insulator at low temperatures, we must first become familiar with the simplest model of a one-dimensional conductor:

$H={H}^{0}+V=\frac{{p}^{2}}{2m}+V(x)$

with ${H}^{0}$ a gas of noninteracting electrons on a line, and $V$ periodic, that is, $V\left(x+a\right)=V\left(x\right),$

the potential from a line of ions spaced $a$ apart. We’ll take the system to have $N$ ions in a total length $L,$ so

$$L=Na$$ and to keep the math simple, we’ll require periodic boundary conditions.

The physics here is that *without*
the potential, the electron eigenstates are *plane
waves*. The effect of the lattice
potential is to partially reflect the waves, like a diffraction grating, generating
components at different wavelengths.
This effect becomes particularly important when the electron wavelength
matches twice the ion spacing. For that
case, the reflected and original waves have the same strength, the electron is
at a standstill. We’ll explore just how
this happens later.

The eigenstates of ${H}^{0}$ are then

${|k\rangle}^{(0)}=\frac{1}{\sqrt{L}}{e}^{ikx},\text{with}{e}^{ikL}=1,\text{so}k=\frac{2\pi n}{L},$

*$n$** *being an integer. The unperturbed energy eigenvalues,

${H}^{0}{|k\rangle}^{(0)}\text{\hspace{0.05em}}\text{\hspace{0.17em}}={E}^{0}{|k\rangle}^{(0)},\text{arejust}{E}^{0}=\frac{{\hslash}^{2}{k}^{2}}{2m}.$

(*Note*: this is to be understood as ${H}^{0}{|{k}_{n}\rangle}^{(0)}\text{\hspace{0.05em}}\text{\hspace{0.17em}}={E}_{n}^{0}{|{k}_{n}\rangle}^{(0)},\text{with}{E}_{n}^{0}=\frac{{\hslash}^{2}{k}_{n}{}^{2}}{2m}\text{and}{k}_{n}=\frac{2\pi n}{L}$.

We are following standard practice here. We shall also write $\sum _{k}f\left(k\right)}\text{meaning}{\displaystyle \sum _{n}f\left({k}_{n}\right)$. )

It’s worth plotting the (*E*,
*k*) curve:

Suppose we have ions with two electrons each to contribute
to this one-dimensional (supposed) conductor.
Assuming they move into these plane wave states, in the system ground
state they will fill up the lowest energy states up to a maximum $k\text{-}$ value denoted by $\pm {k}_{F}$ ($F$ stands for Fermi, this is the *Fermi* momentum.) Where is it?

We know there will be a total of $2N$ electrons.
We also know that the allowed values of $k,$ from the boundary conditions, are ${k}_{n}=2\pi n/L$,
with $n$ an integer. In other words, the allowed $k$ ’s are uniformly spaced $2\pi /L$ apart, meaning they have a *density *of $L/2\pi $ in $k\text{-}$ space, so the total number between $\pm {k}_{F}$ is $L{k}_{F}/\pi $. The $2N$ electrons will have $N$ of each spin, each $k\text{-}$ state can take two electrons (one of each
spin), so $L{k}_{F}/\pi =N=L/a$,
and

${k}_{F}=\pi /a.$

To do perturbation theory, we must find the matrix elements
of $V\left(x\right)$ between eigenstates of *H*^{ 0}:

${}^{(0)}\langle {k}^{\prime}|V{|k\rangle}^{(0)}=\frac{1}{L}{\displaystyle \int {e}^{i(k-k\text{'})x}V(x)dx.}$

This is just the Fourier component ${V}_{k-{k}^{\prime}}$ of $V\left(x\right).$

If $V\left(x\right)$ is periodic with period $a,$

${V}_{k}\ne 0\text{onlyif}k=nK,\text{}n\text{aninteger,}K=2\pi /a.$

In other words, if
a function is periodic with spatial period $a,$ the
only nonzero Fourier components are those having the same spatial period $a.$

Therefore

$V(x)={\displaystyle \sum _{n}{V}_{nK}{e}^{inKx}},\text{and}{V}_{-nK}={V}_{nK}^{*}\text{since}V(x)\text{isreal;}K=2\pi /a.$

The $n=0$ component of $V\left(x\right)$ is of no interest$\u2014$it is just a constant potential, and so can be taken to be zero. Note that this eliminates the trivial first order correction ${E}_{k}^{1}={\text{\hspace{0.17em}}}^{(0)}\langle k|V{|k\rangle}^{(0)}$ to the energy eigenvalues.

We shall consider *only*
the components $n=+1$ and $n=-1$ of $V\left(x\right),$ it turns out that the other components can be
treated in similar fashion. For $n=+1,\text{\hspace{0.17em}}\text{\hspace{0.17em}}n=-1,$ the potential only has nonzero matrix elements
between the plane wave state $k$ and $k+K,\text{\hspace{0.17em}}\text{\hspace{0.17em}}k-K$ respectively.

So, the second order correction to the energy is:

$\begin{array}{c}{E}_{k}^{2}={\displaystyle \sum _{{k}^{\prime}\ne k}\frac{{\left|{}^{(0)}\langle k|V{|{k}^{\prime}\rangle}^{(0)}\right|}^{2}}{{E}_{k}^{0}-{E}_{{k}^{\prime}}^{0}}}\\ =\frac{{\left|{}^{(0)}\langle k|V{|k+K\rangle}^{(0)}\right|}^{2}}{{E}_{k}^{0}-{E}_{k+K}^{0}}+\frac{{\left|{}^{(0)}\langle k|V{|k-K\rangle}^{(0)}\right|}^{2}}{{E}_{k}^{0}-{E}_{k-K}^{0}}\\ =\frac{{\left|{V}_{K}\right|}^{2}}{{E}_{k}^{0}-{E}_{k+K}^{0}}+\frac{{\left|{V}_{-K}\right|}^{2}}{{E}_{k}^{0}-{E}_{k-K}^{0}}.\end{array}$

This result is reasonable provided the terms are small, that is, the energy differences appearing in the denominators are large compared to the relevant Fourier component ${V}_{K}.$ However, this cannot always be true! Notice that the state $k=\pi /a$ has exactly the same unperturbed energy ${E}^{0}$ as the state $k-K=-\pi /a$: in this case, nondegenerate perturbation theory is clearly wrong. In fact, even for states close to $k=\pi /a$, the energy denominator ${E}_{k}^{0}-{E}_{k-K}^{0}$ is small compared with the numerator ${\left|{V}_{-K}\right|}^{2}$, so the series is not converging.

### Quasi-Degenerate Perturbation Theory near the Critical Wavelength

The good news is that, despite the many states near $k=\pi /a$ and $k=-\pi /a$ that are close together in energy, for any *one* state $k$ near $\pi /a$ the potential only has a nonzero matrix
element to *one* other state close in
energy, the state $k-K$,*
*that is, $k-2\pi /a$.

The strategy now is to do what might be called *quasidegenerate*
perturbation theory: to diagonalize the full Hamiltonian in the subspace
spanned by these two states ${|k\rangle}^{(0)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{|k-K\rangle}^{(0)}$. Other states with nonzero matrix elements to
these states are relatively much further away in energy, and can be treated
using ordinary perturbation theory.

The matrix elements of the full Hamiltonian in the subspace spanned by these two states are:

$\left|\begin{array}{cc}{E}_{k}^{0}& {V}_{K}^{*}\\ {V}_{K}& {E}_{k-K}^{0}\end{array}\right|.$

Diagonalizing *within
this subspace* gives energy eigenvalues:

${E}_{\pm}=\frac{1}{2}({E}_{k}^{0}+{E}_{k-K}^{0})\pm \sqrt{{\left(\frac{{E}_{k}^{0}-{E}_{k-K}^{0}}{2}\right)}^{2}+|{V}_{K}{|}^{2}}.$

Notice that, provided $\left|{E}_{k}^{0}-{E}_{k-K}^{0}\right|\gg \left|{V}_{K}\right|$,
to leading order this gives back ${E}_{\pm}={E}_{k}^{0},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{E}_{k-K}^{0},$ the order depending on $k.$ However, as $k$ approaches $\pi /a$, $\left|{E}_{k}^{0}-{E}_{k-K}^{0}\right|$ becomes of order $\left|{V}_{K}\right|$,
and the energies deviate from the unperturbed values. If $k$ is approaching $\pi /a$ from below, ${E}_{k}^{0}{E}_{k-K}^{0}$,
and the lower energy is pushed *downwards*
by the perturbation: ${E}_{k}={E}_{-}{E}_{k}^{0}.$ This is
a common occurrence with almost degenerate states, perturbations cause the
energy levels to “repel” each other.

$\text{For}k=\pi /a,\text{}{E}_{k-K}^{0}={E}_{k-2\pi /a}^{0}={E}_{k}^{0}.$ At this value of $k,$ the unperturbed states are exactly degenerate, and the perturbation lifts the degeneracy to give ${E}_{\pm}={E}_{\pi /a}^{0}\pm \left|{V}_{K}\right|.$

In the graph below, the green (continuous) curve is the unperturbed energy as a function of $k,$ the red curve (with the step) the calculated energy including the leading correction from the periodic potential.

### Energy Gaps and Bands

The energy jump, or gap, of $2\left|{V}_{K}\right|$ at $\left|k\right|=\pi /a$ means that there are no plane wave type
eigenstates with energies in that range$\u2014$attempting
to integrate Schrödinger’s equation in the periodic potential for such an
energy gives exponentially growing and decaying solutions. Such energy gaps in fact are present in real
crystalline solids, the allowed energies are said to be in “bands”. The lowest band for our model is from $k=-\pi /a\text{to}\pi /a.$ Since the allowed values of $k$ are given by $k=2\pi n/L,$ the spacing between adjacent $k$ ’s is $2\pi /L$ and the total number of $k$ ’s in the lowest band is $L/a=N,$ the same as the number of atoms. Since each electron has two spin states, this
implies that a one-dimensional crystal of *divalent*
atoms will *just fill* the lowest band
with electrons. Therefore, any outside
field can only excite an electron to a different state if an energy of at least
$2\left|{V}_{K}\right|$ is supplied$\u2014$for a
small electric field, the filled band of electrons will remain in the ground
state, there will be no current. This
material is an insulator.

On the other hand, if *monovalent*
atoms are used, it is clear that the lowest band is* only half full*, adjacent empty electron states are available. The electrons are free to accelerate if an
external field is applied. Barring the
unexpected, this one-dimensional crystal would be a metal.

Let us now examine how the periodic potential alters the eigenstates. Ignoring the small corrections from plane waves outside the ${|k\rangle}^{(0)},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{|k-K\rangle}^{(0)}$ subspace, the eigenstates to this order have the form

$|k\rangle ={a}_{k}{|k\rangle}^{(0)}+{a}_{k-K}{|k-K\rangle}^{(0)}$

where

$$\frac{{a}_{k-K}}{{a}_{k}}=\frac{{E}_{-}-{E}_{k}^{0}}{{V}_{K}^{*}}$$

from the diagonalization of the $2\times 2$ matrix representing the Hamiltonian in the subspace.

As $k$ increases from 0 towards $\pi /a,$ the plane wave initially proportional to ${e}^{ikx}$ has a gradually increasing admixture of ${e}^{i\left(k-2\pi /a\right)x}$,
until at $k=\pi /a$ the two have equal weight$\u2014$meaning
that the eigenfunction is now a standing wave.
In fact, there are *two* standing wave solutions at $k=\pi /a$,
corresponding to the energies below and above the gap. Taking the atoms to have
an attractive potential, the lower energy wave has a probability distribution
peaking at the atomic positions. The
diffractive scattering that gives a left-moving component to a right moving
wave is known as Bragg scattering. It
also manifests itself in the *group
velocity* of the electronic excitations, ${v}_{\text{group}}=d\omega /dk=\left(1/\hslash \right)dE/dk.$ An electron injected into a one-dimensional
metal would not be a plane wave state, but a wavepacket traveling at the group
velocity. It is evident that for an injected electron with mean value of $k$ close to $\pi /a$,
the electron will move very slowly into the metal. This is to be expected$\u2014$the
eigenstates become standing waves as $k\to \pi /a.$

For three-dimensional crystals, the situation is far more complicated, but many of the same ideas are relevant. Electron waves are now diffracted by whole planes of atoms, and the three-dimensional momentum space is divided into Brillouin zones, with planes having an energy gap across them.

### The Peierls Transition: how Cooling a Conductor Can Give an Insulator

As mentioned in the Introduction, substances very close to
monovalent one-dimensional crystals have been synthesized, and it has been
found$\u2014$surprisingly$\u2014$that at
low temperatures many of them undergo a transition from metallic to *insulating* behavior. *What
happens is that the atoms in the lattice rearrange slightly, moving from an
equally-spaced crystal to one in which the spacing alternates, that is, the atoms
form pairs*. This is called *dimerization*,
and costs some elastic energy, since for identical atoms the lowest state must
be one of equal spacing for any reasonable potential. However, the *electrons* are able to move to a lower energy state by this
maneuver.

Just how this happens can be understood using the
perturbation theory analysis above. For
equally spaced atoms, the electrons *half-fill
*the band, that is, they fill it up (two electrons, one of each spin, per
state) to $\left|k\right|=\pi /2a.$

The crucial point is that if the atoms move together
slightly into pairs, *the crystal has a
new period* $2a$ *instead
of* $a.$ This
means that the *potential* now has a
nonzero component at $K=-\pi /a$,
with a nonzero matrix element between the states $k=\pi /2a$ and $k=-\pi /2a$,
and so on. From this point, we can rerun
the analysis above, except that now the gaps open up at $\left|k\right|=\pi /2a$ instead of at $\left|k\right|=\pi /a.$

The important point is that if the electrons fill all the
states to $\left|k\right|=\pi /2a$,
and none beyond (as would be the case for monovalent atoms) then the opening of
a gap at $\left|k\right|=\pi /2a$ means that *all
the electrons are in states whose energy is lowered*. To find the *total* energy benefit we need to integrate over $k.$

### Calculating the Electronic Energy Gained by Doubling the Lattice Period

It is evident from the above that most of the contribution
comes from fairly close to $k=\pi /2a$ (and of course symmetrically $k=-\pi /2a$ ).
Since we want to find the total lowering in energy, let us study first
the *bare* energy as a function of $k,$ that is, the energy with no potential
present. Of course, there isn’t much to
say: ${E}_{k}^{0}={\hslash}^{2}{k}^{2}/2m$. However, the physics of these one-dimensional
systems concerns only excitations near the “Fermi surface”, the boundary
between filled (low energy) states at zero temperature and empty states. This “Fermi surface” is in fact just two
points in one dimension: $k=\pm \text{\hspace{0.17em}}\pi /2a$.
In the neighborhood of these two Fermi points,
it is an excellent approximation to replace the gently curving ${E}_{k}^{0}={\hslash}^{2}{k}^{2}/2m$ by *straight
line approximations*$\u2014$the slope
being $dE/dk={\hslash}^{2}k/m=\hslash p/m=\hslash v$.

Linearizing in the neighborhood of $k=\pi /2a$, then, we take

${E}_{k}^{0}={E}_{\pi /2a}^{0}+\hslash v\left(k-\pi /2a\right)={E}_{\pi /2a}^{0}+\hslash vq$,

where

$q=k-\pi /2a,$

just $k$ measured from the Fermi point $\pi /2a$.

The variable $q$ is *negative* for the relevant states, since
they are on the lower energy side.

The density of states in $k\text{-}$ space is a constant $2\times L/2\pi =L/\pi $,
remembering the two spin states per *$k\text{-}$** *value.

Recall

${E}_{\pm}=\frac{1}{2}({E}_{k}^{0}+{E}_{k-K}^{0})\pm \sqrt{{\left(\frac{{E}_{k}^{0}-{E}_{k-K}^{0}}{2}\right)}^{2}+|{V}_{K}{|}^{2}}$

but now

$K=\pi /a$

and the lowering of energy of the electrons (counting it as a positive quantity) is:

$2{\displaystyle \underset{0}{\overset{\pi /2a}{\int}}\left({E}_{k}^{0}-{E}_{-}\right)}Ldk/\pi =2{\displaystyle \int \left(\frac{1}{2}({E}_{k}^{0}-{E}_{k-K}^{0})+\sqrt{{\left(\frac{{E}_{k}^{0}-{E}_{k-K}^{0}}{2}\right)}^{2}+|{V}_{K}{|}^{2}}\right)}\text{\hspace{0.05em}}\text{\hspace{0.17em}}Ldk/\pi $

where the extra factor of 2 counts the symmetrical
contribution from the left-hand gap. (In
examining the above expression, recall that for the $k0$ states we are interested in, $k\pi /2a,$ ${E}_{k}^{0}-{E}_{k-K}^{0}$ is *negative*. The integrand on the right-hand side is still
positive, very small for small *k*,
reaching a maximum of $\left|{V}_{K}\right|$ at $k=\pi /2a.$ )

Putting in our linearized energy approximation,

${E}_{k}^{0}={E}_{\pi /2a}^{0}+\hslash v\left(k-\pi /2a\right)={E}_{\pi /2a}^{0}+\hslash vq$,

and remembering that now *$K=\pi /a,$*

${E}_{k-K}^{0}={E}_{-\pi /2a}^{0}-\hslash v\left(k+\pi /2a\right)={E}_{-\pi /2a}^{0}-\hslash vq$.

Since ${E}_{\pi /2a}^{0}={E}_{-\pi /2a}^{0}$,

$({E}_{k}^{0}-{E}_{k-K}^{0})=2v\hslash q.$

Substituting these linearized values in the integral for the total energy lowering:

$\text{\hspace{0.17em}}2{\displaystyle \underset{0}{\overset{\pi /2a}{\int}}\left({E}_{k}^{0}-{E}_{-}\right)}Ldk/\pi =2{\displaystyle \underset{-D}{\overset{0}{\int}}\left(v\hslash q+\sqrt{{(v\hslash q)}^{2}+|{V}_{K}{|}^{2}}\right)Ldq/\pi}$

where in terms of the variable $q$ we have set the lower limit of integration at $-D$: we can safely be vague about this lower limit, as the integral turns out to be logarithmic.

Since the integral is over negative numbers, and we have taken the positive square root, it is zero for zero ${V}_{K},$ as it must be.

The integral can be done exactly, but it is more illuminating to divide the range of integration into $\left|v\hslash q\right|\le \text{\hspace{0.17em}}\left|{V}_{K}\right|$ and $\left|v\hslash q\right|\text{\hspace{0.17em}}\left|{V}_{K}\right|$, then estimate the contributions from these two ranges separately.

First, consider $\left|v\hslash q\right|\le \text{\hspace{0.17em}}\left|{V}_{K}\right|$. Here the integrand is of order $\text{\hspace{0.17em}}\left|{V}_{K}\right|$, and the region $\Delta q$ of integration corresponding to $\left|v\hslash q\right|\le \text{\hspace{0.17em}}\left|{V}_{K}\right|$ is of order $\text{\hspace{0.17em}}\left|{V}_{K}\right|/\hslash v,$ so the integral over this range is of order $\left(L/\hslash v\right){\left|{V}_{K}\right|}^{2}$.

Second, in the region $\left|v\hslash q\right|\text{\hspace{0.17em}}\left|{V}_{K}\right|$, we can write

$\begin{array}{l}2{\displaystyle \int \left(v\hslash q+\sqrt{{(v\hslash q)}^{2}+|{V}_{K}{|}^{2}}\right)}\text{\hspace{0.05em}}\text{\hspace{0.17em}}Ldq/\pi =2{\displaystyle \int \left(v\hslash q+\left|v\hslash q\right|\sqrt{1+\frac{|{V}_{K}{|}^{2}}{{(v\hslash q)}^{2}}}\right)}\text{\hspace{0.05em}}\text{\hspace{0.17em}}Ldq/\pi \\ \end{array}$

and expand the square root term. The leading terms cancel since $q$ is negative, and the main contribution comes from the next term. This gives:

$\Delta E\approx 2|{V}_{K}{|}^{2}{\displaystyle \underset{-D}{\overset{-\left|{V}_{K}\right|}{\int}}\frac{1}{2\hslash v}\frac{Ldq}{\pi \left|q\right|}}=\frac{L|{V}_{K}{|}^{2}}{\hslash v}\mathrm{ln}\frac{\left|{V}_{K}\right|}{D}.$

*The important thing
here is the logarithm*. For
sufficiently small $\left|{V}_{K}\right|$,
this large (negative) term will dominate any term which is just proportional to
${V}_{K}^{2}$. But the elastic energy cost of the lattice
“dimerizing”$\u2014$the atoms
forming pairs, so that the distance between atoms alternates on going along the
chain$\u2014$must be
proportional to ${V}_{K}^{2}$. This leads to the conclusion that some,
probably small, dimerization is always going to happen*$\u2014$**a one-dimensional equally spaced chain with
one electron per ion is unstable*.

This dimerization is known as a *Peierls transition*. Peierls
discovered it in the 1930’s when writing a section on one-dimensional models in
an introductory solid-state textbook. He
put it in the book, but didn’t publish it otherwise. As mentioned in the Introduction, it became
very relevant later when some theories suggested that quasi-one-dimensional
conductors, materials made up of loosely connected chains, each chain having
one electron per atom for a half-filled lowest band, might be high-temperature
superconductors. It was found instead
that many such materials actually became *insulators*
on cooling: the reason was that at high temperatures, the electrons filled
states above and below the point $\pi /2a$ fairly equally, so dimerization did not lower
the overall energy much. On lowering the
temperature, a point was reached where the Peierls transition gave a lower
energy state, and the material became an insulator.

*Reference*: Rudolf
Peierls, *More Surprises in Theoretical Physics*,