We use the analytical potential for water determined by Partridge and Schwenke 51 to obtain SP energies as required by the TD-ADGA algorithm and then we fit these SPs to obtain the PES. We will use the full three-mode representation of the PES and in each of the seven MCs, the PES is fitted to a polynomial function in the displacement coordinates.