$\omega$-dependent solutions are not interesting, since after the substitution $t=x-\omega$
$$0=PV\int_{-\infty}^\infty e^{-\frac{\left(b(x-\omega)+\omega\right)^2}{2}} \frac{f(x)}{i(x-\omega)} \, {\rm d}x \, , \tag{1}$$
you could just set
$$f(x) \equiv e^{\frac{\left(b(x-\omega)+\omega\right)^2}{2}} i(x-\omega) h(x)$$
and you would have
$$0=PV\int_{-\infty}^\infty h(x) \, {\rm d}x$$
which certainly has many solutions for $h$. Therefore we can assume $f$ to not depend on $\omega$ and define
$$f(x) \equiv e^{\frac{b^2x^2}{2}} h(x)$$
after which (1) simplifies to
$$0=PV\int_{-\infty}^\infty e^{\omega x b(b-1)} \frac{h(x)}{i(x-\omega)} \, {\rm d}x$$
where the constant, $\omega$-dependent, factor was removed.

Now we assume $h$ to be entire and use contour integration to get
$$0=\oint_{-\infty}^\infty e^{\omega z b(b-1)} \frac{h(z)}{i(z-\omega)} \, {\rm d}z \\
= PV\int_{-\infty}^\infty e^{\omega x b(b-1)} \frac{h(x)}{i(x-\omega)}\, {\rm d}x - \pi h(\omega) \, e^{\omega^2 b(b-1)} + \lim_{R\rightarrow\infty} \int_0^\pi e^{\omega Re^{it} b(b-1)} \frac{h\left(Re^{it}\right)}{1-\frac{\omega}{R}e^{-it}} \, {\rm d}t \, .$$
The first line follows by Cauchys theorem and the integral is closed along the semi-circle in the upper half-plane. By the assumption of the vanishing PV-integral it follows
$$\pi h(\omega) \, e^{\omega^2 b(b-1)} = \lim_{R\rightarrow\infty} \int_0^\pi e^{\omega Re^{it} b(b-1)} \frac{h\left(Re^{it}\right)}{1-\frac{\omega}{R}e^{-it}} \, {\rm d}t \quad , \quad \forall \omega \in {\mathbb{R}} \, .$$

Now we expand both sides as a power-series in $\omega$. For the LHS we obtain
$$\pi h(\omega) e^{\omega^2 b(b-1)} = \pi \sum_{n,m=0}^\infty c_n \omega^n \frac{\left(\omega^2b(b-1)\right)^m}{m!} = \pi \sum_{k=0}^\infty \omega^k \sum_{\substack{n,m=0 \\ n+2m=k}}^\infty \frac{c_n \left(b(b-1)\right)^m}{m!} \, .$$
Similarly the RHS becomes
$$\lim_{R\rightarrow \infty} \sum_{n,m,k=0}^\infty c_n R^{n-k+m} \omega^{k+m} \frac{\left(b(b-1)\right)^m}{m!} \int_0^\pi e^{it(n-k+m)} \, {\rm d}t \\
=\lim_{R\rightarrow \infty} \sum_{l=0}^\infty \omega^l \sum_{\substack{n,m,k=0 \\ k+m=l}}^\infty c_n R^{n-k+m} \frac{\left(b(b-1)\right)^m}{m!} \\
\times \left( \pi \left[ n-k+m=0 \right] + 0 \left[n-k+m | \text{even} \right] + \frac{2i}{n-k+m} \left[n-k+m | \text{odd} \right] \right) \\
=\pi \sum_{l=0}^\infty \omega^l \sum_{\substack{n,m=0 \\ n+2m=l}}^\infty \frac{c_n\left(b(b-1)\right)^m}{m!} + 2i \lim_{R\rightarrow \infty} \sum_{l=0}^\infty \omega^l \sum_{\substack{n,m,k,p=0 \\ k+m=l \\ n-k+m=2p+1}}^\infty \frac{c_n R^{n-k+m} \left(b(b-1)\right)^m}{(n-k+m)m!} \, .$$
Because of the limit, $p$ starts at $0$. Now equating the expressions found for the LHS and RHS and canceling common terms, it follows that each $\omega^l$-coefficient has to vanish identically
$$0 = \lim_{R\rightarrow \infty} \sum_{\substack{n,m,p=0 \\ n+2m=l+2p+1}}^\infty \frac{c_n R^{2p+1} \left(b(b-1)\right)^m}{(2p+1)m!} \\
= \lim_{R\rightarrow \infty} \sum_{p=0}^\infty \sum_{0\leq 2m \leq l+2p+1} \frac{c_{l+2p+1-2m} R^{2p+1} \left(b(b-1)\right)^m}{(2p+1)m!} \quad , \quad \forall l\in\mathbb{Z}_{\geq 0} \, . \tag{2}$$
If $b=0$ or $b=1$, the $m$-sum disappears except for $m=0$ and so
$$\lim_{R\rightarrow \infty} \sum_{p=0}^\infty \frac{c_{l+2p+1} R^{2p+1}}{(2p+1)} = 0 \quad , \quad \forall l \in \mathbb{Z}_{\geq 0} \, .$$


Since the limit must exist for all $l\geq 0$, this is only possible if $$c_k =0 \quad , \quad \forall k \in \mathbb{N}$$

Note however, that in this case, there is no condition on $c_0$. Therefore $h(x)=\text{ const.}$, depending neither on $x$ nor on $\omega$.

In the case $b\neq 0$ and $b\neq 1$, the requirement (2) also forces $c_0=0$.