Numerical verification of nonlinear response beyond the efficiency-bandwidth limit in time-variant metasurfaces