We study the convergence of the parameter family of series: V α , β ( t ) = p p α exp ( 2 π i p β 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.

trigonometric series lacunary series Hölder continuity fractality random Fourier series prime distribution
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., ). 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 f n ( x ) = sin ( 2 π n x ) and f n + 1 ( x ) = sin ( 2 π ( n + 1 ) x ) for an arbitrary integer 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 π n x ) n 0 a sub-sequence ( 2 π n k ) k 0 , such that the integer sequence grows sufficiently fast, i.e., n k + 1 n k 1 + ρ , ρ > 0 , the quantities f n k ( x ) and f n k + 1 ( x ) behave like independent random variables (see Figure 1 as an example and Section 3).

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

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 π p k ) 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 k sin ( 2 π p k x ) share many of the random properties of lacunary series (see Figure 2 for a first impression or ): for instance, the central limit theorem seems to hold. Unfortunately, this looks difficult to prove (see e.g. , 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 α , β : by introducing appropriate coefficients a k , the walk k a k sin ( π n k x ) can be approximated by a Wiener process, which is an almost everywhere continuous random walk with independent normally-distributed increments (see  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 α , β . Again, we were not able to prove this.

However, we can show that our series 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 a k sin ( 2 k π x ) is in fact famous for these reasons: it belongs to the family of Weierstrass functions F a , b ( x ) = n = 0 a n sin ( b n t ) , which have been extensively studied for its differentiability properties. Under certain conditions on a , b , this function is nowhere differentiable, but Hölder continuous.

Another historical example that is non-differentiable, but multifractal, is the Riemann function R 2 ( x ) = n = 1 inf n 2 sin ( n 2 x ) . Note, that it is not a lacunary series as ( n + 1 ) 2 n 2 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 α , β .

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

2. Convergence and Differentiability

There are basically two factors that influence the smoothness and convergence of a function series k a k exp ( 2 π i n k t ) as ours:

The faster the coefficients a k decrease for k , the smaller is the influence of the higher frequencies. This implies that the series converges better and the resulting function is smoother.

The faster the frequencies n k 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: F a , b ( t ) = n = 0 a n cos ( b n t ) be the family of Weierstrass functions, which have been extensively studied. One knows the following:

([5,6]). 1. If 0 < a b < 1 , a < 1 , b > 1 , then F a , b is differentiable.

2. If 0 < a < 1 < a b , then F a , b is nowhere differentiable. Further, the Hölder exponent is a constant function s = log a log b , i.e., for all t , t 0 , it holds: F a , b ( t ) F a , b ( t 0 ) C t t 0 s .

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 ( n 2 x ) which has the following proven properties:

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

2. If 1 2 < α < 1 , then the series converges in p-norm to a L p -function for p < 2 1 α .

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

4. If α < 5 2 , then R α is not differentiable at any irrational value of x, and its Hausdorff dimension for 3 2 α 5 2 is equal to: dim H ( R α ) = 9 4 α 2 .

If α = 2 , then R 2 is differentiable at x if and only if x = a q where a , q are coprime and four divides q 2 .

5. If α = 2 , the Hölder exponent is discontinuous everywhere. In fact, R 2 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 locally Hölder continuous at x 0 R , if there exist s ( 0 , 1 ] and C , ϵ > 0 , such that: f ( x ) f ( y ) C x y s , for   all x , y B ϵ ( x 0 ) .

We call the supremum of s for which these inequality holds at x 0 the local Hölder exponent.

Let ϕ : R C be a smooth function with compact support supp ( ϕ ) C . We write: ϕ ^ ( u ) = R ϕ ( t ) exp ( i u t ) d t for the Fourier transform of ϕ . Further, let ϕ : R R be given such that the support of its Fourier transform is contained in [ 1 , 1 ] , then the Gabor wavelet transform of a function f : R C is defined by: G ( a , b , λ ) = 1 a R f ( t ) exp ( i λ t ) ϕ t b a d t .

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

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

2.3. Differentiability of <inline-formula> <mml:math id="mm88" display="block"> <mml:semantics> <mml:msub> <mml:mi>V</mml:mi> <mml:mrow> <mml:mi>α</mml:mi> <mml:mo>,</mml:mo> <mml:mi>β</mml:mi> </mml:mrow> </mml:msub> </mml:semantics> </mml:math> </inline-formula>

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 n f ( p ) cos ( 2 π p β t ) , β > 0 , where f is any function of prime numbers. We can state the trivial fact that:

For any β 0 , if 2 f ( x ) ln ( x ) d x < , then the partial sums V β ( n , t ) converge uniformly and absolutely to a continuous function denoted by V β .

We have f ( p ) cos ( 2 π p β t ) f ( p ) for all p. By the Weierstrass M-test, the partial sums V β ( n , t ) converge uniformly and absolutely if p f ( p ) < . Using the Riemann–Stieltjes integral and the prime number theorem, we get: p f ( p ) = 2 f ( x ) d π ( x ) = 2 f ( x ) ln ( x ) d x , where π ( x ) denotes the number of primes x , finishing the proof. ☐

We take now V α , β ( n , t ) = p n p α cos ( 2 π p β t ) and denote with V α , β ( t ) its limit whenever it exists. Then, one can show the following statement:

Let α R and α > 1 .

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

For m 1 , if further α m β > 1 , then the function V α , β ( t ) is C m , i.e., m-times continuously differentiable.

For the first result, we use the properties of the prime zeta function P ( α ) = p p α : it converges absolutely for α > 1 , α R , and diverges for α = 1 (see, e.g., [11,12]). The coefficients p α are an upper bound for the terms p α cos 2 π p β t . Consequently, the Weierstrass M-test implies that for α > 1 and any t [ 0 , 1 ) , V α , β ( n , t ) converges uniformly and absolutely to V α , β ( t ) . As any partial sum 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 , so that one concludes that V α , β ( t ) is continuously differentiable itself with derivative V α , β ( t ) = 2 π p p α + β sin ( 2 π p β t ) . By induction over m, one proves the m-time differentiability of the function. ☐

The result is in accordance with the intuitive smoothness of the series: for fixed α > 1 , the series becomes smoother, the smaller the frequency p β , β 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).

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

First of all, let f : R C be an integrable function and N > 0 , then by using the Riemann–Stieltjes integral (see, e.g., ) and the prime number theorem as above, one knows: p N p α = 2 N x α d π ( x ) 2 N 1 x α ln ( x ) d x .

From this formula and li ( x ) denoting the logarithmic integral function, one deduces (substituting d x by d x 1 α ) for α < 1 : p N p α = ( 1 α ) 1 li N 1 α + O N 1 α exp ( c ln ( N ) .

Approximating the logarithmic integral, this implies: p N p α N 1 α ( 1 α ) ln ( N ) .

If α > 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.,  for a derivation of the formula). We then have by partial summation: p p α = p N p α + n = 1 μ ( n ) n ln ζ ( N , α n ) , with ζ ( N , α ) = ζ ( α ) Π p N 1 p α , where ζ ( α ) = 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 ) n ln ζ ( N , α n ) , with ln ζ ( 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 ) .

For any t , t 0 [ 0 , 1 ) , we choose N = t t 0 1 α . Then, we have with the mean value theorem and using the absolute convergence of the series: V α , β ( t ) V α , β ( t 0 ) p N p α cos ( 2 π p β t ) cos ( 2 π p β t 0 ) + 2 p > N p α p N p α + β t t 0 + 2 p > N p α N α + β + 1 ( β α + 1 ) ln ( N ) t t 0 + 2 C N α C t t 0 2 β + 1 α .

The exponent 1 β < 2 β + 1 α 1 is 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 , which relies on a wavelet transform and the idea of choosing the wavelet transform such that only one frequency of V α , β ( t ) is picked up. Let θ m = min p m β p m 1 β , p m + 1 β p m β and Δ m = p m p m 1 .

We choose a function ϕ whose Fourier transform ϕ ^ has compact support supp ( ϕ ^ ) [ 0 , 1 ] and ϕ ^ ( 0 ) = 1 . We then look at the Gabor-wavelet transform: G m ( θ m 1 , t 0 , p m β ) = θ m k p k α R exp i p k β p m β t ϕ θ m ( t t 0 ) d t = k p k α exp i p k β p m β t 0 R exp i p k β p m β θ m ϕ ( u ) d u , with u = θ m ( t t 0 ) . Substituting ϕ ^ ( y ) = R exp ( i y u ) ϕ ( u ) d u for y = p k β p m β u θ m in the equation, we get: G m ( θ m 1 , t 0 , p m β ) = k p k α exp i p k β p m β t 0 ϕ ^ p k β p m β u θ m .

As the support of ϕ ^ is a subset of the unit interval, it vanishes for any k m , so the expression is reduced to: G m ( θ m 1 , t 0 , p m β ) = p m α

Recall that we have just proven that V α , β is locally Hölder continuous at t 0 R . Further, for all m, it is p m β θ m and θ m 1 ( 0 , 1 ] . Hence, applying Proposition 1, there exists C > 0 such that for all s ( 0 , 1 ) : G m ( θ m 1 , t 0 , p m β ) = p m α C θ m s .

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

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

2.4. Self-Similarity and Fractal Dimension

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

Denote by p k the k-th prime number. We restrict ourselves again to the real part of V 1 , 1 ( n , t ) . The point 1 2 is a global minimum as V 1 , 1 n , 1 2 = 0 and V 1 , 1 n , 1 2 = 1 2 k = 1 n p k 1 as the primes greater than two are odd. Now, consider the point 1 3 : we have V 1 , 1 n , 1 3 = 1 3 1 2 k = 1 , k 2 n p k 1 . More generally, one has: V 1 , 1 n , 1 q = 1 q + l = 1 q 1 cos 2 π l q p k = l mod q n p k 1 , q prime = l = 0 q 1 cos 2 π l q R l , 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 that are congruent to l mod q are approximately the same for all l, that is n Φ ( q ) log ( n ) , where Φ ( q ) denotes the Euler totient function and is equal to q 1 for q prime. Therefore, for any 1 q , q prime, one can use this distribution and the Riemann–Stieltjes integral to show that the difference between the sums p k = l mod q , p k n p k 1 for each l = 1 , , q 1 converges to zero for n , that is: R l , q = p k = l mod q , p k n p k 1 1 q 1 2 n 1 x ln ( x ) d x = 1 q 1 ln ln ( n ) + C .

The factors cos 2 π l q are exactly the prime roots of unity, and the sum l = 0 q 1 cos 2 π l q = 0 . Consequently, one computes: V 1 , 1 n , 1 q 1 q 1 q 1 ln ln ( n ) + C .

As we have V 1 , 1 ( n , 1 ) = p n p 1 ln ln ( n ) + M , one could argue that: V 1 , 1 ( n , t q ) 1 1 q V 1 , 1 ( n , t ) + 1 q , q 3 , prime .

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

Fractal Dimension of <inline-formula> <mml:math id="mm206" display="block"> <mml:semantics> <mml:msub> <mml:mi>V</mml:mi> <mml:mrow> <mml:mi>α</mml:mi> <mml:mo>,</mml:mo> <mml:mi>β</mml:mi> </mml:mrow> </mml:msub> </mml:semantics> </mml:math> </inline-formula>

Further, we compute numerically the box dimension of the graph of V α , β defined in the following way: let A : = [ a , b ] × [ c , d ] be the rectangle such that the graph V α , β ( n , t ) A is contained. We compute then for 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 ] . We denote the number of non-empty intersections by M ( N ) . The box dimension is then given by: dim B V α , β ( n , t ) = lim N ln ( M ( N ) ) ln ( N ) .

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

In recent literature (see, e.g., ), 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 on the natural numbers, which is one if n N is prime, otherwise zero. The fractal dimension is then computed as the average number of ones in randomly-chosen minors n × n of the N × N -matrix S ( n N ) × S ( n N ) . In a recent publication (), 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 α , β depends on its Hölder exponent, which relies on the exponents α , β and on the gaps of frequencies ( 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 p 1 p 2 α .

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 <inline-formula> <mml:math id="mm231" display="block"> <mml:semantics> <mml:msub> <mml:mi mathvariant="bold-italic">V</mml:mi> <mml:mrow> <mml:mi mathvariant="bold-italic">α</mml:mi> <mml:mo>,</mml:mo> <mml:mi mathvariant="bold-italic">β</mml:mi> </mml:mrow> </mml:msub> </mml:semantics> </mml:math> </inline-formula>

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 α , β , 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 π k x ) k and ( cos ( 2 π k x ) k behave like random variables, but strongly dependently. However, if one restricts the sequence of frequencies ( 2 π k ) k 0 to ( 2 π n k ) k 0 where the sequence ( n k ) k 0 has sufficiently fast-growing gaps, i.e., n k + 1 n k 1 + ρ , ρ > 0 ( Hadamard   gap   condition ) , then the sequences sin ( 2 π n k x behave like independent random variables. For example, one has: 1 N k = 1 N sin ( 2 π n k x ) N ( 0 , 1 ) , where 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 ( n k ) k such that the sequence f ( n k ) k 0 for general periodic measurable functions f with vanishing integral exhibits random properties (see the historical overview in ). By introducing weights a k 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  and Weiss in ). Further, it can be shown that the process can be approximated by a standard Brownian motion:

(Philipp-Stout ). Assume the Hadamard gap condition. Assume further that A N : = 1 2 k = 1 N a k 2 and there exists δ > 0 such that lim N a N A N 1 δ = 0 . Then, without changing the distribution of the process: S ( t , x ) = k t a k cos ( 2 π n k x ) , t 0 , it can be redefined on a suitable probability space together with a Wiener process W ( t ) | t 0 such that: S ( t , x ) = W A t + O A t 1 2 ρ , a l m o s t   s u r e l y   f o r   s o m e ρ > 0 .

While the Hadamard growth condition (4) for CLT can be weakened for general sequences ( n k ) k (see ) for coefficients a k = 1 to the optimal growth condition n k + 1 n k 1 + c k k with c k , 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  and were recently sharpened by Berkes, Philipp and Tichy . The difficulties with the prime sequence p k k 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 p k + 1 p k are known. The best results for a lower bound that would be of interest for us do not hold for all k 0 , but only infinitely many. For the upper bound, it is proven by Goldston, Pintz and Yildirim () that lim inf k Δ k log p k = 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  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 ( n k ) , random models were also introduced (see, e.g., ). 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 ( p k x ) or cos ( p k x ) , the base of our series V α , β . Nevertheless, numerical computations strongly suggest that the central limit theorem holds; see Figure 9: we took 10 4 uniformly-distributed points x of the interval π , π and computed the sample average 1 N k = 1 N sin ( p k x ) for N = 78,498, that is, the number of primes 10 6 . 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 ( p k 3 2 x ) , see Figure 10.

4. Concluding Remarks

The properties of the series 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 α , β in this way and do not see a possibility to deduce it from any of them, besides from the trivial fact, that V α , β ( 0 ) is equal to the prime zeta function P ( α ) = p p α . Furthermore, recall that we only consider α R , usually α > 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 Cramér H. On the order of magnitude of the difference between consecutive prime numbers Acta Arith. 1936 2 23 46 10.4064/aa-2-1-23-46 Vartziotis D. Wipper J. The fractal nature of an approximate prime counting function Fractal Fract. 2017 1 10 10.3390/fractalfract1010010 Raseta M. On the strong approximations of partial sums of f(nkx) arXiv 2015 1509.08138 Philipp W. Stout W.F. Almost Sure Invariance Principles for Partial Sums of Weakly Dependent Random Variables American Mathematical Society Providence, RI, USA 1975 Hardy G.H. Weierstrass’s Non-Differentiable Function Trans. Am. Math. Soc. 1916 17 301 325 Jaffard S. Pointwise and directional regularity of nonharmonic Fourier series Appl. Comput. Harmon. Anal. 2010 22 251 266 10.1016/j.acha.2010.02.002 Hardy G.H. Littlewood J.E. Contributions to the arithmetic theory of series Proc. Lond. Math. Soc. 1913 2 411 478 10.1112/plms/s2-11.1.411 Gerver J. The differentiability of the Riemann function at certain rational multiples of π Am. J. Math. 1970 92 33 55 10.2307/2373496 Gerver J. More on the differentiability of the Riemann function Am. J. Math. 1970 93 33 41 10.2307/2373445 Chamizo F. Córdoba A. Differentiability and dimension of some fractal Fourier series Adv. Math. 1996 142 335 354 10.1006/aima.1998.1792 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) 10.1007/BF03014596 Fröberg C.E. On the prime zeta function BIT 1968 8 187 202 10.1007/BF01933420 Rosser J. Schoenfeld L. Approximate formulas for some functions of prime numbers Ill. J. Math. 1962 6 64 94 Cohen H. High Precision Computation of Hardy-Littlewood Constants 2000 Available online: https://www.math.u-bordeaux.fr/~hecohen/ (accessed on 1 January 2018) Cattani C. Fractal patterns in prime numbers distribution 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. Springer Berlin/Heidelberg, Germany 2010 Part II 164 176 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 10.1016/j.physa.2016.05.013 Kahane J.P. A century of interplay between taylor series, Fourier series and brownian motion Bull. Lond. Math. Soc. 1997 29 257 279 10.1112/S0024609396002913 Salem R. Zygmund A. On lacunary trigonometric series Proc. Natl. Acad. Sci. USA 1947 33 333 338 10.1073/pnas.33.11.333 16588761 Salem R. Zygmund A. On lacunary trigonometric series, II Proc. Natl. Acad. Sci. USA 1948 34 54 62 10.1073/pnas.34.2.54 16588787 Erdös P. Gál I. On the law of iterated logarithm I + II Indag. Math. (Proceedings) 1955 58 65 84 10.1016/S1385-7258(55)50010-2 Weiss M. The law of the iterated logarithm for lacunary trigonometric series Trans. Am. Math. Soc. 1959 91 444 469 Erdös P. On trigonometric sums with gaps Magyar Tud. Akad. Mat. Kutató Int. Közl 1962 7 37 42 Gaposhkin V.F. Lacunary series and independent functions Uspekhi Matematicheskikh Nauk 1966 21 3 82 10.1070/RM1966v021n06ABEH001196 Berkes I. Philipp W. Tichy R. Metric discrepancy results for sequences {nkx} and diophantine equations Diophantine Approximation: Festschrift for Wolfgang Schmidt Springer Vienna, Austria 2008 95 105 Goldston D.A. Pintz J. Yildirim C.Y. Primes in tuples I Ann. Math. 2009 170 819 862 10.4007/annals.2009.170.819 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 10.1007/BF00334035 Figures

Graph of (a): sin ( 2 π n k x ) with n k = 2 4 and n k = 2 5 ; and (b): sin ( 2 π n x ) with n = 5 and n = 6 (displayed for x [ 0 , 1 ) at 10 5 points).

Graph of (a): n = 1 10 sin ( 2 n π x ) ; (b): n = 1 1000 sin ( 2 π n x ) ; (c): p 1000 sin ( 2 π p x ) (displayed for x [ 0 , 1 ) at 10 4 points).

Graph of V 1.5 , 2 ( 10 5 , t ) at 5 × 10 4 discrete points in each direction (interpolated).

Graph of V 1.5 , 1.5 ( 10 5 , t ) at 5 × 10 4 discrete points in each direction (interpolated).

Graph of V 1.5 , 1 ( 10 5 , t ) at 5 × 10 4 discrete points in each direction (interpolated).

Graph of V 1 , 1 ( 10 5 , t ) at 5 × 10 4 discrete points.

The graph of the real part of 1 2 V 1 , 1 ( 10 6 , t ) (in black) and V 1 , 1 ( 10 6 , t / 3 ) + 1 3 (in green).

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

Normal distribution of 1 N k = 1 N sin ( p k x ) for x uniformly distributed in π , π .

Normal distribution of 1 N k = 1 N sin ( π p k 3 2 x ) for x uniformly distributed in π , π .