When I perform Born-Oppenheimer Molecular Dynamics (BOMD) simulations of small organic clusters using the CP2K software in NVE condition, I notice that running multiple simulations from the same initial structure produces remarkably similar trajectories. I'm curious about the underlying reason for this behaviour and would like to know how to introduce velocity randomisation in CP2K. I also attached input file of BOMD simulation.