1-Hydroxyethylidene-1,1-diphosphonic acid (HEDP) solutions in the pH range 0.98-13.00 were analysed using FT-Raman spectroscopy and 31 P and 23 Na NMR spectroscopy. Vibrational bands for different protonated species were observed in the Raman spectra, whereas only a single NMR signal that shifted with pH was observed for all samples over the entire pH range. No significant shift in the 23 Na NMR signal was observed, confirming that formation of Na + (aq) complexes did not take place; hence, no interference with the different protonated forms of HEDP occurred. Vibrational bands were assigned using density functional theory(DFT)-calculated spectra of the most likely conformers in solution. Multivariate curve resolution was performed on the Raman spectra in the region containing the PO stretching vibrations to determine the number of protonated species formed over the entire pH range. Chemometric analysis compares very favourably with the experimental species distribution diagram which was generated using the reported log K H values.