Background: Harmonic series approaches to $\lim_{n->\infty}\log(n)$. This gives the Euler–Mascheroni constant. In addition, the harmonic series summed only over the primes approaches to $\lim_{n->\infty}\log(\log(n))$. This gives the Meissel–Mertens constant.
Question: What similar series approaches to $\lim_{n->\infty}\log(\log(\log(n)))$.
My failed attempts: I thought about further filtering out primes to get an even slower-growing series. For example, I only choose the paired primes (twin primes or $k$-twin primes). But then I learned that these series converge using Brun's theorem.
There ought to be someone who has thought about this problem possibly with answers. Can anyone please point me in the right direction? Thanks.
The heuristic argument for why the sum of prime reciprocals diverges as $\log\log n$ is that we may model an integer $n$ as being "prime with probability $\frac1{\log n}$", by the prime number theorem: so its contribution to the sum of prime reciprocals is on average $\frac1{n \log n}$. Further approximating this sum by an integral, we notice that $$\int \frac{1}{x \log x}\,\mathrm dx = \log\log x + C.$$ Similarly, it is true that $$\int \frac{1}{x \log x \log\log x}\,\mathrm dx = \log\log\log x + C,$$ so we want a sum over reciprocals where an integer $n$ has a probability of $O(\frac1{\log n\log\log n})$ of contributing to the sum.
So we could, for example, say: "Take the sum of reciprocals of all primes which have a prime number of digits." Heuristically, a number $n$ has a $\frac1{\log n}$ chance of being prime, but it has $O(\log n)$ digits, and he number of digits has an $O(\frac1{\log \log n})$ chance of being prime.
Now let's try to figure this out carefully. What is the sum of reciprocals of all prime-digit primes less than $n$?
Well, for every prime $p$ such that $10^p < n$, we have the sum of reciprocals of all primes from $10^{p-1}$ to $10^p-1$. By the Meissel–Mertens result, for sufficiently large $p$, this should contribute $(1+o(1))(\log \log 10^p - \log \log 10^{p-1})$, which simplifies to $(1+o(1))\log \frac{p}{p-1})$. As $x \to 0$, $\log(1+x) \sim x$, and we can simplify this to $(1+o(1))\frac1p$.
Therefore for each prime $p$ such that $10^p < n$, the contribution to the sum is $(1+o(1))\frac1p$, which means that the overall sum is asymptotically equal to the sum of $\frac1p$ over all $p$ up to $\log_{10} n$. By the Meissel–Mertens asymptotics again, this should be $\log\log\log n$.