The goal of the protocol presented here is to generate and sample trajectories of configurations of liquid water molecules around catalytic species on a flat transition metal surface. The sampled configurations can be used as starting structures in quantum mechanics-based methods.