Through some properties of von Mangoldt we can extract Riemann zeta zeros alike Fourier transform but without going beyond scale 1 as this hides prime numbers in the process. This is a flat part of one of the possible options:
$\displaystyle f(x)=1-\frac{\ln(x)}{2}+\sum\limits_{n=1}^{M}\left(\Lambda(n)-\frac{n}{M}\right)\frac{\cos(x\ln(n))}{\sqrt{n}}$
M=600
(It is necessary to increase M in order to get a higher accuracy of larger Riemann zeta zeros, as one would expect. Increasing M, however, does not necessarily increase the accuracy of the positions of smaller zeros, as one would expect.)
For small values of M, $f(x)$ is behaving nicer than taking, for example, purely $\Lambda(n)$. If you take $\Lambda(n)-1$ (instead of $\Lambda(n)-\frac{n}{M}$) the function does not behave as smoothly and it is not clear if the error term remains bounded.
This is a great tool as it illustrates nicely the connection between Riemann zeta zeros displayed in this partial Fourier transform and Riemann zeta pole at 1. For example, this is what is the part of the illustration.
We can find from explicit formulae that it is formally:
$\displaystyle \Lambda(n)=\lim_{\epsilon \to 0}\int_{n-\epsilon}^{n+\epsilon}\sum\limits_{\rho}\frac{1}{x^{\rho}} dx$
$\displaystyle \sum\limits_{n=1}^{\infty} \frac{\Lambda(n)}{n^s} =\lim_{\epsilon \to 0}\int_{n-\epsilon}^{n+\epsilon}\sum\limits_{n=1}^{\infty}\sum\limits_{\rho}\frac{1}{n^{s}x^{\rho}} dx$
$\displaystyle \sum\limits_{n=1}^{\infty} \frac{\Lambda(n)}{n^s} =2\epsilon\sum\limits_{n=1}^{\infty}\sum\limits_{\rho}\frac{1}{n^{s}n(\epsilon)^{\rho}} \approx 2\epsilon\sum\limits_{\rho}\zeta(s+\rho)$
If you take $s \to \frac{1}{2}+\rho_{k}$ this last is formally saying that we can expect a sharp drop to the left at each Riemann zeta zero because to the right we have a pole at 1 and Riemann zeta zeros are symmetrical. $f(x)$ in some form might be the answer in what sense the above formal derivation that ignores much of the convergence works.
Another thing is that first Riemann zeta zero can be approximated by simply minimizing
$\displaystyle \frac{\ln(2)\cos(x\ln(2))}{\sqrt{2}}+\frac{\ln(3)\cos(x\ln(3))}{\sqrt{3}}$
around 14, and so on. All this because primes are accessible in the expression.
It seems that $f(x)$ never reaches more than $O(1)$, outside Riemann zeta zeros and region that is initially x < 10 and is slowing moving to the right, (it seems it stays within -5 and 5), for M increased as one is pleased, but I cannot reach any hint of proving it.
Does it look more obvious to someone else?
(In general for the expression above the unbounded oscillations seen to the left are slowly overwhelming the zeros peaks. Still, even for M one million, this will start engulfing only the first zero. The part staying to the right seems bounded all along. Is this visible smoothness anything substantial or should we be happy with $\Lambda(n)-1$ and consider the above useful for small zeros only? That is the question behind the question.)
You need to be more careful with the convergence.
Under the RH $|\psi(x)-x| \le (\log x)^2 x^{1/2}$ so by summation by parts for $\Re(s) > 1/2$
$$\frac{1}{|s|}\ \left|\frac{\zeta'(s)}{ \zeta(s)}+\zeta(s)-\sum_{n=1}^{N+1} \frac{1-\Lambda(n)}{n^s}+\frac{\psi(N\!+\!1)\!-\!N\!-\!1}{N^s}\right| = \left| \sum_{n=N+1}^\infty (n\!-\!\psi(n))\int_n^{n+1} y^{-s-1}dy\right|$$ $$ \textstyle\le \int_N^\infty (\log y)^2y^{1/2-\sigma-1}dy = \frac{d^2}{d^2\sigma} \frac{N^{1/2-\sigma}}{\sigma-1/2}=\frac{N^{1/2-\sigma}}{\sigma-1/2}\left((\log N)^2-\frac{2 (\log N)^2}{\sigma-1/2}+\frac{2}{(\sigma-1/2)^2}\right)$$
Also since $\frac{\zeta'(s)}{\zeta(s)}+\frac{1}{s-1} = \sum_\rho \frac{1}{s-1/2-i\Im(\rho)}- \log 2\pi + 2\frac{\Gamma'(s/2)}{\Gamma'(s/2)}$, you can use the preceding to detect the number of non-trivial zeros in the neighborhood of $s$.