Finite-temperature Wigner phase-space sampling and temperature eﬀects on the excited-state dynamics of 2-nitronaphthalene †

The concept of finite temperature Wigner phase-space sampling allowing the population of vibrationally excited states is introduced and employed to study temperature eﬀects on the absorption spectrum of 2-nitronaphtalene (2NN) and its relaxation dynamics. It is found that, despite the fact that the general deactivation mechanism of 2NN after light irradiation does not change with increasing temperature, i.e. , after excitation to the singlet manifold, 2NN deactivates via internal conversion in less than 100 fs and then undergoes intersystem crossing to the triplet manifold before decaying to the lowest triplet state in less than 150 fs, the intersystem crossing rate increases at higher temperatures. This increase is attributed to the appearance of a new deactivation pathway, which is not operative at lower temperatures. The present example illustrates that appropriate initial conditions beyond the idealized temperature of 0 K are indispensable to obtain reliable excited state mechanisms, which can be then related to experimental conditions.


Introduction
The Wigner function 1 W is the central quantity in the phase-space formulation of quantum mechanics, 2 an alternative to the more commonly known matrix/operator formulation by Heisenberg and Schrödinger or Feynman's path integral formulation. The Wigner function is a representation of the density matrix of the system, and for a quantum system with f degrees of freedom described by the wave function C(q), it reads (1) where q and p denote the coordinates and momenta, respectively, and s is a dummy coordinate variable defining the integration. The Wigner function is used in the simulation of chemical dynamics, e.g., for obtaining an ensemble of molecules with different internal coordinates and momenta as initial conditions in trajectory-based methods. 3 However, one can also perform dynamics simulations completely in the Wigner representation of the density matrix. 4 The Wigner function has the properties ð dpdqWðq; pÞ ¼ 1 (2) ð dpWðq; pÞ ¼ CðqÞ j j 2 ¼ rðqÞ (3) ð dqWðq; pÞ ¼ CðpÞ j j 2 ¼ sðpÞ (4) i.e., it is normalized and upon projection in the coordinate and momentum space, it yields the coordinate and momentum densities r(q) and s(p), respectively, which are appropriate probability distributions being positive semi-definite. In contrast, the Wigner function W itself is often labeled as a ''quasi-probability'' distribution, as it can assume negative values in some regions of the phase space. This feature, however, does not prohibit one from employing the Wigner function in physical problems as has been discussed, e.g., by Feynman. 5,6 Furthermore, for simple cases such as the motion of an atom trapped in a harmonic potential, one can also reconstruct the Wigner function from experimental measurements. 7 This creates a function for the description of the phase space of the system with both positive and negative values. The latter refers to the phase space coordinates that are unavailable to the system in experiment, much like the well-known nodal planes of wave functions. The outcome of a chemical dynamics simulation in general depends strongly on an adequate sampling of initial conditions. 8 In practice, it is common to perform Wigner sampling in the harmonic approximation at zero temperature, 9,10 where one describes a molecule through a sum of uncoupled harmonic oscillators that are restricted to their vibrational ground state. Certainly, both approximations -the harmonic potentials and the zero-temperature formalism -are simplifications because realistic systems are anharmonic and possess a finite temperature. It is, therefore, interesting to study the effects and limitations of such approximations in chemical dynamics simulations.
Investigating the effects of anharmonic contributions to the potentials is a rather difficult task as it requires knowledge of the exact potential. Finite temperature effects, in contrast, can easily be incorporated into the harmonic approximation by allowing population of the vibrationally excited states -as was done by Mitrić, Bonačic-Koutecký, and co-workers, e.g., in ref. 11-17. Barbatti and co-workers 18,19 as well as Du and Lan 20 also mentioned the possibility of including temperature in Wigner sampling in dynamics. Given that initial conditions at different temperatures are available, it is possible to perform simulations for each of the initial conditions generated at different temperatures. Alternatively, it is possible to use the so-called importance sampling approach of Kossoski and Barbatti,21 in which dynamics simulations are explicitly run only at one temperature and then the properties of interest are extrapolated to other temperatures by weighting their probability distribution functions. However, despite its low computational cost, the accuracy of importance sampling can be questionable when the investigated properties are strongly temperature dependent or when they present only negligible values or when they are not present at all at the sampled temperature. In such situations, the explicit simulation of the system at different temperatures provides more accurate results.
Introducing finite temperatures into Wigner sampling can be realized in two approaches. In the first -more general -approach, one can write the temperature-dependent Wigner function W(q,p,T) as a sum of (temperature-independent) Wigner functions W[j i ](q,p) of the individual vibrational states j i weighted by their temperature-dependent populations P i , i.e., Using this formula, one can then step-wise sample the vibrational state and subsequently the coordinates and momenta of the system. This is explained in more detail in Section S1 of the ESI. † For a system described in the harmonic approximation, one can also sample directly from an analytical expression that Wigner et al.
found for a canonical ensemble of harmonic oscillators, i.e., 22 Both approaches present different conceptual advantages. Sampling via eqn (6) avoids dealing with unfamiliar negative quasi-probabilities as W (q,p,T) Z 0 8 (q,p). However, sampling using eqn (6) comes with a loss of information about the vibrational state j n that the system occupies. Although a quantum system can be in a superposition of reference states, this superposition should collapse to a single state once the corresponding property of the system is measured, i.e., once the system is sampled. This information is retained when sampling via eqn (5) is done instead. The discussion of which advantage is favorable when is, despite interesting, beyond the scope of the present work.
In this work, we systematically investigate the effects of Wigner sampling at finite temperature on the excited-state dynamics of the organic chromophore 2-nitronaphthalene (2NN). The deactivation mechanism of 2NN has been previously established by us employing nonadiabatic dynamics simulations, 23 and it is summarized in Fig. 1. Initially, the system is excited to a 1 pp* state of charge-transfer (CT) character [that we label S CT (pp*)], from which it mainly undergoes internal conversion with a time constant of t S B 80 fs to a locally excited 1 np* state [S LE (np*)]. Excited-state population is also transferred from the singlet to the excited triplet manifold via intersystem crossing (ISC) with a time constant of t ISC B 0.7 ps, before the system finally relaxes to the lowest triplet state T 1 with a time constant of t T B 150 fs. The ISC process proceeds via two pathways: a minor fraction (9%) of ISC occurs from the initially excited S CT (pp*) state to a locally excited 3 np* state [T LE (np*)], while the major fraction (91%) of ISC connects the S LE (np*) state to a locally excited 3 pp* state [T LE (pp*)]. The predominance of the S LE (np*) -T LE (pp*) pathway over the S CT (pp*) -T LE (np*) pathway is explained by three key factors. First, the singlet ISC donor state of the major pathway, S LE (np*), is rapidly populated and longer-lived than the other singlet excited states. Second, once the S LE (np*) state is populated, the nuclear configurations of 2NN are very close to the configurations that allow ISC and, thus, only small geometrical arrangements are necessary to transfer to the triplet manifold. Third, at the ISC geometries, the S LE (np*) and T LE (pp*) states share very similar electronic distributions and the transition between both states requires only the change of the angular momentum of one electron -allowing the remaining electronic distribution to remain static.
Previous nonadiabatic dynamics simulations 23 were started from a Wigner ensemble at T = 300 K. For this work, we have additionally performed nonadiabatic dynamics simulations using Wigner ensembles at T = 0 and T = 500 K. As the following results will demonstrate, temperature has a considerable effect on the excited-state dynamics. This is manifested by an increasing rate of ISC at higher temperatures as well as the opening of an additional ISC pathway at higher temperature that explains the increase of ISC rate.

Computational details
The excited-state dynamics of 2NN was simulated as in ref. 23, using the surface hopping including arbitrary couplings (SHARC) approach 24-27 with energies, gradients and couplings calculated with the PBE0 28,29 /DZP 30 level of theory from the ADF2016 program package. 31 The initial conditions required for the nonadiabatic simulations were obtained from Wigner distributions sampled using eqn (5) described above at T = 0, 300, and 500 K, which we shall denote as WIGNER-TX with X = 0, 300, or 500, respectively. For each of the 1000 geometries per ensemble, the ten lowest excited states were calculated and used to simulate an absorption spectrum by convoluting the resulting stick spectra with the Gaussian function using a full-width at half maximum of 0.1 eV (see Fig. 2).
From each WIGNER-TX ensemble, at least 60 geometries were selected stochastically 32 and trajectories were propagated for 500 fs using a nuclear time step of 0.5 fs. The trajectories were started at the corresponding bright excited state located within a 0.5 eV energy range around the maximum of the calculated absorption band (gray areas in Fig. 2), resulting in trajectories starting from the S 1 -S 4 states depending on the initial conditions. Based on the excited states found in the excitation energy range and below, different numbers of states were included in the simulations. Due to convergence problems and hops to inactive states during the simulation time, some trajectories were excluded from the analysis. The total number of trajectories, the number of all initial states, the number of excluded trajectories, and the number of states included in the simulation are listed in Table 1 for all three ensembles.  states for each ensemble. These densities already reveal a temperature effect, as the excited states spread energetically with increasing temperature. This can be seen, e.g., for the S 4 state (middle panel in Fig. 2): its density lies outside of the excitation energy range for T = 0 K, but it reaches inside the excitation energy range for T = 300 K, and reaches even further for T = 500 K. This effect also has the consequence that for the Wigner-T500 ensemble, a number of trajectories were started also from the S 3 and S 4 states, whereas for all other ensembles, all trajectories were started in either the S 1 or S 2 state. Notably, while the density of states is broadened for the Wigner ensembles with increasing temperature, the widths of the absorption bands are hardly affected, with the exception of the lowest-energy absorption band, which is broadened at higher temperature. Notably, however, increasing the temperature leads to a red-shift of the higherenergy absorption bands of 0.1-0.2 eV, though the lowest-energy absorption band remains approximately at the same energy. This shift agrees with a previously suggested general red-shift, which is induced by the vibrational motion of the molecule. 33 For reference, we have presented a table with energies, oscillator strengths, and characters of the low-lying electronic states at the minimum-energy Franck-Condon (FC) geometry in Section S3 in the ESI. †

Time evolution of the excited-state populations
The time evolution of the electronic state populations obtained using the three different ensembles as initial conditions for the nonadiabatic dynamics is shown in Fig. 3. The left panels display the populations of the molecular Coulomb Hamiltonian (MCH) states, 26 i.e., states that are characterized by their spin multiplicity and ordered according to their energy. In this description, the simulations start either on the S 1 or on a higher excited singlet state, S n , with n = 2 for the Wigner-T0/ T300 ensembles and n = 2, 3 or 4 for the Wigner-T500 ensemble. As in the Wigner-T300 ensemble, 23 population is initially transferred from the S n state(s) to S 1 before ISC takes place to the higher-lying triplet states T n (n = 2-6) and subsequently decays to the T 1 state. This general mechanism is summarized as and accounts for almost all processes observed during the first 500 fs. Relaxation from the T 1 state to the S 0 state takes place on a larger time scale and it is not accounted for by our short simulation time.
Alternatively, the dynamics occurring in the singlet manifold can also be described in terms of the so-called spectroscopic states, 26 i.e., states of a particular electronic character. The corresponding time evolution of the state populations is shown in the right-hand side of Fig. 3. Within the spectroscopic representation, the dynamics in the singlet manifold is best described by i.e., after initial excitation to the (bright) S CT (pp*) state, the system relaxes to the S LE (np*) state, regardless of the energetic order of the states. The CT or LE character of the states is discriminated using the dipole moment.   (1) 43 (1) 0 (0) 0 (0) 3 + 2 6 + 2 Wigner-T300 105 (6) 51 (2) 54 (4) 0 (0) 0 (0) 3 + 2 6 + 2 Wigner-T500 60 (2) 12 (0) 36 (2) 1 (0) 11 (0)  states (purple curves in Fig. 3), and it fails for larger simulation times. However, since the S CT (pp*) -S LE (np*) transition takes place on a time scale smaller than ISC (see below), the characterization of the states is valid. The S CT (pp*) -S LE (np*) internal conversion for the Wigner-T300 ensemble was explained in detail in ref. 23. This process also dominates the beginning of the dynamics simulations of the other ensembles, including the Wigner-T500 ensemble, which involves trajectories starting from states S 1 up to S 4 and, thus, an additional charge-transfer S CT (pp*) state. This conclusion is drawn from the analysis of the transition density calculated with the TheoDORE program package [34][35][36][37] at the initially excited states in the trajectories and at the lowest-excited states at the minimumenergy geometry [FC geometry]. Fig. 4 shows the atomic electronhole difference populations of the excited states that are present in the transition with respect to the ground state electron distribution.
As can be seen, the 47 trajectories starting from the S 1 , S 2 and S 3 states in the Wigner-T500 ensemble exhibit an electronhole (e/h) difference population that is similar to that of the S 1 state at the FC geometry (S 1 @FC), which is a bright CT 1 pp* state. For the 11 trajectories starting from the S 4 state, the e/h difference populations resemble more those of the S 3 @FC state -especially the atomic hole population. This state is also a bright singlet CT pp* state [S CT 0 (pp*)]. As the electron population at the nitro group in the hS 4 i states is larger than in the S 3 @FC state, this suggests that there is an admixture of S 1 @FC state character in the hS 4 i states. More importantly for our interpretation of the excited-state dynamics results is to realize that all trajectories start in the CT 1 pp* states (S 1 @FC and S 3 @FC) and not in the LE 1 np* states (S 2 @FC and S 4 @FC).

Temperature effects on intersystem crossing
Having explained the time evolution of the excited-state populations, (Fig. 3) we now proceed to discuss the temperature effects on the dynamics. Table 2 collects the time constants fitted according to the mechanisms proposed in eqn (7) and (8).
As stated before, the initial dynamics is governed by a relaxation within the singlet manifold, where, in terms of spectroscopic states, population transfers from the bright S CT (pp*) state to the S LE (np*) state. This reaction occurs for all three ensembles within similar time constants, between 81 AE 10 fs (Wigner-T300) and 95 AE 12 fs (Wigner-T0). Subsequently, the system undergoes ISC to high-lying triplet states T n , before it relaxes to the T 1 state. Also, the time constants t T describing the relaxation dynamics within the triplet manifolds are similar at the different temperatures (t T = 110-150 fs). Thus, the effect of temperature on the relaxation within both the singlet and the triplet manifolds is quite small (within the error margin), and not worth further discussion.
In contrast, the rate of the ISC process from the singlet to the triplet states shows a clear temperature dependence. The rate of the reaction increases notably with the temperature as t ISC (T 0 ) E t ISC (T 300 ) 4 t ISC (T 500 ). A complete understanding of the reasons behind this temperature dependency of ISC requires an analysis of the electronic structure during the ISC events; as will be shown next, the ISC rate increases at higher temperatures due to an additional ISC pathway that becomes operative at higher temperatures.
Two ISC pathways, which we denoted as major [S LE (np*) -T LE (pp*)] and minor [S CT (pp*) -T LE (np*)], were reported in ref. 23. When characterizing all ISC events only in terms of either of these two pathways, we find in the present work that   (7) and (8) the relative contribution of the minor channel g minor increases at higher temperatures, from 5% at T = 0 K to 20% at T = 500 K, see Table 3. Table 3 also reports the average spin-orbit couplings (SOCs) at the singlet-to-triplet ISC hopping geometries and we find a considerable increase of the SOCs of the minor pathway when going from the Wigner-T0/Wigner-T300 ensembles (hSOCi = 7.9-8.5 cm À1 ) to the Wigner-T500 ensemble (hSOCi = 13.8 cm À1 ). Again, this is due to the additional ISC pathway involving larger SOC, which is a variation of the S CT (pp*) -S LE (np*) minor pathway and becomes effective at higher temperatures (see below). This additional pathway does not exist at T = 0 K or T = 300 K. Fig. 5 shows the natural transition orbitals (NTOs) describing the excited states involved in all the three ISC pathways in the moment of the singlet-triplet transitions. The newly identified ISC pathway is labelled as a delocalized excitation (DE) pathway as it is characterized by a transition from an S CT (rp*) to a T CT (rp*) state. Both singlet and triplet states possess the same p* electron NTO, but differ in their hole NTO that we shall refer to as being of r-type, a hybrid between p-type and n-type orbitals. Both states participating in this ISC transition are CT excited states. However, as the direction of CT is the same in both states, the resulting transition rr between S CT (rp*) and T CT (rp*) does not possess any CT character. Instead, as it is delocalized over the whole chromophore, we denote the pathway containing this transition as DE. Similarly, it is possible to re-label the major and minor ISC pathways by the characters of their ISC transitions as LE and CT pathways, respectively. This is because the pn transition in the major ISC pathway is localized at the nitro group and does not have CT character, while the np transition in the minor ISC pathway is of CT character, shifting electron density from the nitro group to the aromatic ring.
The r NTOs in the DE pathway are a mixture of an n orbital at the nitro group and a p orbital delocalized over the aromatic ring system, similar to the orbitals involved in the CT pathway (Fig. 5). As there is only a small variation of properties between the transitions in the DE and CT pathways -such as the exciton sizes of the participating excited states -it is very difficult to unambiguously assign an ISC hop belonging either to the DE or CT pathway, unless one manually inspects every ISC hopping event, which is an extremely tedious analysis.
Although we did not manually inspect every ISC hopping event, we inspected a small sample of ISC hops, and the analysis suggest that the SOCs between states of the DE pathway are of the order  of 10-20 cm À1 , i.e., larger than the SOCs of the CT pathway (o10 cm À1 ) but smaller than the SOCs of the LE pathway (B40 cm À1 ). The larger SOCs of the DE pathway explain the larger SOCs reported for the minor pathway of the Wigner-T500 ensemble in Table 3, as the minor pathway includes both the CT and DE pathways in this preliminary analysis. The size of the SOCs in the different pathways seems to be correlated with the amount of charge relocalization that is occurring during the ISC hop. The SOCs are largest for the LE pathway, which requires only a change of the angular momentum of one electron and no further charge migration because both NTOs are located at the nitro group (see Fig. 5). In contrast, the SOCs are smallest for the CT pathway, where both NTOs are localized in different regions of the molecule. This fact was already observed in ref. 23. In the present work, we find SOCs of an intermediate size for the DE pathway, as the NTOs involved are delocalized over the complete molecule but show an important overlap, thus limiting the amount of charge redistribution necessary to connect both states. We now examine how the different ISC pathways contribute to the overall ISC rate over time. For this, we plot the time evolution of the average number of ISC hops N Hops (t) during time intervals of Dt = 50 fs for trajectories in singlet states where we differentiate between ISC hops of the LE pathway (previously: 23 major pathway) and ISC hops of either the DE or CT pathway (previously: 23 minor pathway) in Fig. 6. For example, if N Hops (t) = 0.2 at t = 25 fs, this means that on average, 20% of all trajectories that are in a singlet state within the time window of 0-50 fs hop to a triplet state in this time window. Note that these trajectories can subsequently hop back to the singlet manifold.  The time evolution of N Hops can be used as a measurement of the rate of ISC via either the LE pathway or the combination of the DE and CT pathways. Besides the N Hops for the three different ensembles (Fig. 6(a)-(c)), we show the fit functions over time in panel (d).
As can be seen, the number of ISC hops that occur through the LE pathway increases for the first 200-300 fs before it becomes nearly stationary in all three ensembles. This is because the population of the singlet donor state of the LE pathway, the S LE (np*) state, is initially zero and gets populated from the initially excited S CT (pp*) state, and then ISC can occur. In contrast, the majority of ISC hops through either the DE or CT pathway occur at the beginning of the dynamics. It can be noted that these hops are largest for the Wigner-T500 ensemble, both at the beginning as well as throughout the full simulation time, due to the opening of the new DE pathway. This explains the temperature-dependent increase of the ISC rate at higher temperatures as most of the ISC events for the DE + CT pathway occur at short simulation times.

Nuclear motion
In this section, we analyze the nuclear motion involved in the deactivation dynamics of 2NN using a normal mode analysis, as was done in ref. 23. In terms of internal coordinates, we find that the most important motions affect mainly the angle g ONO between the nitro-group atoms, the distances r NO between the nitrogen and oxygen atoms, and the distance r NC a between the nitrogen and its neighboring carbon atom. This was true for the Wigner-T300 ensemble 23 and it also applies to the Wigner-T0 and Wigner-T500 ensembles.
For all three ensembles, the most active normal modes are displayed in Fig. 7(a)-(c), where we show their displacements with respect to the optimized FC geometry at the initial geometries, at the ISC hopping geometries of the LE (''major'') and the DE + CT (''minor'') pathways, and along the dynamics of the trajectories in the 1 np* and 1 pp*/ 1 rp* states (singlet donor states in the ISC), respectively. Due to the difficulties involved in discriminating between the DE and CT pathways, we combine both sets of ISC geometries, and likewise, only distinguish between two types of singlet states, the 1 np* state (defined with m o 8 D) and the 1 pp*/ 1 rp* state (defined with m 4 8 D). The corresponding displacement vectors of the normal modes with respect to the FC geometry are presented in Fig. 7(d).
As expected, the average displacements of the initial geometries (gray histograms in Fig. 7(a)-(c)) are small because the initial geometries are still around the FC geometry. In contrast, the displacements of the ISC geometries (violet and orange) indicate the motions required for the system to undergo ISC, and the displacements of the geometries of the trajectories in the singlet states (blue and green) show in which region of the PES the singlet trajectories spend most of the time during the simulation. Note that for an efficient ISC, it is beneficial that the displacement of the singlet states is similar to that of the ISC geometries, meaning that the system spends time in a region of the PES where ISC is possible.
As can be seen, the normal mode analyses of the Wigner-T0 and Wigner-T300 ensembles are very similar, but different from that of the Wigner-T500 ensemble, particularly at the DE + CT pathway (orange) and its donor 1 pp*/ 1 rp* state (green). These differences can be best understood following the time evolution of the internal coordinates, g ONO , r NO , and r NC a , associated with the most important normal modes, see Fig. 8. This figure also shows the averages of these coordinates at the initial geometries as well as at the ISC hopping geometries of the LE and DE-CT pathways. As can be seen, the average of the latter DE + CT pathway (orange dashed line) varies substantially with temperature. For example, the g ONO angle for the DE + CT pathway is 128.51 (129.01) at T = 0 K (300 K) but 124.71 at T = 500 K. Such a smaller g ONO angle is then more often accessed by the trajectories in the 1 pp* state (green solid line), which are thus more frequently in regions where they can undergo ISC via the DE pathway. This situation is also true for the average r NO and r NC a of the DE + CT pathway at T = 500 K, so that they best overlap with the oscillations of the 1 pp*/ 1 rp* trajectories.
In contrast to the average internal coordinates for the DE + CT pathway, the average internal coordinates of the LE pathway vary less with temperature. This is expected as the electronic states involved in the LE pathway are the same at all temperatures, and, thus, ISC takes place in the same regions of the PES. Only r NC a (violet line) increases notably when going from T = 300 K (1.37 Å) to T = 500 K (1.39 Å). However, accompanying this increase, we find that the oscillations of r NC a of the 1 np* trajectories also increase, thus, probably cancelling any effect that the increase could have on the LE pathway.
Thus, the better overlap of the coordinates of the 1 pp* trajectories with those of the ISC hopping geometries of the DE + CT pathway at higher temperature contributes to the temperature-dependent increase of the ISC rate (see Section 5). Due to the increased energy at T = 500 K, the trajectories in the 1 pp* state are able to visit different regions of the PES where ISC is more accessible than when having less energy. This featurein addition to the larger spin-orbit couplings due to the DE part of the DE + CT pathway -is responsible for the increase of the ISC rate at higher temperatures.

Conclusions
The present work explores the concept of finite-temperature Wigner phase-space sampling and investigates the effects that finite temperature can have on excited-state dynamics simulations. In Wigner sampling, a finite temperature is realized by allowing the population of vibrationally excited states, thereby increasing the vibrational energy by a thermal-energy contribution. Such populations also result in probability distributions of coordinates and momenta that are different from the vibrational ground-state distributions from a hypothetical zero-temperature case. The larger the phase space sampled and the more additional thermal energy available within the finite temperature Wigner sampling, the larger the impact on the excited-state dynamics of certain systems. Despite the trivial nature of this statement, finite temperature effects are rarely taken into account. Therefore, here, we have investigated temperature effects (T = 0, 300 and 500 K) on the absorption spectrum and the excited-state dynamics of the organic chromophore 2NN upon excitation to the lowestenergy absorption band, finding noticeable changes.
In the absorption spectrum of 2NN, increasing the temperature leads to a broadening of the absorption bands. As a consequence, the individual excited states expand over a larger range of energy at T = 500 K, so that setting-up the excited state dynamics requires trajectories starting from a larger number of states than at T = 0 K and T = 300 K, despite using the same excitation-energy range.
The general excited-state decay mechanism of 2NN is found to be the same, regardless of the temperature, i.e., after initial relaxation dynamics in the singlet manifold on a B100 fs time scale, the system undergoes a sub-ps ISC to higher-lying triplet states before decaying to the lowest T 1 state. However, the rate of ISC from the singlet to triplet states increases at higher temperature, especially when going from T = 300 K to T = 500 K. We attributed this increase to the presence of a new ISC channel that involves delocalized singlet-to-triplet transitions. This channel becomes operational only at high temperature, in addition to the two other previously identified pathways involving locally excited and charge transfer transitions.
Clearly, finite temperature Wigner sampling is a simple and straightforward extension of the zero-temperature formalism that yields an improved and more realistic phase-space sampling. As demonstrated here, the effects of the sampling size of the phase space and the energy available can have a non-negligible impact on excited-state dynamics simulations and should be taken into account. Thus, this procedure is well suited to become the new standard to generate vibrational ensembles for the calculation of absorption spectra or to generate initial conditions for excited-state dynamics when anharmonic effects are negligible.

Conflicts of interest
There are no conflicts to declare.