def METHOD_NAME(kigaki_mols):
    # use all of the molecules
    mols = kigaki_mols
    map = {
        "cutoff": 10 * sr.units.angstrom,
        "cutoff_type": "REACTION_FIELD",
        "dielectric": 78.0,
        "platform": "Reference",
    }
    omm = sr.convert.to(mols, "openmm", map=map)
    state = omm.getState(getEnergy=True, enforcePeriodicBox=True)
    energy = state.getPotentialEnergy()
    # get this as a float in kJ mol-1
    energy = energy.value_in_unit(energy.unit)
    # -74975.9 +/- 0.5
    assert mols.energy(map=map).to(sr.units.kJ_per_mol) == pytest.approx(
        energy, abs=0.5
    )