My question: Is there a necessary and/or sufficient condition we can place on suitable continuous $f : [0,\infty) \rightarrow \mathbb{R}$ which allows us to determine whether the process $$X_t := \int_0^t f(t-s) dW_s$$ has differentiable paths?
I was motivated to ask this after looking at the SDE in this thread: \begin{align} dX_t &= Y_t dt \\\\ dY_t &= -X_t dt + b dW_t \end{align} to which I found the solution, when $X_0 = Y_0 = 0$, \begin{align} X_t &= b \int_0^t \sin(t-s) dW_s \\\\ Y_t &= b \int_0^t \cos(t-s) dW_s \\\\ \end{align} Now, from the SDE and the pathwise continuity of solutions, it is clear to see that $X_t$ will be pathwise differentiable (with derivative $Y_t$), whereas $Y_t$ will not be. However this was unintuitive to me just looking at the integrals themselves!
I know from this thread that we cannot hope for differentiable paths when we replace the integrand with just $f(s)$. Moreover, I know that convolutions with white noise of the form \begin{equation} Y_t := \int_{-\infty}^\infty f(t-s) dW_s \end{equation} are pathwise differentiable if the square of the Fourier transform of $f$ decays faster than $|k|^{-1}$ as $|k| \rightarrow \infty$ (the square of the Fourier transform of $f$ is the spectral density of $Y$). However, I don't think this result holds in our case since we aren't integrating over the whole line and, in any case $\sin$ and $\cos$ are not $L^1$.