Higher-order MOM of volume integral equation using hierarchical vector basis functions