We have theoretically elucidated the origin of the isotope effect on phase transition temperature in hydrogen-bonded ferroelectric materials K 3 H(SO 4 ) 2 and K 3 D(SO 4 ) 2 , using the cluster models with multi-component molecular orbital (MC MO) method, which can determine both electronic and nuclear wave functions simultaneously. We have confirmed that the origin of the isotope effect is the localization of the nuclear wave function by deuterization. We have also demonstrated the isotope effect in effective potential energies, geometrical structures, and electronic population on hydrogen nuclei.Keywords: isotope effect on phase transition temperature, hydrogen-bonded ferroelectric material, geometrical isotope effect, K 3 H(SO 4 ) 2 and K 3 D(SO 4 ) 2 , multi-component molecular orbital method, Ubbelohde effect