I wish to estimate the hydrophobic interactions between a protein and ligand for a GROMACS trajectory. While the obvious and tedious way is to calculate distances and angles for each specific interaction, please share any other way to accurately determine the hydrophobic interactions based on geometric criteria.