We analyze a general class of self-adjoint difference operators [Formula: see text] on [Formula: see text], where [Formula: see text] is a multi-well potential and [Formula: see text] is a small parameter. We give a coherent review of our results on tunneling up to new sharp results on the level of complete asymptotic expansions (see [30–35]).Our emphasis is on general ideas and strategy, possibly of interest for a broader range of readers, and less on detailed mathematical proofs.The wells are decoupled by introducing certain Dirichlet operators on regions containing only one potential well. Then the eigenvalue problem for the Hamiltonian [Formula: see text] is treated as a small perturbation of these comparison problems. After constructing a Finslerian distance [Formula: see text] induced by [Formula: see text], we show that Dirichlet eigenfunctions decay exponentially with a rate controlled by this distance to the well. It follows with microlocal techniques that the first [Formula: see text] eigenvalues of [Formula: see text] converge to the first [Formula: see text] eigenvalues of the direct sum of harmonic oscillators on [Formula: see text] located at several wells. In a neighborhood of one well, we construct formal asymptotic expansions of WKB-type for eigenfunctions associated with the low-lying eigenvalues of [Formula: see text]. These are obtained from eigenfunctions or quasimodes for the operator [Formula: see text], acting on [Formula: see text], via restriction to the lattice [Formula: see text].Tunneling is then described by a certain interaction matrix, similar to the analysis for the Schrödinger operator (see [22]), the remainder is exponentially small and roughly quadratic compared with the interaction matrix. We give weighted [Formula: see text]-estimates for the difference of eigenfunctions of Dirichlet-operators in neighborhoods of the different wells and the associated WKB-expansions at the wells. In the last step, we derive full asymptotic expansions for interactions between two “wells” (minima) of the potential energy, in particular for the discrete tunneling effect. Here we essentially use analysis on phase space, complexified in the momentum variable. These results are as sharp as the classical results for the Schrödinger operator in [22].