**Edit:** I simplified the example to a canonical case for clarity.

Given a integral $\int_{\Omega}{g(\mathbf{x})}$ with a well-posed integrand $g(\mathbf{x})$ defined on some multidimensional space $\Omega$, one can integrate it successfully with a Markov Chain Monte Carlo method, in particular using Metropolis-Hastings method. Detailed balance and ergodicity are achieved with any well-posed (samplable) integrand. The value of the desired integral is known to always exist and be finite, even in presence of delta distributions in the integrand. This is given.

In the case of interest, the integrand $f(\mathbf{x})=\delta_{\mathbf{x}_0}(\mathbf{x})+g(\mathbf{x})$ consists of a delta distribution at some unknown location $\mathbf{x}_0\in \Omega$ plus some regular (well-posed) non-zero integrand $g(\mathbf{x})$. This delta distribution cannot be sampled explicitly or using numerical optimization (given). This makes such an integrand unsamplable with random walk or probabilistic sampling. 

I mollify (approximate to the identity) this delta distribution using some mollifier (normalized smooth function $\phi_\epsilon(\mathbf{x})=\epsilon^{-1}\phi\left(\frac{\mathbf{x}}{\epsilon}\right)$ with some bandwidth $\epsilon$). This leads to a tempered integrand $f_\epsilon(\mathbf{x})=\phi_\epsilon(\mathbf{x}-\mathbf{x}_0)+g(\mathbf{x})$. During the integration, at every step $n$, I gradually shrink the parameter $\epsilon_n$ to zero in order to achieve $f_{\epsilon_n} \to f=\delta_{\mathbf{x}_0}(\mathbf{x})+g(\mathbf{x})$ as $n \to \infty$ in spirit of serial tempering. Thus I expect the integration method to be consistent, i.e. to give the proper answer in the limit.

Unfortunately one cannot use the usual parallel or serial tempering here, as the probability of the proposal to descend from a tempered mixture to the original mixture $f$ at the exact location $\mathbf{x}_0$ of the delta distribution is zero.

Thus I have two rather similar questions: 

 1. Would the integral converge to the proper value (which is $\int_{\Omega}{g(\mathbf{x})}+1$ in this case) with the proposed method if the integrand has a delta distribution? What confuses me is that the chain either tends to stuck in the vicinity of the delta distribution or would otherwise just miss it, as the mode of it becomes narrower when $n$ is large. 
 2. What should be the cool-down schedule of the integrand $f_{\epsilon_n}$, i.e. the asymptotic decrease rate of the sequence $\{ \epsilon_n \}$, in order to guarantee the consistent convergence of the MCMC estimate? In other words, what are the conditions for the cool-down schedule that allows the integral to converge before the integrand $f_{\epsilon_n}$ becomes unsamplable?