Suppose $n\geq 1$ is a positive integer. Can we find an asymptotic formula for this product below.
$$\prod_{k=1}^{\infty}\zeta (2kn)=\zeta (2n)\zeta (4n)\zeta (6n) \cdots$$
I tried to use $\zeta (2n)=\dfrac{(-1)^{n+1}B_{2n}(2\pi )^{2n}}{2(2n)!}$ but couldn't get anywhere.
This is not an answer since just based on numerical simulation.
Being stuck with any formal approach I tried, I just used numerical simulation and what I observed is that, if $$P_n=\prod_{k=1}^{\infty}\zeta (2kn)$$ then $\log(P_n-1)$ varies as a linear function of $n$. A basic linear regression $(1\leq n\leq 50)$ gave $$\log(P_n-1)=0.135417-1.390256 n$$ both parameters being statistically significant.
Edit
Inspired by Igor Rivin's and Raymond Manzoni's answers.
Limiting the $\zeta$ function to the first term $(\zeta(s)=1+\frac 1 {2^s})$, it seems that, for large $n$, we can approximate $$P_n\approx 1+\frac 1 {4^n}$$ since $$\log(P_n)=\sum_{k=1}^\infty \log(\zeta(2kn))\approx\sum_{k=1}^\infty \log(1+\frac 1{2^{2kn}})\approx\sum_{k=1}^\infty \frac 1{2^{2kn}}=\frac{1}{4^n-1}$$ $$P_n\approx \exp(\frac{1}{4^n-1})\approx 1+\frac{1}{4^n}$$ $$\log(P_n-1)\approx -n \log(4)$$ ($\log(4)\approx 1.38629$ being quite close to the $1.390256$ from the empirical calculation).