Aquí está mi intento de respuesta, siguiendo la sugerencia de @Lagerbaer. Primero nos subtitute la transformada de Fourier para $\psi_{LP}(k)$,
\begin{equation}
\psi_{LP}(k)=\int dxe^{-ikx}\psi_{LP}(x),
\end{equation}
y obtener
\begin{multline}
\int dxe^{-ikx}i\frac{d}{dt}\psi_{LP}(x)=\int dx\left[\epsilon(k)-i\frac{\gamma(k)}{2}\right]e^{-ikx}\psi_{LP}(x)\\+F_{p}(k)\,\, e^{-i\omega_{p}t}
+\int dx_{3}dx_{1}dx_{2}\\ \times \sum_{q_{1},q_{2}}g_{k,q_{1},q_{2}}\, e^{i(q_{1}+q_{2}-k)x_{3}}\psi_{LP}^{\star}(x_{3})\, e^{-iq_{1}x_{1}}\psi_{LP}(x_{1})\, e^{-iq_{2}x_{2}}\psi_{LP}(x_{2}).
\end{multline}
Con el fin de simplificar esta expresión, tenemos que aplicar la inversa de la transformación a ambos lados (multiplicar por $\int dke^{ikx}$), la obtención de
\begin{multline}
i\frac{d}{dt}\psi_{LP}(x)=\int dx^{\prime}\psi_{LP}(x^{\prime})\int dk\left[\epsilon(k)-i\frac{\gamma(k)}{2}\right]e^{ik(x-x^{\prime})}+F_{p}e^{i(k_{p}x-\omega_{p}t)}\\
+\int dke^{ikx}\int dx_{3}dx_{1}dx_{2}\\
\times\sum_{q_{1},q_{2}}g_{k,q_{1},q_{2}}\, e^{i(q_{1}+q_{2}-k)x_{3}}\psi_{LP}^{\star}(x_{3})\, e^{-iq_{1}x_{1}}\psi_{LP}(x_{1})\, e^{-iq_{2}x_{2}}\psi_{LP}(x_{2}).
\end{multline}
Ahora nos fijamos en una integral de la forma $\int dx^{\prime}\psi_{LP}(x^{\prime})\int dkf(k)e^{ik(x-x^{\prime})}$, para una función genérica $f(k)$. Ampliamos nuestra función genérica en una Serie de Taylor alrededor de $k=0$, y obtener:
\begin{equation}
\sum_{n=0}^{\infty}\frac{f^{(n)}(0)}{n!}\int dx^{\prime}\psi_{LP}(x^{\prime})\int dkk^{n}e^{ik(x-x^{\prime})}=\sum_{n=0}^{\infty}\frac{f^{(n)}(0)}{n!}\left(-i\right)^{n}\int dx^{\prime}\psi_{LP}(x^{\prime})\partial_{x}^{n}\delta(x-x^{\prime})
\end{equation}
donde hemos utilizado la conocida relación $\int dkk{}^{n}e^{ik(x-x^{\prime})}=\left(-i\right)^{n}\partial_{x}^{n}\delta(x-x^{\prime})$.
Mover el operador de la derivada para actuar en $\psi$ (integración por partes), se obtiene finalmente
\begin{equation}
\int dx^{\prime}\psi_{LP}(x^{\prime})\int dkf(k)e^{ik(x-x^{\prime})}=\sum_{n=0}^{\infty}\frac{f^{(n)}(0)}{n!}(i\partial_{x})^{n}\psi_{LP}(x)=f(i\partial_{x})\psi_{LP}(x)
\end{equation}
Por lo tanto, el primer término en el lado derecho de la ecuación se simplifica a $\left[\epsilon(i\partial_{x})-i\frac{\gamma(i\partial_{x})}{2}\right]\psi_{LP}(x)$ y nos queda tratar con el término de interacción.
\begin{equation}
\int dx_{3}\psi_{LP}^{\star}(x_{3})\int dx_{1}\psi_{LP}(x_{1})\int dq_{1}e^{iq_{1}(x_{3}-x_{1})}\int dx_{2}\psi_{LP}(x_{2})\int dq_{2}e^{iq_{2}(x_{3}-x_{2})}\int dkg(k,q_{1},q_{2})e^{ik(x-x_{3})}
\end{equation}
La sustitución de la forma de $g(k,q_{1},q_{2})$, obtenemos
\begin{equation}
g\int dx_{3}\psi_{LP}^{\star}(x_{3})\int dx_{1}\psi_{LP}(x_{1})\int dq_{1}X(q_{1})e^{iq_{1}(x_{3}-x_{1})}\int dx_{2}\psi_{LP}(x_{2})\int dq_{2}X(q_{2})e^{iq_{2}(x_{3}-x_{2})}\int dk\, X^{\star}(k)\, X^{\star}(q_{1}+q_{2}-k)\, e^{ik(x-x_{3})}
\end{equation}
Y aquí es donde me quedé atrapado..