def METHOD_NAME(self, sage_unconstrained, smiles):
    """
    Basic test to see if the amber energy driver is functional.
    Note this test can only use the unconstrained version of Sage because sander applies SHAKE
    constraints in the single-point energy calculation, i.e. uses geometries with constraints
    applied, NOT what is in the coordinate file. See issue
    https://github.com/openforcefield/openff-interchange/issues/323
    """
    mol = Molecule.from_smiles(smiles)
    mol.generate_conformers(n_conformers=1)
    top = mol.to_topology()
    interchange = Interchange.from_smirnoff(sage_unconstrained, top)
    interchange.box = [5, 5, 5]
    interchange.positions = mol.conformers[0]
    omm_energies = get_openmm_energies(interchange, combine_nonbonded_forces=True)
    amb_energies = get_amber_energies(interchange)
    # TODO: More investigation into possible non-bonded energy differences and better reporting.
    #       03/02/2023 manually inspected some files and charges and vdW parameters are
    #       precisely identical. Passing box vectors to prmtop files might not always work.
    omm_energies.energies.pop("Nonbonded")
    amb_energies.energies.pop("vdW")
    amb_energies.energies.pop("Electrostatics")
    omm_energies.compare(amb_energies)