Structure and stability of dications was explored including the inactive (non-breaking) trajectories for 0.01, 1, 2.5 5, 10, 15 and 20 eV, in practice the last steps from BOMD. Geometry optimizations have been carried out at the B3LYP-D3/6-311++G(d,p) level of theory, as well as zero point energy correction were added. All static calculations were performed using the Gaussian16 program .