def make_molecular_hamiltonian(self):
    if self.active_space:
        return self.molecule.get_molecular_hamiltonian(occupied_indices=self.active_space.frozen_reference_orbitals, active_indices=self.active_space.active_orbitals)
    else:
        return self.molecule.get_molecular_hamiltonian()
def make_molecular_hamiltonian(self):
    if self.active_space:
        return self.molecule.get_molecular_hamiltonian(occupied_indices=self.active_space.frozen_reference_orbitals, active_indices=self.active_space.active_orbitals)
    else:
        return self.molecule.get_molecular_hamiltonian()