def METHOD_NAME():
    def add(simulation):
        mc = hoomd.hpmc.integrate.Sphere()
        mc.shape['A'] = dict(diameter=0)
        snapshot = simulation.state.get_snapshot()
        if simulation.device.communicator.rank == 0:
            reference_positions = snapshot.particles.position
            reference_orientations = snapshot.particles.orientation
        else:
            reference_positions = [[0, 0, 0], [0, 0, 0]]
            reference_orientations = [[1, 0, 0, 0], [1, 0, 0, 0]]
        lattice = hoomd.hpmc.external.field.Harmonic(
            reference_positions=reference_positions,
            reference_orientations=reference_orientations,
            k_translational=1.0,
            k_rotational=1.0,
            symmetries=[[1, 0, 0, 0]])
        mc.external_potential = lattice
        simulation.operations.integrator = mc
        return mc, lattice
    return add