A structure optimization of the entire molecule as well as of each fragment, for example by DFT, is followed by a relaxed surface scan along the dissociating bond. Through these data points, a morse function ∆E f it (r) is fitted, to gain a continuous function, describing the electronic energy of dissociation.