HBCalculator: A Tool for Hydrogen Bond Distribution Calculations in Molecular Dynamics Simulations

Hydrogen bonds, crucial noncovalent interactions in molecular systems, significantly impact biological, chemical, and energy-related processes; therefore, characterizing hydrogen bond information is of importance to fundamental studies. This work introduces the HBCalculator, a Tcl-based tool integrated with VMD for calculating 1D and 2D distributions of hydrogen bond density and strength. The tool facilitates spatial analysis, overcoming limitations in existing packages that lack direct spatial distribution output. By employing HBCalculator in MD simulations, three systems of cellulose/water and graphene/water interfaces, were tested to showcase its functionality. The 1D and 2D hydrogen bond distributions reveal insights into interfacial properties, reflecting the impact of material hydrophobicity. The simplicity of usage, along with its potential for diverse molecular systems, positions HBCalculator as a valuable tool for researchers exploring hydrogen bond networks.