When including the environment, the system tends to approach the ground state S 0 slightly faster than when it is isolated, as the decrease in the S 1 population attests for both CT-TSH and TSH-ED. Furthermore, as observed previously, the inclusion of the environment reduces the deviation between the two methods even at long times. Conversely, both trajectory-based methods fail to correctly capture the decrease in the probability of forming the trans conformer at long times, that is expected when comparing reference MCTDH results from Ref. for the 2D and for the 2D + 23BM models. However, CT-TSH reproduces slightly better than TSH-ED the damping of the oscillations of the trans probability over time when the environment is included in the simulation.