In https://academic.oup.com/blms/article-abstract/20/2/121/266256?redirectedFrom=fulltext Vaughan shows the following bounds for the $L^1$-mean of the exponential sum over primes $$\sqrt x\ll \int _0^1\left |\sum _{n\leq x}\Lambda (n)e(\alpha n)\right |d\alpha \ll \sqrt {x\log x}$$ and further says it's likely to be $\asymp \sqrt {x\log x}$. (Note that by Cauchy-Schwarz the upper bound is trivial so the question is whether the $\sqrt {\log x}$ power is really there).
Question: Is there a heuristic reason why this should be so? Or are there any numerical calculations of it?
I initially thought it may not be true based on the following reasoning: Replacing $\Lambda (n)$ by $\omega (n)$, we can get the same $\sqrt x$ lower bound whilst the trivial CS bound is $\sqrt x\log \log x$, however I believe I can show that for $\omega (n)$ the sum really is $\asymp \sqrt x$. I guess too optimistically I thought the removing the $\log \log x$ for $\omega (n)$ may amount to removing the $\sqrt {\log x}$ for $\Lambda (n)$. However, there may really be no link between $\omega (n)$ and $\Lambda (n)$ and there's no reason to expect any similarities in behaviour.
I have tried to do some numerical calculations for the sum with $\Lambda (n)$ and it does seem that the true size is indeed $\sqrt {x\log x}$, but I don't want to completely eliminate the possibility of it being smaller based on my (bad) programming skills either.