# Born-Oppenheimer approximation

The aim of this section is to study the quantum description of the molecular bond. We are going to limit the analysis to diatomic molecules. First of all we consider a molucule composed by two atomic nuclei A and B, whose positions are defined by vectors and . The electrons positions are instead identified by . Both nuclei and electrons positions are refered to the center of mass of the two nuclei. The Hamiltonian we can write for a system like this is:

<dmath>\hat{H}=\hat{T_{N}}+\hat{T_{e}}+U_{eN}+U_{ee}+U_{NN}</dmath>

where and are the kinetic energies associated respectively to the nuclei and the electrons and and represent the Coulomb interactions. Now it's possible to isolate from this Hamiltonian an electron Hamiltonian:

<dmath>\hat{H_e}=\hat{T_e}+U_{eN}+U_{ee}</dmath>

This is a function of the nuclei positions and we can obtain an eigenvalues equation for the electrons:

<dmath>\hat{H_e}\left(\left\lbrace \vec{R}_{i}\right\rbrace\right)\psi^{el}_{n}\left(\left\lbrace \vec{R}_{i}\right\rbrace;\left\lbrace\vec{x}_{i}\right\rbrace\right)= E^{el}_{n}\psi^{el}_{n}\left(\left\lbrace \vec{R}_{i}\right\rbrace;\left\lbrace\vec{x}_{i}\right\rbrace\right)</dmath>

It's possible to demonstrate that solving this equation for the electrons Hamiltonian leads to a solution for the original complete eigenvalue equation which has the form:

<dmath>\psi=\sum_{n}\chi_{n}\left(\left\lbrace\vec{R}_{i}\right\rbrace\right)\psi_{n}\left(\left\lbrace \vec{R}_{i}\right\rbrace;\left\lbrace\vec{x}_{i}\right\rbrace\right)</dmath>

where is a wave function depending only on the nuclei positions. Such an equation is difficult to solve, but it can be semplified through the **Born-Oppenheimer adiabatic approximation**. This model is based on the fact that the mass of the nucleus is much greater than the electron one (). Being so, we can state that electrons move way faster than nuclei. Hence we can assume that electrons follow the motion of the nuclei without changing the energy level they belong to. In other words, there is no energy exchange between nuclei and electrons. In this adiabatic approximation the solution becomes:

<dmath>\psi=\chi_{n}\left(\left\lbrace\vec{R}_{i}\right\rbrace\right)\psi_{n}\left(\left\lbrace \vec{R}_{i}\right\rbrace;\left\lbrace\vec{x}_{i}\right\rbrace\right)</dmath>

The sum over the quantum number disappears, because of the impossibility for an electron to jump to an excited state. Now we can use theese results to write an equation for the total Hamiltonian:

<dmath>\left(\hat{T_N}+U_{NN}+\hat{H_e}\right)\chi_{n}\left(\left\lbrace\vec{R}_{i}\right\rbrace\right)\psi^{el}_{n}\left(\left\lbrace\vec{R}_{i}\right\rbrace;\left\lbrace\vec{x}_{i}\right\rbrace\right)=E_{tot}\chi_{n}\left(\left\lbrace\vec{R}_{i}\right\rbrace\right)\psi^{el}_{n}\left(\left\lbrace\vec{R}_{i}\right\rbrace;\left\lbrace\vec{x}_{i}\right\rbrace\right)</dmath>

At this point we can observe that only acts on , which is actually an eigenstate of , so we can substitute in the equation above with the corresponding eigenvalue . Moreover, again thanks to the greatness of nuclei masses in respect to electrons ones, we can assume that the application of the kinetic term to the component of the total wave function gives a negligible contribution, then we can consider that acts only on . Relying on this last approximation, it's possibile to semplify the term in order to obtain an eigenvalues equation for the sole function:

<dmath>\left[\hat{T_N}+U_{NN}+E^{el}_n\left(\left\lbrace\vec{R}_{i}\right\rbrace\right)\right]\chi\left(\left\lbrace\vec{R}_{i}\right\rbrace\right)=E_{tot}\chi\left(\left\lbrace\vec{R}_{i}\right\rbrace\right)</dmath>

The last equation describes the nuclei dynamics and this means that the equations for and are now uncoupled. depends only on the distance between the nuclei:

<dmath>E_{tot}=E_{tot}\left(\left|\vec{R}_{A}-\vec{R}_{B}\right|\right)</dmath>

Now we can redefine the quantity considering that is just a moltiplicative operator:

<dmath>E^{el}\left(\left\lbrace\vec{R}_i\right\rbrace\right)=E^{el}_n\left(\left\lbrace\vec{R}_i\right\rbrace\right)+U_{NN}</dmath>

In other words, is the eigenvalue for the operator corresponding to the eigenfunction . It's possible to plot in respect to the nuclear distance :

In the graph, the value of the energy in correspondence of the axis is , where and are the energy associated to the idrogenoid Hamiltonian for nucleus and . It can be seen that there is an at which the energy is minimal. This is called **dissociation energy** and it's the energy that has to be given to the molecule in order to break the bond.