To this end, we take advantage of the existing multiple trajectories that are already present in mixed quantumclassical or semiclassical methods. We propose to evolve multiple trajectories and use the trajectory average of the electronic coefficients of these trajectories to compute the reduced density matrix elements at each time step, with the estimator of the reduced density matrix elements evaluated as