Suppose $f(x)$ is a probability density on $\mathbb{R}$. Let $\varphi(t)=\int e^{itx}f(x)dx$ denote the Fourier transform (characteristic function). It is well-known that if $\int |x|^p f(x)dx<\infty$, then $\varphi^{(p)}(t)= i^p \int x^p e^{itx} f(x) dx $. If in addition $\varphi^{(p)}\in L^1$, then through inverse Fourier transform $$ (ix)^pf(x)=\frac{1}{2\pi} \int e^{-ixt} \varphi^{(p)}(t)dt $$ one has the bound $$|f(x)|\le \frac{1}{2\pi|x|^p} \int |\varphi^{(p)}(t)|dt,$$ and hence one obtains the tail decay rate $|f(x)|=O(|x|^{-p})$ as $|x|\rightarrow\infty$.
I wonder if there is a generic way to obtain an estimate better than $O(|x|^{-p})$ by working with suitable conditions on $\varphi(t)$. I have this question because for "nice" $f$ (e.g., $f$ decays like a power function) to satisfy $\int |x|^p f(x)dx<\infty$, one typically needs $|f(x)|=o(|x|^{-p-1})$ in order to have integrability for large $|x|$. Hence there is a gap between the expected order of decay $o(|x|^{-p-1})$ and the order $O(|x|^{-p})$ obtained from $\varphi$.
Another perspective is that Infor the inverse Fourier transform, one has by Riemann Lebesgue Lemma that $\int e^{-ixt} \varphi^{(p)}(t)dt$ tends to zero when $|x|$ tends to infinitelyinfinity. So it seems possible to extract extra decay rate under suitable conditions.