Fractal Curves from Prime Trigonometric Series

Abstract

We study the convergence of the parameter family of series:Vα,β(t)=∑pp−αexp(2πipβt),α,β∈R>0,t∈[0,1)defined over prime numbers p and, subsequently, their differentiability properties. The visible fractal nature of the graphs as a function ofα,βis analyzed in terms of Hölder continuity, self-similarity and fractal dimension, backed with numerical results. Although this series is not a lacunary series, it has properties in common, such that we also discuss the link of this series with random walks and, consequently, explore its random properties numerically.

1. Introduction

The prime numbers are not randomly distributed but, there are random models that capture important properties of the distribution of prime numbers well (e.g., [1]). The random behavior of a deterministic mathematical object can be found elsewhere: there are classical function series that can be approximated by random processes. Let us briefly describe these series: consider the two functions _fn(x)=sin(2πnx) _fn(x)=sin(2πnx) and _fn+1(x)=sin(2π(n+1)x) _fn+1(x)=sin(2π(n+1)x) for an arbitrary integer n∈N n∈N . These behave as strongly dependent random variables if we consider x to be a random real variable uniformly distributed on some interval. However, if one picks from the sequence of frequencies _{(2πnx)n≥0 (2πnx)n≥0} a sub-sequence _{(2πnk)k≥0 (2πnk)k≥0} , such that the integer sequence grows sufficiently fast, i.e., _{nk+1 nk}≥1+ρ,ρ>0 _{nk+1 nk}≥1+ρ,ρ>0 , the quantities _fnk (x) _fnk (x) and _fnk+1 (x) _fnk+1 (x) behave like independent random variables (see Figure 1 as an example and Section 3).

Figure 1. Graph of (a): sin(2πnkx) sin(2πnkx) with nk=24 nk=24 and nk=25; nk=25; and (b): sin(2πnx) sin(2πnx) with n=5 n=5 and n=6 n=6 (displayed for x∈[0,1) x∈[0,1) at 105 105 points).

Now, one can construct a random walk out of these random variables: start at zero. At time k, move _fnk (x) _fnk (x) up. At time N, we find ourselves at S(x,N)=_{∑k=0N fnk} (x) S(x,N)=_{∑k=0N fnk} (x) . This sum is displayed for N=1000 N=1000 in Figure 2 on the left.

Figure 2. Graph of (a): ∑n=110sin(2nπx) ∑n=110sin(2nπx) ; (b): ∑n=11000sin(2πnx) ∑n=11000sin(2πnx) ; (c): ∑p≤1000sin(2πpx) ∑p≤1000sin(2πpx) (displayed for x∈[0,1) x∈[0,1) at 104 104 points).

This example is known as a lacunary Fourier series, that is its frequencies fulfill the growth condition given above. Its random properties are a classical field of research. In the literature, the sequence of prime numbers _{(2πpk)k≥0 (2πpk)k≥0} is often cited as a counterexample for a sequence of frequencies that does not give rise to a lacunary Fourier series: it neither fulfills the growth condition nor alternative conditions on arithmetic patterns that exist in the literature. However, experiments in this article suggest that _∑ksin(2π_pkx) _∑ksin(2π_pkx) share many of the random properties of lacunary series (see Figure 2 for a first impression or [2]): for instance, the central limit theorem seems to hold. Unfortunately, this looks difficult to prove (see e.g. [3], p. 2).

On the other hand, we can look at other manifestations of randomness in lacunary series (e.g., in the example in Figure 2) and try to see if they are also present in our prime series _{Vα,β Vα,β} : by introducing appropriate coefficients _{ak ak} , the walk _{∑k ak}sin(π_nkx) _{∑k ak}sin(π_nkx) can be approximated by a Wiener process, which is an almost everywhere continuous random walk with independent normally-distributed increments (see [4] and Section 3). This implies directly many interesting properties for the series, e.g., the law of iterated logarithm holds. It would be interesting if a similar approximation exists for our series _{Vα,β Vα,β}. Again, we were not able to prove this.

However, we can show that our series _{Vα,β Vα,β} has in fact for specific α,β α,β properties in common with a Wiener process, e.g., its regularity and fractality (see Section 2). The above-mentioned example _{∑k ak}sin(^2kπx) _{∑k ak}sin(^2kπx) is in fact famous for these reasons: it belongs to the family of Weierstrass functions _Fa,b(x)=_∑n=0∞ ^ansin(^bnt) _Fa,b(x)=_∑n=0∞ ^ansin(^bnt) , which have been extensively studied for its differentiability properties. Under certain conditions on a,b a,b, this function is nowhere differentiable, but Hölder continuous.

Another historical example that is non-differentiable, but multifractal, is the Riemann function _R2(x)=_∑n=1inf ⁿ⁻²sin(ⁿ²x) _R2(x)=_∑n=1inf ⁿ⁻²sin(ⁿ²x) . Note, that it is not a lacunary series as ^{(n+1)2 n2}→1 ^{(n+1)2 n2}→1. With our prime series, we place ourselves in between these two historical examples with respect to the growth of its frequencies.

While prime sums are extensively studied in the context of the famous prime conjectures (e.g., for Vinogradov’s theorem and the like), we have not found a treatment of trigonometric series over prime numbers. The reason for this is most probably that these series do not have the necessary form to help to progress in the proofs of the prime conjectures where prime exponential sums play a dominant role. As mentioned above, these series have not been studied in the context of lacunary series as prime numbers neither grow fast enough nor have known arithmetic properties, which are necessary for a straightforward analysis.

By using the results of prime number theory, we are nevertheless able to show conditions on the differentiability and self-similarity of our prime series. Experimentally, we explore also its box dimension in dependence of α,β α,β.

Remark 1.

For most of our questions, we can restrict ourselves, without loss of generality, to the real part _∑p ^p−αcos(2π^pβt) _∑p ^p−αcos(2π^pβt) of the series, which we denote by _Vα,β(t) _Vα,β(t) , as well.

2. Convergence and Differentiability

There are basically two factors that influence the smoothness and convergence of a function series _{∑k ak}exp(2πi_nkt) _{∑k ak}exp(2πi_nkt)as ours:

1. The faster the coefficients _{ak ak} decrease for k→∞ k→∞, the smaller is the influence of the higher frequencies. This implies that the series converges better and the resulting function is smoother.

2. The faster the frequencies _{nk nk}increase or, equivalently, the greater the gaps, the smaller the period of the oscillation becomes, so that one obtains more peaks and sinks in one interval, which increases the fractal character.

2.1. Historical Remarks

The nature of these influences, easily deduced, are also backed by the long history of studies on the following two families of functions (and derived families):

Let:

_Fa,b(t)=∑n=0∞^ancos(^bnt)

be the family of Weierstrass functions, which have been extensively studied. One knows the following:

Theorem 1

([5,6]). 1. If 0<ab<1,a<1,b>1 0<ab<1,a<1,b>1 , then _{Fa,b Fa,b} is differentiable.

2. If 0<a<1<ab 0<a<1<ab , then _{Fa,b Fa,b} is nowhere differentiable. Further, the Hölder exponent is a constant function s=−logalogb s=−logalogb , i.e., for all t,_t0 t,_t0 , it holds:

_Fa,b(t)−_Fa,b(_t0)≤C^t−_t0s.

On the other hand, one has the family of Riemann’s functions (whose authorship by Riemann is apparently only confirmed by Weierstrass) defined by:

_Rα(x)=∑n=1∞^n−αcos(ⁿ²x)

which has the following proven properties:

Theorem 2

([5,7,8,9,10]). 1. If 0<α≤12 0<α≤12 , then the series is not a Fourier series of an ^{L1 L1} -function. If 0<α<12 0<α<12 , then _{Rα Rα} converges at x if and only if x=aq x=aq , where a,q a,q are coprime and four divides q−2 q−2 .

2. If 12<α<1 12<α<1 , then the series converges in p-norm to a ^{Lp Lp} -function for p<21−α p<21−α .

3. If α=1 α=1 , then the series has bounded mean oscillation.

4. If α<52 α<52 , then _{Rα Rα} is not differentiable at any irrational value of x, and its Hausdorff dimension for 32≤α≤52 32≤α≤52 is equal to:

_dimH(_Rα)=94−α2.

If α=2 α=2 , then _{R2 R2} is differentiable at x if and only if x=aq x=aq where a,q a,q are coprime and four divides q−2 q−2 .

5. If α=2 α=2 , the Hölder exponent is discontinuous everywhere. In fact, _{R2 R2} is a function with unbounded variation and multifractal.

In the following, we aim to give a similar description for our function series. Let us start with some preliminary definitions, which are necessary for what follows.

2.2. Preliminary Definitions

We call a function f:R→C f:R→C locally Hölder continuous at _x0∈R _x0∈R , if there exist s∈(0,1] s∈(0,1] and C,ϵ>0 C,ϵ>0, such that:

f(x)−f(y)≤C^x−ys,for allx,y∈_Bϵ(_x0).

We call the supremum of s for which these inequality holds at _{x0 x0}the local Hölder exponent.

Let ϕ:R→C ϕ:R→C be a smooth function with compact support supp(ϕ)⊂C supp(ϕ)⊂C. We write:

ϕ^(u)=_∫Rϕ(t)exp(−iut)dt

for the Fourier transform of ϕ ϕ . Further, let ϕ:R→R ϕ:R→R be given such that the support of its Fourier transform is contained in [−1,1] [−1,1] , then the Gabor wavelet transform of a function f:R→C f:R→Cis defined by:

G(a,b,λ)=1a_∫Rf(t)exp(−iλt)ϕt−badt.

With these notation, we have the following estimation, which is a special case of Proposition 5 in [6]:

Proposition 1 (Jaffard).

Let f:R→C f:R→C be a bounded function. Let G(a,b,λ) G(a,b,λ) be the Gabor wavelet transform of f. If f is locally Hölder continuous at _x0∈R _x0∈R with Hölder coefficient s, then there exists C>0 C>0 such that for all a∈(0,1] a∈(0,1] and for all b∈_B1(_x0)¯ b∈_B1(_x0)¯ and for all λ≥^a−1 λ≥^a−1 , we have:

G(a,b,λ)≤C^{as 1+_x0−bas}.

2.3. Differentiability of _{Vα,β Vα,β}

In the spirit of the results in Section 2.1, we aim to determine which conditions have to be fulfilled by the coefficients and frequencies in our example in order to have a certain degree of differentiability. Firstly, we consider:

_Vβ(n,t)=∑p≤nf(p)cos(2π^pβt),β>0,

where f is any function of prime numbers. We can state the trivial fact that:

Proposition 2.

For any β≥0 β≥0 , if _∫2∞f(x)ln(x)dx<∞ _∫2∞f(x)ln(x)dx<∞ , then the partial sums _Vβ(n,t) _Vβ(n,t) converge uniformly and absolutely to a continuous function denoted by _{Vβ Vβ} .

Proof.

We have f(p)cos(2π^pβt)≤f(p) f(p)cos(2π^pβt)≤f(p) for all p. By the Weierstrass M-test, the partial sums _Vβ(n,t) _Vβ(n,t) converge uniformly and absolutely if _∑pf(p)<∞. _∑pf(p)<∞.Using the Riemann–Stieltjes integral and the prime number theorem, we get:

∑pf(p)=_∫2∞f(x)dπ(x)=_∫2∞f(x)ln(x)dx,

where π(x) π(x) denotes the number of primes ≤x ≤x, finishing the proof. ☐

We take now

_Vα,β(n,t)=∑p≤n^p−αcos(2π^pβt)

and denote with _Vα,β(t) _Vα,β(t)its limit whenever it exists. Then, one can show the following statement:

Theorem 3.

Let α∈R α∈R and α>1 α>1 .

1. Then, the series _Vα,β(n,t) _Vα,β(n,t) converges uniformly and absolutely to a continuous function _Vα,β(t) _Vα,β(t) .

2. For m≥1 m≥1 , if further α−mβ>1 α−mβ>1 , then the function _Vα,β(t) _Vα,β(t) is ^{Cm Cm} , i.e., m-times continuously differentiable.

Proof.

For the first result, we use the properties of the prime zeta function P(α)=_∑p ^p−α P(α)=_∑p ^p−α : it converges absolutely for α>1,α∈R, α>1,α∈R, and diverges for α=1 α=1 (see, e.g., [11,12]). The coefficients ^{p−α p−α} are an upper bound for the terms ^p−αcos2π^pβt ^p−αcos2π^pβt . Consequently, the Weierstrass M-test implies that for α>1 α>1 and any t∈[0,1) t∈[0,1) , _Vα,β(n,t) _Vα,β(n,t) converges uniformly and absolutely to _Vα,β(t) _Vα,β(t) . As any partial sum _Vα,β(n,t) _Vα,β(n,t)is continuous, the limit is a continuous function, as well.

Secondly, for any n and t, we can differentiate the partial sums:

_Vα,β′(n,t)=−2π∑p≤n^p−α+βsin(2π^pβt).

This sequence of derivatives converges uniformly with the same argument as above for α−β>1 α−β>1 , so that one concludes that _Vα,β(t) _Vα,β(t) is continuously differentiable itself with derivative _Vα,β′(t)=−2π_∑p ^p−α+βsin(2π^pβt) _Vα,β′(t)=−2π_∑p ^p−α+βsin(2π^pβt). By induction over m, one proves the m-time differentiability of the function. ☐

Remark 2.

The result is in accordance with the intuitive smoothness of the series: for fixed α>1 α>1 , the series becomes smoother, the smaller the frequency ^{pβ pβ} , β→0 β→0 , or equivalently, the larger the period. Therefore, the peaks and sinks of the oscillation are more and more separated so that the series becomes smoother (see Figure 3, Figure 4 and Figure 5).

Figure 3. Graph of V1.5,2(105,t) V1.5,2(105,t) at 5×104 5×104 discrete points in each direction (interpolated).

Figure 4. Graph of V1.5,1.5(105,t) V1.5,1.5(105,t) at 5×104 5×104 discrete points in each direction (interpolated).

Figure 5. Graph of V1.5,1(105,t) V1.5,1(105,t) at 5×104 5×104 discrete points in each direction (interpolated).

Theorem 4.

If 1<α≤β+1 1<α≤β+1 , then the function is Hölder continuous with Hölder coefficient s≤αβ s≤αβ .

Proof.

First of all, let f:R→C f:R→C be an integrable function and N>0 N>0 , then by using the Riemann–Stieltjes integral (see, e.g., [13]) and the prime number theorem as above, one knows:

∑p≤N^p−α=_∫2N ^x−αdπ(x)∼_∫2N1^xαln(x)dx.

From this formula and li(x) li(x) denoting the logarithmic integral function, one deduces (substituting dx dx by d^x1−α d^x1−α ) for α<1 α<1:

∑p≤N^p−α=^(1−α)−1li^N1−α+O^N1−αexp(−cln(N).

Approximating the logarithmic integral, this implies:

∑p≤N^p−α∼^N1−α(1−α)ln(N).

If α>1 α>1 , we have to use the explicit formula for the prime zeta function to get an estimate for the speed of convergence (see, e.g., [14] for a derivation of the formula). We then have by partial summation:

∑p^p−α=∑p≤N^p−α+∑n=1∞μ(n)nlnζ(N,αn),withζ(N,α)=ζ(α)_Πp≤N1−^p−α,

where ζ(α)=_∑n=1∞ ^n−α ζ(α)=_∑n=1∞ ^n−α denotes the Riemann zeta function and μ μthe Moebius function. Therefore, we get for the tail of the prime zeta function:

∑p>N^p−α=∑n=1∞μ(n)nlnζ(N,αn),withlnζ(N,α)=O^N−αand∑n=1∞μ(n)n=0.

Combining Equations (1) and (2) on the asymptotic of the prime zeta function, we can estimate now the regularity of our function _Vα,β(t) _Vα,β(t).

For any t,_t0∈[0,1) t,_t0∈[0,1) , we choose N=^{t−_t0−1α}. N=^{t−_t0−1α}.Then, we have with the mean value theorem and using the absolute convergence of the series:

_Vα,β(t)−_Vα,β(_t0)≤∑p≤N^p−αcos(2π^pβt)−cos(2π^pβ _t0)+2∑p>N^p−α≤∑p≤N^p−α+βt−_t0+2∑p>N^p−α≤^N−α+β+1(β−α+1)ln(N)t−_t0+2C^N−α≤C^{t−_t02−β+1α}.

The exponent 1−β<2−β+1α≤1 1−β<2−β+1α≤1is not necessarily optimal, but a lower bound. However, it suffices to conclude that the function is Hölder continuous, so that we can derive an upper bound for its Hölder exponent.

For this step, we use a method developed by Jaffard in [6], which relies on a wavelet transform and the idea of choosing the wavelet transform such that only one frequency of _Vα,β(t) _Vα,β(t) is picked up. Let _θm=min_pmβ−_pm−1β,_pm+1β−_{pmβ θm}=min_pmβ−_pm−1β,_pm+1β−_pmβ and _Δm=_pm−_{pm−1 Δm}=_pm−_pm−1.

We choose a function ϕ ϕ whose Fourier transform ϕ^ ϕ^ has compact support supp(ϕ^)⊂[0,1] supp(ϕ^)⊂[0,1] and ϕ^(0)=1 ϕ^(0)=1. We then look at the Gabor-wavelet transform:

_Gm(_θm−1,_t0,_pmβ)=_θm∑k_{pk−α ∫R}expi_pkβ−_pmβtϕ_θm(t−_t0)dt=∑k_pk−αexpi_pkβ−_pmβt0∫Rexpi_pkβ−_pmβθmϕ(u)du,

with u=_θm(t−_t0) u=_θm(t−_t0) . Substituting ϕ^(y)=_∫Rexp(iyu)ϕ(u)du ϕ^(y)=_∫Rexp(iyu)ϕ(u)du for y=_pkβ−_pmβu_θm y=_pkβ−_pmβu_θmin the equation, we get:

_Gm(_θm−1,_t0,_pmβ)=∑k_pk−αexpi_pkβ−_pmβt0ϕ^_pkβ−_pmβu_θm.

As the support of ϕ^ ϕ^ is a subset of the unit interval, it vanishes for any k≠m k≠m, so the expression is reduced to:

_Gm(_θm−1,_t0,_pmβ)=_pm−α

Recall that we have just proven that _{Vα,β Vα,β} is locally Hölder continuous at _t0∈R _t0∈R . Further, for all m, it is _pmβ≥_{θm pmβ}≥_θm and _θm−1∈(0,1] _θm−1∈(0,1] . Hence, applying Proposition 1, there exists C>0 C>0 such that for all s∈(0,1) s∈(0,1):

_Gm(_θm−1,_t0,_pmβ)=_pm−α≤C_θm−s.

The gap _{θm θm} is bounded by _{pmβ pmβ} from above, so that the Hölder coefficient s is bounded by αβ αβfrom above, finishing the proof. ☐

Remark 3.

Let α>1 α>1 be fixed. The bigger the gaps of the frequency, β→∞ β→∞ , the stronger the irregularity of _Vα,β(t) _Vα,β(t) .

2.4. Self-Similarity and Fractal Dimension

The graph of the function _{Vα,β Vα,β} seems to be self-similar for certain α,β α,β . There seems to be an approximate scalar invariance at points ^{q−1 q−1} , where q is prime. Let us make this intuition more precise: look for example at the partial sums _V1,1(n,t)=_∑p≤n ^p−1exp(2πipt) _V1,1(n,t)=_∑p≤n ^p−1exp(2πipt) in Figure 6.

Figure 6. Graph of V1,1(105,t) V1,1(105,t) at 5×104 5×104 discrete points.

Denote by _{pk pk} the k-th prime number. We restrict ourselves again to the real part of _V1,1(n,t) _V1,1(n,t) . The point 12 12 is a global minimum as _V1,1′n,12=0 _V1,1′n,12=0 and _V1,1n,12=12−_{∑k=1n pk−1 V1,1}n,12=12−_{∑k=1n pk−1} as the primes greater than two are odd. Now, consider the point 13 13 : we have _V1,1n,13=13−12_{∑k=1,k≠2n pk−1 V1,1}n,13=13−12_{∑k=1,k≠2n pk−1}. More generally, one has:

_V1,1n,1q=1q+∑l=1q−1cos2πlq∑_pk=lmodqn_pk−1,qprime=∑l=0q−1cos2πlq_Rl,q

That is, we can decompose the partial sum into residue classes of the prime numbers and the roots of unity of cosine. One knows that the number of primes p≤n p≤n that are congruent to lmodq lmodq are approximately the same for all l, that is nΦ(q)log(n) nΦ(q)log(n) , where Φ(q) Φ(q) denotes the Euler totient function and is equal to q−1 q−1 for q prime. Therefore, for any 1q 1q , q prime, one can use this distribution and the Riemann–Stieltjes integral to show that the difference between the sums _{∑pk=lmodq,pk≤n pk−1 ∑pk=lmodq,pk≤n pk−1} for each l=1,…,q−1 l=1,…,q−1 converges to zero for n→∞ n→∞, that is:

_Rl,q=∑_pk=lmodq,_pk≤n_pk−1∼1q−1_∫2n1xln(x)dx=1q−1lnln(n)+C.

The factors cos2πlq cos2πlq are exactly the prime roots of unity, and the sum _∑l=0q−1cos2πlq=0 _∑l=0q−1cos2πlq=0. Consequently, one computes:

_V1,1n,1q∼1q−1q−1lnln(n)+C.

As we have _V1,1(n,1)=_∑p≤n ^p−1∼lnln(n)+M _V1,1(n,1)=_∑p≤n ^p−1∼lnln(n)+M, one could argue that:

_V1,1(n,tq)≈11−q_V1,1(n,t)+1q,q≥3,prime.

See Figure 7. However, keep in mind that these are only asymptotic equivalences, while our partial sum _V1,1(n,t) _V1,1(n,t) does not converge for all t∈[0,1) t∈[0,1) for n→∞ n→∞, so the self-similarity of the graph is certainly not strict.

Figure 7. The graph of the real part of −12V1,1(106,t) −12V1,1(106,t) (in black) and V1,1(106,t/3)+13 V1,1(106,t/3)+13 (in green).

Fractal Dimension of _{Vα,β Vα,β}

Further, we compute numerically the box dimension of the graph of _{Vα,β Vα,β} defined in the following way: let A:=[a,b]×[c,d] A:=[a,b]×[c,d] be the rectangle such that the graph _Vα,β(n,t)⊂A _Vα,β(n,t)⊂A is contained. We compute then for i,j=0,…N−1 i,j=0,…N−1 the intersections _Vα,β(n,t)∩[a+i(b−a)/N,a+(i+1)(b−a)/N]×[c+j(d−c)/N,c+(j+1)(d−c)/N] _Vα,β(n,t)∩[a+i(b−a)/N,a+(i+1)(b−a)/N]×[c+j(d−c)/N,c+(j+1)(d−c)/N] . We denote the number of non-empty intersections by M(N) M(N). The box dimension is then given by:

_dimBVα,β(n,t)=limN→∞ln(M(N))ln(N).

In accordance with our results on the regularity of _{Vα,β Vα,β} , we obtain the following Figure 8 for the (numerically computed) box dimension _{dimB dimB} over the fraction αβ αβ . For α>1 α>1 fixed and β→∞ β→∞ , that is αβ→0 αβ→0 , we expect that the fractal dimension converges to two. On the other hand, for β→0 β→0 , the fractal dimension should converge to one as the graph gets continuously differentiable if α−1β>1+1β α−1β>1+1β.

Figure 8. (a): Box dimension for the graph of Vα,β Vα,β (computed at 105 105 points) in dependence of the fraction of the powers αβ αβ with α∈[1,1.5] α∈[1,1.5] and β∈[0.5,3] β∈[0.5,3] . Remark that Vα,β Vα,β is not convergent for α=1 α=1 . (b): example for α=1.5, β=2 α=1.5, β=2 to show how the box dimension was numerically approximated.

Remark 4.

In recent literature (see, e.g., [15]), the concept of the fractal dimension of prime distribution is studied: it is defined with the help of an indicator function S:N→0,1 S:N→0,1 on the natural numbers, which is one if n∈N n∈N is prime, otherwise zero. The fractal dimension is then computed as the average number of ones in randomly-chosen minors n×n n×n of the N×N N×N -matrix S(n≤N)×S(n≤N) S(n≤N)×S(n≤N) . In a recent publication ([16]), the self-similarity of these images (created by displaying the ones in the matrix in black, zeros in white) is also treated. These concepts are not related to the self-similarity or the fractal dimension of our series: the fractality of _{Vα,β Vα,β} depends on its Hölder exponent, which relies on the exponents α,β α,β and on the gaps of frequencies (^pβ) (^pβ) . In this sense, it certainly depends on the prime distribution, but we do not see a direct application of the concepts cited above. More probable is that in order to analytically compute the fractal dimension, one would have to use the large sieve inequality where on the right side would appear values of _∑p1^p2α _∑p1^p2α .

There is clearly no reverse correlation from the fractal dimension of the graph to the prime distribution as there are other, even multifractal curves as the cited Riemann function whose frequencies are not prime sequences.

3. Random Properties for _{Vα,β Vα,β}

The quite similar behavior of lacunary and random Fourier series allow us to think that it might be possible to capture the random character of the series _{Vα,β Vα,β}, which is the subject of this section. Let us briefly review what is known in the context of lacunary sequences and random variables.

3.1. Lacunary Sequences Behaving as Independent Random Variables: Short Overview

The terms (sin_(2πkx)k (sin_(2πkx)k and (cos_(2πkx)k (cos_(2πkx)k behave like random variables, but strongly dependently. However, if one restricts the sequence of frequencies _{(2πk)k≥0 (2πk)k≥0} to _{(2πnk)k≥0 (2πnk)k≥0} where the sequence _{(nk)k≥0 (nk)k≥0}has sufficiently fast-growing gaps, i.e.,

_{nk+1 nk}≥1+ρ,ρ>0(Hadamard gap condition),

then the sequences sin(2π_nkx sin(2π_nkxbehave like independent random variables. For example, one has:

1N∑k=1Nsin(2π_nkx)→N(0,1),

where N(0,1) N(0,1) is the normal distribution. This was the main observation that has led to study the connections between lacunary and random Fourier series, most importantly the question of which are the optimal growth conditions on the sequence _{(nk)k (nk)k} such that the sequence _{f(nk)k≥0 f(nk)k≥0} for general periodic measurable functions f with vanishing integral exhibits random properties (see the historical overview in [17]). By introducing weights _{ak ak} that obey certain growth conditions themselves, one can recover several limit theorems in complete analogy with random variables. In particular, the Central Limit Theorem (CLT) and the Law of Iterated Logarithm (LIL) are true (see the results by Salem-Zygmund in [18,19], Erdös-Gál in [20] and Weiss in [21]). Further, it can be shown that the process can be approximated by a standard Brownian motion:

Theorem 5

(Philipp-Stout [4]). Assume the Hadamard gap condition. Assume further that _AN:=12_{∑k=1N ak2}→∞ _AN:=12_{∑k=1N ak2}→∞ and there exists δ>0 δ>0 such that _{limN→∞aN AN1−δ}=0 _{limN→∞aN AN1−δ}=0 . Then, without changing the distribution of the process:

S(t,x)=∑k≤t_akcos(2π_nkx),t≥0,

it can be redefined on a suitable probability space together with a Wiener process W(t)|t≥0 W(t)|t≥0 such that:

S(t,x)=W_At+O_At12−ρ,almost surely for someρ>0.

While the Hadamard growth condition (4) for CLT can be weakened for general sequences _{(nk)k (nk)k} (see [22]) for coefficients _ak=1 _ak=1 to the optimal growth condition _{nk+1 nk}≥1+_ckk _{nk+1 nk}≥1+_ckk with _ck→∞ _ck→∞ , one has observed that sequences with much slower growth can nevertheless satisfy the CLT if they fulfill certain arithmetic conditions, more precisely bounds on the number of solutions for the diophantine equation. Results in this direction started with Gaposhkin [23] and were recently sharpened by Berkes, Philipp and Tichy [24]. The difficulties with the prime sequence _{pkk≥0 pkk≥0} are on both sides: Firstly, while it is sure that the prime sequence is not a Hadamard sequence, neither precise lower, nor upper bounds for the prime gap _pk+1−_{pk pk+1}−_pk are known. The best results for a lower bound that would be of interest for us do not hold for all k≥0 k≥0 , but only infinitely many. For the upper bound, it is proven by Goldston, Pintz and Yildirim ([25]) that lim_{infk→∞Δk}log_pk=0. lim_{infk→∞Δk}log_pk=0. Secondly, there is no building law for prime numbers known, and the infinite recurrence of certain patterns like twin primes is only conjectured, but not completely proven. On the other hand, the random character of prime numbers is often invoked without being analytically established anywhere although the random model by Cramér [1] is widely used and reproduces some results very efficiently (but fails in other aspects, e.g., in forecasting the size of the prime gap). For the question on the convergence of functions f(_nk) f(_nk) , random models were also introduced (see, e.g., [26]). Obviously, this is a broad and intensively-studied mathematical subject where we do not dare to make contributions. Therefore, we stay more closely to our studied series.

3.2. The Central Limit Theorem

Because of the reason mentioned above, we have not been able to show the central limit theorem for the random variables sin(_pkx) sin(_pkx) or cos(_pkx) cos(_pkx) , the base of our series _{Vα,β Vα,β} . Nevertheless, numerical computations strongly suggest that the central limit theorem holds; see Figure 9: we took ^{104 104} uniformly-distributed points x of the interval −π,π −π,π and computed the sample average 1N_∑k=1Nsin(_pkx) 1N_∑k=1Nsin(_pkx) for N= N= 78,498, that is, the number of primes ≤¹⁰⁶ ≤¹⁰⁶ . We computed the histogram for the values of the sample average, which experimentally tends to a normal distribution as the size of the sample tends to infinity. To confirm this observation, we did the same computation for sin(_pk32x) sin(_pk32x) , see Figure 10.

Figure 9. Normal distribution of 1N∑k=1Nsin(pkx) 1N∑k=1Nsin(pkx) for x uniformly distributed in −π,π −π,π .

Figure 10. Normal distribution of 1N∑k=1Nsin(πpk32x) 1N∑k=1Nsin(πpk32x) for x uniformly distributed in −π,π −π,π .

4. Concluding Remarks

The properties of the series _{Vα,β Vα,β}that we have discussed in this article are intimately related to the distribution of prime numbers, and this was mostly due to the unanswered questions on prime numbers for which the analytical access to our series is limited. Therefore, knowledge about the distribution and bounds for the gaps of prime numbers would imply more or less directly the properties where we were restricted to a numerical approach.

Although the series might be reminiscent of the Riemann zeta function or other number-theoretical functions, we did not construct _{Vα,β Vα,β} in this way and do not see a possibility to deduce it from any of them, besides from the trivial fact, that _Vα,β(0) _Vα,β(0) is equal to the prime zeta function P(α)=_∑p ^p−α P(α)=_∑p ^p−α . Furthermore, recall that we only consider α∈R α∈R , usually α>1 α>1, so that no zeros of the Riemann zeta function come into the play for the scope of this article.

Acknowledgments

We would like to thank Juri Merger and Florian Pausinger for his helpful hints and questions.

Author Contributions

D.V. conceived of and designed the experiments. D.V. and D.B. performed the experiments. D.V. and D.B. analyzed the data. D.B. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

1. Cramér, H. On the order of magnitude of the difference between consecutive prime numbers. Acta Arith. 1936, 2, 23–46.

2. Vartziotis, D.; Wipper, J. The fractal nature of an approximate prime counting function. Fractal Fract. 2017, 1, 10.

3. Raseta, M. On the strong approximations of partial sums of f(nkx). arXiv, 2015.

4. Philipp, W.; Stout, W.F. Almost Sure Invariance Principles for Partial Sums of Weakly Dependent Random Variables; American Mathematical Society: Providence, RI, USA, 1975.

5. Hardy, G.H. Weierstrass’s Non-Differentiable Function. Trans. Am. Math. Soc. 1916, 17, 301–325.

6. Jaffard, S. Pointwise and directional regularity of nonharmonic Fourier series. Appl. Comput. Harmon. Anal. 2010, 22, 251–266.

7. Hardy, G.H.; Littlewood, J.E. Contributions to the arithmetic theory of series. Proc. Lond. Math. Soc. 1913, 2, 411–478.

8. Gerver, J. The differentiability of the Riemann function at certain rational multiples of π. Am. J. Math. 1970, 92, 33–55.

9. Gerver, J. More on the differentiability of the Riemann function. Am. J. Math. 1970, 93, 33–41.

10. Chamizo, F.; Córdoba, A. Differentiability and dimension of some fractal Fourier series. Adv. Math. 1996, 142, 335–354.

11. Landau, E.; Walfisz, A. Über die nichfortsetzbarkeit einiger durch dirichletsche reihen definierter funktionen. Rendiconti del Circolo Matematico di Palermo 1920, 44, 82–86. (In German)

12. Fröberg, C.E. On the prime zeta function. BIT 1968, 8, 187–202.

13. Rosser, J.; Schoenfeld, L. Approximate formulas for some functions of prime numbers. Ill. J. Math. 1962, 6, 64–94.

14. Cohen, H. High Precision Computation of Hardy-Littlewood Constants. 2000. Available online: https://www.math.u-bordeaux.fr/~hecohen/ (accessed on 1 January 2018).

15. Cattani, C. Fractal patterns in prime numbers distribution. In Lecture Notes in Computer Science, Proceedings of the International Conference on Computational Science and Its Applications (ICCSA 2010), Fukuoka, Japan, 23–26 March 2010; Taniar, D., Gervasi, O., Murgante, B., Pardede, E., Apduhan, B.O., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; Part II, pp. 164–176.

16. Cattani, C.; Ciancio, A. On the fractal distribution of primes and prime-indexed primes by the binary image analysis. Phys. A Stat. Mech. Its Appl. 2016, 460, 222–229.

17. Kahane, J.P. A century of interplay between taylor series, Fourier series and brownian motion. Bull. Lond. Math. Soc. 1997, 29, 257–279.

18. Salem, R.; Zygmund, A. On lacunary trigonometric series. Proc. Natl. Acad. Sci. USA 1947, 33, 333–338.

19. Salem, R.; Zygmund, A. On lacunary trigonometric series, II. Proc. Natl. Acad. Sci. USA 1948, 34, 54–62.

20. Erdös, P.; Gál, I. On the law of iterated logarithm I + II. Indag. Math. (Proceedings) 1955, 58, 65–84.

21. Weiss, M. The law of the iterated logarithm for lacunary trigonometric series. Trans. Am. Math. Soc. 1959, 91, 444–469.

22. Erdös, P. On trigonometric sums with gaps. Magyar Tud. Akad. Mat. Kutató Int. Közl 1962, 7, 37–42.

23. Gaposhkin, V.F. Lacunary series and independent functions. Uspekhi Matematicheskikh Nauk 1966, 21, 3–82.

24. Berkes, I.; Philipp, W.; Tichy, R. Metric discrepancy results for sequences {nkx} and diophantine equations. In Diophantine Approximation: Festschrift for Wolfgang Schmidt; Springer: Vienna, Austria, 2008; pp. 95–105.

25. Goldston, D.A.; Pintz, J.; Yildirim, C.Y. Primes in tuples I. Ann. Math. 2009, 170, 819–862.

26. Schatte, P. On a law of iterated logarithm for sums mod 1 with applications to Benford’s law. Probab. Theory Relat. Fields 1988, 77, 167–178.