Unit of harmonic constant in atomic confining potentials method


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

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!

It’s N/m. The value 450 in the sample job that ships with Q-Chem (harm-conf.in) coincides with 450 N/m here,
I will update the manual to make this explicit.

Thank you very much for the quick reply and the explanation!