Hi,
add
$nacme
to the control file (and optionally activate weight derivatives in $dft to improve accuracy of the quadrature).
Then run excited state gradient calculation with egrad as usual. egrad will then run in the non-adiabatic mode, the resulting gradients and nonadiabatic coupling matrix are printed to the output and in the control file as $couplingvector.
Regards,
Uwe