We split the Hamiltonian into two parts $H=H_n+H_{~n}$, where $H_n$ describes a particle in a single delta-function potential well, and $H_\overline{n}$ is the perturbation by the other delta functions:
We split the Hamiltonian into two parts $H=H_n+H_{\hat{n}}$, where $H_n$ describes a particle in a single delta-function potential well, and $H_\hat{n}$ is the perturbation by the other delta functions:
such that $H_n|n\rangle = -\epsilon_0|n\rangle = -\hbar^2\kappa^2/2m |n\rangle$ with $\kappa=mV_0/\hbar^2$. We can now calculate
$$
\langle n | H |n \rangle = \epsilon_0 + \langle n |H_\overline{n}|n\rangle
\langle n | H |n \rangle = \epsilon_0 + \langle n |H_\hat{n}|n\rangle
$$
Note that the last term represents the change in energy of the wavefunction $|n\rangle$ that is centered at the $n$-th delta function caused by the presence of the other delta functions. This term yields