We present an approach for carrying out non-adiabatic molecular dynamics simulations of systems in which non-adiabatic transitions arise from the coupling between the classical atomic motions and a quasi-continuum of electronic quantum states. Such conditions occur in many research areas, including chemistry at metal surfaces, radiation damage of materials, and warm dense matter physics. The classical atomic motions are governed by stochastic Langevin-like equations, while the quantum electron dynamics is described by a master equation for the populations of the electronic states. These working equations are obtained from a first-principle derivation. Remarkably, unlike the widely used Ehrenfest and surface-hopping methods, the approach naturally satisfies the principle of detailed balance at equilibrium and, therefore, can describe the evolution to thermal equilibrium from an arbitrary initial state. In addition, unlike other schemes, there is no need to explicitly propagate wave functions in time.