Submitted:
26 October 2023
Posted:
27 October 2023
You are already at the latest version
Abstract
We examine a linear q-difference differential equation which is singular in complex time t at the origin. Its coefficients are polynomial in time and bounded holomorphic on horizontal strips in one complex space variable. The equation under study represents a q-analog of a singular partial differential equation, recently investigated by the author, which comprises Fuchsian operators and entails a forcing term that combines polynomial and logarithmic type functions in time. A sectorial holomorphic solution to the equation is constructed as a double complete Laplace transform in both time t and its complex logarithm log t and Fourier inverse integral in space. For a particular choice of the forcing term, this solution turns out to solve some specific nonlinear q-difference differential equation with polynomial coefficients in some positive rational power of t. Asymptotic expansions of the solution relatively to time t are investigated. A Gevrey type expansion is exhibited in a logarithmic scale. Furthermore, a formal asymptotic expansion in power scale is displayed, revealing a new fine structure involving remainders with both Gevrey and q-Gevrey type growth.
Keywords:
Asymptotic expansion
; Borel-Laplace transform
; Fourier transform
; initial value problem
; formal power series
MSC: 35C10; 35C20
1. Introduction
In this work, we draw attention to a family of singular linear difference differential equations modelled as
for vanishing initial data , where from the leading term of (1) are integers, , represent polynomials with complex coefficients and stands for a polynomial in its arguments with holomorphic coefficients relatively to the space variable z on a horizontal strip in designed as , for some prescribed width . The forcing term is a logarithmic type function represented as a polynomial in both complex time variable t and inverse of its complex logarithm with coefficients that are bounded holomorphic on the strip .
This paper is a natural sequel of the recent study [9] by the author. Indeed, in [9], we focused on the next singularly perturbed linear partial differential equations shaped as
for given initial data , for integers appearing in the principal term of (2), complex polynomials , as above and where represents a polynomial in whose coefficients are bounded holomorphic w.r.t z on the strip and relatively to a complex parameter on some fixed disc centered at 0 for some radius . The forcing term comprises coefficients that rely polynomially on complex time t, analytically in on and holomorphically in z on . This term also entails logarithmic type functions stated as truncated Laplace transforms along a fixed segment for some radius that involve the inverse complex logarithm . When the radius is taken large, the expression of the forcing term h becomes proximate to maps that are similar to the forcing term of (1) described above, namely a polynomial in both and with bounded holomorphic coefficients on .
Under suitable constraints set on the profile of (2), we were able to construct a set of genuine bounded holomorphic solutions to (2), for p in a finite subset of the natural numbers , expressed as a complete Laplace transform of integer order in the monomial , a truncated Laplace transform of order 1 in the inverse and inverse Fourier integral in the space variable z,
where the so-called Borel-Fourier map is
- –
- analytic near and relatively to and and has (at most) exponential growth of order along some well chosen unbounded sector centered at 0 and containing the halfline for , with respect to .
- –
- continuous and subjected to exponential decay in phase .
As a result, these functions define bounded holomorphic maps on domains , for well selected bounded sector edged at 0 and where is an appropriate set of bounded sectors centered at 0. At this point, it is crucial to notice that these solutions cannot be represented as complete Laplace transform in the map . It turns out that the radii are related by a rule of the form , for some suitable constant and positive integers .
Besides, asymptotic features of these solutions have been examined in [9]. It appears that the family owns asymptotic expansions of Gevrey type in two distinguished scales of functions. Indeed, for each , the partial map holds a generalized asymptotic formal expansion (in a sense defined in the classical textbooks [5] and [14])
on the domain , in the scale of logarithmic functions with bounded holomorphic coefficients on . These asymptotic expansions are revealed to be of Gevrey 1 on , giving rise to constants for which the error bounds
occur for all integers , provided that , , where stands for the Gamma function in x. On the other hand, all the partial maps , , share a common generalized asymptotic formal expansion
on , in the scale of monomial with bounded holomorphic coefficients on , for some open domain containing all the sectors , . Moreover, these asymptotic expansions happen to be of Gevrey order on each sector, meaning that constants can be pinpointed for which the error estimates
hold for all integers , whenever . At last, we proved in [9] that the coefficients and of both formal expansions and solve explicit differential recursion relations with respect to that might be handy for effective computations.
In the present investigation of the problem (1), we plan to follow a similar roadmap as in [9]. Namely, we plan to build up genuine sectorial solutions to (1) and describe their asymptotic expansions as time t borders the origin, instead of a perturbation parameter which does not appear in (1). We notice that our main problem (1) can be viewed as a analog of (2) where the Fuchsian operator is substituted by the discret dilation operator . This terminology stems from the plain observation that the quotient neighbors the derivative as q tends to 1. The problem (2) involves at first sight only powers of the basic differential operator of Fuchsian type . However, the conditions imposed on (2) allows to express it also by means of powers of the basic differential operator of so-called irregular type . The same fact is acknowledged for the problem (1) under study for which difference operators of the form where appear, see (22). These operators are labeled of irregular type in the literature by analogy with the differential case. We quote the classical textbooks [2] and [3] as references for analytic aspects of differential equations wih irregular type and the book [15] for analytic and algebraic features of difference equations with irregular type. This suggests that in the building process of the solutions to (1), the classical Laplace transform of order ought to be supplanted by a Laplace transform of order similarly to our earlier work [11] where some related initial value difference differential problem was handled.
We now describe a little more precisely the main statements of this paper achieved in Theorem 1 and Theorem 2. Namely, under fitting restrictions on the shape of (1) listed in Subsection 2.2 and complemented in the statement of Theorem 1 in Subsection 4.3, we can establish the existence of a bounded holomorphic solution to (1) on a domain , for some small radius , where stands for an open sector centered at 0 with large opening that does not contain the halfline , see (19), for thoroughly chosen directions . In addition, the map is modelled through a triple integral which entails a Fourier inverse, a Laplace and a complete Laplace transforms
where the Borel-Fourier map is
- –
- analytic on a unbounded sector centered at 0 containing the halfline with respect to where it has (at most) exponential growth of order .
- –
- analytic relatively to on some open halfstripwith small width and on a small disc .
- –
- continuous and submitted to exponential decay in phase .
At this stage, we emphasize that the geometry of the Borel space in the variable for the map differs significantly from the one of the Borel-Fourier map in (3). Indeed, the map is in general not analytic near while possesses this property. As we will realize later on, this will be the root of the dissemblances observed between the asymptotic properties of the solutions of (2) and the solution u of (1). Besides, the partial map is only holomorphic on some fixed disc but is analytic on a full halfstrip which allows the solution to be expressed as a complete Laplace transform in in direction while is represented as a truncated Laplace transform along the segment . A direct by-product of this observation is that the forcing term of (1) can be presented as an exact polynomial in both time t and inverse complex logarithm while the forcing term has to be only considered as proximate to such a polynomial in t and . Some interesting aftermath is reached when is chosen a mere monomial in t and since in that case solves an explicit nonlinear ordinary differential equation with polynomial coefficients in some positive rational power , , displayed in (34). As a result, turns out to be an exact holomorphic solution to some specific nonlinear difference differential equation with bounded holomorphic coefficients with respect to z on and polynomial in , stated in (36). Contrastingly, the equation (2) becomes close to some nonlinear partial differential equation as but no information can be extracted about the existence of an exact genuine solution to the limit nonlinear problem.
It is worthwhile noting that in the recent years much attention has been drawn on on nonlinear difference equations and especially on those related to the so-called Painlevé equations. For a comprehensive overview on major studies for Painlevé equations and more generally for integrable discrete dynamical systems, we refer to the book [7]. In this trend of research we quote the novel paper [6] where the authors construct convergent generalized power series with complex exponents on sectors that are solutions to nonlinear algebraic difference equations. In the context of nonlinear difference differential equations we mention an important result by H. Yamazawa obtained in [17]. Indeed, he considers equations with the shape
for , , for some integers , some real number , where F is a well prepared analytic function in its arguments. Under non resonance conditions of the so-called characteristic exponent associated to (9) at , he has constructed convergent logarithmic type solutions of the form
where the coefficients and are holomorphic on a common disc and where stands for a analog of the characteristic exponent .
In the second part of Theorem 1, we exhibit for the solution of (1) a generalized asymptotic expansion of Gevrey type in a logarithmic scale for t in the vicinity of 0. The statement is similar to the one reached in [9] for the solutions of (2). Indeed, the partial map is shown to possess a generalized formal series
with bounded holomorphic coefficients on the domain as asymptotic expansion of Gevrey order 1 with respect to t on , leading to estimates of the form
for some constants , for all integers , whenever . Furthermore, in Section 4.4, Proposition 6, we provide explicit and simple difference and differential recursion relations displayed in (154) and (155) for the coefficients , , intended for practical use. The existence of such a formal expression (10) is shown in a comparable way as (4) for the partial maps in the problem (2). Namely, it is based on sharp estimates of some exponential decay for the differences of neighboring analytic solutions , disclosed in (120), to some related difference differential equation which comprises an homography action, see (116) and (118) in Proposition 4. In the process, we use a classical result known as the Ramis-Sibuya theorem (see Theorem (R.S.) in the subsection 4.2) which ensures the existence of a common Gevrey asymptotic expansion for families of sectorial holomorphic functions.
In the second main result of this paper, stated in Theorem 2, a generalized asymptotic expansion of the solution is established in the scale of monomials . This statement differs notably from the one obtained for the partial maps in the problem (2). Namely, the holomorphic solution to (1) can be split into a sum where
- the map owns a formal expressionwith bounded holomorphic coefficients on the domain as generalized asymptotic expansion of so-called Gevrey order . It means that two constants can be found with the error boundsfor all integers , all .
- the map has the null formal series as asymptotic expansion of order 1 in a logarithmic scale as t tends to 0. Indeed, two constants can be sorted with the estimatesfor all integers , provided that .
At this point, we stress the fact that the generalized expansion of Gevrey type (7) obtained for the solutions of (2) in the monomial scale are obtained by means of the Ramis-Sibuya theorem (see Theorem (R.S.) in Subsection 4.2) through precise estimates of some exponential decay for the differences of the consecutive maps relatively to on the intersections . These estimates were achieved according to the fact that the Borel-Fourier maps are analytic at in (3). In contrast, for the problem (1) under study, as observed earlier in this introduction, any of the partial Borel-Fourier map appearing in (8) for any admissible direction is not analytic near , only on sectors centered at 0. Therefore, no bounds for differences of solutions to (1) for different directions can be rooted out and the recipe using the Ramis-Sibuya theorem fails to be applied. Instead we introduce a new approach based on a specific splitting of the triple integral (8) defining and on the observation that the partial map can be analytically continued near provided that remains on the small disc , see Proposition 10. Besides, whereas explicit differential recursions could be provided for the coefficients , of the formal expansions (6), no such relations are achieved for the coefficients , of (12). However, explicit formulas (displayed in (207)) for , , can be presented as double truncated Laplace, Laplace transforms and inverse Fourier integral of derivatives of the partial Borel-Fourier map at the origin.
2. Setup of the main initial value problem and an associated set of difference-differential problems with homography action
2.1. Accounts on Laplace transforms of order k and Fourier inverse maps
This concise subsection presents the basic material about Laplace transforms and Fourier inverse maps that will be handled to built up the solution of our main problem under study.
Let be an integer and set a positive real number. We present the definition of a Laplace transform of order k as described in our former work [10]. In the construction of this analog of the classical Laplace transform of order k, the Jacobi Theta function of order k defined as the Laurent series
for any plays a prominent role.
We remind that the set of zeros of this analytic function is given by and is contained on the real line . The next lower bounds for the Jacobi Theta function attesting its so-called exponential growth of order k on a domain bypassing this set of zeros are essential. Let . A constant relying on and independent of can be chosen such that
provided that with for all .
Definition 1.
Let be a disc of some radius centered at 0 and be an open unbounded sector edged at 0 with bisecting direction in . Let us consider a holomorphic function assumed to be continuous up to the closure and subjected to the bounds
for all , for some given positive constants , and some integer . We select some direction such that . The Laplace transform of order k of f in direction γ is assigned as
where stands for a halfline in direction γ.
Let be some fixed real number. The integral transform represents a bounded holomorphic function on the domain , for any radius constrained by
and where
In the special case is an entire function with Taylor expansion conforming to the bounds (16), its Laplace transform of order k, (17) does not depend on the direction and defines a bounded holomorphic function on under the restriction (18) which possesses a Taylor expansion given by the convergent series .
The next Banach space of continuous function on with exponential decay was introduced in [4].
Definition 2.
Let be real numbers. We denote the vector space of continuous functions such that
is finite. The space endowed with the norm becomes a Banach space.
We recall the definition of the inverse Fourier transform acting on the space .
Definition 3.
Let with , . The inverse Fourier transform of f is given by
for all . The function extends to an analytic bounded function on the strips
for all given .
The next lemma described how the inverse Fourier integral is transformed under the action differential operators and products.
Lemma 1.
a) Let f be an element of for , . Define the function which belongs to the space . Then, the next identity
occurs for all , for any .b) Take and set the convolution product of f and g
Then, ψ belongs to and moreover, the next equality
holds for all , provided that .
2.2. Layout of the main problem
Throughout this subsection, we unveil the principal initial value problem under investigation in this work. It is shaped as follows,
for vanishing data , where stands for the difference operator acting on t by means of for some given real number .
The set I represents a finite subset of and are positive integers that are subjected to the next list of technical constraints:
- The inequalityholds for all .
- The restrictionsare required for all .
The maps , and for are polynomial required to fufill the next features:
- The degrees of Q and of are constrained by the relationfor all .
- We assume the existence of an open sectorial domain with inner radius (resp. outer radius ) given byfor some opening , which satisfies the next inclusionfor all . Furthermore, the inner and outer radii together with the aperture of will be suitably constrained later on in the work.
The coefficients , , are built up through the next procedure. For , we consider maps that belong to the Banach space , for given real numbers and constrained to
for all . We introduce the constants
for all on which restrictions will be set in due course of the paper. We define the coefficient as the inverse Fourier transform
for all , provided that . According to Definition 3, the maps stand for bounded holomorphic functions on the strips for any prescribed .
The forcing term is described in term of the next construction. Let be finite subsets of the positive natural numbers . For all , , we deem some maps which appertain to for and given above. We introduce the next polynomial
in the variables , with coefficients in . We bring in the map
where stands for a halfline in some given direction and is the negative real axis.
Owing to Definition 1 and Definition 3, this map is well defined provided that
- the variable belongs to , for any fixed and radius subjected to (18) where , for any given ,
- the variable is not vanishing and obeys the constraint , for some ,
- the variable z is kept in the strip for any .
However, we can further simplify the expression of . Taking heed of Definition 1, we notice that turns out to be a polynomial in ,
whose coefficients are expressed through sums over of Laplace transforms in direction ,
with bounded holomorphic coefficients on . Besides, according to the definition of the Gamma function and Cauchy’s theorem, we acknowledge that
provided that with . On that account, it follows that can be expanded as a polynomial in both variables and with bounded coefficients on , for . Namely, we get
where we define
for all and . At last, we configure the forcing term as the logarithmic type function
Here stands for the principal value of logarithm, namely provided that . Furthermore, we observe that
for some , whenever and close enough to 0.
In the particular case and for some positive integers , we make the noteworthy remark that the solution of the linear main equation (21) actually solves a special nonlinear difference-differential equation with polynomial coefficients in some positive rational power of time t stated in (36). Indeed, let the forcing term have the particular shape
where is given by the expression (31). By direct computation, we check that the forcing term satisfies the next nonlinear ordinary differential equation with polynomial coefficients in
Let us recast the main equation (21) in the form
where the difference-differential operator P is polynomial in t, with bounded holomorphic coefficients in z on the strip , for . The combination of (34) and (35) gives rise to the next nonlinear equation
2.3. A set of related difference-differential equations with an homography action
In this subsection, the main problem is embedded in a set of auxiliary problems which comprise three independent complex variables which will be the object of study in the forthcoming sections.
We seek for solutions to (21) for prescribed vanishing initial data at of the form
for some expression in the three independent variables and z.
The next computations hold for any given rational number ,
where
- the dilation acts on relatively to through ,
- the homography is applied on with respect to the variable by means of
As a result, it follows that the expression (formally) solves the equation (21) under the condition if the expression fulfills the next equation
under the constraint . Later on, we plan to build a genuine solution to (21) and in order to investigate its asymptotic expansion in some particular scale described in Subsection 4.3, we are required to complement the above single equation (37) by a family of auxiliary problems stated underneath.
For any given direction with (modulo ) and given positive real radius , we define a new forcing term
where stands for a segment of length in direction and is the halfline appearing in the formula (29). Owing to Definition 1, we notice that the map does not rely on the direction . However, it hinges on the direction and radius . We display the next problem
for given vanishing initial data .
3. Analytic solutions to the associated set of difference and differential problems under homography action
In this section, we intend to exhibit analytic solutions to the problems (37) and (39) we came up with in Subsection 2.3.
3.1. Profile of the analytic solutions and joint convolution difference equations
We search for a solution to (37) (resp. (39) for modulo ) in the form of a double Laplace, Laplace transform and inverse Fourier integral with the shape
(resp.
Here it is assumed (this fact will be justified later on in the work) that the so-called Borel-Fourier map appertains to a Banach space labelled which consists in functions with so-called exponential growth of order w.r.t , exponential growth in and exponential decay relatively to the mode m. This space is described in the next
Definition 4.
We consider the constants as prescribed in Section 2. Let and , be real numbers. We set as an unbounded sector edged at 0 with bisecting direction . We introduce the open half strip
for some given real width . We denote the vector space of all valued continuous maps on the domain , holomorphic w.r.t on the product , such that the norm
is finite. The vector space equipped with the norm represents a Banach space.
Our main objective is to establish some convolution difference equation that the Borel-Fourier map is asked to obey. On the way, we need some additional features on the Laplace transforms under multiplication by a monomial and action of difference operators. These properties have already been discussed in our past work [10]. Besides, we describe the action of the homography relatively to the variable on both expressions (40) and (41).
Lemma 2.
Let the map supposed to belong to the Banach space . Then, the next identities hold.
In line with the above technical lemma together with Lemma 1, the next statement follows.
3.2. Solving the convolution difference equation (48) on unbounded sectors and half strips
In the course of this subsection, we prove the existence and unicity of a solution to the convolution difference reached in Lemma 3.
Our scheme consists in reorganizing the equation (48) as a fixed point equation (displayed later on in (98)). On the way, we are asked to divide our equation by the next Fourier mode depending map with two complex variables
provided that and . An essential factorisation of the above map is provided in the next lemma.
Lemma 4.
For a convenient choice of the inner radius , outer radius and aperture of set up in (25), one can distinguish an unbounded sector edged at 0 with suitable bisecting direction along with an appropriate strip and a small radius ρ for which the next splitting of the map holds. Let written in the factorized form
for some radius and complex number with . Let us take . Then, one can decompose in the form
for some well chosen complex number (depending on and m and which remains bounded relatively to m), for some , close to 0 and some (for some fixed constant ). With the above factorisations (50), (51), one can express the map in the form of a non vanishing product
for given and .
Proof.
We choose appropriately the sectorial domain given in Subsection 2.2 and select an unbounded sector edged at 0 with bisecting direction chosen in a way that the next constraint
holds for all , all for some small positive numbers . Let be given. We can factorize it in the form (50). We set
Notice that remains bounded and penned in a small domain we denote which is located at some small positive distance of the real axis, when m spans the real numbers according to the condition (25) imposed.
In the next step, we select the strip and the disc in a way that
As a result, when one takes some element , we can write it in the form (51) for some real number for some and some real number that can be chosen close to 0.
By construction of , we get in particular that
In consequence of the combined factorisations (50) and (51) together with the above identity (56), the next computations hold
which is exactly the announced expression (52). In particular, this product is non vanishing since , for all owing to (25) and considering that but close to the origin, the piece enclosed by brackets in (52) cannot vanish. □
Let us consider the next linear map
In the next proposition it is stated that the map stands for a shrinking map on some fittingly chosen ball of the Banach space discussed in Definition 4.
Proposition 1.
We select the sectorial domain , the unbounded sector together with the strip and the disc as in Lemma 4. Then, provided that the constants displayed in (27) are small enough, for , an adequate radius can be chosen for which the map enjoys the next two properties
- The inclusionis granted, where denotes the closed ball of radius ϖ centered at 0 in the space .
- The Lipschitz conditionholds for all .
Proof.
We first aim our attention to the inclusion (59). Let us prescribe some real number and take some element of subjected to the condition
We plan to disclose norm estimates for each piece of the map . We first focus on norm upper bounds for the elements involved in the sum over I. The next technical lemma is crucial in this respect.
Lemma 5.
Proof.
According to Definition 4, we next upper bounds hold for the element ,
for all , all , . We deduce first upper bounds
whenever , and . The resulting bounds (61) will be reached after several steps of computations. Namely,
1) We provide upper bounds for the function
for . Since are polynomials, we get a constant with
for all . Besides, owing to the assumption (25), a constant can be pinpointed with the lower bounds
for all and from the definition of the constants , we know that
for all . The collection of bounds (65), (66) and (67) together with the triangular inequality enable the next estimates
At last, according to Lemma 2.2 from [4] or Lemma 4 from [12], we call to mind that the quantity
is finite under the assumption (24) and (26). Therefore, a constant can be singled out with
for all .
2) We focus on upper estimates for the quantity
provided that with for some fixed real number . We need to perform the next expansions
together with
Owing to the freshman classical limit and equivalence relation as x tends to 0, we reach a two constants with
provided that . Furthermore, since and are both increasing maps on , we observe the inequalities
whenever . From the two expansions (71), (72) and the bounds (73), (74) together with the assumption (22), we arrive at the next estimates
provided that with , for some constant .
3) We supply upper bounds for the quantity
for with and where , , close to 0 and for some constant , according to the decompositions (50) and (51). We recast in the form
Taking heed of our assumption (23), we obtain a constant with
for all and , for close to 0. Furthermore, a constant can be singled out with
as long as and , for close to 0. From (78) and (79), we deduce a constant such that
whenever with and all .
4) We establish bounds for the quantity displayed in (70) provided that with , where has been fixed in 2). A mere observation yields a constant with
for all with .
5) We present bounds for the piece
for with and where , , close to 0 and for some constant , according to the decompositions (50) and (51). We rearrange as follows
Bearing in mind the condition (23), we get a constant with
provided that and , for close to 0. On the other hand, a constant can be set with
as long as and , for close to 0. Due to (84) and (85), a constant can be picked out such that
for all with and all .
6) As a consequence of the list of estimates (75), (80), (81) and (86), we obtain a constant with
for and where , , close to 0 and for some constant , according to the decompositions (50) and (51).
In conclusion, on the basis of the factorization (52) for the map together with the bounds (69) and (87) combined with the bounds (63), we arrive at the next inequality
for all , all . Notice that this last inequality is tantamount to the awaited bounds (61) for the constant
□
We need control on the norm of the last term of related to the forcing term of the equation (48).
Lemma 6.
There exists a constant such that
Proof.
In view of the factorization (52) and the definition (28) of , we notice that
for all and for which the splittings (50) and (51) hold, and all . Besides, by Definition of , constants can be found such that
for all . Furthermore, we can pinpoint a constant for which
hold where and with the decompositions (50) and (51). As a result, combining (90), (91) and (92) gives rise to the next upper estimates
for all and , where
keeping in mind that , for . At last, it remains to notice that the due inequality (89) results from (93) by taking heed of Definition 4. □
We select the constant suitably together with the constants , for , taken close enough to 0 in a way that the next inequality
holds where appears in Lemma 5 and stems from Lemma 6. Eventually, the expected inclusion (59) prompts from the bounds (61) and (89) under the restriction (94).
We discuss the second item addressing the shrinking feature (60). We take two elements in the closed ball from whose radius has been prescribed in the first item (59). According to Lemma 5, under the conditions (22), (23), (24), (25) and (26) listed in Subsection 2.2, the next inequality
holds for the constant introduced in Lemma 5. We set the constants , for , small enough allowing the next inequality
to hold. The Lipschitz property (60) is a straight consequence of (95) under the requirement (96).
In the forthcoming proposition, we provide a solution to the convolution difference equation (48) established in Lemma 3.
Proposition 2.
Let us prescribe the sectorial domain , the unbounded sector together with the strip and the disc as in Lemma 4. Then, the constants defined in (27) and a constant can be fittingly chosen in a manner that a unique solution to the convolution difference equation (48) can be built up in the space under the condition
Proof.
We select as in Proposition 1. We mind the closed ball in the Banach space which represents a complete metric space for the distance deduced from the norm. The proposition 1 states that induces a contractive map from the metric space into itself. According to the classical Banach fixed point theorem, it follows that owns a unique fixed point inside the ball , we denote . It means that
for all , and . By transfering the term
from the right to the left handside of (48) and dividing the resulting equation by the map displayed in (49), we observe that (48) can be rearranged into the fixed point equation (98). On that account, the unique fixed point obtained in precisely solve (48), which yields Proposition 2. □
3.3. Analytic solutions to the auxiliary equations (37) and (39)
In the next proposition, we craft analytic solutions to the associated set of difference and differential problems under the action of homographic maps established in Subsection 2.3.
Proposition 3.
The sectorial domain , the unbounded sector together with the strip and the disc are prescribed as in Lemma 4.
-
We define the mapwhere the Borel-Fourier map is built up in Proposition 2 and solves the convolution difference equation (48). The map (99) boasts the next two qualities
- -
- It defines a bounded holomorphic function on the product for some given , where stands for the set (19) and is a disc centered at 0 with radius subjected to the constraintand . Besides, represents a bounded sector edged at 0 with bisecting direction π with radius , submitted to the next condition: there exists some real number withfor all , where , for fixed in Definition 4.
- -
- It solves the auxiliary equation (37) for prescribed initial data .
-
For a direction (modulo ), we shape the mapwhere is the Borel-Fourier map mentioned in the above item. The map (102) enjoys the next two properties
- -
- It represents a bounded holomorphic function on the product , for the domain , disc and constant given in the first item. Furthermore, stands for a bounded sector centered at 0 with bisecting direction and with radius chosen as in the first item and subjected to the next restriction: some positive real number can be found withfor all .
- -
- It obeys the auxiliary equation (39) for given vanishing initial data .
Proof.
We discuss the first item. We parametrize and in the form and for . Then, owing to (15) and (97), we get
for all with for all and . In order to provide upper bounds for the right handside of (104), we propose the next alternative.
Assume that with as above under the constraint . Then, one can single out a constant such that
for all with .
Assume that and for a radius under the constraint (100). The next three expansions are useful. Namely,
together with
and
Since holds as x is close to 0 and owing to the classical limit , we get from (107) and (108) two constants with
for all . As a result, we get from the computation (106) and bounds (109) that
At last, from the assumption (100) and requirement , the next bounds
are deduced from (110).
On the other hand, taking heed of (101), we observe that
provided that , where and according to the claim that .
As a consequence of the above bounds (105) along with (111) and (112), we deduce that the map is well defined and represents a bounded holomorphic function on the product under the above requirements (100) and (101).
Recall that the Borel-Fourier map has been constructed as a solution of the associated convolution difference equation (48) in Proposition 2. From Lemma 3, we deduce that obeys the auxiliary equation (37) on the domain for prescribed initial data .
We turn to the second item. Let and be parametrized as follows and with , . Bearing in mind (15) and (97), we obtain a constant such that the next inequality
holds provided that and . According to (103), we notice that
under the restriction . By dint of the upper bounds (105) in a row with (111), (112) and (114), we acknowledge the fact that is bounded and stands for a holomorphic map on the product under the assumptions (100) and (103). Since the Borel-Fourier map solves the convolution difference equation (48) as shown in Proposition 2, we deduce from Lemma 3 that conforms the auxiliary equation (39) on the domain for given vanishing initial data . □
4. Construction of a holomorphic solution to the main initial value problem (21) and its Gevrey asymptotic expansion relatively to complex time t in logarithmic scale.
4.1. A finite set of genuine solutions to related initial value problems.
We restate the definition of a good covering in as described in the textbook [8], Section XI-2.
Definition 5.
Let be an integer. A set of bounded sectors edged at 0 is deemed with the next three attributes
- Any two consecutive sectors and have non empty intersection , for , where the convention is assumed.
- The intersection of any three sectors is reduced to the empty set for all distinct non negative integers less than .
- The union covers some punctured neighborhood of 0 in .
Such a set is tagged a good covering in .
A notion of fitting set of sectors is discussed in the next definition.
Definition 6.
Let be an integer. A finite set of bounded sectors is minded with the next three constraints.
- For each , the sector is edged at 0, with bisecting direction and is subjected to the condition that some real number can be singled out withfor all .
- There exists an index with . All the sectors , have the same radius which obeys the restrictionwhere is introduced in the above item and is declared in Definition 4.
- The set forms a good covering in in the sense of Definition 5.
A set endowed with the above three features is called a fitting set of sectors.
In the oncoming proposition, we exhibit analytic solutions to the auxiliary problems (37) and (39) where the directions span the set of bisecting directions of some fitting set of sectors. Furthermore, sharp estimates of their consecutive differences are provided which are essential in the study of their asymptotic expansions in the variable that will be described in the next Subsection 4.2.
Proposition 4.
Let the sectorial domain , the unbounded sector together with the strip and the disc be arranged as in Lemma 4. Consider a fitting set of sectors and assign a radius a with . Then, provided that the constants are taken close enough to 0 in accordance with the requirements of Proposition 2, the properties described in the forthcoming three items hold.
- For each (where stems from Definition 6 2.) the equationwhere the forcing term is given by the triple integral formula (38), possesses a bounded holomorphic solution on the domain , where stands for the set (19), for a radius fulfilling (100), which observes the condition . Furthermore, the map is embodied in a Fourier inverse and a double Laplace, Laplace transformwhere the Borel-Fourier map belongs to the Banach space (introduced in Definition 4) constrained to the bounds (97).
- The equationwith forcing term is displayed in (29) and expressed as a polynomial in (30), holds a bounded holomorphic solution on the domain where the set and radius are given in the above item, under the vanishing condition . In addition, the map is expressed through a Fourier inverse and a double Laplace, Laplace transformwhere the Borel-Fourier map is described in the former item.
- The neighboring differences of the maps are controlled by the next bounds. For all , two constants can be found such thatfor all , all , provided that for a well chosen radius . Here we adopt the convention that .
Proof.
The first two items are direct corollaries of the statement of Proposition 3 and the definition of a fitting set of sectors chosen at the onset of Proposition 4.
We focus on the third item which demands more labor and hinges on paths deformations arguments. We distinguish two different situations.
Case 1.
Let or . We discuss only the subcase since the other alternative can be treated in a similar manner. By construction, we notice that (modulo ). According to Proposition 2, for any prescribed and , the partial map is analytic on the union . As a result, the oriented path can be bent into the union of
- –
- The halfline
- –
- The arc of circle
and the classical Cauchy’s theorem enables the difference to be reorganized as a sum of two contributions. Namely,
for all , all and . We need to control the first piece of (121)
Drew on the bounds (104), (105) together with (111) and (112), we split the halfline in the union of two segments and and we are reduced to provide bounds for the next two quantities and for
where
and
Indeed,
and
Now, we set for some real number . Hence,
provided that . As a result of (124), (125) and (126), we deduce from the splitting (123) that
for all , all and , where
In the next step, we display bounds for the second piece of (121)
According to Definition 6 1. of fitting set of sectors, we notice that the lower bounds
for all whenever the angle belongs to . By breaking up the halfine into the segments and , similar computations as above yield the bounds
for all , all , as long as .
Case 2.
Assume that . We observe that both directions and are not equal to modulo . Owing to Proposition 2, for any fixed and , the partial map is analytic on the disc . On these grounds, we can deform the oriented path into a single arc of circle
and rewrite the difference as a single triple path integral
for all , all and . Upper bounds are asked for the quantity
The definition 6 of fitting sets of sectors allows the next lower bounds
to hold for all whenever the angle is taken in . Using the partition of the halfline in two segments and , comparable estimates as the ones performed in the case 1. give rise to the next bounds
for all , all , provided that , where is given by the expression (128).
4.2. Gevrey asymptotic expansions for the bounded holomorphic solutions to the family of auxiliary problems (116) and (118).
In the next proposition, asymptotic expansions of Gevrey type are achieved for the maps , that are displayed in Proposition 4, relatively to the variable .
Proposition 5.
For the constants and fixed in Proposition 4, we denote the Banach space of valued bounded holomorphic functions on the product endowed with the sup norm. Then, for all , the partial maps , viewed as bounded holomorphic maps from the bounded sector into , share a common formal power series
with coefficients , , that belong to , as Gevrey asymptotic expansion of order 1 on . It means that, for each , two constants can be chosen in a way that the next error bounds
hold for all integers , all , whenever and .
Proof.
In the proof, we apply the next result known as the Ramis-Sibuya theorem that we rephrase for the sake of completeness and clarity for the reader (see Lemma XI-2-6 in [8]).
Theorem (R.S.)
Let be a Banach space over the field of complex numbers and let be a good covering in as outlined in Definition 5. For all , we consider holomorphic functions that enjoy the next two features
- The maps are bounded on for all .
- The difference stands for a holomorphic map on the intersection which is exponentially flat of order k, for some integer , meaning that one can select two constants for whichholds provided that , for all . By convention, we set and .
Then, a formal power series with coefficients belonging to can be singled out, which is the common Gevrey asymptotic expansion of order relatively to u on for all the maps , for . It attests that two constants can be chosen with the result that the error bounds
hold for all integers , all , all .
For each , we introduce the map set as
In view of Proposition 4, we acknowledge that
- –
- The set of sectors forms a good covering in owing to Definition 6 3.
- –
- For each , the map is bounded holomorphic on the sector .
- –
- For each , the difference suffers the boundsfor the constants and displayed in (120), provided that .
Thereupon, the claims 1. and 2. of Theorem (R.S) are matched for the family of maps with the constant . The existence of the formal power series (136) which represents the collective Gevrey asymptotic expansion of order 1 relatively to on for all the maps , follows. As a result, the error bounds (137) are warranted. □
4.3. Statement of the first main result.
In this subsection, a bounded holomorphic solution to our main initial value problem (21) is shaped. This solution is favored with an asymptotic expansion in some logarithmic scale that reveals to be of Gevrey type. The next theorem represents the first main achievement of our work.
Theorem 1.
Let the sectorial domain , the unbounded sector together with the strip and the disc be duly prescribed as in Lemma 4. Then, assuming that the constants are in the vicinity of 0 as specified by the requirements of Proposition 2 and that the radius is close enough to 0, the equation
has a bounded holomorphic solution on the domain for vanishing initial data . In addition, the map can be expressed as a triple integral comprising a Fourier inverse, a Laplace and Laplace transforms
where the Borel-Fourier map originates from the Banach space (see Definition 4) and is restrained to the bounds (97).
The function enjoys a generalized asymptotic expansion of Gevrey type in a logarithmic scale as t tends to 0. More precisely, one can single out a formal series
with bounded holomorphic coefficients on the domain , which stands for an asymptotic expansion of Gevrey order 1 in the scale of logarithmic functions of the map with respect to t on the domain . In other words, two constants can be found with the aim that the next error bounds
hold for all integers , all , provided that .
Proof.
We select a fitting set of sectors and we take the index for which according to Definition 6 2. By definition of the principal value of the logarithm , for , whenever , we check that
as long as , provided that we take sufficiently close to 0, where has been disclosed in the third item of Proposition 4. We define
where the map is described in the second item of Proposition 4. By construction of , we ascertain that represents a bounded holomorphic function on the product .
Besides, according to the second item of Proposition 4, we know that the map stands for a solution to the equation (118) on the domain . On the basis of the computations made in Subsection 2.3, we deduce that the map solves the main equation (21) on the domain , constrained to the initial value condition .
4.4. Computational features related to the formal power series (136).
In this subsection, we establish that the formal series (136) which represent the asymptotic expansion of Gevrey type for the holomorphic maps actually solve some functional partial differential equation. On the journey, we notice that its coefficients , fulfill some handy recursion relations that might be of interest for concrete applications.
Proposition 6.
Proof.
We depart from the equation (118) recast in the form
provided that , and . We remind the reader the next useful classical result which relates the coefficients of an asymptotic expansion of a holomorphic map f to its high order derivatives.
Proposition
([2], Proposition 8, p. 66) Let be a holomorphic map from a bounded open sector G centered at 0 into a complex Banach space endowed with a norm . The following two statements are equivalent
- There exists a formal power series with coefficients in subjected to the next feature. For all closed subsector S of G centered at 0, there exists a sequence of positive real numbers such thatfor all , all integers .
- All derivatives of order n, are continuous at the origin and there exists a sequence of elements in such thatfor all integers .
As a result of the above proposition, we deduce from the asymptotic expansion (137) in the particular case (meaning that ) the next limits
for all integers . On the basis of the above limits, in order to reach recursion relations for the coefficients , , our strategy consists in searching for recursion relations for the related th derivatives of the map relatively to . On the way, we need to take the th derivative with respect to of the left and right handside of the equation (147). However, the equation (147) involves composition of with explicit homographic maps and we are asked to explicitely compute their higher order derivatives. In order to overcome this difficulty, we will apply a rule to evaluate high order derivatives of compositions of functions which has been introduced in [16] and is suitable for Gevrey estimates. This rather new identity allows us to avoid computations with the cumbersome combinatorial classical Faa-Di-Bruno formula and enables us to present very practical recursions relations. Indeed, we recall this higher order chain rule (Theorem 2.1 in [16]) under stronger assumptions (which will be sufficient for our scope) as stated in the previous work of the author [13].
Lemma 7.
Let be open sets in . Let and be holomorphic functions. Then, the n-th order derivative of the composite function is given by the formula
for all integers and .
In the next lemma, we perform an auxiliary computation which entails the homographic maps appearing in the main equation (147).
Lemma 8.
For any integer , we set
Then, for all integers with , the next identity
holds for all , with the convention that when .
Proof.
Direct computations show that
and hence
for all . We deduce that
for all and all integers . It follows from (150) that
which coincides with the formula (149) in the case under the convention that . On the other hand, when , we deduce from (150) that
which yields the awaited identity (149) by setting in the formula (152). □
On the ground of the above lemmas and based on equation (147), we can derive some recursion relation on the sequence of th derivatives of with respect to . Namely,
for all , all , all and all .
In the next step, we let tend to 0 on the sector in both identities (147) and (153). According to the limits (148) and bearing in mind that the maps and are holomorphic relatively to , we get the next recursion relations for the coefficients , . Namely,
together with
for all , all and .
In the last part of the proof, we show that the formal power series (145) obey the functional equation (146). Our approach hinges on the next technical lemma where the Taylor expansion of the composition of the formal series (145) with some homographic map is explicitely computed.
Lemma 9.
Let be an integer. The next formal Taylor expansion
holds, for all and .
Proof.
By mere composition, we notice that
On the other hand, the geometric series allows to write
and taking its derivative of order with respect to yields the expansion
with the notation if and if , for all . From (159), for all integers , we deduce the next identity
As a result of (157) and (160), we deduce that
Besides, by straight calculus, we observe that
for all and . Eventually, the combination of (161) and (162) yields the awaited formal Taylor expansion (156). □
According to the fact observed in (30) that the map defines a polynomial in the variable , it follows that its Taylor expansion
is convergent (and actually a finite sum) near the origin with respect to , for all and .
5. Fine structure of Gevrey/Gevrey asymptotic expansions in combined power and logarithmic scales for the holomorphic solution to the initial value problem (21).
5.1. Solving the convolution difference equation (48) on some neighborhood of the origin
In order to study the equation (48) in the Borel space near the origin in and Fourier space on , we introduce the next Banach space.
Definition 7.
Let be real numbers. For a given real number , we denote the vector space of all continuous valued functions on , holomorphic with respect to on , such that the norm
is finite. The vector space endowed with the norm is a Banach space.
We plan to solve the next convolution difference equation
provided that , and , with some function in the Banach space .
In preparation for achieving our goal, we rearrange the equation (165) as a fixed point equation (disclosed later on in (188)). Along the road, we need to divide our equation by the map displayed in (49) whenever , and the mode m belongs to . Lower bounds for the map are provided in the next lemma.
Lemma 10.
Let the inner radius , outer radius and aperture of introduced in Subsection 2.2 be chosen as in Lemma 4. Let be the radius fixed in Lemma 4. Then, for a proper choice of radius , taken close enough to 0, one can find a constant with
for all , all , all .
Proof.
Take a fixed . We introduce the complex number
Observe that remains bounded and parked in a small domain we denote which is located at some small positive distance of the origin, when m varies within the real numbers, owing to the requirement (25). We select the radius accordingly to the condition
Now, let us take an arbitrary complex number . We decompose it in the form
for some real numbers close to 0 for given above. By construction of in (167), the next identity
holds. Select some arbitrary . We split it in a factorized form
for some angle and radius with the constraint . The combined splitting (169) and (171) together with the identity (170) enables the factorisation of the map
Besides, provided that the radius is chosen in the vicinity of the origin, we can find a constant with
for all , all , all close to 0. At last, the factorization (172) and the lower bounds (173) give rise to (166). □
In the ongoing proposition, we check that the map introduced in (58) represents a shrinking map on some appropriately selected ball in the Banach space examined in Definition 7.
Proposition 7.
We fix the sectorial domain and the radius ρ, b as in Lemma 10. Let be real numbers fixed as in Subsection 2.2. Then, assuming that the constants presented in (27) are small enough, for , for all radius chosen large enough, the map given by (58) is favoured with the next two features
- The inclusionis granted, where denotes the closed ball of radius centered at 0 in the space .
- The Lipschitz conditionholds for all .
In particular, since the radius can be taken arbitrarily large, we observe that the map turns out to be well defined on the whole space where the shrinking property (175) holds true.
Proof.
Let us focus on the first item of the proposition. We first provide bounds for the forcing term of disclosed in the next
Lemma 11.
There exists a constant such that
Proof.
In the next lemma, we come up with bounds for the linear part of the map .
Lemma 12.
One can find a constant such that
for all .
Proof.
Let us take . We provide bounds for the function
By definition of the space , we notice that
for all , all and all . Owing to the assumption (22), we notice that provided that . Hence,
whenever , and . Then, according to the lower bounds (166) together with (180), we deduce that
and bearing in mind the estimates (69) where the map is introduced in (64), we reach
for all , and . At last, we arrive at some constant for which the norm bounds
holds. □
Now, we select the constants , for , small enough and take a radius large enough in a way that the next inequality
holds where the constant appears in Lemma 12 and shows up in Lemma 11. Eventually, the bounds (176) along with (178) under the restriction (184) trigger the expected inclusion (174).
In the second part of the proof, we address the shrinking property (175). Let us choose two arbitrary elements , in the closed ball whose radius has been prescribed in the first item (174). Owing to Lemma 12, the following inequality
holds for the constant stemming from Lemma 12. We prescribe the constants , for , small enough allowing the next inequality
to hold. The Lipschitz property (175) is a straight consequence of (185) under the requirement (186).
The next proposition provides a solution to the convolution difference equation (165) inside the space .
Proposition 8.
We prescribe the sectorial domain together with the radius ρ, b as in Lemma 10. Let be real numbers fixed as in Subsection 2.2. Assume that the constants , , are chosen small enough in a suitable way as in Proposition 7. Then, for all radius large enough, a unique solution to the convolution difference equation (165) can be constructed in the space under the requirement
Proof.
Select a radius as in Proposition 7. The closed ball stands for a complete metric space for the distance . The proposition 7 claims that the map induces a contractive map from the metric space into itself. The classical Banach fixed point theorem allows the map to possess a unique fixed point located inside de ball that we denote . As a result, the next identity
holds provided , , for all . At last, under the conditions imposed, we observe that the convolution difference equation (165) can exactly be rearranged after a division by the map as (188). As a consequence, the unique fixed point obtained in fully solves (165). This yields Proposition 8. □
5.2. Link between the solutions and to the convolution difference equation (48), (165).
In order to unveil the analytic relation between the two solutions and to the same convolution difference equation considered in Subsection 3.2 and Subsection 5.1, we introduce a new auxiliary Banach space.
Definition 8.
Let be given positive real numbers and let be an unbounded sector edged at 0 with bisecting direction . We denote the vector space of all continuous maps on the product , holomorphic relatively to the couple on the domain , for which the norm
is a finite quantity. The vector space equipped with the norm is a Banach space.
In the next proposition, we claim that the map displayed in (58) is well defined on the space where it boasts a Lipschitz feature.
Proposition 9.
We prescribe the sectorial domain and the radius b, ρ as in Lemma 10. We set the constants as in Subsection 2.2. We select an unbounded sector as in Lemma 4. Then, assuming that the constants introduced in (27) are close enough to 0, for all , the map declared in (58) is well defined on the whole space and is subjected to the next Lipschitz condition
for all , belonging to .
Proof.
The proof of Proposition 9 mirrors in the very details the one of Proposition 7 and will not be presented in this work in order to avoid redundancy. □
The following proposition establish the awaited analytical connection between and .
Proposition 10.
Let the sectorial domain and the radius b, ρ be prescribed as in Lemma 10. The constants are set as in Subsection 2.2 and the unbounded sector is chosen as in Lemma 4. Then, provided that the constants given by (27) are taken in the vicinity of the origin for all , the next identity
holds for all , all , all . In particular, for given and , the partial map is the analytic continuation of the partial map on the full disc .
Proof.
According to Proposition 2, we know that the map belongs to the Banach space . According to Definition 8 it follows that the restricted map , for , and belongs to . On the other hand, we know from Proposition 8 that the map belongs to the space . As a result, the restricted map on also belongs to . Furthermore, according to (98) and to (188), we observe in particular that the next two idendities
holds as functions provided that , and . At last, if one sets and in the inequality (190), it follows from (192) that
It implies that , from which the expected identity (191) follows. □
5.3. Statement of the second main result.
In this subsection, we exhibit a fine structure for the asymptotic expansion of Gevrey/Gevrey type for the solution to the equation (139) which combines both a logarithmic scale and a power scale. The next statement represents the second deed of our work.
Theorem 2.
We consider the function displayed in (140) which solves our main initial value problem (139) for vanishing initial data built up in Theorem 1. Then, the map can be broken up as a sum of two functions
where
- –
- the map is bounded holomorphic on the domain and possesses a generalized asymptotic expansion of so-called Gevrey type in a power scale as t tends to 0. It means that one can distinguish a formal power serieswith bounded coefficients on the domain which represents a generalized asymptotic expansion of Gevrey order in the scale of monomials of the map with respect to t on the domain . Namely, two constants can be singled out for which the next error boundshold for all integers , all , provided that .
- –
- the map is bounded holomorphic on the domain and carries the null formal series as asymptotic expansion of Gevrey order 1 in a logarithmic scale as t tends to 0. In other words, two constants can be identified in order that the following error boundshold for all integers , all , as long as .
Proof.
Our idea consists in the splitting of the triple integral representation of given by (140) into three specific contributions
where
and
in a row with
where the integration paths are stated as follows
along with
where the positive real numbers are prescribed in Lemma 10.
In the next first main proposition, we provide asymptotic expansions for the first piece relatively to t.
Proposition 11.
There exists a sequence of maps , , that are well defined and bounded holomorphic relatively to on the product which are submitted to the bounds
for some well selected constants and , where and are the two constants arising in (15), for all integers , provided that and . For any given natural number , the next decomposition
holds for on , where the remainder term stands for a bounded holomorphic function on and is monitored by means of the bounds
for the constants appearing in (201) and for a suitable small radius , as long as and .
Proof.
Let be fixed as in Lemma 10. Owing to Proposition 8, we know in particular that the partial map is bounded and analytic on the disc for any prescribed and . As a result, we can apply the Taylor formula with integral remainder of some fixed order to that function and get the next expansion
provided that , and . According to Proposition 10, we know that the function coincides with the map for , all and . Hence, from the identity (204), we deduce the next development
for all , all and all . This last formula (205) enable the expansion of the map in the form
where
for .
In the next step, we provide upper bounds for the maps , . We first need to remind the reader the next formula
for all which has been applied in our recent work [10], see Lemma 3 therein, from which we deduce the splitting
for all . As a result, one can further break up the term as follows
where
In the next lemma, we focus on bounds for the function .
Lemma 13.
For all , the map is well defined and bounded holomorphic with respect to on . Furthermore, there exists two constants and such that
for all on , provided that .
Proof.
We remind from Proposition 8 that the map belongs to the space and that a constant can be pinpointed with the bounds
provided that , and . Besides, from the classical Cauchy’s formula, we know the next integral representation
to hold for and , where the integration is realized along any positively oriented circle centered at 0 with radius R subjected to . On account of (213) and the bounds (212), we reach the estimates
for all , all . As a result of (214), we arrive at
where
for all , all with . □
In the next lemma, bounds for the second piece of (209) are determined.
Lemma 14.
For all , the map
is well defined and stand for a bounded holomorphic function relatively to on . In addition, the next upper bounds
hold for all , , where the constants and are prescribed in Lemma 13 and where and are the two constants appearing in (15).
Proof.
The technical estimates displayed in the next lemma are crucial.
Lemma 15.
The next inequality
holds for all , all integers , where and are the two constants appearing in (15).
Proof.
Owing to (15), we first observe that
for all and . Based on (220), we deduce that
where the quantity is derived by performing the change of variable in the integral along the segment above and stands for
In the next step, we reach upper bounds for . By coarse upper estimates, we first get
Then, at last, we show that the constant can be computed in an exact manner. Indeed, we make the change of variable
in the integral , which gives rise to
On the other hand, we recall the Gaussian identity
which is valid for any given real number , that has been already used in our former paper [11] and stems from the book [1], Chapter 10, p. 498. This last identity enables the straight computation of (224) as follows
for any integer .
In the next lemma, we address bounds for the remainder part of the expansion (206) for .
Lemma 16.
Proof.
We first need to upper bound the next quantity
relatively to t and N. Indeed, owing to (220), we deduce
where the element is obtained by applying the change of variable in the integral along the segment overhead and stands for
In the next step, we merely observe that
where is given in the inequality (223). According to the computation made in (225), we notice that
At last, with the combination of (230), (231), (232) and (233) we arrive at
provided that , for any given integer .
Besides, owing to the classical Cauchy’s formula, the next integral representation
holds for all , , and , where the integration is performed along a positively oriented circle centered at with small radius chosen in a way that . From (235) together with (212), we deduce the useful bounds
for all , , and . Eventually, the gathering of (234) and (236) gives rise to
for all and where the constant is defined in (216). □
In the second main proposition, we show that the second piece has the null formal series as Gevrey asymptotic expansion of order 1 in a logarithmic scale with respect to t.
Proposition 12.
The map is well defined and bounded holomorphic relatively to on the product . Furthermore, for some well chosen constants and any given integer , the next error bounds
hold provided that and , where and are the two constants stemming from (15).
Proof.
According to (97) in Proposition 2, one can find a constant for which the map is subjected to the next upper bounds
provided that , and . Besides, the bounds (234) for the quantity (229) in the special case yields the next estimates
for all . Furthermore, a constant can be singled out with
as long as , for chosen small enough. In the next step, we remind the reader the following technical estimates that are taken from Lemma 14 of [10]. Namely, for any given real number , one can select a constant such that
for all integers , all real numbers . Based on (242) for the constant and specific value , we deduce from (241) that
for all integers , provided that .
In the last principal proposition, the third piece is shown to have the null formal series as asymptotic expansion of Gevrey order in the scale of monomials relatively to t.
Proposition 13.
The map is well defined and bounded holomorphic relatively to on the product . In addition, for some suitable constants , , and any given integer , the next error bounds
hold provided that and , where and are the two constants appearing in (15).
Proof.
We further break up the integral in two parts
where
with and for
along the segment .
As stated in (97) in Proposition 2, the map is subjected to the next upper bounds
provided that , and . On the other hand, we need the next technical upper bounds.
Lemma 17.
One can single out two constants such that
provided that and .
Proof.
In the next lemma, we exhibit Gevrey type estimates on both segments and .
Lemma 18.
The next two Gevrey type estimates hold.
- On the segment , we get thatholds for all , all , for all integers .
- On the segment , we arrive atprovided that , all , for all integers .
Proof.
1) Consider and . In particular, we notice that and . It follows that
As a result, the inequality (249) becomes
2) Let us take . In particular . We select small enough and fulfilling (100) in a way that
for all . The inequality (249) is then changed into
The next estimates have been presented in Lemma 12 of our recent work [10]. Namely, for any prescribed real number , the next inequality
occurs for all integers , all positive real numbers . In particular, the next upper bounds
along with
hold for all , for all integers .
We return to the proof of Theorem 2. On the ground of the decomposition (197), we set
According to Proposition 11 and Proposition 13, we observe that represents a bounded holomorphic map on the domain . Moreover, is submitted to error bounds of the form (195) for the sequence of functions , given by , which represent bounded holomorphic maps on the domain , owing to the upper bounds (201).
On the other hand, we assign
As claimed by Proposition 12, we check that stands for a bounded holomorphic function on . Furthermore, is subjected to error bounds shaped in (238). Theorem 2 is established. □
References
- G. E. Andrews, R. Askey, R. Roy, Special functions. Encyclopedia of Mathematics and Applications, 71. Cambridge University Press, Cambridge, 1999.
- W. Balser, Formal power series and linear systems of meromorphic ordinary differential equations. Universitext. Springer-Verlag, New York, 2000. xviii+299 pp.
- O. Costin, Asymptotics and Borel summability, Chapman and amp; Hall/CRC Monogr. Surv. Pure Appl. Math., 141 CRC Press, Boca Raton, FL, 2009, xiv+250 pp.
- O. Costin, S. Tanveer, Short time existence and Borel summability in the Navier-Stokes equation in R3, Comm. Partial Differential Equations 34 (2009), no. 7-9, 785–817.
- A. Erdélyi, Asymptotic expansions. Dover Publications, Inc., New York, 1956. vi+108 pp.
- R. Gontsov, I. Goryuchkina, A. Lastra, On the convergence of generalized power series solutions of q-difference equations, Aequationes Math. 96 (2022), no. 3, 579–597.
- J. Hietarinta, N. Joshi, F. Nijhoff, Discrete systems and integrability. Cambridge Texts in Applied Mathematics. Cambridge University Press, Cambridge, 2016. xiii+445 pp.
- P. Hsieh, Y. Sibuya, Basic theory of ordinary differential equations. Universitext. Springer-Verlag, New York, 1999.
- S. Malek, Double-scale Gevrey asymptotics for logarithmic type solutions to singularly perturbed linear initial value problems, Results Math. 77 (2022), no. 5, Paper No. 198, 61 pp.
- S. Malek, On a partial q-analog of a singularly perturbed problem with Fuchsian and irregular time singularities, Abstr. Appl. Anal. 2020, Article ID 7985298, 32 p. (2020).
- S. Malek, On parametric Gevrey asymptotics for a q-analog of some linear initial value problem, Funkc. Ekvacioj, Ser. Int. 60, No. 1, 21–63 (2017).
- S. Malek, On Gevrey asymptotics for some nonlinear integro-differential equations. J. Dyn. Control Syst. 16 (2010), no. 3, 377–406.
- S. Malek, On functional linear partial differential equations in Gevrey spaces of holomorphic functions, Ann. Fac. Sci. Toulouse, Math. (6) 16, No. 2, 285-302 (2007).
- F. Olver, Asymptotics and special functions. Computer Science and Applied Mathematics. Academic Press, New York-London, 1974. xvi+572 pp.
- J-P. Ramis, J. Sauloy, C. Zhang, Local analytic classification of q-difference equations, Astérisque (2013), no. 355, vi+151 pp.
- T. Yamanaka, A new higher order chain rule and Gevrey class. Ann. Global Anal. Geom. 7 (1989), no. 3. 179–203.
- H. Yamazawa, Holomorphic and singular solutions of q-difference-differential equations of Briot-Bouquet type, Funkcial. Ekvac. 59 (2016), no. 2, 185–197.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.