Warning: I am probably going to mess up factors of $\pm i$, $\pi$, etc. So take what I say modulo these constants.
We know that if $f(x)$ is sufficiently nice, then integration by parts yields $\mathcal{F}\{x^n f(x)\} = i^n \frac{d^n}{dk^n}\hat{f}(k)$. Since the Fourier transform is linear, if $p(x)$ is a polynomial, then we can express $\mathcal{F}(p(x)f(x))$ as a polynomial in the derivatives of $\hat{f}(k)$. In your case, the function $f(x)$ is a Gaussian, hence has Gaussian Fourier transform, so the question is reduced to:
If $p$ is a polynomial in $i d/dk$, is there a closed form expression for $p^n e^{-k^2}$? Let $h_n(k) = e^{k^2} p^n e^{-k^2}$ Note that $h_n(k)$ is a polynomial in $k$, and we can find an explicit recurrence relation via \begin{align} h_{n+1}(k) &= e^{k^2} p p^n e^{-k^2} \\ &= e^{k^2} p ( e^{-k^2} h_n(k) ) \end{align}
For example, in the linear case we have $p = id/dk + a$, and we obtain \begin{align} h_{n+1}(k) &= e^{k^2} (id/dk + a)(e^{-k^2} h_n(k)) \\ &= e^{k^2}(-2ki e^{-k^2} h_n(k) + i e^{-k^2} h_n'(k) + a e^{-k^2} h_n(k) \\ &= -2ki h_n(k) + i h_n'(k) + ah_n(k) \\ &= (-2ki + a) h_n(k) + i h_n'(k). \end{align} Note that since the $h_n(k)$ are polynomials, all of this can be computed symbolically. Also note that setting $a=0$ recovers (up to constants and conventions) the usual Hermite polynomials. Actually I expect that in this case it should be easy to work out an explicit relation to the Hermite polynomials (by shifting $k$), but I'll leave that to you.
The quadratic case will certainly be much harder as you'll get a more complicated recurrence relation. It might be possible to find a generating function (as in the case of Hermite polynomials) and to use it to find a closed-form expression.