configuration = nlread("methanol_box.hdf5")[0]

builder = OPLSPotentialBuilder()
calculator = SoftMatterCalculator(builder)
configuration.setCalculator(calculator)

energy = TotalEnergy(configuration)
print(f"The total energy of the configuration is {energy.evaluate()}")

