Hello!

I am running constrained geometry optimization using the atomic confining potential method described here:

https://manual.q-chem.com/latest/subsec_HarmonicConfiner.html

To find appropriate harmonic constant in the computations I previously ran a harmonic vibrational calculation and used the force constant of a certain mode as the harmonic constant. I noticed that the unit of force constants in the output of the harmonic vibrational calculation is mDYN/Angstrom, but I have not found any information about the unit used for the harmonic constant in the atomic confining potential. Based on the value used in the example of the link above I assume it is N/m but I am not sure. Can someone help me where I can find information about that?

Thank you very much in advance!