In quantum chemistry, the potential energy surfaces are obtaine within the adiabatic or Born-Oppenheimer approximation. This corresponds to a representation of the molecular wave function where the variables corresponding to the molecular geometry and the electronic degrees of freedom are separated. The non separable terms are due to the nuclear kinetic energy terms in the molecular Hamiltonian and are said to couple the potential energy surfaces. In the neighbourhood of an avoided crossing or conical intersection, these terms cannot be neglected. One therefore usually performs one unitary transformation from the adiabatic representation to the so-called diabatic representation in which the nuclear kinetic energy operator is diagonal. In this representation, the coupling is due to the electronic energy and is a scalar quantity which is much more easy to estimate numerically.