The first full-dimensional accurate potential energy surface was developed for the CO + H2O system based on ca. 102 000 points calculated at the CCSD(T)-F12a/AVTZ level using a permutation invariant polynomial-neural network (PIP-NN) method.