We have performed high-level quantum mechanical calculation for multiple proton transfer in HOCl + HCl and $H_2O$ + $ClONO_2$ on water clusters, which can be used as a model of the reactions on ice surface in stratospheric clouds. Multiple proton transfer on ice surface plays crucial role in these reactions. The structures of the clusters with 0-3 water molecules and the transition state structures for the multiple proton transfer have been calculated. The energies and barrier heights of the proton transfer were calculated at various levels of theory including multi-coefficient correlated quantum mechanical methods (MCCM) that have recently been developed. The transition state structures and the predicted reaction mechanism depend very much on the level of theory. In particular, the HF level can not correctly predict the TS structure and barrier heights, so the electron correlation should be considered appropriately.