All calculations were performed using Gaussian 09 program
package [18]. The geometry of each compound (B6AOH), radical
(B6AO), radical cation (B6AOH+) or anion (B6AO) was optimized
using Density Functional Theory (DFT) method with B3LYP [19,20],
M05-2X [21] and M06-2X [22] functionals without any constraints
(energy cut-off of 105 kJ mol1, final RMS energy gradient under
0.01 kJ mol1 Å1). Calculations were performed in 6-311++G⁄⁄
basis set [23]. Solvent (benzene and water) contribution to the