NOTE: I've edited the question one last time, to be much simpler, in the hopes of getting more focused and preciseresponses.
SETUP: Let $p_n$ denote the $n$th prime, let $p_x = p_{\lceil x \rceil}$ for all $x > 0$, let $\pi(x)$ denote the prime counting function, and let $\operatorname{li}(x) = \int_0^x \frac{dt}{\log t}$ denote the logarithmic integral function.
QUESTION: Does there exist a $t \in \mathbb{R}$ such that $ \pi(x)-\frac{x}{p_{\operatorname{li}(x)}}\operatorname{li}(x) \neq O(\sqrt{x} (\log x)^t)$?
MOTIVATION:
The function $\frac{x}{p_{\operatorname{li}(x)}}\operatorname{li}(x)$ seems to approximate $\pi(x)$ much better than $\operatorname{li}(x)$ does, at least for small $x$. The spirit of this post is to inquire whether or not $\frac{x}{p_{\operatorname{li}(x)}}\operatorname{li}(x)$ is a genuine improvement over $\operatorname{li}(x)$ in approximating $\pi(x)$ for large $x$.
Note that $$\pi(x)-\frac{x}{p_{\operatorname{li}(x)}}\operatorname{li}(x) = E(x)+H(x),$$ where $E(x) = \pi(x)-\operatorname{li}(x)$ and $H(x) = \frac{\operatorname{li}(x)}{p_{\operatorname{li}(x)}}(p_{\operatorname{li}(x)}-x)\sim \frac{1}{\log(x)}(p_{\operatorname{li}(x)}-x)$, and one has $E(x)\geq 0$ if and only if $H(x) \leq 0$ (so all "interference" in $E(x)+H(x)$ is destructive).
QUESTION: Let Let $\Theta \in [1/2,1]$$\Theta$ denote the supremum of the real parts of the zeros of the Riemann zeta function. The infimum of all $t > 0$ such that $E(x) = o(x^t)$ is $\Theta$. Likewise, the infimum of all $t > 0$ such Note that $H(x) = o(x^t)$ is alsoboth $\Theta$. Clearly, then, the infimum of all$E(x)$ and $t > 0$ such that $\pi(x)-\frac{x}{p_{\operatorname{li}(x)}}\operatorname{li}(x) = o(x^t)$ is less than or equal to$H(x)$ are $\Theta$. Does equality hold? If not, that answers my question fully. If so$O(x^\Theta \log x)$, let's assume that equality holds. Let $\Theta_1$ be the infimum of all $t \in \mathbb{R}$ such that $E(x) = o(x^\Theta (\log x)^t)$. Then $\Theta_1$ is also the infimum of all $t \in \mathbb{R}$ such that $H(x) = o(x^\Theta (\log x)^t)$. Thus the infimum all $t \in \mathbb{R}$ such that $\pi(x)-\frac{x}{p_{\operatorname{li}(x)}}\operatorname{li}(x)$ is $o(x^\Theta (\log x)^t)$and neither is less than or equal to $\Theta_1$$O(\sqrt{x}/\log x)$. Does equality hold?
MOTIVATION: Note that $x \geq p_y$ if and only if $\pi(x) \geq y$, for all $x, y > 0$. (This defines a monotone Galois connection from $\mathbb{R}_{>0}$ to itself, which expresses an adjoint relationship in the corresponding poset category). Thus $x \geq p_{\operatorname{li(y)}}$$x \geq p_{\operatorname{li}(y)}$ if and only if $\pi(x) \geq \operatorname{li}(y)$. The quantity $p_{\operatorname{li(x)}}\pi(x) - x \operatorname{li}(x)$ is thus natural to consider---discrete on one side, continuous on the other, and involving "both sides" of the duality, not just one. If the answer to the QUESTION above is NO, then $p_{\operatorname{li(x)}}\pi(x) - x \operatorname{li}(x)$ were$p_{\operatorname{li}(x)}\pi(x) - x \operatorname{li}(x)$ is smaller than $x(\pi(x)-\operatorname{li}(x))$ in order of growth (in the sense asked in the QUESTION above), then that would be interesting. Itwhich would showimply that the product $p_{\operatorname{li(x)}}\pi(x)$$p_{\operatorname{li}(x)}\pi(x)$ of the right-left adjoint pair is better approximated by $x \operatorname{li(x)}$ than either adjoint is by $x$ and $\operatorname{li(x)}$, respectively.