We used surface hopping for the nonadiabatic dynamics to build the nuclear ensembles. An issue not directly related to the spectrum simulations but crucial to the results' accuracy is the velocity adjustment after hopping between states. The conventional way of making such adjustments yields artifacts when applied to large molecules, such as pyrene. We discuss in this subsection how this problem was solved.