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:

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:

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

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:

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:

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:

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:
![{\displaystyle \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)}](//restbase.wikitolearn.org/en.wikitolearn.org/v1/media/math/render/svg/cd34c6f457772f56a52d1f204d9db8405aca9ba6)
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:

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

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.