# Electronic dynamics¶

Electron nuclear dynamics within the Ehrenfest *ansatz* is based on the
integration of the following equation of motion (EOM) for the reduced one body
density matrix:

As noted above the EOM contains a purely electronic driving term in which the dynamics may be generated, either from an external time dependent driving field in \(H\) or from the fact that the density matrix might be in a non-equilibrium state and therefore does nos commute with \(H\).

The second term in equation (1) is part of the Ehrenfest *ansatz*
and includes the non adiabatic coupling matrix \(D\), the elements of which
are defined as \(D_{\kappa \nu} = \langle \phi_\kappa | \dot{\phi_\nu}
\rangle\), where \(| \phi_\kappa \rangle\) and \(| \phi_\nu \rangle\) are
a pair of atomic orbitals and, in a localised basis set such as the one used in
DFTB+, the time derivative of an atomic orbital with respect to time is a
function of the velocity of the nuclei to which it is attached. The
non-adiabatic term introduces a mechanism by which nuclear motion can induce
electronic transitions.

The forces used to propagate the nuclei are calculated from the instantaneous expectation value of the force operator adding the corresponding self consistent and repulsive terms. This force comes from a non equilibrium density matrix and is the form in which an excited density matrix may induce nuclear motion.

The EOM for the density matrix is integrated using a Leapfrog scheme in which the density matrix at time \(t_{i+i}\) is obtained from its value a t time \(t_{i-1}\) and its derivative at time \(t_i\):

\(\rho_{i+1}=\rho_{i-1}+2\Delta t \dot{\rho}_i.\)

The nuclear motion is integrated using a velocity Verlet algorithm.

The propagation of the electronic and nuclear dynamics may be used for the calculation of absorption spectra, to study the response to constant or pulsed excitation and even the simulation of pump-probe spectroscopy.

All input related to electronic dynamics is located within the
*ElectronDynamics* block of the input file. Electron dynamics is run after
ground state self consistency is achieved, with the possibility of reading
previously converged charges for the ground state or the state of the
time-evolving density matrix to continue time propagation.