First hyperpolarizabilities of a large number of push−pull substituted conjugated systems with heteroaromatic spacers have been calculated. The static, nonresonant components were computed at the ab initio level (4-31G basis) using the coupled perturbed Hartree−Fock approach and at the AM1 level employing the finite field method. Sum-over-states procedure has also been used with the AM1/CI method to compute β0 and β at an excitation energy of 1.17 eV. The computed β values at the various levels are reasonably similar and exhibit the same trends. The largest values are obtained with a donor on pyrrole and an acceptor on thiophene or thiazole. The variations do not always inversely follow the order of delocalization energies of the heterocyclic rings. The trends in the dipole moment changes and transition energies between the ground and first excited charge-transfer state primarily determine the variations in the computed β values.