Help: Error encountered in gmx hbond hydrogen bond calculation

GROMACS version:2023.3
GROMACS modification: Yes/No
Hello everyone, I currently want to use the gmx hbond command to calculate hydrogen bonds between specified protein atom groups. First, I generated the specified atom groups using the command gmx make_ndx -f result.gro -o index.ndx. Then, I processed the trajectory file with gmx trjconv -f result.xtc -s md.tpr -o 40-50ns.xtc -b 40000 -e 50000 -pbc mol -center -ur compact to center the protein and save a 10ns trajectory file. Finally, I used the command gmx hbond -f 40-50ns.xtc -s md.tpr -n index.ndx -num -hbn -hbm to generate three files: hbmap.xpm, hbnum.xvg, and hbond.ndx. The hbmap.xpm and hbnum.xvg files are shown in the figures below. However, in the [ hbonds_Protein ] section of the hbond.ndx file, the specific atoms forming the hydrogen bonds are not correctly displayed, and the middle atom always shows an index of -1. Could anyone explain what might be causing this issue?
hbond.ndx:


hbnum.xvg:

Could you try the gmx hbond tool in the 2025 release? There has been a complete rewrite of the code and the issue is likely fixed there.

Okay, thank you for the reminder. I’ll give it a try