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 N-inversion barriers and phenyl rotations about the N-C bonds. The state-specific corrections developed in Ref. 61 were not used, since they were derived specifically for monomeric adiabatic states.