Scaling active spaces in simulations of surface reactions through sample-based quantum diagonalization