An efficient numerical algorithm is employed which enables one to perf
orm multidimensional integrations of complicated integrands. Temperatu
re dependence of the second virial coefficient for water is reproduced
using the Matsuoka-Clementi-Yoshimine intermolecular water-water pote
ntial. Metastable states are shown to occupy significant domain in the
water dimer phase space.