Ehrenfest (EH) dynamics is a MQC approach for propagating the coupled electron-photon-nuclear dynamics. We choose to treat the electronic and photonic DOFs as fully quantum mechanical while treating the nuclear DOFs as classical ones. In this way, we can write the total wavefunction (Eq. 38) in an approximate form, including only the electronic and photonic DOFs explicitly and the nuclear DOFs as parameters as,