The electronic states of a singly ionized on-axis double-donor complex (D2+) confined in two identical vertically coupled, axially symmetrical quantum dots in a threading magnetic field are calculated. The solutions of the Schrödinger equation are obtained by a variational separation of variables in the adiabatic limit. Numerical results are shown for bonding and antibonding lowest-lying artificial molecule states corresponding to different quantum dot morphologies, dimensions, separation between them, thicknesses of the wetting layers, and magnetic field strength.