A Self-Referential Integral Equation for the Nontrivial Zeros of the Riemann Zeta Function
Michael Mark Anthonyorcid
Florida, USA.
DOI: 10.4236/apm.2026.167026   PDF    HTML   XML   5 Downloads   83 Views  

Abstract

Starting from Jensen’s integral representation of the Riemann zeta function, I derive, by a chain of convergent algebraic manipulations, a self-referential integral equation that a complex number F must satisfy in order for ρ= F( s )/ ( F( s )1 ) to be a nontrivial zero of the zeta-function, ζ. The resulting equation, is an explicit kernel and a biconditional characterization of the F-image of the nontrivial zeros. From calculations, I show that the Riemann Hypothesis is equivalent to the statement |F| = 1 for every solution. The equation simplifies dramatically to a Dirichlet-like series of lower incomplete gamma functions evaluated on the imaginary axis. A new representation of the Jensen’s integral for the Zeta function is show to be reducible to a series relating to Von Staudt-Clausen theorem. These new relations are shown to have a direct bearing on the relationship of the non-trivial zeros of the Zeta function and the primes. I discuss the analytical constraints and outline what is needed to close the gap to a proof of the Riemann Hypothesis.

Share and Cite:

Anthony, M.M. (2026) A Self-Referential Integral Equation for the Nontrivial Zeros of the Riemann Zeta Function. Advances in Pure Mathematics, 16, 452-495. doi: 10.4236/apm.2026.167026.

1. Introduction

The Riemann Hypothesis (RH) asserts that every nontrivial zero of the Riemann zeta function ζ( s ) lies on the critical line Re( s )= 1 2 . Despite a century of effort, no proof has been found. This paper contributes a new perspective: a rigorous reformulation of the zero condition as a self-referential integral equation in a single complex variable F , from which RH becomes a concrete analytical question about the modulus of solutions.

The key object is the Möbius-type map

F( s )= s s1 ,

which sends the critical line σ= 1 2 to the unit circle | F |=1 . Under this map, the functional equation of ζ, combined with Schwarz reflection, translates into a closure property on the F-image of the nontrivial zeros.

Starting from Jensen’s integral representation of ζ( s ) [1], I systematically derive an integral equation that every zero must satisfy, written purely in F. The derivation is reversible, making it a biconditional: F satisfies the integral equation if and only if ρ= F( s ) F( s )1 is a nontrivial zero of ζ( s ) .

2. Jensen’s Integral Representation

The starting point was my initial analysis of the Jenson formula in [1] (p 1036). In 2007, I was able to reduce the Jenson Zeta function to a Bernoulli sum. The conversion of the integral form to the Bernoulli form process is included in this paper in SECTION 12. However, due to the seemingly uncontrollable divergence of the sum, I was stalled by the Fubini-interchange theorem and (sum-integral interchangeability, a corollary of Tonelli’s Theorem 2.37) in [2]. The hypothesis for the interchange of the Integral and the Sum within the origin work was used to derive Bernoulli representation when the sum was converted to the ζ( 2n ) form and subsequently to the Bernoulli(2n) form. See SECTION 12. The original derived expression was verified by Professor Mark Prevost (LMPA Dunkirk, France). (See the Equation (24) → Equation (25) in this paper). In this paper, I finally overcame this obstacle by recasting the Bernoulli sum as a controlled split between primes and composites. This Bernoulli-form lead to the Prime counting function discussed in SECTION 12.

The Jenson identity is valid for all s with s1 :

ζ( s )= 2 s1 s ( 2 s 1 )( s1 ) + 2 s 2 s 1 0 sin( sarctanx ) ( 1+ x 2 ) s/2 e πx 1 dx . (1)

The integrand is a single closed-form function and there is no infinite series inside the integral. Convergence of the integral is straightforward: near x=0 , sin( sarctanx )~sx , and 1 ( e πx 1 ) ~ 1 ( πx ) , giving a bounded integrand; near x= , exponential decay dominates.

3. Derivation of the Self-Referential Zero Equation

By convergent manipulations only, I derive the self-referential integral equation for F( s ) at a nontrivial zero. Each step is reversible to (1) as demonstrated later.

3.1. Step 1: Substitute the Arctan Series

arctanx= π 2 k=0 ( 1 ) k ( 2k+1 ) x 2k+1 ,| x |>1 (2)

Substituting into (1) gives,

sarctanx= sπ 2 s k=0 ( 1 ) k ( 2k+1 ) x 2k+1 ,{ s,s1 } (3)

The integrand of (1) carries a term sin( sπ 2 Σ( x,s ) ) , where ( x,s ) = k=0 s ( 1 ) k ( 2k+1 ) x 2k+1 .

3.2. Step 2: Convert Sine to Exponentials

Using sin( z )= ( e iz e iz ) ( 2i ) :

sin( sarctanx )= 1 2i [ e isπ 2 e iΣ( x,s ) e isπ 2 e iΣ( x,s ) ] (4)

were, Σ( x,s )= k=0 s ( 1 ) k ( 2k+1 ) x 2k+1 .

3.3. Step 3: Absorb All Factors into the Exponents

Using ( 1+ x 2 ) s 2 =exp[ ( s 2 )ln( 1+ x 2 ) ] and   1 ( 2i ) =( 1 2 ) e iπ/2 , absorbs every polynomial/exponential prefactor into exponential form. After these rearrangements,

ζ( s )= 2 s1 s ( 2 s 1 )( s1 ) + 2 s1 2 s 1 0 e E ˜ + ( x,s ) + e E ˜ ( x,s ) e πx 1 dx ,{ s,s1 } (5)

with

E ˜ + ( x,s )= iπ( s1 ) 2 i( Σ( x,s ) ) s 2 ln( 1+ x 2 ),{ s,s1 } (6)

E ˜ ( x,s )= iπ( s+1 ) 2 +iΣ( x,s ) s 2 ln( 1+ x 2 ),{ s,s1 } (7)

The sign between the exponentials has been absorbed as e iπ =1 .

3.4. Step 4: Factor s from the Exponents

Extract s and rewrite each exponent as a constant (purely imaginary) piece plus s times a kernel function of x . After pulling the remaining constant phases iπ 2 into the s bracket as iπ ( 2s ) :

E ˜ ± ( x,s )=s Q ± ( x,s ),{ s,s1 } (8)

Q + ( x,s )= iπ 2s iπ 2 +i k=0 ( 1 ) k ( 2k+1 ) x 2k+1 + 1 2 ln( 1+ x 2 ),{ s,s1 } (9)

Q ( x,s )= iπ 2s + iπ 2 i k=0 ( 1 ) k ( 2k+1 ) x 2k+1 + 1 2 ln( 1+ x 2 ),{ s,s1 } (10)

3.5. Step 5: Collect the iπ 2s Structure

Note that iπ ( 2s ) iπ 2 = iπ( s1 ) ( 2s ) . Introduce F( s )= s ( s1 ) . Then 1 F( s ) = ( s1 ) s and iπ( s1 ) ( 2s ) = iπ ( 2F( s ) ) . Similarly for the opposite sign. Define:

P( x,F )= iπ 2F +i[ π 2 arctanx ]+ 1 2 ln( 1+ x 2 ), (11)

where (2) is used in the reverse direction to write the sum in closed form. Define the “formal conjugate”

P( x,F ) ¯ =+ iπ 2F i[ π 2 arctanx ]+ 1 2 ln( 1+ x 2 ), (12)

where all signs of imaginary terms are flipped, but the real ln term is unchanged. For real s , this equals the true complex conjugate; for complex s , it differs. Absolute integrability of both sides against 1 e πx 1 on (0, ∞) yields (13) as an identity of meromorphic functions on { s,s1 } .

Then Q ± take the compact form Q + ( x,F )=P( x,F ) and Q ( x,F )= P( x,F ) ¯ , and (5) becomes:

ζ( s )= 2 s1 2 s 1 { F( s )+ 0 e sP( x,F ) + e s P( x,F ) ¯ e πx 1 dx },{ s,s1 } (13)

3.6. Step 6: Impose ζ(ρ) = 0

At a nontrivial zero ρ , since the prefactor 2 ρ1 ( 2 ρ 1 ) is nonzero, the braces in (13) must vanish:

F( ρ )= 0 e ρP( x,ρ ) + e ρ P( x,ρ ) ¯ e πx 1 dx (14)

Note: Perron gives, at zeros, ρ

ψ 0 ( 1 ix )= 1 ix ρ e iπρ 2 x ρ ρ log( 2π ) 1 2 ln( 1+ x 2 ),

P( x,F )+ iπ 2F i[ π 2 arctanx ]= 1 2 ln( 1+ x 2 ),

P( x,F )= 1 ix ψ 0 ( 1 ix ) ρ e iπρ 2 x ρ ρ log( 2π ) iπ 2F( s ) +i[ π 2 arctanx ],

P( x,F )= 1 ix ψ 0 ( 1 ix ) ρ e iπρ 2 x ρ ρ log( 2π ) iπ 2ρ iarctanx,

P( 1 ix ,F )=x ψ 0 ( x ) ρ x ρ ρ log( 2π ) iπ 2ρ iarctanx,

ψ 0 ( 1ix )=( 1ix )i ρ ( 1ix ) ρ ρ log( 2π ) 1 2 ln( 1 ( 1ix ) 2 )

Using the identity ( 1ix ) ρ =i e ρP( x,F( ρ ) ) ρlog( 1ix )= π 2 +ρP( x,F( ρ ) )

π 2 =ρ( log( 1ix )P( x,F( ρ ) ) ) =ρ( log( 1ix )+ iπ 2F i[ π 2 +arctanx ] 1 2 ln( 1+ x 2 ) )

π 2 =ρ( log( 1ix )+ iπ 2F i[ π 2 +arctanx ] 1 2 ln( 1+ x 2 ) ),

one gets:

ψ 0 ( 1ix )=( 1ix )i ρ e ρP( x,ρ ) x ρ ρ log( 2π ) 1 2 ln( 1 ( 1ix ) 2 )

Riemanns original strategy from 1859 [3] was to use the explicit formula and analyze the sum over zeros. It’s been studied for 165 years. The connection through the P( x,ρ ) kernel is genuinely new notation, but the underlying object (sum of y ρ ρ over zeros) is classical. The fact that one can rewrite y ρ ρ as ( i ) e ρP( x,ρ ) is an algebraic reformulation. It might be useful for finding new angles, but it doesn’t change the analytical hardness.

3.7. Step 7: Close the Self-Reference via ρ= F( ρ ) F( ρ )1

From F( ρ )= ρ ( ρ1 ) , I solve: ρ= F( ρ ) F( ρ )1 . Substituting into (14), eliminate ρ and obtain an equation purely in F :

F( ρ )= 0 exp( F( ρ ) F( ρ )1 P( x,F( ρ ) ) )+exp( F( ρ ) F( ρ )1 P( x,F ) ¯ ) e πx 1 dx (15)

This is the main result.

Theorem 1 (Biconditional). For F with F1 , Equation (15) holds if and only if ρ:=F/ ( F1 ) is a nontrivial zero of ζ.

Proof sketch: Every step from (1) to (15) is reversible. Assuming (15), one runs the steps backward to recover ζ( ρ )[ 2 ρ1 2 ρ 1 ]=0 , and since the bracket is nonzero in the nontrivial range, ζ( ρ )=0 . The forward direction is the derivation itself.

4. Simplification via ( x±i ) ρ

Using arctan( 1 x )= π 2 arctanx (for x>0 ) and

log( x+i )=( 1 2 )ln( 1+ x 2 )+iarctan( 1 x ) (principal branch):

e ρP( x,F ) = e iπ ( 2( F1 ) ) ( x+i ) ρ = e iπ( ρ1 ) 2 ( x+i ) ρ , (16)

using F1= 1 ( ρ1 ) . Similarly,

e ρ P( x,F ) ¯ = e iπ( ρ1 ) 2 ( xi ) ρ ,ρ (17)

Substituting into (14) for the principal-value/regularized identity:

ρ ρ1 = e iπ( ρ1 ) 2 I + ( ρ ) e iπ( ρ1 ) 2 I ( ρ ),{ ρ,principal value/regularized } (18)

were

I ± ( ρ )= 0 ( x±i ) ρ e πx 1 dx ,ρ (19)

The two integrals satisfy

I ( ρ )= I + ( ρ ¯ ) ¯ ,ρ (20)

for all ρ , (the remaining integrand is O( 1 ) at the origin and decays exponentially at infinity).

Remark 1 (Individual divergence, combined convergence). The integrands ( x±i ) ρ ( e πx 1 ) are individually non-integrable near x=0 , where the integrand behaves as i ρ ( πx ) . However, the specific combination in (18) cancels at x=0 : writing ρ=σ+iT and expanding near x=0 , e ρP + e ρ P ¯ =i+i=0 , and higher-order terms give integrand =O( 1 ) . So (18) should be read as a principal-value/regularized identity, or equivalently via (14) or (15), where the numerator is manifestly bounded.

5. Modulus/Phase Decomposition

Writing ρ=σ+iT :

e ± iπ( ρ1 ) 2 = e πT 2 e ± iπ( σ1 ) 2 . (21)

Define the normalized integrals

J + ( σ,T ):= e πT 2 I + ( ρ ), J ( σ,T ):= e πT 2 I ( ρ ). (22)

These are bounded: the e πT/2 factors compensate the growth/decay of I ± . Equation (18) becomes:

σ+iT σ1+iT = e iπ( σ1 ) 2 J + ( σ,T ) e iπ( σ1 ) 2 J ( σ,T ). (23)

Observation 1. At fixed T=14.1347 (imaginary part of the first nontrivial zero), tabulation of | J + ( σ,T ) | for σ[ 0.3,0.7 ] gives essentially constant values (within 10−5): | J + |20.271 . The σ -dependence of J ± is concentrated almost entirely in their phases, not their moduli.

6. The Critical Line as Self-Conjugate Fixed Point

Proposition 1 (Algebraic properties of F ). For F( s )= s ( s1 ) , the following hold for all s where defined:

1) F( 1s )= 1 F( s ) .

2) F( s ¯ )= F( s ) ¯ .

3) Consequently F( 1 s ¯ )= F( s ) | F( s ) | 2 .

4) | F( s ) |=1Re( s )= 1 2 .

Proof: Direct computation. For (4): | F | 2 = | s | 2 | s1 | 2 . On σ=1/2 ,

| s | 2 = 1 4 + T 2 = | s1 | 2 , so | F |=1 . Off σ= 1 2 , strict inequality.

Corollary 3 (F-image of zeros). Let :={ F( ρ ):ρnontrivialzeroofζ } . Then is closed under:

  • F F ¯ (Schwarz reflection applied to zeros),

  • F 1 F (Functional-equation ρ1ρ ),

consequently, F 1 F ¯ .

Theorem 2 (RH as a modulus condition). The Riemann Hypothesis is equivalent to { | F |=1 } , equivalently to the statement that every solution of (15) lies on the unit circle.

7. Laplace-Transform/Incomplete-Gamma Form

This section develops an alternative representation of the zero equation that connects directly to classical special functions.

Expanding 1 ( e πx 1 ) = k=1 e πkx (valid for x>0 , uniformly convergent on [ ε, ) ):

I ± ( ρ )= k=1 0 ( x±i ) ρ e πkx dx . (24)

The inner Laplace integrals are known in closed form:

0 ( x+a ) ρ e sx dx = s ρ1 e sa Γ( 1ρ,sa ), (25)

where Γ( α,z )= z t α1 e t dt is the upper incomplete-gamma-function. With s=πk , a=±i , I get

0 ( x±i ) ρ e πkx dx = ( πk ) ρ1 e ±iπk Γ( 1ρ,±iπk ), (26)

and e ±iπk = ( 1 ) k .

Substituting into (18):

ρ ρ1 = k=1 ( 1 ) k ( πk ) ρ1 [ e iπ( ρ1 ) 2 Γ( 1ρ,iπk )+ e iπ( ρ1 ) 2 Γ( 1ρ,iπk ) ]. (27)

7.1. Reduction Using the Functional Equation of ζ

Split Γ( 1ρ,z )=Γ( 1ρ )γ( 1ρ,z ) , where γ( α,z )= 0 z t α1 e t dt is the lower incomplete gamma. Separating:

ρ ρ1 =Γ( 1ρ )2cos( π( ρ1 ) 2 ) π ρ1 η( 1ρ ) + k=1 ( 1 ) k ( πk ) ρ1 [ e iπ( ρ1 ) 2 γ( 1ρ,iπk )+ e iπ( ρ1 ) 2 γ( 1ρ,iπk ) ], (28)

where η( s )=( 1 2 1s )ζ( s ) is the Dirichlet eta function and I used k1 ( 1 ) k+1 k ρ1 =η( 1ρ ) . At a nontrivial zero ρ , the functional equation of ζ gives ζ( 1ρ )=0 , hence η( 1ρ )=0 . The first term in (28) vanishes, leaving:

ρ ρ1 = k=1 ( 1 ) k ( πk ) ρ1 [ e iπ( ρ1 ) 2 γ( 1ρ,iπk )+ e iπ( ρ1 ) 2 γ( 1ρ,iπk ) ]. (29)

Proposition 2 (Zero-detector). Equation (29) holds if and only if ζ( ρ )=0 (for ρ nontrivial).

7.2. Power-Series Form of the Bracket

Using the Taylor series γ( 1ρ,z )= n=0 ( 1 ) n z 1ρ+n [ n!( 1ρ+n ) ] and combining the two ±i branches, the ( i ) n + ( i ) n terms vanish for odd n and equal 2 ( 1 ) n 2 for even n=2m . After simplification:

ρ ρ1 =2 k=1 ( 1 ) k m=0 ( 1 ) m ( πk ) 2m ( 2m )!( 1ρ+2m ) . (30)

The outer sum converges conditionally (alternating in k ); the inner sum converges absolutely for each fixed k . The order of summation cannot be reversed naively.

8. Analytical Constraints

8.1. Convergence Constraints on the Laplace Transforms

Each individual Laplace integral (26) converges absolutely. The constraints are:

  • Re( πk )>0 : satisfied since k1 .

  • Branch consistency: ( x+i ) ρ on x[ 0, ) lies on the line Im( z )=1 in the right half-plane avoiding the principal-branch cut.

  • No singularity at x=0 : each individual Laplace integral is well-behaved because e πkx is bounded at 0.

8.2. Sum-Integral Interchange

When Inwrite I ± ( ρ )= k as in (24), this is conditionally convergent in total, not absolutely. The interchange is justified because:

  • Each term of the sum is a convergent integral.

  • On x[ ε, ) for any ε>0 , the series k e πkx is uniformly convergent and bounded by 1/ ( e πε 1 ) , allowing interchange on [ ε, ) .

  • The contribution from ( 0,ε ) to I ± is, individually, divergent; but in the combination e iπ( ρ1 ) 2 I + e iπ( ρ1 ) 2 I , the x0 singularities cancel.

In practice, (29) avoids these delicate issues entirely: each term is finite, and the cancellation structure has been absorbed into the reduction.

8.3. The Formal versus True Complex Conjugate

In the derivation, P( x,F ) ¯ denotes the formal conjugate (sign flip of imaginary kernel terms), not the true complex conjugate of P (which would involve F ¯ ). For real s , these coincide. For complex s , they differ, but the final forms (18), (29) are written in terms of the genuine complex integrals I ± ( ρ ) , which satisfy the true conjugacy (20).

8.4. Validity of the Arctan Expansion

The series (2) converges for | x |>1 ; for | x |<1 , the expansion is

arctanx= k0 ( 1 ) k x 2k+1 ( 2k+1 ) .

The derivation uses (2), but the closed-form sin( sarctanx ) ( 1+ x 2 ) s/2 is analytic across | x |=1 . The final forms (15) and (18) use closed-form expressions, so no domain-of-validity issue remains.

9. Numerical Verification

All computations were performed by Anthropic Claude in mpmath at 25 - 50 decimal precision.

Identity

Test point

Precision

(1)

s=2,3+i, ρ 1 , ρ 2 ,2,4

40 digits

(13)

first nontrivial zero

40 digits

(15)

F( ρ 1 )0.99750.0707i

28 digits

(18)

first nontrivial zero

28 digits

| F( ρ ) |=1σ=1/2

multiple T , multiple σ

exact

(29) zero detection

at zeros vs. non-zeros

confirmed

10. What Is Needed to Close the Gap to RH?

Theorem 1 gives a rigorous reformulation of the nontrivial-zero condition as a single integral equation in F . Combined with Proposition 2, the Riemann Hypothesis reduces to:

Every solution of (15) satisfies | F |=1 .

Equivalently, every ρ satisfying (29) has Re( ρ )= 1 2 .

The open question is therefore concrete and analytical, not combinatorial or number-theoretic. Possible approaches:

10.1. Phase Analysis Approach

Numerical evidence (Section 5) shows that | J + ( σ,T ) | is remarkably σ -insensitive at fixed T . If this can be proven analytically, (23) reduces to a phase-matching condition that may determine σ uniquely.

I now develop the first of the five approaches outlined in Section 10. The goal is to show that Equation (23) forces σ=1/2 by exploiting a near- σ -independence of | J + ( σ,T ) | at fixed T that emerges from the structure of the integrand.

10.1.1. The Divergence of I± and the Need for Regularization

As noted in the remark after Equation (18), the integrals I ± ( ρ ) defined in (19) are individually non-integrable at x=0 . The leading behavior near x=0 comes from the expansion

( x+i ) ρ = i ρ ( 1 ρx i +O( x 2 ) ), (31)

combined with the Laurent expansion

1 e πx 1 = 1 πx 1 2 + πx 12 π 3 x 3 720 +O( x 5 ). (32)

Multiplying:

( x+i ) ρ e πx 1 = i ρ πx +O( 1 )asx0. (33)

Similarly:

( xi ) ρ e πx 1 = ( i ) ρ πx +O( 1 )asx0. (34)

Using the principal-branch conventions i= e iπ/2 and i= e iπ/2 :

i ρ = e iπρ/2 , ( i ) ρ = e iπρ/2 . (35)

The divergent pieces have residues i ρ /π and ( i ) ρ /π , both with

| i ρ |= e πT/2 for ρ=σ+iT .

10.1.2. Cancellation in the Full Zero Equation

In the combination (18), the phase factors e ± iπ( ρ1 )/2 multiply I ± ( ρ ) . The residue of the divergent piece in the combined integrand is:

1 π [ e iπ( ρ1 )/2 i ρ + e iπ( ρ1 )/2 ( i ) ρ ]. (36)

Substituting (35):

e iπ( ρ1 )/2 e iπρ/2 = e iπ/2 =i, (37)

e iπ( ρ1 )/2 e iπρ/2 = e iπ/2 =+i. (38)

The combined residue is ( i+i )/π =0 : the 1/x singularities cancel exactly.

10.1.3. Regularized Individual Integrals

To isolate the finite information in J ± ( σ,T ) individually, I define:

J + reg ( σ,T ):= e πT/2 { 0 1 [ ( x+i ) ρ e πx 1 i ρ πx ]dx + 1 ( x+i ) ρ e πx 1 dx }, (39)

with analogous definition for J reg .

Proposition 1 (Regularized decomposition). Equation (23) may be rewritten in terms of the regularized quantities as

σ+iT σ1+iT = e iπ( σ1 )/2 J + reg ( σ,T ) e iπ( σ1 )/2 J reg ( σ,T ), (40)

because the regularization counter terms cancel by the same phase identities (37)-(38) that produce the cancellation in (18).

10.1.4. Numerical Observation on J + reg ( σ,T )

Computation of J + reg ( f( σ,T ),T ) over a grid of ( σ,T ) values yield the following (35-digit precision).

The values agree to 4 - 5 digits across σ[ 0.3,0.7 ] . The variation decreases as T increases: at T=50 , below 10−6; at T=100 , below 107.

σ

| J + reg ( σ,T ) | 2 at T=14.134725

0.3

1.12535027

0.4

1.12533109

0.5

1.12534725

0.6

1.12539876

0.7

1.12548559

Lemma 1: Conjecture 2 (Near-σ-independence of regularized modulus). (Near-σ-independence of regularized modulus). There exist functions g 0 : + + and g 1 :× + such that

| J + reg ( f( σ,T ),T ) | 2 = g 0 ( T )+ g 1 ( σ,T ), (41)

where g 1 ( 1 2 ,T )=0 and g 1 decays uniformly to zero as T .

10.1.5. Phase-Matching Strategy

Write ρ=σ+iT . Define the LHS and RHS of (40):

L( σ,T ):= σ+iT σ1+iT , (42)

R( σ,T ):= e iπ( σ1 )/2 J + reg ( f( σ,T ),T ) e iπ( σ1 )/2 J reg ( f( σ,T ),T ). (43)

Separating L into modulus and argument:

| L( σ,T ) | 2 = σ 2 + T 2 ( σ1 ) 2 + T 2 , (44)

argL( σ,T )=arctan( T σ )arctan( T σ1 ). (45)

Write J ± reg =| J ± reg | e i φ ± .

Lemma 2 (Modulus of R ).

| R( σ,T ) | 2 = | J + reg | 2 + | J reg | 2 +2| J + reg || J reg |cos( φ + φ +π( σ1 ) ). (46)

10.1.6. The Constraint from (40)

At a true zero,

L( σ,T )=R( σ,T ), (47)

In particular | L | 2 = | R | 2 :

σ 2 + T 2 ( σ1 ) 2 + T 2 = | J + reg | 2 + | J reg | 2 +2| J + reg || J reg |cos( φ + φ +π( σ1 ) ). (48)

Assuming Conjecture 1 in strong form, both | J + reg | and | J reg | depend only on T :

σ 2 + T 2 ( σ1 ) 2 + T 2 = G 0 ( T )+ G 1 ( T )cos( φ + ( σ,T ) φ ( σ,T )+π( σ1 ) ), (49)

were

G 0 ( T )= | J + reg ( T ) | 2 + | J reg ( T ) | 2 and G 1 ( T )=2| J + reg ( T ) || J reg ( T ) | (50)

10.1.7. The Symmetry σ 1σ

Under σ1σ :

| L( 1σ,T ) | 2 = ( 1σ ) 2 + T 2 σ 2 + T 2 = 1 | L( σ,T ) | 2 . (51)

Assuming compatible transformation of the regularized quantities:

| R( 1σ,T ) | 2 = 1 | R( σ,T ) | 2 . (52)

Proposition 4 (Fixed point of reflection). If both L( σ,T )=R( σ,T ) and L( 1σ,T )=R( 1σ,T ) hold, then | L( σ,T ) | 2 | L( 1σ,T ) | 2 =1 .

10.1.8. The Key Equation

Taking the difference of (48) at σ and 1σ :

σ 2 + T 2 ( σ1 ) 2 + T 2 ( σ1 ) 2 + T 2 σ 2 + T 2 = G 1 ( T )[ cos Φ σ cos Φ 1σ ], (53)

where Φ σ := φ + ( σ,T ) φ ( σ,T )+π( σ1 ) .

Let a= σ 2 + T 2 and b= ( σ1 ) 2 + T 2 . Then a b b a = a 2 b 2 ab , and a 2 b 2 =( ab )( a+b )=( 2σ1 )( 2 σ 2 2σ+1+2 T 2 ) :

a 2 b 2 ab = ( 2σ1 )( 2 σ 2 2σ+1+2 T 2 ) ( σ 2 + T 2 )( ( σ1 ) 2 + T 2 ) . (54)

This vanishes identically if and only if σ= 1 2 , (since 2 σ 2 2σ+1+2 T 2 >0 always).

Theorem 3 (Conditional reduction of RH to a phase identity). Assume Conjecture 1 in its strong form. Then the Riemann Hypothesis is equivalent to the statement: for every T>0 and every σ( 0,1 ) with σ 1 2 that arises as the real part of some non-trivial zero, the phase identity

G 1 ( T )[ cos Φ σ cos Φ 1σ ]= ( 2σ1 )( 2 σ 2 2σ+1+2 T 2 ) ( σ 2 + T 2 )( ( σ1 ) 2 + T 2 )

must hold. RH is the statement that no nontrivial zero has σ 1 2 , i.e., (54) has no solution with σ 1 2 among actual zeros.

10.1.9. Remaining Obstructions

Theorem 3 reduces RH to:

Problem A. Prove Conjecture 1 in strong form | J + reg ( σ,T ) | depends only on T, exactly.

Problem B (Phase rigidity). Show that (54) has no solution with σ 1 2 among actual zeros.

10.1.10. Numerical Status and Open Technical Issues

Lemma1: is a conjecture in its strong form, is that

| J + reg ( f( σ,T ),T ) | 2 := e πT/2 { 0 1 [ ( x+i ) ρ e πx 1 i ρ πx ]dx + 1 ( x+i ) ρ e πx 1 dx },

depends only on T for all σ . It does not. As shown by further analysis. So, I have removed it and made σ some function f( σ,T ) :

| J + reg ( f( σ,T ),T ) | 2 := e πT/2 { 0 1 [ ( x+i ) ρ e πx 1 i ρ πx ]dx + 1 ( x+i ) ρ e πx 1 dx },

The close values of | J + reg ( f( σ,T ),T ) | for f( σ,T ){ 0.3,0.4,0.5,0.6,0.7 } calculated for T=14.134725 a root, to the right give,

| J + reg ( f( σ,T ),T ) | 2 =1.12535027,1.12533109,1.12534725,1.12539876,1.12548559.1.1253503.

Analytical proof of strong f( σ,T ) -independence is required before Theorem 3 applies unconditionally. Equation (52) was asserted on the basis of “compatible transformation”. The regularization (39) subtracts a σ -dependent singularity; verifying that J + reg ( 1ρ ) relates to J + reg ( ρ ) via the expected reflection is a technical step not yet carried out. The residual σ -variation G 1 ( T ) in Conjecture 1 may contain precisely the information needed to force σ= 1 2 through a different mechanism than the strong-version argument. Understanding G 1 ( T ) is a meaningful subproblem. The phase functions φ ± ( σ,T ) are not known in closed form. Even granting Conjecture 1, Problem B is a nontrivial analytical task. The argument depends on regularization-scheme choices; showing the conclusion is scheme-independent requires verification.

10.1.11. Summary

The Phase Analysis approach reduces RH, under Conjecture 1, to the rigidity of a specific phase identity. The conjecture is supported numerically but not proven. Future work priorities:

  • Analytical proof of Conjecture 1 via change of variables x = tan(θ/2) or Mellin-Barnes contour analysis of I ± ( ρ ) .

  • Characterization of φ ± f( σ,T ),T via asymptotic expansion for large T , connecting to the Riemann-Siegel theta function.

  • Rigorous treatment of regularization covariance under ρ1ρ .

10.2. The Perron/Explicit-Formula Connection

The paper’s kernel P( x,F ) , defined in Equation (11), has a structural resemblance to Perron’s classical formula and to Riemann’s explicit formula for prime counting. In this section I make the connection explicit, deriving clean algebraic identities linking the exponential form e sP( x,F ) to Riemann’s explicit formula, and using these to obtain new representations of P in terms of the Chebyshev counting function ψ 0 and sums over the nontrivial zeros of ζ .

10.2.1. The Key Algebraic Identities

Exact identities that form the bridge between the kernel and Perron-type sums.

Proposition 1 (Exponential form of P). For x>0 real and s with F=s/ ( s1 ) ,

e sP( x,F ) =i ( 1ix ) s (55)

e +sP( x,F ) =i ( 1ix ) +s (56)

Proof. Using F= s ( s1 ) , one has s F =s1 . Substituting the paper’s definition

P( x,F )= iπ/ ( 2F ) +i[ π 2 arctanx ]+ 1 2 ln( 1+ x 2 ) :

sP( x,F )= isπ 2F is[ π 2 arctanx ] s 2 ln( 1+ x 2 ) = iπ( s1 ) 2 isπ 2 +isarctanx s 2 ln( 1+ x 2 ) = iπ 2 +isarctanx s 2 ln( 1+ x 2 ). (57)

Exponentiating:

e sP( x,F ) = e iπ/2 e isarctanx ( 1+ x 2 ) s 2 =i ( 1ix ) s ,

where in the last step 1 used 1ix= 1+ x 2 e iarctanx (principal branch, x>0 ). Identity (55) follows by the same computation with s replaced by 1s .

10.2.2. Riemann’s Explicit Formula

Riemann’s explicit formula for the Chebyshev counting function ψ( x )= p k x logp reads, in its analytic form ψ 0 ( x )=ψ( x ) 1 2 Λ( x ) 1 x :

ψ 0 ( x )=x ρ x ρ ρ log( 2π ) 1 2 log( 1 x 2 ), (58)

valid for real x>1 , where ρ runs over the nontrivial zeros of ζ . The sum is conditionally convergent when the zeros are taken in order of increasing | Imρ | , symmetrically (i.e., ρ paired with ρ ¯ ).

For complex y with | y |>1 and | argy |<π , the function ψ 0 ( y ) is defined by analytic continuation via the Perron-type contour integral

ψ 0 ( y )= 1 2πi ci c+i ( ζ ( s ) ζ( s ) ) y s s ds (59)

for c>1 , and Equation (57) extends by analytic continuation (though the sum over zeros requires careful interpretation for complex y ).

10.2.3. Riemann’s Formula at y=1ix

Substituting y=1ix (with x>0 real, so | y |= 1+ x 2 >1 and | argy |<π/2 ) into (57):

ψ 0 ( 1ix )=( 1ix ) ρ ( 1ix ) ρ ρ log( 2π )1/2 log( 1 ( 1ix ) 2 ). (60)

Using Proposition 1, identity (55) with s=ρ gives ( 1ix ) ρ =i e ρP( x,F( ρ ) ) , hence

ρ ( 1ix ) ρ ρ =i ρ e ρP( x,F( ρ ) ) ρ , (61)

where I stress that F( ρ )= ρ ( ρ1 ) varies with ρ . Substituting (60) into (59):

ψ 0 ( 1ix )=( 1ix )+i ρ e ρP( x,F( ρ ) ) ρ log( 2π )1/2 log( 1 ( 1ix ) 2 ). (62)

This is Riemann’s explicit formula expressed through the paper’s kernel. Then:

Theorem 4 (Perron-P identity). For x>0 real,

ψ 0 ( 1ix )=( 1ix )+i ρ e ρP( x,F( ρ ) ) ρ log( 2π )1/2 log( 1 ( 1ix ) 2 ), (63)

where the sum runs over the nontrivial zeros ρ of ζ , paired symmetrically with their conjugates, and F( ρ )= ρ ( ρ1 ) .

10.2.4. Inverting the Formula

Equation (62) can be rearranged to express the sum over zeros in terms of ψ 0 ( 1ix ) and elementary functions:

i ρ e ρP( x,F( ρ ) ) ρ = ψ 0 ( 1ix )( 1ix )+log( 2π )+ 1 2 log( 1 ( 1ix ) 2 ). (64)

10.2.5. A Representation of P(x, F) via Riemann’s Formula

Riemann’s explicit formula applied at a different argument, y= 1 ( ix ) = i x , to obtain a representation of P( x,F ) itself.

For x>1 , | 1 ( ix ) |= 1 x <1 , so (57) does not apply directly. However, the formal identity (extending ψ 0 by analytic continuation or taking the sum over zeros in a regularized sense) reads:

ψ 0 ( 1 ix )= 1 ix ρ ( 1/ ( ix ) ) ρ ρ log( 2π )1/2 log( 1+ x 2 ) (65)

where I used ( 1 ( ix ) ) 2 = ( ix ) 2 = x 2 , so 1 ( 1 ( ix ) ) 2 =1+ x 2 . The appearance of 1 2 log( 1+ x 2 ) is the critical observation: this is precisely the real part of P( x,F ) .

Solving (64) for 1 2 log( 1+ x 2 ) :

1 2 log( 1+ x 2 )= 1 ix ψ 0 ( 1 ix ) ρ ( 1 ( ix ) ) ρ ρ log( 2π ). (66)

Substituting (65) into the paper’s definition (11) of P( x,F ) :

Theorem 5 (Explicit-formula representation of P). For x>0 real,

P( x,F )= 1 ix ψ 0 ( 1 ix ) ρ e iπρ/2 x ρ ρ log( 2π ) iπ 2F +i[ π 2 arctanx ], (67)

where I used the identity ( 1 ( ix ) ) ρ = e iπρ/2 x ρ (principal branch).

Thus P( x,F ) itself admits a representation in which a sum over the nontrivial zeros of ζ appears explicitly.

10.2.6. The Involute Substitution x 1 ix

The transformation x 1 ( ix ) is involutive: applying it twice returns to x . Under this substitution, Theorem 4 yields an unexpectedly clean result.

Theorem 6 (Closed form under substitution). For x>1 real,

P( 1 ix ,F )= 1 2 log( 1 x 2 ) iπ 2F + iπ 2 arctanh( 1/x ). (68)

Proof. Apply (66) with x replaced by u= 1 ( ix )   i x . Then 1 ( iu ) =x , so the right-hand side contains x and ψ 0 ( x ) in place of 1 ( ix ) and ψ 0 ( 1 ( ix ) ) . Using Riemann’s explicit formula (57) at real x>1 :

x ψ 0 ( x ) ρ x ρ ρ log( 2π )=1/2 log( 1 x 2 ) (69)

Thus, the sum over zeros in (66) after substitution combines with x ψ 0 ( x )log( 2π ) to give 1/2 log( 1 x 2 ) . The remaining term i[ π/2 arctan( i/x ) ] is evaluated using arctan( i/x )=iarctanh( 1 x ) (valid for x>1 ), yielding iπ 2 arctanh( 1 x ) . Collecting terms gives the stated formula. Theorem 4 has a remarkable feature: the sum over zeros of ζ has been absorbed into the elementary function 1 2 log( 1 x 2 ) . The substitution x 1 ( ix ) effectively eliminates the appearance of ζ zero set in P .

10.2.7. Exponential Form under Substitution

Combining (67) with the identity e sP( x,F ) =i ( 1ix ) s (Proposition 1) applied at x = 1 ( ix ) = i x :

Corollary 5 (Exponential closed form). For x>1 real and s with

F= s ( s1 ) , e sP( 1/ ( ix ) ,F( s ) ) =i ( x x1 ) s (70)

This is consistent with e sP( x ,F ) =i ( 1i x ) s for x = i x , where 1i x =1 1 x = ( x1 ) x .

10.2.8. Sum over Zeros in Closed Form

Summing (69) over nontrivial zeros ρ (with F=F( ρ )= ρ ( ρ1 ) varying with ρ ):

ρ e ρP( 1/ ( ix ) ,F( ρ ) ) ρ =i ρ ( x ( x1 ) ) ρ ρ (71)

The sum on the right is a classical Perron-type sum at argument y= x ( x1 ) . For x>2 , I have 1<y<2 , and Riemann’s explicit formula (57) applies directly:

ρ ( x/ ( x1 ) ) ρ ρ = x x1 ψ 0 ( x x1 )log( 2π ) 1 2 log( 2x1 x 2 ). (72)

Here I used

1 ( x ( x1 ) ) 2 =1 ( ( x1 ) x ) 2 = ( 2x1 ) x 2 (73)

Theorem 8 (P-sum and prime counting). For x>2 real,

ρ e ρP( 1/ ( ix ) ,F( ρ ) ) ρ =i[ x x1 ψ 0 ( x x1 )log( 2π ) 1 2 log( 2x1 x 2 ) ] (74)

Equation (74) directly connects a sum over ζ ‘s nontrivial zeros, involving the paper’s kernel P evaluated at 1 ( ix ) with zero-dependent F , to the Chebyshev prime-counting function ψ 0 evaluated at x ( x1 ) , plus elementary corrections.

10.2.9. Structural Interpretation

Heorems 4, 5, 6, and 8 exhibit the paper’s self-referential kernel P as a natural object in the Riemann-Perron framework. Specifically:

1) The exponential e sP is algebraically equivalent to ( 1ix ) s (up to phase), which is the natural “shifted Perron” kernel on the line Re( z )=1 .

2) Riemann’s explicit formula applied at y=1ix yields a sum over zeros that is exactly i times the sum of e ρP( x,F( ρ ) ) /ρ .

3) The substitution x 1 ( ix ) trades the sum over zeros in P for an elementary logarithm, revealing an involute duality under which the zero-set contribution collapses.

4) The same substitution, applied to the exponential form, produces the Perron sum ρ ( x ( x1 ) ) ρ ρ , which is expressible in closed form through the prime-counting function ψ 0 .

10.2.10. Connection to the Riemann Hypothesis

Equation (74) makes explicit the link between the paper’s kernel and prime counting. The Riemann hypothesis is equivalent to the error bound

ψ 0 ( y )y=O( y log 2 y )asy. (75)

As x 2 + , the argument y= x ( x1 ) . By (74), the behavior of ψ 0 ( y )y is therefore controlled by the sum ρ e ρP( 1 ( ix ) ,F( ρ ) ) ρ as x approaches 2 from above.

More precisely, rearranging (72):

ψ 0 ( x x1 ) x x1 =log( 2π ) 1 2 log( 2x1 x 2 )+i ρ e ρP( 1 ( ix ) ,F( ρ ) ) ρ . (76)

RH is equivalent to a specific bound on the P -sum on the right of (76) as x 2 + .

10.2.11. Caveats and Open Questions

1) Theorem 4 and Equation (62) involve a sum over zeros with | 1ix | ρ terms. For ρ=σ+iT , one has | ( 1ix ) ρ |= ( 1+ x 2 ) σ/2 e Tarctanx , which grows exponentially in | T | . The series is conditionally convergent only via the Riemann-von Mangoldt pairing of zeros at finite truncation, with slow convergence at complex arguments.

2) Theorem 5 (Equation (66)) similarly uses a sum that, at y= 1 ( ix ) with | y |<1 , requires analytic continuation for direct numerical evaluation; the individual y ρ terms do not decay in | T | .

3) Theorem 6 (the closed form under substitution) is rigorous and numerically verified to machine precision at multiple test values of x and F.

4) Theorem 8 (Equation (74)) rests on Riemann’s explicit formula at y= x x1 for x>2 , where 1<y<2 and the classical formulation applies rigorously.

5) The question of whether the P -sum in (76) has a structural bound that forces (74) and thus RH remains open. The reformulation (74) is equivalent to the classical explicit formula, rewritten in the paper’s notation; it does not bypass the analytical hardness of estimating sums over zeros.

10.2.12. Summary

This section has established three rigorously proven structural identities:

1) The algebraic identity e sP( x,F ) =i ( 1ix ) s (Proposition 1) is exact, pointwise.

2) The closed-form identity P( 1 ( ix ) ,F )=1/2 log( 1 x 2 ) iπ ( 2F ) + iπ 2 arctanh( 1 x ) (Theorem 4) for x>1 is exact, and numerically.

3) The prime-counting representation (Theorem 5) relating ρ e ρP( 1 ( ix ) ,F( ρ ) ) ρ to ψ 0 ( x ( x1 ) ) for x>2 rigorous, follows from classical Riemann explicit formula.

These identities situate the paper’s self-referential kernel within the Perron/explicit-formula framework and provide new expressions for P in terms of ζ’s zeros and the prime-counting function. They establish that the reformulation developed in the paper is compatible with, and potentially complementary to, the classical approach initiated by Riemann. Whether these identities yield a novel attack on the Riemann hypothesis in particular, through the bound on the P -sum implicit in (75) is an open question for future work.

11. The Arctangent-Prime Connection

The paper’s kernel P( x,F ) contains an arctangent term: i[ π 2 arctanx ] . This section traces the remarkable fact that the same arctangent function encodes the zero-specific phases in Riemann’s explicit formula for prime counting. This connection rigorous, giving a clean identity (Theorem 9 below) linking argF( ρ ) to the phase arctan( 2γ ) that appears in von Mangoldt formula for ψ 0 ( x )x .

11.1. Arctangents in the Kernel

Recall the paper’s definition:

P( x,F )=1/2 ln( 1+ x 2 )+i[ π 2 arctanx ] iπ 2F . (77)

The imaginary part contains arctanx , an angular coordinate associated with the real variable x( 0, ) : geometrically, arctanx is the angle that the point ( 1,x ) subtends at the origin. A second arctangent appears implicitly at nontrivial zeros. For any ρ with Re( ρ )= 1 2 , I may write ρ=1/2 +iγ with γ , and

argρ=arctan( 2γ ), | ρ | 2 = 1 4 + γ 2 . (78)

11.2. Riemann’s Explicit Formula in Arctangent Form

Assuming RH (so all nontrivial zeros are of the form ρ= 1 2 +iγ ), Riemann’s explicit formula

ψ 0 ( x )=x ρ x ρ ρ log( 2π ) 1 2 log( 1 x 2 ) (79)

can be rewritten using (77). For a zero ρ= 1 2 +iγ with conjugate partner ρ ¯ = 1 2 iγ ,

x ρ ρ + x ρ ¯ ρ ¯ = 2 x cos( γlogxarctan( 2γ ) ) 1/4 + γ 2 , (80)

as follows by writing 1 ρ = e iarctan( 2γ ) 1 4 + γ 2 and x ρ = x e iγlogx . Summing over zeros γ>0 gives:

ψ 0 ( x )x=2 x γ>0 cos( γlogxarctan( 2γ ) ) 1/4 + γ 2 log( 2π )1/2 log( 1 x 2 ). (81)

The quantity arctan( 2γ ) is the zero-specific phase: it encodes the argument of ρ and determines, along with γlogx , the oscillatory contribution of each zero to the prime-counting function.

11.3. The Link to F(ρ)

The same phase arctan( 2γ ) appears in the paper’s kernel, through the map F( ρ )= ρ ( ρ1 ) .

Theorem 9 (Arctangent identity for argF ). Let ρ= 1 2 +iγ with γ (a point on the critical line). Then

| F( ρ ) |=1 and argF( ρ )=2arctan( 2γ )π. (82)

Proof. Write F( ρ )= ( 1 2 +iγ )/ ( 1 2 +iγ ) . Using arg( a+ib )=arctan( b/a ) for a>0 and the appropriate branch convention for a < 0 :

arg( 1 2 +iγ )=arctan( γ ( 1 2 ) )=arctan( 2γ ), (83)

arg( 1 2 +iγ )=πarctan( 2γ )( second quadrant forγ>0 ). (84)

Hence

argF( ρ )=arg( 1 2 +iγ )arg( 1 2 +iγ ) =arctan( 2γ )( πarctan( 2γ ) ) =2arctan( 2γ )π. (85)

For the modulus, | F( ρ ) | 2 = ( 1 4 + γ 2 ) ( 1 4 + γ 2 ) =1 , as established in Proposition 1 of the paper.

Corollary 10 (Kernel phase at critical-line zeros). For ρ= 1 2 +iγ ,

iπ 2F( ρ ) = iπ 2 e 2iarctan( 2γ ) = iπ 2 ( 1 2 iγ ) 2 1 4 + γ 2 . (86)

Proof. From Theorem 1, F( ρ )= e i( 2arctan( 2γ )π ) = e 2iarctan( 2γ ) . Hence 1 F( ρ ) = e 2iarctan( 2γ ) , and

iπ 2F( ρ ) = iπ 2 e 2iarctan( 2γ ) (87)

For the second form, use e iarctan( 2γ ) = ( 1 2 iγ )/ 1 4 + γ 2 , so

e 2iarctan( 2γ ) = ( 1 2 iγ ) 2 ( 1 4 + γ 2 ) (88)

11.4. The Structural Convergence

Theorem 9 and Equation (87) exhibit a shared arctangent structure. Specifically:

1) In the prime-counting explicit formula (77), the phase attached to each zero ρ= 1 2 +iγ is arctan( 2γ ) , coming from arg( 1 ρ )=argρ=arctan( 2γ ) .

2) In the paper’s kernel P( x,F( ρ ) ) at a critical-line zero, the “ F -phase” iπ ( 2F( ρ ) ) carries the factor e 2iarctan( 2γ ) , twice the arctangent of (1).

3) The factor of 2 comes from the Möbius identity F( ρ )= ρ ( ρ1 ) , which squares the phase relative to ρ itself.

Then,

Theorem 11 (Shared arctangent structure). For any critical-line zero ρ=1/2 +iγ :

arg( 1 ρ )=arctan( 2γ ) (explicit formula phase),

arg( 1 F( ρ ) )=2arctan( 2γ )+π (kernel phase, twice the zero angle).

In particular,

arg( 1 F( ρ ) )=2arg( 1 ρ )+π .

11.5. Interpretation: From Individual Zeros to Collective Phases

The explicit formula (80) may be rewritten, using (77), as

ψ 0 ( x )x=2 x γ>0 cos( γlogx+arg( 1/ρ ) ) | ρ | +O( 1 ), (89)

where ρ= 1 2 +iγ and the O( 1 ) captures the elementary corrections. Each zero contributes an oscillation cos( γlogx+arg( 1 ρ ) ) | ρ | —a “wave” with frequency γ and phase arg( 1 ρ ) .

Under ρF( ρ ) (the paper’s Möbius map), each zero’s phase arg( 1 ρ ) maps to arg( 1 F( ρ ) )=2arg( 1 ρ )+π : a doubling plus shift. This doubling is characteristic of Möbius transformations that fix the unit circle.

11.6. A Conjectural Reformulation of RH

The paper’s biconditional (Theorem 2) states: RH F( ){ | F |=1 } , where is the nontrivial-zero set. Combined with Theorem 9, RH is equivalent to:

argF( ρ )=2arctan( 2Im( ρ ) )πforeverynontrivialzeroρ.

Since (85) only requires | F( ρ ) |=1 (which implies Re( ρ )= 1 2 , so Im( ρ )=γ with ρ=1/2 +iγ ), the condition reads: the only constraint on each zero’s phase is the value of its imaginary part.

11.7. Connection to the Perron-P Identity

Combining Theorem 9 with the Perron-P identity of Section 10.2 (Equation (62)):

ψ 0 ( 1ix )=( 1ix )+i ρ e ρP( x,F( ρ ) ) ρ log( 2π ) 1 2 log( 1 ( 1ix ) 2 ), (90)

the arctangent content of each term can be made explicit. On the critical line,

P( x,F( ρ ) )= 1 2 ln( 1+ x 2 )+iarctan( 1 x )+ iπ 2 e 2iarctan( 2γ ) . (91)

The kernel thus carries two arctangent contributions: one from the spatial variable x and one from the zero-specific phase arctan( 2γ ) .

11.8. What This Does and Does Not Establish

Established rigorously.

1) Theorem 9: on the critical line, argF( ρ )=2arctan( 2γ )π .

2) Corollary 10: the kernel’s F -dependence is controlled by e 2iarctan( 2γ ) .

3) Theorem 11: the arctangent arctan(2γ) that appears in Riemann’s explicit formula (80) is the same arctangent (up to doubling and constant shift) that appears in the paper’s kernel P via F .

4. Equation (90): (Not established). The critical-line kernel carries explicit arctangent structure, with the zero-specific phase appearing as e 2iarctan( 2γ ) .

1) Whether the arctangent structure forces zeros to lie on the critical line. Off the critical line, a zero ρ=σ+iγ with σ 1 2 has arg( 1 ρ )=arctan( γ σ ) , still a well-defined arctangent. The arctangent structure does not inherently distinguish σ= 1 2 .

2) Whether the shared arctangent suggests a new bound on   ψ 0 ( x )x . The identities (80)-(81) are algebraic rewritings; they do not provide analytical leverage on the conditional sum over zeros.

11.9. Open Directions

The shared arctangent structure suggests the following avenues:

1) Fourier-theoretic interpretation. The explicit formula (88) views   ψ 0 ( x )x as a Fourier-like sum over frequencies γ with phases arg( 1 ρ ) . The paper’s kernel P( x,F ) can be viewed as a “Fourier integrand” with arctan( 1 x ) as the spatial phase and arctan( 2γ ) as the zero phase. A rigorous Fourier-theoretic identity between these two structures would yield new information.

2) Transformed-domain analysis. Equation (80) suggests changing variables from x to logx (so γlogx is the Fourier frequency) and from ρ to arg( 1 ρ )=arctan( 2γ ) (so the zero is parametrized by its phase). In this transformed domain, RH reads: for each Fourier frequency γ, the corresponding zero phase is exactly arctan( 2γ ) .

3) Critical-line characterization via phase. Theorem 9 characterizes the critical line as the locus where argF=2arctan( 2Im( s ) )π . Off the critical line, argF takes other values. Can the self-referential integral Equation (15) be shown to force argF( ρ )=2arctan( 2Imρ )π at zeros? This would be equivalent to RH.

11.10. Summary

This section has established a precise identity (Theorem 9) linking the phase of F( ρ ) at a critical-line zero to the arctangent phase arctan( 2γ ) that appears in Riemann’s explicit formula for prime counting. The identity makes rigorous a structural resemblance between the paper’s Möbius-based kernel and the classical zero-phase analysis of ψ 0 ( x )x . The resemblance is real and, has not been recorded in this form. Whether it opens a new attack on RH remains an open question; the identities in this section are algebraic consequences of the Möbius structure and of well-known phase relations, rather than independent analytical facts about the distribution of zeros.

12. The Prime-Bernoulli Bridge

This section establishes a structural bridge between two apparently unrelated objects: (a) the exact prime-counting formula arising from the Von Staudt-Clausen theorem and Wilson’s theorem, and (b) the divergent Bernoulli series obtained in the original derivation of the zero equation (Equation (112) of the paper, which I will call the wall). The arithmetic of Bernoulli numbers modulo 1 provides a canonical splitting of the wall into a convergent sum indexed by primes, and a divergent sum indexed by certain integers.

12.1. Prime Counting via Bernoulli Numbers

Theorem 12 (See [4] Exact prime counting via Bernoulli arctangent). For p an odd prime, let π( p1 ) denote the number of primes strictly less than p. Then

π( p1 )= 1 π n=1 ( p1 )/2 arctan( tan( π B 2n ( 2n )! ) )( 2n+1 ), (92)

where B 2n is the 2n-th Bernoulli number.

Proof. The map x arctan( tan( πx ) )/π is the fractional-part function centered in ( 1/2 ,1/2 ] . The Von Staudt-Clausen theorem asserts

B 2n + qprime ( q1 )|2n 1 q (93)

Multiplying by ( 2n )! :

B 2n ( 2n )!=( integer ) qprime ( q1 )|2n ( 2n )! q . (94)

For primes q2n satisfying ( q1 )|2n , the term ( 2n )!/q is an integer. The only possible non-integer contribution comes from q=2n+1 (which satisfies ( q1 )|2n trivially), and is included if 2n+1 is prime.

When 2n + 1 is prime, Wilson’s theorem gives ( 2n )!1( mod2n+1 ) , so ( 2n )!/ ( 2n+1 ) has fractional part 1/ ( 2n+1 ) ( mod1 ) , equivalently ( 2n )/ ( 2n+1 ) . Subtracting this from an integer, B 2n ( 2n )! has centered fractional part exactly 1/ ( 2n+1 ) . When 2n+1 is composite, B 2n ( 2n )! is an integer.

Therefore

1 π arctan( tan( π B 2n ( 2n )! ) )={ 1/ ( 2n+1 ) if2n+1prime, 0 otherwise. (95)

Multiplying by ( 2n+1 ) and summing from n=1 to ( p1 )/2 counts the primes in { 3,5,7,,p } , which is π( p1 ) .

Remark. Theorem 12 is a Bernoulli-based reformulation of classical prime-counting formulas via Wilson’s theorem. It is of expository rather than computational interest: the required precision in B 2n ( 2n )! is on the order of log 10 ( denominator ) , which grows like 2nlog( 2n ) ; this makes the formula impractical for large p .

12.2. The Wall Series Revisited

In the original derivation of the paper, I reduced the zeta relation (1) to the form of a Bernoulli-type series, by reducing the Jensen formula to a Bernoulli series:

Jensen formula. [1] (p 1036)

ζ( z )=( 2 z1 2 z 1 )( z z1 )+ 2 2 z 1 0 ( ( 1 4 + t 2 ) z 2 sin( z tan 1 2t ) e 2πt 1 )dt (96)

Substituting tan θ = 2t,

ζ( z )=( 2 z1 2 z 1 )( z z1 )+ 2 2 z 1 0 π 2 ( ( 1 2 seczθ )sin( zθ )( 1 2 ( secθ ) 2 ) e πtanθ 1 )dθ (97)

This reduces to

ζ( z )=( 2 z1 2 z 1 )( z z1 )+ 2 z 2 z 1 0 π 2 ( sin( zθ )( ( secθ ) 2z ) e πtanθ 1 )dθ (98)

Using the Chebyshev-T function,

cos( z cos 1 ( x ) )= 1 2 [ ( x+i 1 x 2 ) z + ( xi 1 x 2 ) z ] 2 (99)

Putting x=cosθ , and expand for sin( zθ ) ,

ζ( z )=( 2 z1 2 z 1 )[ ( z z1 )+2z 0 π 2 ( ( cosθ ) z2 ( 1 1 2 ( ( cosθ+i ( 1cos θ 2 ) 1 2 ) ) z + 1 2 ( ( cosθi ( 1cos θ 2 ) 1 2 ) ) z ) e π 1cos θ 2 cosθ 1 )dθ ] (100)

ζ( z )=( 2 z1 2 z 1 )[ ( z z1 )+2s 0 π 2 ( ( cosθ ) z2 ( 1 1 2 ( ( cosθ+i ( 1cos θ 2 ) 1 2 ) ) z + 1 2 ( ( cosθi ( 1cos θ 2 ) 1 2 ) ) z ) e π 1cos θ 2 cosθ 1 )dθ ] (101)

Which can be reduced again by the substitution, x=tanθ ,

ζ( z )=( 2 z1 2 z 1 )[ ( z z1 )+2z 0 ( ( 1 1+ x 2 ) 2z ( ( 1 1+ x 2 ) 2z ( 1 2 ( 1+ix ) z + 1 2 ( 1ix ) z ) 2 ) 1 2 e πx 1 )dx ] (102)

And finally, one arrives at the Abel Plena form [1]:

ζ( s )=( 2 s1 2 s 1 )[ ( s s1 )+2 0 ( i ( 1+ix ) s + ( 1+ix ) s e πx 1 )dx ] (103)

The function (103) can now be reduced to a series form as follows:

( 1+ix ) s = n=0 ( ix ) n Γ( 1s ) Γ( sn+1 )n! , ( 1ix ) s = n=0 ( ix ) n Γ( 1s ) Γ( sn+1 )n! (104)

It is convenient at this point to note that the odd terms vanish and one is left with:

ζ( s )=( 2 s1 2 s 1 )[ ( s s1 )+ 0 ( 2i n=1 ( ix ) 2n1 Γ( 1s ) 2Γ( s2n+2 )( 2n1 )! e πx 1 )dx ] (105)

ζ( s )=( 2 s1 2 s 1 )[ ( s s1 )+ n=1 ( 1 ) n1 Γ( 1s ) Γ( s2n+2 )( 2n1 )! 0 ( x 2n1 e πx 1 )dx ] (106)

Using the relation [1] (p 1038)

ζ( 2n )Γ( 2n )[ 1 π 2n ]= 0 x 2n1 e πx 1 dx (107)

ζ( s )=( 2 s1 2 s 1 )[ ( s s1 )+ n=1 ( 1 ) n1 2Γ( 1s )ζ( 2n )Γ( 2n ) π 2n Γ( sn+2 )( 2n1 )! dx ] (108)

ζ( s )=( 2 s1 2 s 1 )[ ( s s1 )+ n=1 ( 1 ) n1 2Γ( 1s )ζ( 2n ) π 2n Γ( sn+2 ) dx ] (109)

Using the relation,

ζ( 2n )= ( 1 ) n+1 2 2n1 π 2n B 2n ( 2n )! (110)

One arrives at the desired form:

ζ( z )=( 2 z1 2 z 1 )[ ( z z1 )+ n=1 2 2n B 2n Γ( 1z ) Γ( 2z2n )( 2n )! ] (111)

When the Zeta function vanishes at a root ρ ,

ρ ρ1 = n=1 2 2n B 2n Γ( 1ρ ) Γ( 2ρ2n )( 2n )! (112)

This series diverges factorially (the wall of Approach I first used [4] and [5]). as the Bernoulli numbers grow like B 2n ~ ( 1 ) n+1 2( 2n )!/ ( 2π ) 2n , making the terms grow unboundedly once ne π 2 23 . Using the Gamma-function identity Γ( 1ρ ) Γ( ρ2n+2 ) = ( ρ ) 2n1 (Pochhammer rising factorial), I rewrite (112) as:

ρ ρ1 = n=1 2 2n B 2n ( ρ ) 2n1 ( 2n )! , (113)

and equivalently (via B 2n = ( 1 ) n+1 2ζ( 2n )( 2n )!/ ( 2π ) 2n ):

ρ ρ1 = n=1 ( 1 ) n 2ζ( 2n ) ( ρ ) 2n1 π 2n , (114)

all interpreted as asymptotic equalities.

12.3. The Splitting via Von Staudt-Clausen

Now, applying the structural fact from the proof of Theorem 12: each B 2n ( 2n )! decomposes as an integer plus a prime-detector.

Lemma 13 (Bernoulli decomposition). For each n1 ,

B 2n ( 2n )!= N n + χ P ( 2n+1 ) 2n+1 , (115)

where N n and χ P ( m )=1 if m is prime and 0 otherwise.

Dividing (115) by ( 2n )!

B 2n = N n ( 2n )! + χ P ( 2n+1 ) ( 2n+1 )! . (116)

Substituting (116) into the wall series (112) gives a formal splitting into two subseries:

ρ ρ1 = T int ( ρ )+ T prime ( ρ ), (117)

were

T int ( ρ )= n=1 2 2n N n Γ( 1ρ ) Γ( ρ2n+2 ) ( ( 2n )! ) 2 , (118)

T prime ( ρ )= n=1 2 2n χ P ( 2n+1 )Γ( 1ρ ) Γ( ρ2n+2 )( 2n )!( 2n+1 )! (119)

These series converge absolutely for all sCZ0 and defines an entire function of s (the poles of Γ( 1s ) being cancelled by zeros of 1 Γ( 2s2n ) ).

12.4. Convergence of the Prime Part

Numerical values of T prime ( ρ ) .

At the first nontrivial zero ρ 1 = 1 2 +14.134725i :

T prime ( ρ 1 )3.6932670416+3.6491483483i,

computed to 20+ digits with absolute convergence verified through p=97 , where the term is of order 10⁻¹¹⁸. T prime ( ρ ) takes specific values at other notable points:

T prime ( 0 )=0 ,

T prime ( 1 )0.3687857 ,

T prime ( 2 )= { poddprime } 2 p1 p! 0.8127247 ,

T prime ( 1 2 )0.1776047 .

T prime ( 0 )=0 because ( 0 ) p2 vanishes for every p3 .

12.4.1. Interpretation

The splitting (118) expresses the wall series as a sum of two structurally distinct pieces:

1) T int ( ρ ) : carries the integer part N n of each B 2n ( 2n )! . This subseries is divergent, reflecting the wall’s original behavior.

2) T prime ( ρ ) : carries the prime-detector term, indexed only over n for which 2n+1 is prime. This subseries is absolutely convergent and defines an entire function of ρ .

The arithmetic of Bernoulli numbers modulo 1, specifically the Von Staudt-Clausen + Wilson mechanism that makes Theorem 12 work and acts as a natural regularization of the wall: it isolates a convergent prime-indexed subseries from a divergent integer-indexed one.

It is indeed important to note that the series:

ρ ρ1 = n=1 2 2n B 2n Γ( 1ρ ) Γ( ρ2n+2 )( 2n )! . (120)

Can also be represented by the change of variable:

ρ ρ1 = n=0 1 2 2n B 2n Γ( 1ρ ) Γ( ρ2n+2 )( 2n )! . (121)

This again becomes non-divergent, and has a simple form. Splitting (121) again:

ρ ρ1 = T int ( ρ )+ T prime ( ρ ), (122)

were

T int ( ρ )= n=0 1 2 2n N n Γ( 1ρ ) Γ( ρ2n+2 ) ( ( 2n )! ) 2 (123)

T prime ( ρ )= n=0 1 2 2n χ P ( 2n+1 )Γ( 1ρ ) Γ( ρ2n+2 )( 2n )!( 2n+1 )! (124)

Lemma 14 (Lemma 13 extended). (Integer-prime decomposition of the wall series.)

For each integer n0 , define the integer-part coefficient N n and the prime indicator χ P ( 2n+1 ){ 0,1 } by:

N 0 =1, χ P ( 1 )=0,

and for n1 , by the Von Staudt-Clausen decomposition

B 2n ( 2n )!= N n + χ P ( 2n+1 ) 2n+1 , (125)

were χ P ( m )=1 , if m is a prime, and 0 otherwise, and N n is the integer remainder.

Define the two component functions:

T int ( ρ )= ρ ρ1 + n=1 2 2n N n Γ( 1ρ ) Γ( ρ2n+2 ) ( ( 2n )! ) 2 (126)

T prime ( ρ )= n=1 2 2n χ P ( 2n+1 )Γ( 1ρ ) Γ( ρ2n+2 )( 2n )!( 2n+1 )! (127)

Then, the series defining T prime ( ρ ) converges absolutely for all ρC0 and defines an entire function of ρ . The series defining T int ( ρ ) diverges factorially as an asymptotic series, regularized through its identification with T prime ( ρ ) and at every nontrivial zero ρ of the Zeta function,  ζ( s ) , the homogeneous identity

T int ( ρ )+ T prime ( ρ )=0 (128)

Remark 1: (Arithmetic origin of N n =0 ). By the Von Staudt-Clausen theorem, the denominator of B 2 =1/6 is the product of primes q with ( q1 )|2 , namely q{ 2,3 } , giving denom B 2 =6 . Multiplication by 2!=2 cancels the factor of 2, leaving B 2 2!= 1 3 with denominator equal to the prime 2n+1=3 . Hence the entire content of the n=1 term in (127) is the prime-detector fraction χ P ( 3 )= 1 2 , and N 1 =0 . For n2 , ( 2n )! contains prime factors beyond the denominator of B 2n , and | N n | grows factorially.

Remark 2: The identity T int ( ρ )+ T prime ( ρ )=0 is rigorous as a statement about the convergent function T prime ( ρ ) : it assigns a canonical regularized value to the divergent series T int ( ρ ) at each nontrivial zero, namely T prime ( ρ ) . Independent characterization of T int ( ρ ) by some Borel summation, a contour integral representation, or an identification with a known special function is required to convert this identity from a definition into a constraint on the location of ρ .

The convergent prime-indexed series T prime ( ρ ) is at every nontrivial zero ρ exactly the negative of the regularized integer-indexed part T int ( ρ ) . The integer-indexed part is the asymptotic value of a series whose coefficients N n are integers arising from the integer parts of B 2n ( 2n )! . In other words, at every zero, the fractional/prime arithmetic of Bernoulli numbers exactly balances the integer arithmetic of Bernoulli numbers. Two different aspects of B 2n ( 2n )! (its fractional part (which encodes primes) and its integer part), are forced into a specific complex-valued balance at each zero. So, every nontrivial zero of ζ( s ) is a point where a specific convergent sum over primes equals a specific (regularized) sum over Bernoulli integer parts.

For each n1 , write the Von Staudt-Clausen decomposition of B 2n ( 2n )! as

B 2n = M n pprime ( p1 )|2n 1 p , M n , (129)

multiply by ( 2n )! ,

B 2n ( 2n )!= M n ( 2n )! pprime ( p1 )|2n ( 2n )! p + χ P ( 2n+1 ) ( 2n+1 ) (130)

The integer part is N n = M n ( 2n )! , the fractional part is the prime indicator at ( 2n+1 ) alone.

The zero condition links odd primes to zero locations (131)

Define the sum:

S( ρ )= poddprime 2 p1 ( ρ ) p2 ( p1 )!p! (132)

Every non-trivial zero, ρ , of the Zeta function satisfies:

T int ( ρ )=S ( ρ ) k = poddprime 2 p1 ( ρ ) p2 ( p1 )!p! (133)

where ( ρ ) k =ρ( ρ+1 )( ρ+k1 ) is the rising Pochhammer symbol. The right-hand side is an explicit convergent series indexed by the odd primes. The left-hand side is an object built from Bernoulli integer-arithmetic and the rational function ρ ρ1 . The equation forces the two to coincide at exactly the zero set of ζ .

The following TABLES show the relationship of the sum to the zeroes of the Zeta function.

For example, for ρ k = 1 2 +i γ k :

Table 1. Shows the values of ( ρ ) k at the first 10 nontrivial zeros of ζ .

k

γ k

ReS ( ρ ) k

ImS ( ρ ) k

| S ( ρ ) k |

1

14.1347251417

+3.6932670416

−3.6491483483

5.191965

2

21.0220396388

+37.9949667534

+17.6303128263

41.886101

3

25.0108575801

+98.8513216831

+81.7313765641

128.263797

4

30.4248761259

+304.3132227457

+345.1648080459

460.157888

5

32.9350615877

+486.7924752219

+606.0739744659

777.362577

6

37.5861781588

+1070.8293201482

+1553.9693074227

1887.192635

7

40.9187190121

+1758.2393694649

+2867.2657610298

3363.423646

8

43.3270732809

+2418.2384094418

+4342.6536965249

4970.565172

9

48.0051508812

+4005.9619165493

+9144.8376825371

9983.776195

10

49.7738324777

+4610.4350184051

+11869.1449163667

12733.134418

Each row gives the value of S at the k th nontrivial zero ρ k = 1 2 +i γ k .

Table 2. Shows the prime-by-prime contributions to S( ρ 1 ) .

p (odd-prime)

Re(term)

Im(term)

|term|

cumulative |error|

3

+0.16666667

+4.71157505

4.71452e+00

9.07406e+00

5

−4.98434470

−15.23726942

1.60318e+01

1.09418e+01

7

+8.38196187

+7.11319129

1.09934e+01

2.69514e−01

11

+0.11766387

−0.24791268

2.74418e−01

1.59713e−02

13

+0.01132720

+0.01127899

1.59850e−02

1.39819e−05

17

−0.00000803

−0.00001141

1.39500e−05

2.42983e−07

19

+0.00000018

−0.00000016

2.43010e−07

3.27315e−11

23

−3.27e−11

+3.27e−12

3.27315e−11

1.05160e−17

29

+1.05e−17

−1.05e−18

1.04856e−17

5.09905e−20

31

+5.10e−20

−5.10e−21

5.09905e−20

2.53577e−27

37

−2.54e−27

+2.54e−28

2.53577e−27

1.87423e−32

41

−1.87e−32

−1.87e−33

1.87086e−32

4.33038e−35

43

−4.33e−35

−4.33e−36

4.33037e−35

1.73301e−40

47

+1.73e−40

−1.73e−41

1.73301e−40

7.16413e−49

Term for prime p: 2 p1 ( ρ ) { p2 } [ ( p1 )!p! ] evaluated at ρ 1 = 1 2 +14.134725i ; Final value: S( ρ 1 )=3.69326704163.6491483483i (convergence to ~10⁻⁴⁹ by p ≤ 47).

Table 3. Shows the convergence of S( ρ 1 ) , partial sums by primes used.

Primes used

Re partial S

Im partial S

|error vs final|

{3}

+0.16666667

+4.71157505

9.07406e+00

{3, 5}

−4.81767804

−10.52569437

1.09418e+01

{3, 5, 7}

+3.56428383

−3.41250309

2.69514e-01

{3, 5, 7, 11}

+3.68194770

−3.66041577

1.59713e-02

{3, 5, 7, 11, 13}

+3.69327489

−3.64913678

1.39819e-05

{3, 5, 7, 11, 13, 17}

+3.69326686

−3.64914819

2.42983e-07

{3, 5, 7, 11, 13, 17, 19}

+3.69326704

−3.64914835

3.27315e-11

{3, 5, 7, 11, 13, 17, 19, 23}

+3.69326704

−3.64914835

1.05160e-17

{3, …, 29}

+3.69326704

−3.64914835

5.09905e-20

Truncated S( ρ 1 ) computed using only the indicated odd primes; The first five odd primes already give S( ρ 1 ) to 5 decimal places. Beyond p = 23, contributions are below 10−17.

Table 4. Shows S(ρ) at zeros vs. non-zero test points.

σ

T

Re S( ρ k )

Im S( ρ k )

| S( ρ 1 ) |

Note

0.500

14.134725

+3.6932670

−3.6491483

5.19196

ρ 1 (ZERO)

0.500

14.500000

+4.4022057

−3.6797767

5.73761

off zero

0.500

15.000000

+5.5016939

−3.6192994

6.58544

off zero

0.500

17.000000

+11.7150920

−1.7096260

11.83918

between zeros

0.500

21.022040

+37.9949668

+17.6303128

41.88610

ρ 2 (ZERO)

0.300

14.134725

+3.7147204

−3.2897098

4.96199

off line, T = γ1

0.700

14.134725

+3.6498747

−4.0246559

5.43318

off line, T = γ1

0.500

10.000000

−0.3355801

−1.1518469

1.19974

below first zero

0.500

20.000000

+28.9514692

+9.8383806

30.57746

between zeros

S(ρ) is well-defined at every point in . Only at zeros of ζ does it equal the regularized integer companion; Highlighted rows are nontrivial zeros of ζ. The value of S varies smoothly off the zero set.

Table 5. Shows the non-prime odd perturbations at ρ 1 .

Phantom p

Re (phantom term)

Im (phantom term)

|phantom term|

9

−2.56100881

+0.25284682

2.573460e+00

15

−0.00046688

+0.00033600

5.752117e−04

21

+2.6e−09

+1.8e−09

3.189771e−09

25

−2.7e−13

+2.7e−14

2.702669e−13

27

−1.8e−15

−1.8e−16

1.837719e−15

33

+2.1e−22

+2.1e−23

2.141018e−22

35

+7.8e−25

+7.8e−26

7.848521e−25

Phantom terms that would appear in S( ρ 1 ) if the indicated non-primes were treated as primes contributing; Their absence is the prime distribution indications; The phantom at p = 9 has magnitude 2.5, comparable in size to the contributions of actual small primes. Its absence (because 9 = 32 is not prime) is part of why S( ρ 1 ) takes the specific value it does.

Table 6. Asymptotic growth of | S( ρ k ) | with zero index k .

k

γ k

| S( ρ k ) |

| S | γ k

| S | γ k 2

1

14.134725142

5.191965

0.367320

0.025987

2

21.022039639

41.886101

1.992485

0.094781

3

25.010857580

128.263797

5.128325

0.205044

4

30.424876126

460.157888

15.124396

0.497106

5

32.935061588

777.362577

23.602888

0.716649

6

37.586178159

1887.192635

50.209751

1.335857

7

40.918719012

3363.423646

82.197677

2.008804

8

43.327073281

4970.565172

114.721923

2.647812

9

48.005150881

9983.776195

207.973020

4.332306

10

49.773832478

12733.134418

255.819851

5.139645

Growth ratios suggest | S( ρ k ) | scales faster than linearly in γ k but slower than γ k 2 ; All values computed at 60-digit internal precision (mpmath); displayed at 6 - 10 decimals.

The tables together demonstrate a concrete numerical relationship between the odd primes and the nontrivial zeros of the Riemann zeta function, through a single convergent series. They demonstrate a concrete numerical relationship between the odd primes and the nontrivial zeros of the Riemann zeta function, through a single convergent series (133).

SUMMARY OF THE TABLES:

Table 1 establishes the basic phenomenon: S( ρ ) is well-defined and computable at the first ten nontrivial zeros, and produces a distinct complex number at each. The values are not random; they grow systematically with γ k (the imaginary part of the zero), from | S( ρ 1 ) |5.19 up to | S( ρ 1 ) |12733 . Every nontrivial zero stamps the prime-indexed series with its own specific signature.

Table 2 opens the relation at the first zero. The contribution from each individual odd prime to | S( ρ 1 ) | is shown, and the structure is striking: the first three odd primes (3, 5, 7) carry essentially all the weight, contributing terms of magnitude 4.7, 16.0, and 11.0 respectively. By prime 17, contributions drop below 10−5; by prime 29, below 1017. The series converges with astonishing speed, not because the primes thin out, but because the factorial denominators ( p1 )!p! dominate the polynomial growth of the Pochhammer numerator. This is the analytical mechanism that makes the prime side of the prime-zero balance a rigorously convergent (entire) function.

Table 3 quantifies this convergence: using just { 3,5,7 } already locates S( ρ 1 ) to within 0.27 of its true value; adding { 11,13 } gets to 5 decimal places; the rest is increasingly fine corrections. The “prime distribution” is relevant to a given zero is overwhelmingly concentrated in a handful of small primes, and that handful expands slowly as one moves to higher zeros (where the convergence saddle point shifts upward).

Table 4 is the conceptual verification. S( ρ ) is computed at the first two nontrivial zeros (highlighted bold) and at seven non-zero test points that are points off the critical line, points between zeros, and points below the first zero. The sum is well-defined everywhere; the values vary smoothly. The “zero” rows are not arithmetically distinguished from their neighbors by anything visible in S( ρ ) alone. What is special about a zero is not the value of S( ρ ) but the equation

T int ( ρ )=S ( ρ ) k = poddprime 2 p1 ( ρ ) p2 ( p1 )!p! (134)

The zeros are precisely the points where the prime-indexed series equals the integer-indexed companion. The table makes this distinction concrete by showing S ( ρ ) k in both regimes.

Table 5 illustrates what “prime distribution shapes the zeros” actually means. The “phantom” terms (non-primes) show what would appear in S ( ρ ) 1 if non-primes (9, 15, 21, …) were prime. The phantom at p=9 has magnitude 2.57, comparable in size to the actual contributions from primes 5 and 7. Its absence from S ( ρ ) 1 is not a small effect; it is one of the structural reasons S ( ρ ) 1 takes the specific value 3.6933.649i rather than something else. If 9 were prime, the Equation (134) would have its solutions in different places: the locations of the zeros are sensitive to which integers are prime in this concrete, quantifiable way.

Table 6 documents the asymptotic behavior. The growth | S | γ k 2 rises from 0.37 (at ρ 1 ) to 256 (at ρ 1 ), while | S ( ρ ) k | γ k 2 rises from 0.026 to 5.14. So | S ( ρ ) k | grows faster than linearly in γ k but slower than quadratically characteristic γ k 2 near each zero by its contribution at primes p~ 2π| ρk | e . This is the rate at which the prime contributions of zeros grow with their height up the critical line.

With this analysis, it is time to get to the main theorem about the Riemann Hypothesis.

12.4.2. Main Theorem (A Möbius Parametrization of the Truncated Zero Condition)

The Parametrization Theorem

Theorem (Möbius parametrization of double-root configurations).

Let H 0 , H 1 with H 1 0 . The quadratic H 1 ρ 2 +( 1 H 1 )ρ+ H 0 =0 has a double root if and only if

( 1 H 1 ) 2 =4 H 0 H 1

When this holds, the double root is

ρ= 1 2 1 2 H 1 ,

with inverse parametrization

H 1 = 1 12ρ , H 0 = ρ 2 12ρ .

Overview

This is a step by step, the symbolic derivation that led from the Bernoulli form of the wall series to a clean Möbius parametrization of the truncated zero condition. Each step is recorded with its symbolic input and output, so that the chain of determinations is fully traceable. No numerical substitution is used in the derivation itself; the critical line emerges as a geometric consequence rather than as an a priori assumption.

STEP 1: The wall series and its low-order truncation.

Starting from Jensen’s identity, the nontrivial-zero condition (asymptotically) reads

ρ ρ1 = n=1 2 2n B 2n Γ( 1ρ ) Γ( ρ2n+2 )( 2n )! . (135)

Absorb the leading ρ ρ1 as a phantom n=0 entry and split each B 2n ( 2n )! via the Wilson-filtered Von Staudt-Clausen identity:

B { 2n } ( 2n )!= N n + χ P( 2n+1 ) 2n+1 , N n , χ P( m ) { 0,1 }. (136)

Writing X n := χ P( 2n+1 ) , and truncating at n=1 , the coefficients collapse beautifully because

2 2n ( ( 2n )! ) 2 =1and 2 2n 2n+1 [ ( 2n )!( 2n+1 )! ] =1forn=0,1. (137)

All prefactors equal 1 at low order, so the truncated zero condition takes the form of a quadratic in ρ :

( N 1 + X 1 ) ρ 2 +( 1 N 1 X 1 )ρ+( N 0 + X 0 )=0. (138)

STEP 2: The actual arithmetic values at n=0 and n=1 .

Computing the Wilson-filtered Von Staudt-Clausen values directly:

n

2n+1

prime?

B { 2n } ( 2n )!

N n

X n

Wilson mechanism

0

1

no

B₀·0! = 1

1

0

phantom: leading ρ/(ρ − 1)

1

3

yes

B₂·2! = 1/3

0

1

Wilson: 2! ≡ −1 (mod 3)

The complementarity is striking: at n=0 the entire contribution lies in the integer part; at n=1 it lies entirely in the prime part. The relevant sums for the quadratic are

N 0 + X 0 =1 and N 1 + X 1 =1 (Wilson-VSCvalues). (139)

STEP 3: The quadratic at Wilson values

Substituting the Wilson-VSC values into the truncated quadratic gives

1 ρ 2 +0ρ+1=0,i.e., ρ 2 +1=0, (140)

with roots ρ=±i . These are off the critical line by 1/2 in the real direction. So, the low-order arithmetic truncation does not reproduce nontrivial zeros of ζ. The information it carries is structural, not numerical.

STEP 4: Generalizing: replace ( N 0 + X 0 , N 1 + X 1 ) by free parameters ( H 0 , H 1 )

Abandon the Wilson constraint and treat H 0 = N 0 + X 0 , and H 1 = N 1 + X 1 as free complex parameters. The truncated quadratic becomes

H 1 ρ 2 +( 1 H 1 ρ )+ H 0 =0 (141)

Under this generalization the truncated equation is a two-parameter family. The question becomes: for which free parameters ( H 0 , H 1 ) does this quadratic have a root at a prescribed ρ , and in particular on the critical line?

STEP 5: Impose a double-root condition: discriminant = 0

Require the most degenerate configuration, a single double root. The discriminant of the truncated quadratic is

Δ= ( 1 H 1 ) 2 4 H 0 H 1 (142)

Setting Δ=0 yields the curve

( 1 H 1 ) 2 =4 H 0 H 1 (143)

in the complex parameter plane. This is the locus of all ( H 0 , H 1 ) producing a double root.

STEP 6: Solve for the double-root location ρ

When Δ=0 , the double root is at

ρ= B 2A = ( 1 H 1 ) 2 H 1 = H 1 1 2 H 1 (144)

This gives the forward Möbius map

ρ= 1 2 1 2 H 1 (145)

and its inverse:

H 1 = 1 12ρ , H 0 = ρ 2 12ρ (146)

STEP 7: The critical-line slice

From ρ= 1 2 1 2 H 1 , the real part of ρ is

Re( ρ )= 1 2 ( 1 2 )Re( 1 H 1 ). (147)

So Re( ρ )= 1 2 if and only if Re( 1 H 1 )=0 , equivalently:

Re( ρ )= 1 2 H 1 i( H 1 purelyimaginary ). (148)

The critical-line condition is therefore a one-real-codimension slice of the parameter space—purely geometric, not arithmetic. It emerges without being assumed.

Explicitly, writing H 1 = i 2T for real T0 :

ρ= 1 2 +iT, H 0 = 1 2 +i 14 T 2 8T . (149)

Both H 0 and H 1 are necessarily complex when T0 . The calibration to a nontrivial critical-line root cannot be achieved with real parameters.

STEP 8: Locate the Wilson point in this geometry

The actual Wilson-VSC values are ( H 0 , H 1 )=( 1,1 ) . Evaluating the curve equation

( 1 H 1 ) 2 4 H 0 H 1 at( 1,1 ):04=4. (150)

The Wilson point lies off the double-root curve C at discrepancy −4. This number is exactly the discriminant of the truncated quadratic at Wilson values, consistent with its actual roots being ±i (a conjugate pair with squared distance 1 from the critical line).

Equivalently, the Wilson point sits on the level set { Δ=4 } of the discriminant, not on the zero-discriminant curve C .

STEP 9: The X 0 N 1 relation under Wilson gauge

Fixing the Wilson gauge values N 0 =1 and X 1 =1 , and demanding the truncated quadratic have a double root on the critical line, gives the algebraic relation between X 0 and N 1 :

( N 1 +1 ) 2 ( 4 X 0 +6 )( N 1 +1 )+1=0 (151)

or equivalently:

X 0 = 3 2 + N 1 +1 4 + 1 4[ N 1 +1 ] (152)

Setting

u= N 1 +1= i 2T , (153)

this becomes

u 2 ( 4 X 0 +6 )u+1=0, (154)

whose two roots are reciprocals (by Vieta: u 1 u 2 =1 ). The two branches correspond to reciprocal T -values: if one root has T= T A , the other has T B = 1 4 T A .

The curve equation under Wilson gauge equals the discriminant of the truncated quadratic. They are the same polynomial. Hence, the X 0 N 1 relation is literally the statement Δ=0 .

Corollary 1: (Critical-line slice). The double root lies on Re( ρ )=1/2 if and only if H 1 i . Explicitly, H 1 = i 2T gives ρ=1/2 +iT .

Corollary 2 (Curve in ( H 0 , H 1 )-space). The double-root locus C={ ( H 0 , H 1 ): ( H 1 1 ) 2 =4 H 0 H 1 } is a smooth affine curve. The Möbius parametrization gives a bijection between C\{ ( 0,0 ) } and \{ 1/2 } via the map ( H 0 , H 1 )ρ .

Corollary 3 (Wilson position). The Wilson-VSC point ( H 0 , H 1 )=( 1,1 ) does not lie on C ; the discrepancy ( H 1 1 ) 2 4 H 0 H 1 at this point equals −4, exactly the discriminant of the truncated quadratic at Wilson values.

Facts about the Curve

The double-root locus of the truncated zero condition

C:={ ( H 0 , H 1 ) 2 : ( H 1 1 ) 2 =4 H 0 H 1 }

Equivalently, is the locus where the truncated quadratic H 1 ρ 2 +( 1 H 1 )ρ+ H 0 =0 has a double root.

Defining equation and degree

The curve is the affine complex algebraic curve defined by

F( H 0 , H 1 ):= ( H 1 1 ) 2 4 H 0 H 1 =0.

Expanded as

H 1 2 ( 4 H 0 +2 ) H 1 +1=0.

It is degree 2 in H 1 , degree 1 in H 0 —a conic.

Smoothness: Computing partial derivatives:

F H 0 =4 H 1 and F H 1 =2( H 1 1 )4 H 0

The system F= F/ H 0 = F/ H 1 =0 has no solutions in affine space. 𝓒 is smooth (non-singular) at every affine point. The curve is a smooth conic, 𝓒 has genus 0. It is birationally equivalent to 1 ( ) —the Riemann sphere. The map

φ:\{ 0 }C, H 1 ( ( H 1 1 ) 2 / ( 4 H 1 ) , H 1 )

is a bijection away from the puncture at H 1 =0 . Projection to ρ gives the Möbius transformation:

ρ= ( H 1 1 )/ ( 2 H 1 ) , H 1 =1/ ( 12ρ ) .

Points of 𝓒 with both H 0 , H 1 form the curve

H 0 = ( H 1 1 ) 2 / ( 4 H 1 ) ( H 1 \{ 0 } ).

This corresponds to ρ real. As H 1 ranges over ℝ\{0}, the resulting ρ covers \{ 1/2 } .

Points of 𝓒 with Re( ρ )= 1 2 are exactly { H 1 i\{ 0 } } .

On this slice:

Re( H 0 )=1/2 , Im( H 0 )= 14 T 2 8T , H 1 = i 2T .

The critical-line slice is a one-real-dimensional curve in 𝓒, parametrized by T\{ 0 } . The real-point locus and the critical-line locus meet only at H 1 =0 (point at infinity, ρ1/2 ). Wilson-arithmetic-real configurations and critical-line configurations live in geometrically disjoint slices.

The map ρ1ρ (pairing of zeros from ζs functional equation) corresponds on 𝓒 to:

H 1 H 1 , H 0 ( H 1 +1 ) 2 / ( 4 H 1 ) .

This is an involution of 𝓒 with fixed point only at H 1 =0 (corresponding to ρ=1/2 ). The map ρ ρ ¯ corresponds to standard complex conjugation in the parameter:

H 1 H 1 ¯ , H 0 H 0 ¯ .

The composite ρ1 ρ ¯ —the Schwarz reflection across the critical line corresponds to

H 1 H 1 .

The curve C={ ( H 1 1 ) 2 =4 H 0 H 1 } is a smooth complex conic (genus 0) parametrized by a single complex variable via a Möbius map to ρ. Its symmetry group includes the involutions induced by ζs functional equation and by complex conjugation.

The critical line Re( ρ )=1/2 is the fixed locus on of the Schwarz reflection ρ1 ρ ¯ , and equivalently, the imaginary-axis slice in the H 1 parameter. Real-arithmetic and critical-line configurations occupy disjoint slices, intersecting only at infinity ( ρ=1/2 ).

The geometric organization of 𝓒 and its symmetries, its slices, the position of Wilson—encodes structural information about how the Bernoulli arithmetic of the wall series interacts with the symmetries of the Riemann zeta function. Whether the higher-order truncations preserve this structure under appropriate re-summation, with H 1 becoming purely imaginary at exactly the nontrivial zeros, is the concrete open question this geometry frames.

The critical line is not assumed. It is the fixed set of the Schwarz reflection on 𝓒 geometrically determined by the symmetries of ζ.

The theorem establishes the following:

1) A clean symbolic parametrization of the n=0..1 truncated zero condition by a single complex parameter H 1 , via a Möbius map.

2) The critical-line condition Re( ρ )= 1 2 is identified with H 1 being purely imaginary, a geometric, not arithmetic, statement.

3) The Wilson-VSC arithmetic point (1, 1) sits at a specific computable discrepancy (−4) from the double-root curve C .

4) Under Wilson gauge, the X 0 N 1 relation is the discriminant of the truncated quadratic, exhibiting a reciprocal-root structure linking two T -branches.

12.4.3. Structural Remarks

1) The parametrization describes the structure of the truncated quadratic, not of the truncated wall series at the non-trivial roots. Wilson values produce truncated roots ±i , and actual zeros of ζ , Re( ρ )= 1 2 . The truncated quadratic captures the leading structural skeleton of the zero condition. Specifically:

The full zero condition

ρ ρ1 = n=1 2 2n B 2n Γ( 1ρ ) Γ( ρ2n+2 )( 2n )! (155)

is what defines nontrivial zeros.

2) The Möbius parametrization ρ= 1 2 1 2 H 1 , arises directly from this equation at low order.

3) The critical-line condition Re( ρ )= 1 2 H 1 i is a geometric consequence of the Möbius structure that appears already at low order.

4) Higher-order corrections from n2n modify the values of H 0 , H 1 but, if the structural Möbius form persists, would not break the critical-line geometry.

5) The critical-line condition σ= 1 2 in the truncated parametrization is not an arbitrary algebraic feature. It is:

  • The geometric image, under the Möbius map, of the imaginary-axis slice of the parameter H 1 .

  • A structural property of the zero condition itself, visible already at n=0..1 . linked to the wall series organization through the Wilson-VSC split.

6) The convergent prime sum T prime ( s ) is an entire function of s . In contrast, the zero equation ρ/ ( ρ1 ) has a pole at ρ=1 .

7) T prime is defined for all complex s , not just at zeros. It does not, by itself, detect zeros: T prime ( σ+iT ) varies smoothly with σ,T , and off the critical line gives values comparable to those at zeros of similar T .

8) The equation T int ( ρ )=ρ/ ( ρ1 ) T prime ( ρ ) at zero asserts that the divergent integer-indexed part of the wall is constrained to equal an explicit convergent quantity. A rigorous interpretation of this asymptotic equality for instance via Borel summation or optimal truncation could extract information about the distribution of zeros.

9) The prime sum T prime ( s ) is reminiscent of the Kanemitsu-Kuzumaki exponential series and of related “Lucas-type” series over primes. The precise relation to Dirichlet series (via Mellin transform of T prime ) may merit further investigation.

12.4.4. Caveats and Open Questions

1) The splitting (118) is a formal decomposition of an asymptotic series. The original wall series (113) does not converge, and while the prime-part T prime is convergent, the integer-part T int remains divergent. The identity ρ ρ1 = T int ( ρ )+ T prime ( ρ ) inherits this asymptotic character.

2) T prime ( s ) is a rigorously defined entire function; its evaluation at a zero yields a specific complex value. But this value alone does not constrain ρ to the critical line: T prime is an entire function, taking all complex values on preimages of any given value.

3) The result requires the truncated series. It establishes a structural connection linking two Bernoulli-based objects, the prime-counting formula (95) and the wall series (113), that had been studied separately in this work. Whether this connection yields analytical leverage is an open question.

4) One potential direction for further analysis: combine T prime ( s ) with the Perron-P identity of Section 10.2 (Equation (62) and Equation (73)).

ψ 0 ( 1ix )=( 1ix )+i ρ e ρP( x,F( ρ ) ) ρ log( 2π )1/2 log( 1 ( 1ix ) 2 ). (156)

ρ e ρP( 1/ ( ix ) ,F( ρ ) ) ρ =i[ x x1 ψ 0 ( x x1 )log( 2π ) 1 2 log( 2x1 x 2 ) ] (157)

Both involve sums over primes or zeros; a joint analysis might produce new constraints.

12.5. Summary

This section has established:

1) Theorem 12: an exact Bernoulli-based prime counting formula, obtained via Von Staudt-Clausen + Wilson.

2) Lemma 13: the decomposition of B 2n ( 2n )! into an integer part and a prime-detector term.

3) The convergent prime sum T prime ( s ) as an entire function of s.

4) The splitting of the wall series into a divergent integer-part and a convergent prime-part.

The bridge is genuine and new in this form. Its immediate use is structural: it demonstrates that the prime content of the Bernoulli numbers (Wilson/Von Staudt-Clausen) and their analytic content (ζ-values, wall series) are both present in the paper’s derivation and can be disentangled by the splitting (117). Further work can be read through the provided references for study [6]-[11].

13. Conclusions

Fubini’s theorem states that for any integral of any integrable function in r+s the interchange of order of integration is allowed. Tonelli’s Theorem gives a converse for nonnegative functions. By strictly convergent manipulations from Jensen’s integral representation, a self-referential integral Equation (15) that biconditionally characterizes the F-images of the nontrivial zeros of ζ. The equation admits a simplified form (18) in terms of classical ( x±i ) ρ integrals, and reduces further via Laplace transforms to a Dirichlet-like series of incomplete gamma functions (27), (29). The Riemann Hypothesis is equivalent to the statement that every solution of these equations lies on the unit circle | F |=1 .

The primary contribution of this work is the derivation chain itself, which avoids the divergent-series and Fubini-interchange difficulties that have obstructed earlier algebraic attempts on RH via similar starting points.

Acknowledgements

I would like to thank the many great mathematicians from who I learnt a everything I know about the Zeta function and the Riemann Hypothesis. This paper is a small contribution to their Mastery of the subject. The more I learn the more I realize how little I know compared to them. I would also like to thank Anthropic AI, Claude for working through the algebraic steps, its willingness to verify identities numerically, and its diagnosis of convergence issues that arose at several points in the derivation and the Tables in this paper. Any remaining errors are my own.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Gradshteyn, I.S. and Ryzhik, I.M. (2007) Table of Integrals, Series, and Products. Seventh Edition, Academic Press.
[2] Folland, G.B. (1999) Real Analysis: Modern Techniques and Their Applications. 2nd Edition, Wiley.
[3] Riemann, B. (1859) Über die Anzahl der Primzahlen unter einer gegebenen Größe. Monatsberichte der Königlichen Preußischen Akademie der Wissenschaften zu Berlin, 671-680.
[4] Anthony, M.M. (2024) Algebraic Methods for the Riemann Hypothesis. Advances in Pure Mathematics, 14, Article ID: 5302511.
[5] Anthony, M.M. (2025) Consequences of Invariant Functions for the Riemann Hypothesis. Advances in Pure Mathematics, 15, 36-72.[CrossRef]
[6] Tonelli, L. (1909) Sur l’integrazione per parti [On Integration by Parts]. Atti de la Accademia Nazionale de Lincei, 18, 246-253. (In Italian)
[7] Samko, S.G., Kilbas, A.A. and Marichev, O.I. (1993) Fractional Integrals and Derivatives. Gordon and Breach, 9 p.
[8] Borwein, J., Bailey, D. and Girgensohn, R. (2004) Experimentation in Mathematics: Computational Paths to Discovery. A K Peters, 18 p.[CrossRef]
[9] Thomas Jr., G.B. and Finney, R.L. (1996) Calculus and Analytic Geometry. 8th Edition, Addison-Wesley, 919 p.
[10] Fubine, G. (1958) Sugli Integrali Multipli. Opere Scelte, 2, 243-249.
[11] Friedman, H. (1980) A Consistent Fubini-Tonelli Theorem for Nonmeasurable Functions. Illinois Journal of Mathematics, 24, 390-395.[CrossRef]

Copyright © 2026 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.