Towards a low-order density matrix functional theory (motivated by exact solution of model aroms) with two and four electrons