The Bennett acceptance ratio method is implemented in modern molecular dynamics systems, such as Gromacs.
Python-based code for MBAR and BAR is available for download at [https://github.com/choderalab/pymbar].