Semi-empirical Hartree-Fock calculations for KNbO3 and KTaO3

As a first step in modeling the electronic structure of Perovskite-type ferroelectric mixed crystals K(Nb,Ta)O3, semiempirical calculations for pure KNbO3 and KTaO3 are performed with the intermediate neglect of the differential overlap (INDO) quantum chemical method. The calculations are mostly done for 40-atom supercells. The choice of the INDO parameters based on the comparison of results with ab initio and experimental data is discussed. INDO results for the equilibrium geometry and the (Gamma) -TO phonon frequencies are given. The results show that the accuracy of the INDO method is sufficient for reliably reproducing the energy differences on the order of 1 mRy (per formula unit) scale which are typical when dealing with ferroelectric instabilities.