An adequate preparation of the initial state of the system, point (i) above, depends on the process under investigation or on the experiment that one wants to simulate, and, thus, a general "recipe" for the correct procedure does not exist. Our group addressed this issue in earlier work focusing independently either on the nuclear 177 or on the electronic 142 aspects within the exact-factorization framework, but the problem has been discussed in the literature by other authors . The choice of initial conditions, namely the nuclear positions and momenta together with the electronic coefficients/populations at the initial time, can affect the calculated dynamics and observables, like photoproducts, lifetimes, and quantum yields, in many different ways. The dependence on the initial nuclear distribution might result more severe for the exactfactorization-based coupled-trajectory methods than for other approaches, since properly reconstructing the nuclear density along the dynamics is crucial to compute the quantum momentum that enters the expression of the classical force, in CTMQC ( ), and of the evolution equation of the electronic coefficients, in CTMQC (15) and in CTTSH (17).