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.