Electron-phonon interactions in the monoanions of B, N-substituted acenes such as B3N3F6 (1f) and B5N5F8 (2f) are studied, and compared with those in the monoanions of B3N3H6 (1h) and B5N5H8 (2h), and B3N3D6 (1d) and B5N5D8 (2d). The low frequency modes around 500 cm-1 as well as the frequency modes higher than 1000 cm-1 strongly couple to the lowest unoccupied molecular orbitals (LUMO) in 1f and 2f. The total electron-phonon coupling constants (l(LUMO)) are estimated to be 2.710 and 2.054 eV for 1f and 2f, respectively, and those are estimated to be 0.342 and 0.235 eV for 1d and 2d, respectively, while those were estimated to be 0.340 and 0.237 eV for 1h and 2h, respectively. That is, the l(LUMO) value increases much more significantly by H-F substitution than by H-D substitution in B, N-substituted acenes. The larger displacements of B and N atoms in the vibronic active modes in 1f and 2f than those in 1d and 2d due to larger atomic mass of fluorine than that of deuterium, and the phase patterns difference between the LUMO in 1f and 2f, in which the atomic orbitals between N and its neighboring F atoms form strong sigma-antibonding interactions, and that in 1d and 2d, in which the atomic orbitals between two neighboring B and N atoms form weak pi-bonding and pi-antibonding interactions, are the main reason why the l(LUMO) value increases much more significantly by H-F substitution than by H-D substitution. The reorganization energies between the neutral molecules and the corresponding monoanions are estimated to be 0.122, 0.063, 0.733, and 0.830 eV for 1h, 2h, 1f, and 2f, respectively. Therefore, the estimated reorganization energies between the neutral molecules and the corresponding monoanions for 1f and 2f are much larger than those for 1h and 2h. (C) 2004 American Institute of Physics.