Skip to main content
11 of 13
deleted 24 characters in body
martin
  • 1.9k
  • 11
  • 25

Zeros of the semiprimes

Let $P$ be the prime zeta function

\begin{aligned} {\displaystyle P (s) = \sum _ {p\, \in\mathrm {\, primes}} {\frac {1} {p^{s}}} = {\frac {1} {2^{s}}} + {\frac {1} {3^{s}}} + {\frac {1} {5^{s}}} + {\frac {1} {7^{s}}} + {\frac {1} {11^{s}}} + \cdots } \end{aligned}

and define the semiprime zeta function as

\begin{aligned} & \zeta_{\Omega_ 2} (s) = \exp \sum _k^\infty \frac{P (k s)^2+P(2 k s)}{2 k} \end{aligned}

When these zeros are plugged into

\begin{aligned} &\sum _{\gamma} -\frac{\cos (\gamma \log t)}{\ \sqrt{\gamma}} \end{aligned}

where $\gamma$ is the imaginary part of the zeros of $\zeta_{\Omega_ 2} (\frac{1}{2}+it)$, the sum presumably then should peak at the semiprimes, (just as peaks are formed at the primes when summed over the imaginary parts of the Riemann Zeta function).

But the zeros are hard to find. I calculated the near-zeros along the line $1/2 +it$ with

func[t_] := 1/Abs[Exp[(PrimeZetaP[(1/2 + I t)]^2 + PrimeZetaP[2 (1/2 + I t)])/2]];
data = Table[{t, func[t]}, {t, .01, #, .01}] &[(**)10^4];
peaks = Pick[data, PeakDetect[data[[;; , 2]], .01, .001], 1];

in Mathematica, which takes a long time to evaluate, so I uploaded the data here:

peaks = ToExpression[Import["https://raw.githubusercontent.com/martinq321/peaks/main/omega2"]];

Below is a plot of the absolute value of the first term of $\zeta_{\Omega_ 2} (s)$

$\left| \exp \left(P \left (1/2+ i t \right)^2/2+P \left(1+ 2i t \right)/2\right)\right|$ with the red lines marking near-zero points:

enter image description here

Given the above, the plot below shows

\begin{aligned} &\sum _{k \in \text{peaks}} -\frac{\cos (k \log t)}{\ \sqrt{k}}, \quad 200 \geq t \geq 250 \end{aligned}

where, despite the noise, the peaks show at the semiprimes and the primes:

enter image description here

and with the noise cleared up a bit:

enter image description here

li2 = Select[peaks, #[[2]] > 5 &][[All, 1]]; x1 = 200; x2 = 250;
data1 = Table[(Sum[-Cos[li2[[p]] Log[t]]/Sqrt@p, 
{p, Length@li2}]), {t, x1, x2, .1}];
ListLinePlot[BilateralFilter[data1, 1., 1, MaxIterations -> 25], 
PlotStyle -> {AbsoluteThickness[3]}, DataRange -> {x1, x2}, 
ImageSize -> 700, Frame -> True, AspectRatio -> 1/8]

What is going on here? I presume the peaks also showing at the primes because many of the 'zeros' are very close to the zeros of the zeta function. Is there a way to filter them out? Is there a better way to find the zeros of $\zeta_{\Omega_ 2} (s)$?

martin
  • 1.9k
  • 11
  • 25