In this work, we use the Born-Oppenheimer approximation, where the potential between atoms can be approximated as a function of distance between the two nuclei, to study the four-quark bound states. By this approximation, Heitler and London calculated the spectrum of the hydrogen molecule, which includes two protons (heavy) and two electrons (light). Generally, the observed exotic mesons Zb(10610), Zb(10650), Zc(3900) and Zc(4020) (Zc(4025)) may be molecular states made of two physical mesons and/or diquark-anti-diquark structures. Analogous to the Heitler-London method for calculating the mass of the hydrogen molecule, we investigate whether there exist energy minima for these two structures. Contrary to the hydrogen molecule case where on...