Non-adiabatic coupling matrix elements for large CI wavefunctions