For a given bond, I would like to perform an NBO calculation for each structure while scanning the bond lengths.
Can I use Opt=ModRedundant in combination with pop=nbo? Simply specifying both keywords would only perform the NBO calculation for the last structure.
Is there an easy way to retrieve the structures obtained from the scan as separate input files (e.g. .xyz or .gjf)?
I can use Gaussian 16 and Gaussview 5 or Gaussview 6.