Look into our groups (CATCO group) work (which would directly answer your question):
Decoding chemical information from vibrational spectroscopy data: Local vibrational mode theory, DOI:10.1002/wcms.1480
The local mode program enables users to gather information on intrinsic bond strength for all types of bonds including hydrogen bonds through the use a local vibrational force constants - providing the most adequate measure of bond strength. In the article above the software is introduced - I believe our group has released the software or is near to releasing the software for public use.
There are articles are group has published based upon determining the strength of hydrogen bonds in various systems - I have listed a few below:
A Comprehensive Analysis of Hydrogen Bond Interactions Based on Local Vibrational Modes DOI: 10.1002/qua.24118
Different Ways of Hydrogen Bonding in Water - Why Does Warm Water Freeze Faster than Cold Water? DOI: 10.1021/acs.jctc.6b00735
Article Strengthening of hydrogen bonding with the push-pull effect
You can find more related articles from our groups and download copies of the publications for free here:
It may be useful for you to check our latest article (DOI:10.1039/D0RA05802E), in which we estimate intramolecular H-bond strength with different methods using Gaussian and Multiwfn packages.
A reliable and very easy way of evaluating H-bond energy is using the relationship based on AIM theory proposed in my recent work J. Comput. Chem., 40, 2868–2881 (2019) DOI: 10.1002/jcc.26068. This analysis can be easily and very quickly realized using Multiwfn (http://sobereva.com/multiwfn). For Gaussian user, one only need to generate .wfn or .wfx or .fch first (see beginning of Chapter 4 of Multiwfn manual on how to generate them), then load the file into Multiwfn, enter topology analysis module, search critical point, then read electron density corresponding to the H-bond you are interested, and finally substitute this value into the formula proposed in the J. Comput. Chem. paper to predict H-bond energy. Please carefully check Section 4.2.1 of Multiwfn manual, this is a very detailed example showing how to use this method to predict H-bond energy between 2-pyridoxine and 2-aminopyridine in their dimer.
By the way, Multiwfn is also able to perform this analysis based on output file of GAMESS-US (you need to change the suffix to .gms), and based on .molden file generated by Molpro, see beginning of Chapter 4 of Multiwfn manual for detail.