“…Using these considerations, a total of 10 rate constants (with k d,i,j,k,l,m,n and k ib,i,j,k , which could not have been obtained from literature [27][28][29][30][31][32]) for all potential substituents had to be determined, that is k iha,i,j,k,l,m,n , k iar,i,j,k,l,m,n,o , k it1,i,j,k,l,m,n , k it2,i,j,k,l,m,n , k ha,i,j,k,l,m,n , k ar,i,j,k,l,m,n,o , k t1,i,j,k,l,m,n , k t2,i,j,k,l,m,n , k b,i,j,k , and k e,i,j,k,l,m,n . The rate constants and the associated activation energies and pre-exponential factors were determined utilizing molecular modeling [22,23,33]. This was achieved by the semi-empirical Parameterized Model number 3 method calculations, using unrestricted Hartree-Fock spin pairing and self-consistent field controls of 0.001 convergence limit and 1000 iterations limit.…”