The time-dependent potentials emerging from the exact factorization can be constructed from the output of a vibronic wavepacket quantum simulation. If the electron-nuclear state of the system |Ψ(ϕ, q, t) is expanded in the diabatic basis with coefficients χ (d) k (ϕ, q, t) (k = 0, 1), the nuclear wavefunction in the chosen gauge is χ(ϕ, q, t) = |χ