Critical points in condensed matter illuminate universality

Every person is unique. No two people are identical. We differ in physical appearance, personality, fingerprints, heartbeat, gait, and DNA. Such differences are used to identify criminals and in video surveillance of citizens by nation states. Yet in other ways all humans are the same. We all have brains, hearts, and lungs. All our bodies use the same biochemistry to stay alive: whether to breathe oxygen, digest food, or fight infections. On some level we have common aspirations: to survive, to be loved, to be happy, and to find meaning and purpose. Yet these aspirations find many expressions. Humans have certain universal qualities and properties, yet at a finer level of detail there is a particularity of each of these properties. They are at one level the same but are not the same at another level. 

All academic disciplines search for universals; they develop categories, concepts, and theories that overarch particularities. Biologists classify species of plants and animals and types of cells and viruses. All biological systems use the same molecules (DNA, RNA, and proteins) and chemical reactions. The same genetic code uses the information encoded in a piece of DNA to make proteins with specific functions. Anthropologists study the immense diversity of human cultures and societies. This diversity can be described in terms of universal concepts such as kinship, family, ritual, community, economics, law, and morality. Linguists study the common structures and grammars of the thousands of different human languages.  Although the world we live in is diverse, disciplines have each discovered some universals.

Condensed matter physicists study diverse states of matter and the transitions between them. A surprising discovery is that there is much more universality than might be expected, particularly given the chemical and structural diversity of materials. In this chapter, I will discuss the nature of this universality, how it emerges, and the length scales associated with transitions between different states of matter. Landau’s great insight was that many of the chemical and structural details of materials are irrelevant to understanding phase transitions. Furthermore, a precise classification of different types of phase transitions, into what are called universality classes, can be made. For example, superconducting, superfluid, and a subset of magnetic transitions are in the same class. The determinants of the universality classes are the symmetry of the state and the spatial dimensionality of the system. None of the other details matter.

Many phase diagrams (such as the Figure above) include a critical point, located at the end of a boundary between two different states of matter. A common example is the critical point that occurs at a specific temperature and pressure for a transition between a liquid and a gas. Understanding the physical properties of a material close to its critical point was a great challenge for theoretical physics, lasting a hundred years, and was only solved in the 1970s. The powerful theoretical ideas and techniques that were developed provide a quantitative way to relate the properties of a system at one length scale to properties at a different length scale. These techniques also have application to a wide range of other problems and fields including elementary particle physics, chaos theory, fractals, polymers, and machine learning. New insights were gained into universality and emergent phenomena.

An extract from "The Critical Point," chapter 6, Condensed Matter Physics: A Very Short Introduction

Quantum justification for classical discussions of potential energy surfaces in chemistry

 In computational quantum chemistry, the Born-Oppenheimer approximation (BOA) is used to determine potential energy surfaces (PES) for electronic states of molecules. It is standard practise to identify local minima on a PES with molecular structures. Molecular binding energies are identified with the difference in energy between these minima and the energies of the isolated atoms of which the molecule is composed.

 A chemical reaction between two molecules A and B to produce C can be understood in terms of the PES for the composite system consisting of all the atoms in A and B. The dynamics of the chemical reaction can be described in terms of a path on the PES that goes from the local minimum associated with A and B infinitely far apart to the minimum associated with the structure C. The path will pass through a saddle point on the PES and this is identified with a transition state in the chemical reaction and its energy determines the activation energy for the chemical reaction. The local curvature of the PES near a minima can be used to determine force constants for harmonic motion and the associated vibrational frequencies of a molecule.

This picture is a completely classical one and so may motivate a claim that chemists mix classical and quantum concepts and calculations in an ad-hoc manner. (This kind of argument is often used by philosophers to claim that chemistry cannot be reduced to physics). This is unfair because a quantum description of the nuclear dynamics can be given in terms of quantum wave-packets, consistent with the Heisenberg uncertainty principle, and whose dynamics is defined by the quantum equation for nuclear motion that is given by the BOA. Furthermore, the dynamics of the centre of a wave-packet is given by classical equations of motion (Ehrenfest’s theorem). Thus, the classical language used by chemists can be viewed as a justified and compact version of a quantum description.

Structural isomers. Isomers are associated with different local minima on the electronic ground state potential energy surface for a given combination of atoms. To understand a quantum description of isomers, consider a reaction coordinate associated with an isomerisation reaction (i.e. conversion of one isomer to the other). There are three energy scales of relevance: the energy difference between the ground state energy of the two isomers, the magnitude of the barrier height (activation energy), and the quantum zero-point energy associated with vibrations in the direction of the reaction coordinate. Denote these energies as dE, Eb , and Ezp, respectively. If dE ~ Ezp << Eb then the nuclear probability density rho(R) for the vibrational ground state will have two local maxima, corresponding to the geometries of the two isomers. If dE >> Ezp then the nuclear probability density rho(R) for the vibrational ground state will have only one local maxima, corresponding to the lower energy isomer geometry. However, the geometry of the higher energy isomer can be found as a local maximum in the nuclear probability density rho(R) for one of the excited vibrational ground states. 

 Figure. Potential energy surface associated with the two structural isomers of HOCO. TS_n denote different transition states associated with the chemical reaction OH + CO -> H + CO2. Taken from Bui et al.

Condensed matter physics in flatland

Adventures in Flatland

In everyday life we think of most objects as having three dimensions. But what would life be like in a two-dimensional world? For one thing, it would be harder to move around. We could no longer step over things but would have to move around them. In 1884 Edwin Abbott published Flatland: A Romance of Many Dimensions, under the pseudonym, A. Square, a satirical novella about social life in Victorian England. People are represented by geometrical objects. Men are represented by shapes such as triangles and hexagons. Women are represented by lines. The social status of men increases with the number sides that their shape has and how many of the sides are of the same length. Abbott’s book created limited interest and was largely forgotten by the 1920s. Interest revived when theoretical physicists started to think about worlds in different dimensions. This interest was stimulated by Albert Einstein’s theories of relativity, that proposed that we live in a four-dimensional world, not a three-dimensional one. Time is the fourth dimension, and there is an intimate and concrete connection between time and space. Attempts to unify gravity with other fundamental forces has led to physicists proposing and studying theories with more than four dimensions.

Changing the number of spatial dimensions leads to different physics because it changes what is mathematically possible. In three dimensions, there were only five highly symmetrical shapes known as Platonic solids (tetrahedron, cube, octahedron, icosahedron, and dodecahedron). In contrast, in two dimensions it is possible to make an infinite number of symmetrical shapes, known as regular polygons, shapes made of straight lines of equal length such as squares or hexagons. Similarly, the number of Bravais lattices differ in two and three dimensions. Changing the number of spatial dimensions changes both what is mathematically possible and what is physically possible.

What would condensed matter physics be like in Flatland? This question received limited attention before the 1970s. Occasionally, theoretical physicists would investigate mathematical models of crystals or magnets in one or two dimensions just because the mathematics was simpler and more tractable than in three dimensions. The goal was to obtain insight into physics in three dimensions. We will consider a famous example, the Ising model. 

In the 1970s, several surprising developments led to significant interest in condensed matter physics in spatial dimensions different from the usual three. First, it became possible to make a wide range of material systems that were two-dimensional. Secondly, theoretical work showed that states of matter, and phase transitions between them, can be qualitatively different in one, two, and three spatial dimensions. And thirdly, considering different numbers of spatial dimensions turned out to be very fruitful for theory, particularly for understanding phase transitions near critical points. 

An extract from "Adventures in Flatland," chapter 5 in Condensed Matter Physics: A Very Short Introduction

What does the Born-Oppenheimer approximation mean for emergence?

Most philosophical debates about the emergence of molecular structure centre around the issue of irreducibility. Specifically, can the existence of structures be predicted from quantum theory without assuming their existence or invoking classical concepts? I will argue that the answer is yes, contrary to much of the philosophical literature, which relies heavily on the widespread use of the Born-Oppenheimer approximation (BOA) in quantum chemistry calculations. However, the fact that these arguments for irreducibility are weak does not mean that emergence (defined in terms of novelty) is not central to chemistry.

In a previous post, I discussed recent work showing how the BOA is not necessary for quantum chemistry and that molecular structure can be defined independently of it.

However, since the BOA plays a central role in the philosophical arguments, it is worth reviewing what it is and what it does and does not assume or mean.

In 1927, Born and Oppenheimer introduced an approximation to allow the solution of the full quantum equations for electrons interacting with charged nuclei. Without the BOA, much of theoretical chemistry and solid-state physics would be incredibly difficult in practice. The approximation is based on the separation of time and energy scales associated with electronic and nuclear motion. It leads to the concept of potential energy surfaces for electronic states. They define an effective theory for the dynamics of the atomic nuclei in a molecule or solid.

The full Hamiltonian (given earlier) can be denoted by

 

where the first term is the kinetic energy operator for the nuclei. In the Born-Oppenheimer approximation (BOA) the full wavefunction is written as a product of a nuclear wavefunction and an electronic wavefunction.

Substituting this in the eigenvalue equation for the full Hamiltonian leads to separate eigenvalue equations for the electronic and nuclear wavefunctions, assuming terms depending on gradients with respect to R of the electronic terms can be neglected.

 

In the first equation, the nuclear co-ordinates appear as parameters not as operators. This is central to the philosophical debates.

The second equation can be viewed as an effective Hamiltonian for the nuclear degrees of freedom. The function E_e(R) defines the potential energy surface of the molecule. 

In the BOA the nuclear probability distribution defined above is

 

As discussed in the earlier post, the structure of many molecules can be defined in terms of the value of R at which the probability is maximum. 

I make four points about the BOA that are relevant to philosophical debates about whether molecular structure is predictable in a logically consistent manner from quantum theory.

1. The BOA does treat the nuclear degrees of freedom quantum mechanically. They are described by the nuclear wavefunction Phi(R), which is determined by the second eigenvalue equation. Consequently, the BOA does not violate Heisenberg’s uncertainty principle, contrary to some claims in the philosophy literature.

2. The BOA is not ad hoc. Corrections to it can be calculated and have been for many molecules. These corrections are typically small, being of order (me/Mi)^1/2. Exceptions, such as near conical intersections (where the potential energy surfaces for two electronic states touch) are well-known and well-studied.

3. For most small molecules, the results of BOA calculations compare favourably with wave-functions obtained from solutions of the full quantum Hamiltonian. When there are differences, they are largely small quantitative differences. When the differences are qualitative, they have largely been anticipated from knowledge of the limitations of the BOA.

4. The Born-Oppenheimer approximation is an example of a general approach to quantum mechanics problems, discussed by Migdal. Consider a system composed of two subsystems that have dynamics on two vastly different time scales, termed fast and slow. The effects of the fast system on the slow system can be treated by adding a potential energy term to the Hamiltonian operator of the slow system. 

In forthcoming posts, I will discuss quantum justifications for classical descriptions of nuclear dynamics on potential energy surfaces and then discuss philosophers' views about the BOA and molecular structure.

Condensed matter physics is about how order emerges from disorder

 The order of things

Life and the world around us sometimes appears chaotic and random. We may feel this way about traffic, weather, economics, social change, politics, or our personal relationships. Perhaps that is why many yearn for regularity, predictability, order, and stability. Science is a search for patterns and order in the natural world. Condensed matter physics is about how order emerges from disorder.

This chapter explores how different states of matter are associated with different types of ordering of the atoms in the material. The symmetry of the state reflects the type of ordering, i.e., the patterns associated with the state. There is also a rigidity associated with the ordering and the rigidity determines the nature of the deviations from perfect ordering and results in entities such as vortices that are central to the physical properties of the state of matter.

The association of a state of matter with a specific type of ordering is illustrated in Figure 15 by an analogue with the dodgem bumper cars at an amusement park. A quiet day at the park is not much fun as collisions between cars are rare. In other words, there is little correlation between the relative locations and speeds of the cars. In comparison, on a busy day at the park the spatial separation of the cars is small, and their positions and speeds are more correlated with one another than on a quiet day. But, in both cases, there is no ordered arrangement of the cars. In contrast, after the park closes the cars are parked and arranged in an orderly manner. There is a rigidity associated with their spatial arrangement. One car cannot be moved without moving others. These three states of the dodgem cars are an analogue of three states of matter: gas, liquid, and crystal. 

Figure 15. A dodgem car analogue for the three states of matter: crystal, liquid, and gas. The only ordered arrangement is for the crystal (car park after hours) and this is associated with a specific symmetry and rigidity. The liquid and gas (busy and quiet day) only differ in density and the amount of correlation between the positions of the different atoms (dodgem cars).

In the dodgem car analogue, there are other possible types of ordering. In some amusement parks there is a track, and the cars are meant to all go in the same direction. The symmetry between clockwise and anti-clockwise of the track is then broken.  In the car park, Figure 15 shows cars that are symmetrical with respect to front and back. However, real cars have a front and back, and so can be parked either front first or back first. Hence, several types of ordering are possible: all cars park back first, all cars park front first, cars are front first or back first at random, alternating patterns of front first and back first as one goes along a row, alternating rows of front first and back first, and so on. These different types of ordering in the car park all have analogues in different solid states of matter.

Liquid crystals involve unique types of ordering. These materials are composed of elongated organic molecules, such as those shown in Figure 16. At high temperatures the material is in a liquid state and the orientations and positions of the molecules are random. The liquid has both continuous translational and rotational symmetry. At low temperatures the molecules form a solid crystal without the continuous translational and rotational symmetry of the liquid state. As the crystal is heated the temperature increases and there is a phase transition to the liquid crystal state, in which all the molecules point in the same direction, but their positions are random. Hence, the liquid crystal state has the continuous translational symmetry of the liquid, but not its continuous rotational symmetry, like the crystal. As the temperature increases further there is a transition to the liquid state (Figure 16). In terms of the dodgem car analogue the liquid crystal state is similar to when cars park in a field all pointing in the same direction but there are no grid lines, and their positions are then random.

The existence of a state in between a liquid and crystal was first proposed in 1888 by botanist and chemist Friedrich Reinitzer who was doing research on cholesterol at the Institute for Plant Physiology in Prague. He performed a heating experiment similar to that described in Figure 4. Instead of one melting transition he observed transitions at two distinct temperatures. 

Figure 16. Liquid crystals. (a) An example of the type of elongate organic molecule found in these materials. Each molecule can be represented by an oval shape. (b) In the nematic liquid crystal state, the molecules tend to point in the same direction, but their positions are random. 

There are multiple alternative orderings for liquid crystals with names such as nematic, smectic, chiral nematic, discotic, and chlorestic. In the smectic phase molecules form layers of oriented molecules. The character of the liquid crystal state can be detected by shining polarised light on the material. Liquid crystal displays (LCDs) in electronic devices use the property that an electric field can orient the molecules, and this changes the interaction of the material with polarised light.

For solid crystals the nature of the ordering and the symmetry associated with a specific crystal structure is clear once the spatial arrangements of the atoms in the crystal are determined, such as by X-ray diffraction. For other states of matter, such as superconductors, superfluids, and antiferromagnets, the nature of the ordering and the symmetry is often not apparent and has only been determined with significant scientific insight. 

An extract from "The order of things," chapter 4 in Condensed Matter Physics: A Very Short Introduction.

The emergence of molecular structure from quantum theory

Most debates about the emergence of molecular structure centre around the issue of irreducibility. Specifically, can the existence of molecular structures be predicted from quantum theory without assuming their existence or invoking classical concepts?

Consider a molecule that contains Ne electrons and Nn atomic nuclei (ions). The full quantum-mechanical Hamiltonian for the system is 

where e is the electronic charge, rj is the position of the j-th electron, Zi  and Mi are the charge and mass, respectively, of the i’th ion with position co-ordinate Rj. This is the Hamiltonian that Laughlin and Pines dubbed “The Theory of Everything” because if the solution (i.e., eigenstates and eigenvalues of the Hamiltonian operator) could be found it would describe almost all of chemistry and materials science.

This Hamiltonian treats the electrons and nuclei on an equal footing. 

For isomers, the Hamiltonian is identical. However, as will be discussed in a later post, that does not preclude solutions to the Hamiltonian that can describe isomers. 

The Hamiltonian has global translational and rotational symmetry, where all the particles undergo the same rotation or translation. In contrast, molecular structures may have discrete rotational symmetries. However, this is not necessarily a problem, as an eigenstate of a quantum problem can transform according to a non-trivial irreducible representation of the symmetry. For example, except the s-orbitals all the orbitals of the hydrogen atom are spatially anisotropic.

The electrons are identical particles and so have permutation symmetry. They are fermions with spin-1/2 and so any eigenstate must be antisymmetric under the exchange of two electrons. The energies associated with this exchange are crucial to the formation of chemical bonds and the stability of molecular structures.

If two or more atoms in the molecule are identical, then any exact eigenstate must be consistent with permutation symmetry. If a nucleus is composed of an even (odd) number of nucleons, then it is a boson (fermion) with integer (half-integer) spin, and eigenstates must be symmetric (antisymmetric) under exchange of identical nuclei. However, the corresponding exchange energies are relatively small (because the quantum delocalisation of the nuclei is small) and consequently most practical calculations of the eigenstates do not make this requirement of the eigenstates. Nevertheless, if the electrons and nuclei are treated on equal footing, this should be done. Although this is challenging, it has been done recently, as discussed below. 

Full quantum solutions of the Hamiltonian

In most computational quantum chemistry, the Hamiltonian is solved in the Born-Oppenheimer approximation, which will be introduced and discussed later. This is a source of some confusion and contention in philosophical discussions about the emergence of molecular structure.

Due to advances in methodology and computational power over the past few decades, it has become possible in practise to solve the full quantum Hamiltonian for small molecules. There are three levels of complication associated with this: quantum nuclear motion, rotational symmetry, and some nuclei being identical particles. There are also two challenges: first, finding the eigenstates and second, deducing the molecular structure from the eigenstates.

To begin, I consider the simplest case and ignore the complications associated with rotational symmetry or identical nuclei. This provides some insight and undermines some objections in the philosophical literature.

The ground state eigenfunction can be written as

where r and R are 3Ne and 3Na -dimensional vectors, respectively. Note that this function will have a complicated structure as it will depend on the spin states of all the electrons, denoted by s.

A probability distribution (reduced density matrix) for the positions of the nuclei is given by

where the sum is over all the electron spin degrees of freedom.

For many molecules, but not all, this probability distribution will have a unique global maximum at the coordinates R_0. This set of coordinates defines the geometry of the molecular structure. The physics underlying the existence of well-defined maxima is that the mass of the nuclei is much larger than the mass of the electrons, and as a result, the zero-point motions of the nuclei are much smaller than the separation of the nuclei in the molecular structure.

Note that the nuclear probability distribution is regularly measured in scattering experiments (using X-rays, neutrons, or electrons), and its maxima are used to determine the structures of molecules and crystals. The Debye-Waller factor is a measure of the width of the probability distribution. At low temperatures, it is determined by quantum zero-point motion. In other words, it is well established experimentally that classical molecular structures are an approximation to a fluctuating quantum structure.

Not every molecule will have a probability distribution with a unique maximum. An example is ammonia. As discussed further below, it has two maxima; each represents an umbrella geometry, and they are related by an inversion symmetry. The ground state wavefunction of the whole system is a superposition of two quantum states, each being associated with one of the two umbrella geometries, and the electronic and nuclear degrees of freedom are entangled with one another.

A general quantum definition of molecular structure

Lang et al. have recently overcome the challenges mentioned above to determine molecular structure in a manner that treats the electrons and nuclei on an equal footing with regard to quantum theory. They have considered both rotational symmetry and nuclear permutation symmetry and given a general definition of molecular structure involving nuclear probability densities calculated from the full wavefunction. They have explicitly performed these calculations for D3+, (where D is deuterium). The result is that the molecule has the same triangular structure that is observed experimentally and calculated using the Born-Oppenheimer approximation. This work is significant because it explicitly shows that molecular structure can be predicted in practice, not just in principle, from quantum theory.

In a forthcoming post, I will discuss the Born-Oppenheimer approximation and some of the confusion associated with it.

Symmetry matters in condensed matter physics

 Snowflakes form incredibly diverse structures, seen when they condense onto a plate of glass. Every snowflake is different. On the other hand, every snowflake is the same. They are all composed of ice, a solid state of water. Every snowflake is composed of units that have a six-fold symmetry (Figure 8). Every snowflake is composed solely of water molecules. This paradox of the particular and the universal is at the heart of condensed matter physics. Although diversity prevails anything is not possible. No snowflake has five-fold symmetry. Snowflakes have enchanted scientists for a long time. The astronomer Johannes Kepler studied them and in 1611 wrote a small book about them as a gift for his patron. Kepler suggested snowflakes provided clues to deeper questions about the composition of matter. Today, Kenneth Libbrecht, a physicist at Caltech, has spent most of his career studying snowflakes and has produced beautiful volumes of photographs of them.

Figure 8. A snowflake shows a six-fold symmetry, just like a hexagon. The snowflake appears identical when it is rotated by an angle of sixty degrees about an axis passing through its centre and perpendicular to the page.

Condensed matter physicists ask several questions about snowflakes. What is the reason for the six-fold symmetry of the snowflake? What is the connection between the macroscopic properties of snowflakes and the properties of the underlying microscopic constituents, molecules of H2O? How is the diversity of snowflake shapes possible? Is there a phase diagram that defines the external conditions under which the different shapes form?

There is a long history in art, architecture, philosophy, and science, of associating symmetry with beauty and perfection. The ancient Greek philosopher Plato was a proponent of this view. He studied a particular class of solid shapes: cube, tetrahedron, octahedron, icosahedron, and dodecahedron. Plato identified the first four shapes with the four “elements”: earth, wind, fire, and water, respectively, and the fifth with the heavens. Each of these solid shapes is highly symmetric. Every face of a Platonic solid is the same shape (square, triangle, pentagon,...) and each of those shapes has edges of equal length. 

Like Plato, Kepler believed that “God is a geometer” and that God’s creation should reflect the perfection of God. These convictions led Kepler to propose in 1597 that the orbits of the planets around the Sun were circular and that the Platonic solids determined the relative size of the orbits. Later this model for the solar system was shown not to be true. In fact, Kepler himself became famous because he showed that the planets moved in elliptical, not circular orbits. Nevertheless, Kepler’s model was the beginning of a long history of successfully relating physical laws to symmetry and geometry.

A key discovery in physics from the past century is that symmetry is central to understanding a wide range of physical phenomena, whether colliding billiard balls, the allowed energies of an atom, the fundamental forces of nature, or different states of matter. Symmetries determine what is physically possible. For example, that energy cannot be created or destroyed is a consequence of the fact that physical laws do not change with time.

In this Chapter I explore three key ideas. First, transitions between different states of matter are associated with changes in symmetry. Thus, symmetry provides a criterion for specifying the qualitative difference between distinct states of matter. Second, for a specific state of matter the relevant symmetry constrains what is physically possible. Third, symmetry is central to making connections between the macroscopic and microscopic properties of a state of matter. The next chapter will explore how symmetry is associated with the type of ordering that occurs in a state of matter.

An extract from Chapter 3, Condensed Matter Physics: A Very Short Introduction