Computing the product of p/(p - 2) over the odd primes

Lets do this as explicitly as possible (I would like to find the constant and error term)

1 The Original Problem First, consider the identity $\log\left(1-\frac{2}{p}+\frac{1}{p^2}\right) +\log\left(1- \frac{1}{(p-1)^2} \right)=\log\left(1- \frac{2}{p} \right)$. From this, it follows that

$$\sum_{2<p\leq n}\log\left(1-\frac{2}{p}+\frac{1}{p^2}\right) +\sum_{2<p\leq n} \log\left(1- \frac{1}{(p-1)^2} \right) = \sum_{2<p\leq n}\log\left(1-\frac{2}{p}\right)$$

Multiplying by negative one and exponentiating both sides yields

$$\left(\prod_{2<p\leq n}1-\frac{1}{p}\right)^{-2} \cdot \prod_{2<n<p} \left(1- \frac{1}{(p-1)^2} \right)^{-1} = \left(\prod_{2<p\leq n}1-\frac{2}{p}\right)^{-1} $$

Recall $ \Pi_2=\prod_{2<p} \left(1- \frac{1}{(p-1)^2} \right)$ is the Twin Prime Constant, and that the one product on the left hand side converges to the reciprocal of this. It is then by one of Mertens formulas that we know $$\left(\prod_{2<p\leq n}1-\frac{1}{p}\right)^{-1}=\frac{1}{2}e^\gamma \log n + O(1)$$ where $\gamma$ is the Euler Mascheroni Constant. (Specifically this is Theorem 2.7 (e) in Montgomery Multiplicative Number Theory I. Classical Theory) Upon squaring this asymptotic result, we are able to conclude:

$$\left(\prod_{2<p\leq n}1-\frac{2}{p}\right)^{-1} = \frac{1}{4}e^{2\gamma}\Pi_2^{-1} \log^2n + O(\log n)$$

Hope that helps,

Note: The reason I substituted the twin prime constant in for that other product is because it converges very fast in comparison with the error term. I can give more details if desired, but I'll leave it as an exercise.

2 What is best possible? Can the error term be made better? Yes. It turns out we can make that error term a lot better. By using the Prime Number Theorem we find $$\prod_{2<p\leq n} \left( 1-\frac{1}{p}\right)^{-1}=\frac{1}{2}e^\gamma \log n + e^{-c\sqrt{\log n}}$$ where $c$ is the constant used in the proof of the Zero Free Region. Since $e^{-\sqrt{\log n}}$ decreases faster than any power of $\log$, we obtain a much better result upon squaring this estimate. Precisely we have:

$$\left(\prod_{2<p\leq n}1-\frac{2}{p}\right)^{-1} = \frac{1}{4}e^{2\gamma}\Pi_2^{-1} \log^2n + O\left( e^{-c\sqrt{\log n}} \right)$$

(again the convergence to $\Pi_2$ is much to rapid to interfere with the error term)

I would be willing to bet that this is the best we can do, and that any better would imply stronger results regarding the error term for $\pi(x)$, the prime counting function.

3 Numerics Just for fun, the exact constant in front of the $\log^2n$ is: 1.201303. How close is this? Well for:

$n=10$ we get an error of $0.630811$

$n=50$ we get an error of $1.22144$

$n=100$ we get an error of $0.63493$

$n=1000$ we get an error of $0.438602$

$n=10^4$ we get an error of $0.250181$

$n=10^5$ we get an error of $0.096783$

$n=10^6$ we get an error of $0.017807$

Where each time the error is positive. That is the product seems to be slightly larger than the asymptotic (but converging fairly rapidly). However, my intuition tells me it is almost certain (I will not prove it here) that the error term oscillates between negative and positive infinitely often.

4 Under Riemann Hypothesis

If we assume the Riemann Hypothesis, the error term is bounded by $$\frac{C\log^2 x}{\sqrt{x}}$$ for some constant C. By analyzing the above data with numerical methods, the error seems to be best fitted by $\frac{C\log x}{\sqrt {x}}$.

Not sure if this is what you were looking for,

Taking logs,

We get $\sum (\log p - \log (p-2)) = \sum (\frac{2}{p} + \mathcal{O}(\frac{1}{p^2}))$

Now it is well known that $ \sum_{p \le n} \frac{1}{p} = \log \log n + A + \mathcal{O}(\frac{1}{\log^2n})$

And thus

Your product is asymptotically

$$K \log^2 n$$

for some constant $K$.

As you noticed in your question, you can compute the $K$ using the prime-zeta function (and knowing the $A$ above).

Do you know this paper by Rosser and Schoenfeld?

It contains a number of bounds that might be useful for you, including (3.29 in the paper) $$\prod_{p\le x} \frac{p}{p-1} < e^C (\log x)\left(1+\frac{1}{2 \log^2 x} \right)$$ for $286 \le x$, where $C$ is Euler's constant.

Their proof of this result might be useful, and/or other results in the paper might be.