Phosphorescence can also have significant contributions from vibronic couplings arising from symmetry reduction during vibrational motion. Different methods are available to incorporate vibronic couplings into the simulations, such as the nuclear ensemble approach (NEA) , direct vibronic coupling , adiabatic Hessian , and vertical gradient . Torres and co-authors recently compared them, and the interested reader may know more about the pros and cons of each one in their paper . We particularly favor NEA, which we have been developing and applying for many years to simulate different spectrum types and generate initial conditions for dynamics .