The value of the harmonic restraint has been chosen to maintain the integrity of the cluster without creating artefacts from the pulling force. The onset distance of the flat-bottom potential was varied according to the size of the clusters so that the density of the cluster was close to the experimental density of the (bulk) liquid.