AbstractA finite‐field method for the calculation of polarizabilities and hyperpolarizabilities is developed based on both an energy expansion and a dipole moment expansion. This procedure is implemented in the MOPAC semiempirical program. Values and components of the dipole moment (μ), polarizability (α), first hyperpolarizability (β), and second hyperpolarizability (γ) are calculated as an extension of the usual MOPAC run. Applications to benzene and substituted benzenes are shown as test cases utilizing both MNDO and AM1 Hamilto