When CalculationMode=td, PCMTDLevel it sets the way the time-depenendent solvent polarization is propagated.
If PCMTDLevel=eq, the solvent is always in equilibrium with the solute or the external field, i.e.,
the solvent polarization follows instantaneously the changes in solute density or in the external field.
PCMTDLevel=neq and PCMTDLevel=eom are both nonequilibrium runs.
If PCMTDLevel=neq, solvent polarization charges are splitted in two terms:
one that follows instantaneously the changes in the solute density or in the external field (dynamical polarization charges),
and another that lag behind in the evolution w.r.t. the solute density or the external field (inertial polarization charges).
If PCMTDLevel=eom, solvent polarization charges evolves following an equation of motion, generalizing 'neq' propagation.
The equation of motion used here depends on the value of PCMEpsilonModel.