anchor.diff:              end do
           end associate
   
  -    end function slater_condon0_mol
  +    end function slater_condon0_mol_orb_list
   
       pure function slater_condon1_mol(sys, occ_list, i, a, perm) result(hmatel)
   
  
positive.code: 
      end function slater_condon2_mol_excit
  
  end module hamiltonian_molecular
  
negative.code:         integer, intent(in) :: i, j, a, b
          logical, intent(in) :: perm
  
          ! < D | H | D_{ij}^{ab} > = < ij || ab >
  
          hmatel = get_two_body_int_mol(sys%read_in%coulomb_integrals, i, j, a, b, sys%basis%basis_fns) &
  