There are very deep reasons for this - the Born-Oppenheimer approximation, a central simplifying assumption, imposes some constraints on the system that make it easier to model with (semi)classical dynamics, even when you marry it with modern techniques such as DFT.