Riesz potentials (also called Riesz fractional derivatives) and their Hilbert transforms are computed for the Korteweg-de Vries soliton. They are expressed in terms of the full-range Hurwitz Zeta functions and . It is proved that these Riesz potentials and their Hilbert transforms are linearly independent solutions of a Sturm-Liouville problem. Various new properties are established for this family of functions. The fact that the Wronskian of the system is positive leads to a new inequality for the Hurwitz Zeta functions.

1. Introduction

In recent years the theory of fractional derivatives and integrals called Fractional Calculus has been steadily gaining importance for applications. Ordinary and partial differential equations of fractional order have been widely used for modeling various processes in physics, chemistry, and engineering (see, e.g., [13] and the references therein). Recent theoretical developments shed new light on the interpretation and properties of fractional derivatives. Having written the latter in the form of Stieltjes integrals, Podlubny [4] found new physical and geometric interpretation of these structures relating them to inhomogeneity of time. Extension of the classical maximum principle to the case of a time-fractional diffusion equation appeared in the recent work of Luchko [5]. In the present paper we are concerned with Riesz fractional derivatives (also called Riesz potentials; see [6, page 88], and [7, page 117]) that are defined as fractional powers of the Laplacian with . They are well known for their role in investigating the solvability of nonlinear partial differential equations, and the Korteweg-de Vries equation (KdV henceforth) in particular (see, e.g., [711] and the references therein). In the current work, Riesz potentials of KdV solitons are computed and their relation to ordinary differential equations is established.

We continue the study of Riesz fractional derivatives of solutions to Korteweg-de-Vries-type equations started in [12]. After appropriate rescaling, KdV can be written in the form It is well known that the fundamental solution of the Cauchy problem for the linearized KdV is expressed in terms of the Airy function of the first kind and its Hilbert transform (conjugate) in terms of the Scorer function . The papers [1214] were devoted to the study of fractional properties of the Airy functions and their conjugates and to the establishing of related properties for KdV-type equations.

Although there exists extensive literature on solitons, as far as we know, a study of their fractional properties is still missing. A preliminary investigation of Riesz potentials for a KdV soliton was carried out in [15]. In this paper the emphasis was put on the issue of whether solitons inherit fractional properties of fundamental solutions. Riesz potentials of a soliton, , where , , and , and their Hilbert transforms, , were obtained in terms of the Hurwitz Zeta function of a complex argument, with and . It was proved in [15] that the zero mean properties hold for both and with . This confirmed the predictions based on the properties of fundamental solutions of the linearized Cauchy problems.

The goal of the current paper is to go further and to study Riesz potentials of solitons as solutions of differential equations. We intend to show that these functions and their Hilbert transforms form linearly independent systems of solutions for a second-order ordinary differential equation in a self-adjoint form. This fact may be helpful in understanding the issue of using these structures as intrinsic mode functions in signal processing (see [16, 17] and the references therein), that is, in using Riesz potentials for expansions. In this context it is interesting to point out that the graphs of the functions reveal a striking similarity to those of the Airy wavelets generated by the function (see [18, page 34], and Figure 3 below).

For the analysis to follow; we employ the full-range Hurwitz Zeta functions: and (symmetric and antisymmetric combinations of and ), recently introduced in [19] for . We prove that the functions , are solutions of the Sturm-Liouville problem for . Here , , and is a real function. The essential point consists in proving that the Wronskian of and is positive for all and . It allows one to prove that and are positive and to estimate the number of zeros of and on any bounded interval.

The fact that this Wronskian is positive also leads to a new inequality for the Hurwitz Zeta functions where , , and the bar over the letter denotes complex conjugation. As far as we know, there are no results in the literature on the arguments of the Hurwitz Zeta functions. However, (1.4) provides some information on this issue. Indeed, setting and in (1.4) allows one to deduce the relation for , , and .

The paper is organized as follows. In Section 2, we provide the necessary information on the special functions involved. Section 3 is devoted to the study of Riesz potentials for KdV solitons and their Hilbert transforms. In Section 3.1, the main properties of these functions are summarized. Sturm-Liouville problem (1.2) is derived in Section 3.2. Section 3.3 deals with the properties of the Wronskian . Zeros of the functions and are studied in Section 3.4. In Section 4, the inequality (1.4) is discussed.

2. Preliminaries

Introduce the Fourier transform of the function by and the inverse Fourier transform by

For real and , define Riesz potentials of a function by the formula (see [7, page 117]) provided that the integral in the right-hand side exists. Define derivatives of with respect to by provided that these integrals exist.

Introduce the Hilbert transform of the function by (see [20, page 120]) where and denotes the Cauchy principal value of the integral. According to our choice of the Fourier transform , one can see that on , where is the identity operator. Also, and , where the operator is defined by (2.3).

Next, introduce the Trigamma function by (see [21, page 260, 6.4.1]) Notice that (see [21], page 260, 6.4.7) Also, the following asymptotic expansion holds:

The Hurwitz (generalized) Zeta function is defined by (see [22, page 88]) This implies that This function has the integral representation where is the Gamma function. In two particular cases, we have (see [22]) that where is the Riemann Zeta function.

The singularity of as is given by the relation The asymptotic expansion of for large is (see [23]) where are the Bernoulli numbers.

Introduce the full-range forms of Hurwitz Zeta functions (see [19]) Representations (2.16) and (2.17) imply that Hence, follow the symmetric and antisymmetric properties of the functions, It follows from (2.16) that is a periodic function of , with the unit period. It is even with respect to and . The function is odd about and . These functions satisfy the functional differential equations

Consider now that and denote by the complex conjugate of . Then . It implies that, for , In a similar way,

Denote by the Wronskian of the functions and , that is, . For reader's convenience, we present here [24, Theorem ].

Theorem 2.1. Let , be real valued and continuous for . Let and be real valued solutions of the equation satisfying the condition Let be the number of zeros of on . Then

In conclusion of this section, we would like to quote an interesting result concerning integrals over the real axis (see [25]).

Theorem 2.2. For any integrable function and with , Moreover, the above formula holds true if where is any sequence of positive constants, are any real constants, and the series is convergent.

3. Fractional Derivatives of A KdV Soliton and Their Conjugates

In this section, we consider Riesz fractional derivatives of a KdV soliton and their Hilbert transforms and establish their properties. We notice that all the graphs were obtained with the Mathematica 6 software.

We take the soliton solution of (1.1) in the form with (see Figure 1) and introduce the function where Notice that the functions and form a conjugate pair (see [20, page 120]) since

The next statement was proved in [15]. Using the functions (2.16) and (2.17), we rewrite it in a more convenient form.

Theorem 3.1. The functions and have the following representations for and : where and are the full-range Hurwitz Zeta functions (see (2.16) and (2.17)) and is the Gamma function.

3.1. Properties of the Functions and

In this subsection, we collect the properties of the functions and . Some of them were established in [15] and some are given for the first time as follows.

Properties of the Functions and
The functions and satisfy the functional differential equations where , , and the prime denotes differentiation with respect to . This follows from (3.3) and the relation .The functions are even and the functions are odd on (see [15] and Figures 14).The function is periodic with the period = . It follows from the periodicity of with the unit period.(4)For all (see [15]), These properties are reflected on the graphs (see Figures 3 and 4).(5)For all and , These relations follow from the differentiation of the identities in (3.7).(6)For all and ,
Moreover, where is any sequence of positive constants, are any real constants, and the series converges. These relations follow from (3.7) and Theorem 2.2.(7)The functional sequence of Riesz potentials converges pointwise to the soliton , and the functional sequence converges pointwise to the conjugate soliton for (see [15]). Notice that (see Figures 1 and 2). Here the conjugate soliton is given by Equation (3.12) can be recovered from (3.5) with thanks to (2.13).

Remark 3.2. The conjugate soliton (3.12) is an algebraic solitary wave for extended KdV: obtained by applying the Hilbert transform to (1.1) and setting . The term “algebraic solitary wave” is explained by the fact that has a decay for large . More precisely (see [15]),

(8)For , the functions and are the elements of . Moreover, they are orthogonal in the principal value sense, namely, for all ,

Equations (3.12) and (3.14) imply that . The fact that and with are the elements of follows from their asymptotics obtained in [15], namely, Orthogonality of and follows from the fact that all with are even functions and all with are odd functions of (see Property  2 and Figures 14).(9)At the point , one has for all where is the Riemann Zeta function (see [15]).

3.2. Sturm-Liouville Problem

It is convenient to represent in the exponential form, namely, where Exponential representation (3.19) allows one to deduce the boundary value problem for the functions . It turns out that for all the functions solve the equation where Here is the Wronskian of and . Below we shall use a shorter notation .

Taking into account the behavior of for large (see [15]), we can restate the obtained results in another form. Indeed, are solutions of the Sturm-Liouville problem corresponding to . Without loss of generality, we can choose the constant in (3.22) to be positive. In the next subsection, we shall prove that . It implies that and in (3.22).

It follows from (3.7) that for the functions also satisfy the zero mean condition This reflects the oscillatory behavior of for (see Figures 3 and 4).

The graph of the arctangent function is shown in Figure 5. conveniently serves as a zero counter for both functions: and . It possesses zeros at the points where has zeros and has jumps at the points where has zeros.

Remark 3.3. Observe that a general solution of (3.24) can be written in the form where , .

Remark 3.4. We would like to point out that the differential equation given by (3.21) can be factored in the following way (see [26, page 269]):

3.3. Wronskian of and

Lemma 3.5. The following properties hold for the Wronskian for and all :

Proof. We start with (3.29). Taking into account relations in (3.3) and the fact that , we can write Since the functions are even and are odd with respect to , is even for . Differentiation of (3.31) with the help of (3.6) and (3.18) yields . Next, we turn to the proof of (3.30). Since the functions and are linearly independent solutions of the equation (3.24), 0 for all and . It remains to establish the sign of the Wronskian. In view of (3.17) and (3.18), By Abel's formula, for all , where the function is continuous on . After some simplification, we have that This representation yields (3.30). The lemma is proved.

Three-dimensional graph of is given in Figure 6.

Remark 3.6. What does the positivity of the Wronskian yield for the soliton and its conjugate? For , (3.30) simplifies to read We notice that for , for , and (see Figure 2). Integrating the inequality over the interval for and the inequality over for yields the estimate

3.4. Zeros of the Functions and

This subsection is devoted to the estimates of the number of zeros for the functions in question. By a strictly monotone change of variable Equation (3.21) can be reduced to the equation where and . Therefore, any nontrivial solution of (3.39) can have not more than a finite number of zeros on any bounded interval ([24], page 323).

Theorem 3.7. Let be the number of zeros of the function on the interval , where . Then the following inequality holds: where

Proof. Since we chose in (3.22), Therefore, Therefore, by Theorem 2.1 of Section 2, for the interval we have the estimate where Here we have used the fact that for .

Theorem 3.8. The zeros of separate and are separated by those of .

Proof. This follows from Sturm's Separation Theorem (see [24], page 335).

4. Inequality for Hurwitz Zeta Functions

Here we discuss a new inequality for the Hurwitz Zeta functions which follows from Lemma 3.5. The next statement is a corollary of this lemma.

Corollary 4.1. For and with , the following inequality holds: This inequality can also be written in another form:

Proof. Dropping positive terms in front of the full-range Hurwitz Zeta functions in (3.4) and (3.5) and using (3.30) lead to

Remark 4.2. Setting and , we can rewrite (4.1) in the form It shows that, for , , Introduce the scalar product for the complex-valued functions and by the formula Then for and ,

Remark 4.3. The proof of (4.1) becomes quite difficult when one approaches it from the point of view of special functions. For example, the use of integral representations (2.21) and (2.22) yields The change of variables leads to It is not clear at all that the integral (4.9) is positive for all and . However, (4.1) shows that it is. Multiplication of the series representations (2.16) and (2.17) does not make the proof any easier.