The photoinduced dynamics in fulvene upon excitation to S 1 and decay to S 0 as well as the relevant CIs have been studied in the literature quite abundantly . It has been found that the relaxation dynamics is driven by two CIs: (i) the first one is a sloped CI seam and is associated to a stretching of the ethenylic bond; (ii) the second one is a peaked CI and is associated to the twist of the methylene group . The LVC model of fulvene contains the intersection seam (i), but does not reproduce the peaked intersection (ii). We assume, that (ii) is not important in the simulated short time dynamics as argued previously . In Fig. we show the time trace of the populations of the diabatic state 1 (upper panel) and of the adiabatic state S 1 (lower panel). Starting in S 1 , the system experiences a nearly complete population transfer to S 0 within the first 15 fs but the population of S 1 rapidly increases again reaching a plateau value of 40 -50% depending on the simulation method. At around 30 fs, we observe again a fast decrease in the population of S 1 before reaching a second plateau, lower than the previous one. These oscillations repeat until the end of the simulated dynamics, with decreasing amplitude, such that a stepwise population transfer from S 1 to S 0 emerges, reminiscent of the reflection dynamics of the original one-dimensional Tully model #3. CTMQC, CCTTSH and TSH-ED agree with each other fairly well, while MTE population does not de-cay as much as the other methods. Once again, CTTSH does not manifest internal consistency: while CTTSH (F) shows the lowest population at most times, CTTSH (P) behaves similar to MTE, but this incorrect behavior is corrected by CCTTSH. The stepwise population transfer is also present in the diabatic representation, as shown in the upper panel of Fig. . There, comparison with MCTDH is reported for all methods except CTTSH. CT-MQC, CCTTSH and TSH-ED yield excellent agreement with the reference, while MTE shows, as in the adiabatic representation, a less complete population transfer after 30 fs. Since the population traces are similar to each other and, therefore, do not reveal much information about the different dynamics generated by the trajectory-based methods, we analyze once again the vibrational dynamics. We focus in this case on the vibrational mode q 24 , as this is the stretching mode involving the phenyl ring and the CH 2 moiety, as shown in the top right panel of Fig. using a ball-and-stick representation of fulvene. This mode is associated to the decay via the sloped CI discussed above.