Quantum justification for classical discussions of potential energy surfaces

 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.