We have been working with crosslinking thermosetting polymers and ultimately will be exploring interfaces between systems with different levels of crosslinking. The chemists in our research group are asking us to characterize the types of bonds that exist in a given volume and I am looking for a script that can start us down the path. An existing script that generates a bond length distribution would be a good start to work from. Is there such a script available (where do I look for such scripts in general?). The idea would be to characterize distances between atoms types and for interatomic distances less than the sum of vdW radii, characterize the type of non-covalent interaction further. I realize this is a relatively deep subject, but for a start a simple bond distribution function would be something to work from. In particular, we are interested in categories such as O-H (hydrogen bonding), vdW interactions, H-H interactions, and possible benzene ring Pi bond interactions with nearby H. Any suggestions would be gratefully received.