A protocol that uses enhanced QM/MM method to investigate the isotopic effect on the double proton transfer process in porphycene is presented here.