I am interesting in bounding the arithmetic sum
$$ \sum_{n \leq x} \frac{\mu(n)^2}{\varphi(n)}$$
(The motivation is that this is a sum that comes up a lot in sieving primes, in particular in the Selberg sieve.) It is not too hard to show that this is $\geq \sum_{n \leq x} \frac{1}{n} \approx \log x$, and in fact the difference be described succinctly as follows: if $s(m)$ denotes the largest squarefree factor of $m$, then
$$\sum_{n \leq x} \frac{\mu(n)^2}{\varphi(n)} = \sum_{s(m) \leq x} \frac{1}{m}$$
In particular, one sees that the difference
$$ \sum_{n \leq x} \frac{\mu(n)^2}{\varphi(n)} - \sum_{n \leq x} \frac{1}{n} = \sum_{m > x, s(m) \leq x} \frac{1}{m}. $$
I would like to be able to bound how small this actually is. Montgomery-Vaughan cites an unpublished result of De Bruijn that is is actually $O(1)$, so I would certainly like to see that.
The closeness of these sums is rather unintuitive to me (apart from the vague comment that $\varphi(n) \approx n$ unless $n$ has a lot of distinct prime factors, which should be relatively rare), so I would also appreciate any insight into why we could expect this.