Multiple‐band theory of dynamics for interacting adsorbates coupled to phonons. I. Variationally optimized Hamiltonian