<A HREF="https://www.projecteuclid.org/journals/communications-in-mathematical-physics/volume-75/issue-3/Double-wells/cmp/1103908148.full">Double wells,</A> by Evans M. Harrell, Comm. Math. Phys. **75**, 239 (1980), should be sufficiently rigorous. The result for the tunnel splitting, theorem 2.22, applies to a symmetric double-well potential in arbitrary number $n$ of dimensions, constructed as the sum $V(x-f)+V(x+f)$ of two single-well potentials $V(x)$, bounded and of compact support, after translation by $\pm f$.

The tunnel splitting is given for large $f$ in terms of the eigenfunction $\Phi(x)$ for the lowest eigenstate in the potential $V(x)$,
$$E_1-E_2=\left(\int\Phi(x-f)[V(x-f)+V(x+f)]\Phi(x+f)\,dx\right)(1+o(f^{-n}).$$