I'm struggling to derive the second-quantised expression of the Hamiltonian operator of a many-particle system. Starting from \begin{equation*} \renewcommand{\vec}[1]{\mathbf{#1}} \newcommand{\ket}[1]{\lvert{#1}\rangle} \newcommand{\bra}[1]{\langle{#1}\rvert} \newcommand{\braket}[2]{\langle{#1}\vert{#2}\rangle} H=\sum_{i=1}^n\frac{\vec p_i^2}{2m}, \end{equation*} I want to derive \begin{equation*} H=-\frac{\hbar^2}{2m}\iint\Psi^\dagger(\vec{x})\triangle\Psi(\vec{x})\,\mathrm{d}\vec{x}\,\mathrm{d}\vec{x}'. \end{equation*}
I know that, by definition, the second quantised form of $H$ is written as \begin{equation*} H=\iint\Psi^\dagger(\vec{x})\bra{\vec{x}}H\ket{\vec{x}'}\Psi(\vec{x}')\,\mathrm{d}\vec{x}\,\mathrm{d}\vec{x}'. \end{equation*} I don't know how to tackle the term \begin{equation*} \bra{\vec{x}}\vec{p}_i^2\ket{\vec{x}'} \end{equation*} that appears in the derivation. I tried simplifying the problem down to two bodies only, i.e. \begin{equation*} H=-\frac1{2m}(\vec{p}_1^2+\vec{p}_2^2)= -\frac1{2m}(\vec{p}^2\otimes 1+1\otimes\vec{p}^2), \end{equation*} for which I have calculated \begin{equation*} \begin{split} \bra{\vec{x}_1,\vec{x}_2}\vec{p}^2\otimes 1+1\otimes\vec{p}^2\ket{\vec{x}'_1,\vec{x}'_2}&= \bra{\vec{x}_1}\vec{p}^2\ket{\vec{x}'_1}\braket{\vec{x}_2}{\vec{x}'_2}+ \braket{\vec{x}_1}{\vec{x}'_1}\bra{\vec{x}_2}\vec{p}^2\ket{\vec{x}'_2}=\\ &= -\hbar^2\delta(\vec{x}_1-\vec{x}'_1)\delta(\vec{x}_2-\vec{x}'_2)\triangle_1-\hbar^2\delta(\vec{x}_1-\vec{x}'_1)\delta(\vec{x}_2-\vec{x}'_2)\triangle_2=\\ &= -\hbar^2\delta(\vec{x}_1-\vec{x}'_1)\delta(\vec{x}_2-\vec{x}'_2)(\triangle_1+\triangle_2), \end{split} \end{equation*} but (assuming I did it right) I don't know how to link it to the result I want to prove. In this last equation I have $n$ Laplacian operators and $n$ coordinate vectors, while in the previous one I had only one...