In all MD simulations we used the so-called added potential (added to each of the four molecules of the tetramer) which corrects the low AM1 inversion barriers and phenyl rotations about the N-C bonds. We note, however, that the state-specific corrections developed in Ref. 49 were not used, since they were derived specifically for monomeric adiabatic states.