EN
We determine the states of a pair of bound magnetic polarons in a diluted magnetic semiconductor by taking into account both the molecular nature of the two-electron wave function in the Heitler-London approximation, as well as the thermodynamic fluctuations of magnetization due to localized 3d spins. We extend the theory of a single bound magnetic polaron to molecular case.