The geometries of NH n (CH 3 ) 3Àn , NH n (CH 3 ) þ 3Àn (n ¼ 0-3) and their corresponding coupling complexes [NH n (CH 3 ) 3Àn Á Á ÁNH n (CH 3 ) þ 3Àn (n ¼ 0-3)] were determined using density functional theory (DFT) and ab initio methods at the 6-311þG* basis set level, and the relative stability is predicted to decrease in the orderfor the four stable encounter complexes. The inapplicability is also discussed for the DFT methods in predicting the dissociation energy curves especially with long contact distance in which the DFT methods give abnormal behavior for the dissociation of the complexes due to the 'inverse symmetry breaking' problem. The contact distance dependences of the activation energy, the coupling matrix element and the electron-transfer (ET) rate were determined with the MP2/6-311þG* or MP2/6-31G* method. The results show that ET reactions occur chiefly only over a small range of contact distances where the favorable ones are 1.93 coupling systems. The most optimum contact distances for the above four coupling systems are 2.4, 3.0, 2.7 and 3.7 Å where the corresponding maximum ET rates are 2.53 Â 10 5 s À1 (H 3 , respectively, and the corresponding maximum ET rate order is H 3 NÁ Á ÁNH þ 3 > CH 3 H 2 NÁ Á Á NH 2 CH þ 3 > (CH 3 ) 2 HNÁ Á ÁNH(CH 3 ) þ 2 > (CH 3 ) 3 N Á Á ÁN(CH 3 ) þ 3 . Increasing substituents on the active N centers may significantly change the ET rate and other kinetic parameters.