where ρ IJ = c * I c J is the electronic density matrix element. Since the probability should not be negative, we set f IJ to 0 if f IJ < 0. Then, the non-adiabatic transition, i.e. stochastic switch from the current occupied state |I to another state |K , occurs if the following expression is met