The age-old Riemann Hypothesis is a conjecture that the Riemann Zeta function has its zeros only at the negative even integers and complex numbers with real part and is firstly stated in the essay by Bernard Riemann .
The fact that the (warped) Riemann Zeta function under certain conditions will arise as the solution of a second order ordinary differential equation on the time domain while the Riemann Zeta function may arise as the same solution but on the warped frequency domain (domain of the Dirichlet transform (93) image functions defined originally on a time domain) gave rise to possible presumptions that the two mentioned functions may be a Dirichlet Transform of each other. Since the extended form of the Riemann Zeta function is found to be an eigenfunction of the Fourier Transform operator, a composite Mellin Transform is found that projects into preserving scale properties and is de facto a line-invariant transform between and ,since there exists a well-known functional relationship between and ,the so-called functional equation of the extended Riemann Zeta function
A fact that confirms the age-old conjecture of Bernard Riemann.
The -form of the eigenfunctions of the Fourier transform and the Dilation operator are thoroughly studied in . The discrete and collapsed forms of these hyperbolic (eigen)functions are used as analytical input signals in      for an analytical signal analyzing affine model.
In this manuscript, the discrete super-position of the mentioned (warped) eigenfunctions are used to find a new eigenfunction of the Fourier Transform and to build solutions of the known and studied second order Ordinary Differential Equations (ODEs) on both time and warped frequency domain.
The solutions of the momentum equation of the Navier Stokes for both 2d and 3d versions are related to Dirichlet series (95), on the warped frequency domain. An appropriate inverse Dirichlet transform may be used to retrieve the solution on the time domain.
A general solution of the Navier Stokes 3d equation for the viscid incompressible flows is constructed from the begin conditions in the warped frequency domain, using geometric properties of the 2d version. This solution is called a minimal solution to the 3d Navier Stokes equation for the viscid incompressible flows. Those types of functions may contain domain singularities, which will cause the blow-up time when the numerical pivots strike a singularity. Please do see Appendix A for a study of the Navier Stokes d.e.
2. The Euler Differential Equation and a Think Experiment
Consider the solutions for the Euler differential equation
that may be deducted from the second order polynomial characteristic equation
[With substitution of to eliminate the non-constant terms does transpose to and followed by the substitution of delivers the found characteristic polynomial Equation (3).]
The solutions of the Euler differential equation  is given by
and D the discriminant equivalent with
If we think the variable t as restricted to the set of natural numbers ,and and as constants and equal to one and take improper super-position of all identified solutions, then the Riemann Zeta Function will arise as the family bundle of solutions of (2) as follows:
with the zeros of the characteristic polynomial (3) defined as for
and as when
In our think experiment we have got both the Riemann Zeta function and the warped Riemann Zeta function and the derivative of the Riemann Zeta function as solutions to the Equation (2).
The Riemann Zeta function arises also as a special case of the super-position of the general solutions (see TableB1 for examples in Appendix B) of other second order ordinary differential equations like the confluent hyper geometric differential equation called hyper geometric function of the first and second kinds, generalized Laguerre d.e., the Bessel d.e. etc. and the incompressible Navier-stokes equation
,but in the warped frequency domain. We would have expected then that the Riemann Zeta
warped frequency domain (as an image of the Dirichlet Transform operator). This is very fundamental in the sense that with proper assumptions and a proper time-frequency transformation these types of functions can co-exist both on the time and warped frequency domain as we are going to see next. Therefore it was very important for us to analyze the behavior of such functions on the warped frequency domain. Appendix B consists of a study of the relationship between discrete solutions of selected 2nd order ODEs and the Riemann Zeta Function. Eigenfunctions of the Fourier Transform and the Dilation operators are thoroughly discussed in the thesis of Garas in  and the denomination “warped frequency domain” that is used in this manuscript is adapted from .
The Mellin transform operator   does consist of the warping operator , p. 122 and further, followed by a Fourier transform. Its image will therefore reside on the warped frequency domain, when the initial function does reside on the time domain.
3. The Extended Riemann Zeta Function as an Eigenfunction of the Fourier Transform
Consider the monomial represented by and its defined unitary Fourier transform (with ordinary frequency) . If we take the improper sum of both and its Fourier Transform, a procedure also known as the Poisson Summation formula, then will the extended form of the Riemann Zeta function in the complex plane emerge as follows,
using the functional equation for the extended Riemann Zeta function (1). In other words the extended Riemann Zeta function is invariant under the Fourier Transform.
This is a very important find because we can now derive a general formula for the Mellin Transform with fractional derivatives of the Fourier Transform and use it to calculate the Mellin Transform of the Riemann Zeta function.
 if all frequency components of are positive (or if is an analytical signal), can be seen as a Fourier Transform with transform variable and with
and using the formula for the Fourier Transform of the μ-th derivative of a Fourier integrable function f,
with i the imaginary unit and extended from the set of natural numbers to the set of complex numbers and with
If we take the Riemann Zeta function as the analytical signal
will above formula yield for
with using the fact that the fractional derivative of
for and . The identified scaling and dilative effect (or Doppler effect) makes this transform ideal for analytical signal reconstruction and signal analysis purposes applied in radar technology.
4. The Proof of the Age-Old Riemann Hypothesis
For normalization purposes and to be coherent with the used Fourier transform we will use instead of .
An affine or a line-preserving map between and may be formed by the composite Mellin transform operator ,with
and the warping operator on variable n and the Fourier transform. ( ,similarly with the discrete Mellin Transform used in ,entails a manipulation with the summation index variable to calculate the result). Using the found composite Mellin transform on will yield
Q.e.d. The “quod erat demonstrandum” abbreviation refers to the confirmation of the age-old Riemann Hypothesis, knowing that there exists a functional equation between and .
The popularized expression for the Riemann Zeta function in the complex plane on the so-called critical line, ,is also a solution of higher order
ODEs when the pair of complex conjugate solutions of the corresponding characteristic polynomials are repeating as we have seen in our thought experiment. Its warped version (thus warping of a warped function with ) can be represented as an even characteristic polynomial, since the Riemann Zeta function is an eigenfunction of the Fourier transform. Hence
and with ,with defined as in the On-Line Encyclopedia of Integer Sequences as OAEIS A053117 ,that is comparable with a characteristic polynomial of a higher order ODE for some . Here are and respectively the Chebyshev polynomials of resp. the first and second kinds . Since the index v has to be taken even in the cumulative polynomial gives the opportunity to form repeating pairs of conjugate complex roots as the resulting characteristic polynomial may then be decomposed in product of lower order parabolic expressions which is a conditio sine qua non for the warped Riemann Zeta function to exist as a solution of the discussed ODE.
The expression is equivalent with the Fourier series of the zero-th Bernoulli polynomials,
The last definition can be found in the public accessible on-line Wikipedia.
If we do interchange the variable x and n with each other, a logic way to calculate the Dirichlet Transform of an Ordinary Generating Function resulting in a Dirichlet series, yields:
what equivalent is with for and is called the Hurwitz formula for the generalized extended zeta function and with the well-known functional equation for the Riemann Zeta function will emerge.
The extended generalized zeta function is invariant under the Dirichlet Transform (and so under the Mellin Transform), when transforming only the domain of its arguments will undergo a variable swap as we have seen. The Hurwitz Zeta function is therefore an eigenfunction under both the Dirichlet and the Mellin transforms.
 p. 177.
5. Conclusions and Future Work
Using the fact that the extended Riemann Zeta function is an eigenfunction under the Fourier Transform, it is proved that there exists a line-invariant composite Mellin transform operator that projects into namely
what confirms the validation of the age-old Riemann Hypothesis, knowing that there exists a well-known functional equation between and .
De facto they are conjugate pairs of each other using repetitively a fractional derivative Fourier Transform operator or the so-called Mellin Transform operator (12). The Mellin transform will expand the transform variable for from the set of natural numbers to the set of real numbers and will expand the transform variable from the set of real numbers to the set of (negative) complex numbers .
Our future work will entail the study of the full spectrum of the eigenspaces belonging to the (fractional derivative) Fourier Transform operator and the behavior of the so-called conjugate pairs on their respectively domains with regard to known solution procedures for partial differential equations.
Minimal solution for the 3d viscid incompressible Navier Stokes equation: an application to physics
The solution of the Laplace equation in cylindrical coordinates is of paramount importance in many applications and this is not an exception for finding a minimal solution of the 3d Navier Stokes equations for the viscid incompressible media. The viscid flow arises when the viscous forces are prominent present in the equations. The solution is called minimal because a 2d geometric solution is extended to the 3d case, where domain singularities might be inherited.
If we extend the solutions to the complex space and take the super-positioned result of all positive discretized solutions will yield a polynomial or hyperbolic progression, which is not surprising as we are going to see.
Hence all second order ordinary differential equations that we have studied in this manuscript do have a polynomial progression as a solution as we are going the see in the next paragraphs.
The solution of the Mellin Transformed Laplace equation in cylindrical coordinates on the radius variable
If we discretize ,then take the superposition of all discrete solutions, is still a general solution of above equation, namely
and using instead of s and taking and (to simplify the expression), then will the solution become
since is invariant under the Fourier Transform or
if we would have taken summation over the variable x.
If we take now the Dirichlet Transform of the solution (33) for ,will this result in with if we would have used a general term . Please note that for ,yields
with the Heaviside step function .
Above, the so-called functional equation of the Riemann Zeta function (1) is used, and
The Mellin Transform of (32), for and ,may be represented as
the series expansion of the begin condition f around the value 0. Equation (38) is de facto a solution on the warped frequency domain to the eigenvalue problem Equation (30) or the Helmholtz equation that typically arises in physical problems involving partial differential equations in both space and time (e.g. the viscid Navier Stokes incompressible d.e.).
When the begin conditions differ per independent solution legs, the same method may be applied but now instead of using ,its real part (with regard to the term with ) must be used and respectively its imaginary part (with regard to the term with ) must be used to calculate (38).
The time domain equivalent of (38) may be calculated by the inverse Dirichlet transform and is defined by
The Dirichlet Transform of a function f is defined as
what is equivalent with
The above solution (38) can also be used for the reconstruction of the general solution (for the velocity profile, and thus also for the pressure warped frequency distribution) of the incompressible 2d-Navier Stokes equation in cylinder coordinates, taking
since the Laplace equation in cylindrical coordinates is congruent with the 2d momentum equation for an incompressible flow. The solution for the (warped) velocity distribution function is therefore:
With the begin function f equivalent with 1 does Equation (46) yield the Dirichlet series taking :
Figure A1. warped frequency domain.
Figure A2. time domain.
The 3d Navier Stokes system of equations for the incompressible viscid flow  is reducible to an extent that its general and exact solutions are determined by the (general and exact) solutions of the 2d Laplace equation in cylindrical coordinates (as above) as follows. We will call this solution a minimal solution.
Consider the NVS 3d systems of equations for the incompressible viscid flow:
Equation (48) per directional component in cylindrical coordinates:
To solve the above 3d NVS system of equation we will have to reduce the r.h.s. of each equation to the sum of the pressure term and the external force and equating the Laplace expression in cylindrical coordinates to zero (solve up the Laplacian to get the velocity distribution), and then expand this using time convolution to with the appropriate Gaussian diffusion Green’s function to solve the enclosed diffusion equation
and then use it to calculate the pressure distribution as follows
To simplify this we will expand the super-position of discrete solutions Equation (32) over the whole set of integers ,then
The last function is rescaled along both axes compared to its time-domain version, and is called a “compressed” eigenfunction of the Mellin Transform operator, the corresponding transform pair on the time domain is equivalent with
In two dimensions it will be congruent with
The interested reader may try to derived it, the warped frequency 3d pressure distribution in Cartesian coordinates would be then, taking
out of Equation (54) and disregarding the terms for and the Heaviside time restriction in the different directions,
based on an i, j, k indexed classical Hamilton’ian unit vectors with magnitude where and
, , , , , and where Erf(x) the Error function is, defined as
On the edge of the domain the following is true for the normal vectors
The total energy of each wave system is equal on both the time and frequency domain  . In symbols,
with above the expression for the Wigner-Ville Distribution function ( a class of dual energy distributions), and below, the integrand [is equivalent to the Ambiguity Function] as defined in   
is the analytical signal on the frequency domain and its Fourier Transform inverse on the time domain and the total (kinetic) energy of the system and r is a warp parameter that may be set to yield a factor
in the integrand and what may be calculated with the Residue theorem. Expression (65) will yield an unbounded wave energy-level of for every component with .
The Wigner-Ville time-frequency distribution function is a unitary affine function that suffices the following equation
The total wave energy does inherit the known singularities within the Dirichlet series from the components of the distribution function of what
Figure A3. Contour of 3d warped pressure distribution.
makes it also unbounded.
The last maximum value of the total energy is also supported by the fact that the Wigner-Ville Distribution function definition also does have a factor of the hyperbolic cosecant with half argument in its integrand (defined on ,with poles at what may lead to variants of the improper Riemann contour integral 
Below schematically, the application of the presented formulas for solving up the Laplace equation in cylindrical coordinates in an analytic signal analyzing model used by . The constructed transformation is scale invariance:
An analytic signal analyzing model
Out of the general solution is it clear that the distribution functions for and on the warped frequency domain do inherit the known singularities of the Dirichlet series, naturally depending on the begin-conditions which may cause shifts in singularities at ,removable singularities at
with p a prime number as a consequence of the Euler denomination of the Dirichlet series ( for the Riemann Zeta function) in the form of the product formula
only for is totally multiplicative . Solving up the zeros of the denominator in the last expression
does result in the mentioned complex removable singularities. So may be said that the distribution functions and are defined on for each component.
The begin-conditions may also dictate the periodicity of the general solution. There exists no general solution that is independent of the begin-solutions, since the solution is constructed out of the begin-conditions. The begin-conditions may cause singularities or may annihilate them, depending on the initial problem and its basic solution. The basic solution of the Laplace equation in cylindrical coordinates, which is the solution core for both the velocity and pressure frequency distributions (in the incompressible case) for the Navier Stokes equations in 2d and to certain extent also for the 3d version, does contain and consequently inherits the singularities observed in the cited improper Riemann contour integral in particular or from the Dirichlet series in general. The input signal
of the used model may manifest or cause singularities on both time and warped frequency domains with both the summation parameter indexed by f (hyperbolas) and by (polynomials). On the time domain the signal may experience singularities due to the vast powers that are used that may be unpractical for most commercial available micro-processors to compute [for example when f is chosen as the series of orthogonal polynomials like the Hermite polynomials or due to the known singularities of the Zeta function. On the warped frequency domain the singularities are due to the properties of the formed Dirichlet series output signals for both cases.
Interesting to note is that the so-called blow-up time within the numeric schemes may manifest whenever the numerical pivoting method strikes a known singularity of the general solution.
When the analytical signal condition is dropped or when negative discrete solutions of the Laplace equation (read the momentum equation) are allowed, the motion solution of the 3d Navier Stokes equations for the viscid incompressible flow will become a (multi-dimensional) Gaussian function with no singularities on neither the time domain nor the warped frequency domain.
The Riemann Zeta Function represented as a Dirichlet Transform of discrete super-position of second order ODEs solutions.
The confluent hyper geometric functions are solutions to the Weber differential equations, obtained after separation of variables when decomposing the Laplacian in cylindrical coordinates;
The substitution converts the weber equation to
with solution ,where is the Hermite polynomial.
In general the solutions of the second order ordinary confluent hypergeometric d.e.
are of the form
and are called the confluent hypergeometric function of the first and second kind respectively. The first kind version is also denoted as or with
and the Pochhammer symbol defined as
An integral representation for is
The confluent hypergeometric function of the second kind, can be transposed to the improper contour integral
by choosing ,and . Namely
Taking an infinitely summation at both sides of the latter expression and substituting for n n*n and for s will emerge as super-position of the solutions, namely
If we take instead of both sides as argument then the following expression will arise
is a solution to the Weber d.e. with p an integer and equal to null. is the other linear independent solution of Weber d.e.
With substitution converts the Weber d.e. to with solution .
the Hermite polynomial and the modified Hermite polynomial.
The modified Hermite polynomial is defined as
If we denote by and compare it to the found expression of i.e.
and substituting s by -s yields,
If we interchange the parameters z and n with each other for ,and superposition all the elements with each other will the expression of emerge as above. The generalized Laguerre polynomials can be represented by
and are related to the Riemann Zeta function as described above. A known transform operator  that can transpose polynomial expressions into Dirichlet series is the Dirichlet transform
that transforms the ordinary generating functions
(OGF) to (95)
This transform is discussed in the chapter of combinatorics of  p. 177 and deals primarily with size conditions of of the permutations with auxiliary parameter variable k.
The Dirichlet series with
where it may be used in conjunction with the Laplace-Borel Transform first, as intermediary aid, to first convert the Exponential Generating Functions (EGF)
into the OGF form what may be the case for the confluent hyper-geometric functions of the first kind. The Laplace-Borel Transform is defined as
In  Vojta mentioned about the cannibalistic feature of the Mellin Transform with regard to the Laplace Transform and “many other functions” what may validate the two step approach to calculate the Dirichlet Transform result for the generalized Laguerre polynomials.
TableB1 shows the study of several second order ODEs with regard to their solutions on the time domain, and the improper summation of their Dirichlet transforms on the warped frequency domain. The transform of the solution double argument functions entails a variable swap between the domain variables and .
is one of the independent solutions of the confluent hypergeometric differential equation. The row elements in the TableB2 are enclosed in the integer sequence OIES A113025  and the summation of the diagonal elements are the Dirichlet coefficients of .
Table B1. Second order ODEs (Ordinary Differential Equations) and their solutions on time and warped frequency domains.
Table B2. Progression of the modified Bessel function of the first kind .
 Ovarlez, J.P. (1994) La Transformation de Mellin Transform et L’analyse des signaux Large Bande. Coll. Temps-Frequence, Ondelettes et Multiresolution: Theorie, Modeles et Applications, Vol. 9, 1-8.
 Wikipedia Contributors (2017) Chebyshev Polynomials.