Best-of-both-worlds computational approaches to difficult-to-model dissociation reactions on metal surfaces