We provide the exact
analytical form of diatomic molecular orbitals,
as given by the solutions of a single-electron diatomic molecule with
arbitrary nuclear charges, using our recently developed method for
solving Schrödinger equations. We claim that the best representation
of the wave function is a factorized form including a power prefactor,
an exponentially decaying term, a modulator function on the exponential,
and additional factors accounting for nodal surfaces and the magnetic
quantum number. Applying our method, we have identified unexpected
extreme points along the potential energy curves, hence revealing
the limitations of the well-known concepts of bonding and antibonding.