High order numerical simulations for the polymer self-consistent field theory using the adaptive virtual element and spectral deferred correction methods.