def METHOD_NAME(self):
    super().METHOD_NAME()
    algorithm_globals.random_seed = 8
    self.reference_energies = [
        1889.95738428,
        3294.21806197,
        4287.26821341,
        5819.76975784,
    ]
    self.qubit_mapper = DirectMapper()
    import sparse as sp  # pylint: disable=import-error
    watson = WatsonHamiltonian(
        quadratic_force_constants=sp.as_coo(
            {
                (0, 0): 605.3643675,
                (1, 1): 340.5950575,
            },
            shape=(2, 2),
        ),
        cubic_force_constants=sp.as_coo(
            {
                (1, 0, 0): -89.09086530649508,
                (1, 1, 1): -15.590557244410897,
            },
            shape=(2, 2, 2),
        ),
        quartic_force_constants=sp.as_coo(
            {
                (0, 0, 0, 0): 1.6512647916666667,
                (1, 1, 0, 0): 5.03965375,
                (1, 1, 1, 1): 0.43840625000000005,
            },
            shape=(2, 2, 2, 2),
        ),
        kinetic_coefficients=sp.as_coo(
            {
                (0, 0): -605.3643675,
                (1, 1): -340.5950575,
            },
            shape=(2, 2),
        ),
    )
    basis = HarmonicBasis([2, 2])
    self.vibrational_problem = watson_to_problem(watson, basis)