In Gaussian, there is a restart option for the PES scan. Is there such an option in QChem?
This would be relatively easy to do by recognizing where the previous job stopped and setting the boundaries of the scan appropriately, using as input the final geometry from the previous job. Example in the manual include scanning a bond length,
$scan
stre 1 2 1.0 2.0 0.5
$end
or a torsion angle and a bond length,
$scan
tors 1 2 3 4 -180 180 30
stre 2 3 1.5 1.6 0.05
$end
In the latter case you would need to determine which is the outer loop, bond length or angle.
I see. Thank you for the clarification!