Understanding the current-induced vibrational dynamics in molecular nanojunctions is critical for gaining insight into the stability of such systems. While it is well known that Joule heating at higher bias voltages plays an important role for the stability of the nanojunction, a different mechanism caused by current-induced nonconservative forces has been reported to cause vibrational instabilities already at much lower voltages. In this paper, we apply a mixed quantum-classical approach based on electronic friction and Langevin dynamics to a model system for which vibrational instabilities have previously been reported. The electronic friction and other electronic forces can be calculated exactly using nonequilibrium Green's functions. We perform Langevin simulations of the vibrational dynamics to calculate the steady-state vibrational excitation and compare the results to a simplified eigenmode analysis, identifying the parameter regimes in which runaway modes appear.
Mäck et al. (Wed,) studied this question.