Hello all. Is there a simple way of performing a constrained minimization where one property (pair distance between two atoms, say) is held fixed and the rest of the degrees of freedom are relaxed? In Gaussian this is called a "bond length scan". I've been shown a way to do it in NWChem, but it's not nearly as concise as it seems like it should be.