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.
pop(always,nbo6)
. Never ever use the built in nbo module of Gaussian. This is incredibly old, has no support and plenty of know issues. This is only ever suitable for reproduction of values and showing that they are wrong. Also, part 2 of the question: chemistry.stackexchange.com/q/62530/4945 $\endgroup$