Thomas Dreyfus

Alberto Lastra
email: alberto.lastra@uah.es

S T Éphane Malek
email: stephane.malek@math.univ-lille1.fr

ON THE MULTIPLE-SCALE ANALYSIS FOR SOME LINEAR PARTIAL q-DIFFERENCE AND DIFFERENTIAL EQUATIONS WITH HOLOMORPHIC COEFFICIENTS

Keywords: Mathematics Subject Classification. 35C10, 35C20 Asymptotic expansion, Borel-Laplace transform, Fourier transform, Formal power series, Singular perturbation, q-difference-differential equation

, where the first author makes distinction among the different q-Gevrey asymptotic levels by successive applications of two q-Borel-Laplace transforms of different orders both to the same initial problem and which can be described by means of a Newton polygon.

Introduction

This work is devoted to the study of a family of linear q-difference-differential problems in the complex domain. It can be arranged into a series of works dedicated to the asymptotic study of holomorphic solutions to different kinds of q-difference-differential problems involving irregular singularities such as [START_REF] Lastra | On q-Gevrey asymptotics for singularly perturbed q-difference-differential problems with an irregular singularity[END_REF], [START_REF]On parametric Gevrey asymptotics for singularly perturbed partial differential equations with delays[END_REF], [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF], [START_REF] Lastra | On q-asymptotics for linear q-difference-differential equations with Fuchsian and irregular singularities[END_REF], and [START_REF]On singularly perturbed q-difference-differential equations with irregular singularity[END_REF]. The study of q-difference and q-differencedifferential equations in the complex domain is a promising and fruitful domain of research. In the literature, one may find other interesting approaches to these problems. We refer to [START_REF] Ramis | About the growth of entire functions solutions of linear algebraic q-difference equations[END_REF] as a reference, and contributions in the framework of nonlinear q-analogs of Briot-Bouquet type partial differential equations in [START_REF] Tahara | q-Analogues of Laplace and Borel transforms by means of q-exponentials[END_REF]. We provide [?, ?] as novel studies in this direction.

The study of q-difference equations has also been under study in different applications in the last years. Some advances in this respect are [?, ?, ?], and the references therein.

The main aim of this work is to study a family of q-difference-differential equations of the form (1.1) Q(∂ z )σ q,t u(t, z, )

= ( t) d D 1 σ d D 1 k 1 +1 q,t R D 1 (∂ z )u(t, z, ) + ( t) d D 2 σ d D 2 k 2 +1 q,t R D 2 (∂ z )u(t, z, ) + D-1 =1 ∆ t d σ δ q,t (c (t, z, )R (∂ z )u(t, z, )) + σ q,t f (t, z, ),
where D, D 1 , D 2 are integer numbers larger than 3, Q, R D 1 , R D 2 and R for = 1, . . . , D -1 are polynomials of complex coefficients, and ∆ ≥ 0, δ , d ≥ 1 are nonnegative integers for every 1 ≤ ≤ D -1. The numbers d D 1 and d D 2 are positive integers. q stands for a real number with q > 1.

We consider the dilation operator σ q,t acting on variable t, i.e. σ q,t (t → f (t)) := f (qt), and the generalization of its composition given by σ γ q,t (t → f (t)) := f (q γ t), for any γ ∈ R. We also fix positive integer numbers k 1 and k 2 with

1 ≤ k 1 < k 2 .
As in the former work [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF] of the third author, the coefficients c (t, z, ) and the forcing term f (t, z, ) represent bounded holomorphic functions in the vicinity of the origin in C 2 w.r.t (t, ) and on a horizontal strip H β = {z ∈ C/|Im(z)| < β} of width 2β > 0 relatively to the space variable z. However, a new additional constraint is required on the growth of the Taylor expansion of each c according to the mixed variable t , see (4.6). It implies that the functions c (t, z, ) can be extended as entire functions in the monomial t in the whole plane C with so-called q-exponential growth of some order related to k 1 and k 2 (this terminology will be explained later in the paper).

Two singularly perturbed terms on the right hand side of equation (1.1) are distinguished. This makes a crucial difference with respect to the previous work [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF] in which only one term appears, whilst the multi-level q-Gevrey asymptotic behavior comes from the forcing term. More precisely, in that previous work we focused on families of q-difference-differential equations that can be factorized as a product of two operators in so-called normal forms each enjoying one single level of q-Gevrey asymptotics. In the present work, the appearance of these two terms would cause a multilevel q-Gevrey phenomenon in the study of the asymptotic solution of (1.1) regarding the perturbation parameter. Our approach is to follow a two-step procedure of summation of the formal solution, which makes the two q-Gevrey asymptotic orders emerge.

Another important difference compared to our previous contribution [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF] is that we are now able to handle holomorphic coefficients in time t whilst only polynomial coefficients were considered in [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF]. This fact relies on new technical bounds for a q-analog of the convolution of order k displayed in Proposition 2. [START_REF] Hsieh | Basic theory of ordinary differential equations[END_REF].

The point of view we use here is similar to the one performed in the work of the first author, see [START_REF] Dreyfus | Building meromorphic solutions of q-difference equations using a Borel-Laplace summation[END_REF], and is related to direct constraints on the shape of the main equation via a possible description by a Newton polygon. It is important to stress that this approach is specific to the q-difference case. Namely, such a direct procedure for producing two different Gevrey levels in the differential case for the problem stated in the work [START_REF]On parametric Gevrey asymptotics for singularly perturbed partial differential equations with delays[END_REF] is impossible due to very strong restrictions related to a formula used in the proof and appearing in [START_REF] Pravica | Solutions of a class of multiplicatively advanced differential equations[END_REF], see formula (8.7) p. 3630. In that case, only a proposal via factoring the main equation did actually work, as performed in our joint work [START_REF]On parametric multisummable formal solutions to some nonlinear initial value Cauchy problems[END_REF].

Let us briefly review the steps followed in order to achieve our main results in the present work.

Let 0 ≤ p ≤ ς -1. First, we apply q-Borel transformation of order k 1 to equation (1.1) in order to obtain our first auxiliary problem in a Borel plane, problem (4.12). A fixed point result in a complex Banach space of functions under an appropriate growth at infinity lead us to an analytic solution, w dp k 1 (τ, m, ) of (4.12). More precisely, w dp k 1 (τ, m, ) defines a continuous function defined in (U dp ∪ D(0, ρ)) × R × D(0, 0 ), where U dp is an infinite sector of bisecting direction d p , and holomorphic with respect to the variables τ and in (U dp ∪ D(0, ρ)) and D(0, 0 ), respectively. In addition to that, it holds that this function admits q-exponential growth of order κ at infinity with respect to τ in U dp .This result is described in detail in Proposition 4.2.

A second auxiliary problem in the Borel plane is constructed by applying the formal q-Borel transformation of order k 2 to the main problem (1.1). The second auxiliary equation is stated in (4.32). A second fixed point result in another appropriate Banach space of functions allow us to guarantee the existence of an actual solution of the second auxiliary problem, w dp k 2 (τ, m, ), defined in S dp × R × D(0, 0 ) and holomorphic with respect to τ and in S dp and D(0, 0 ), respectively. Here, S dp stands for an infinite sector with vertex at the origin and bisecting direction d p . Moreover, this function suffers q-exponential growth of order k 2 at infinity w.r.t. τ on S dp . This statement is proved in Proposition 4. [START_REF] Dreyfus | Building meromorphic solutions of q-difference equations using a Borel-Laplace summation[END_REF].

As a matter of fact, the key point in our reasoning is the link between the q-Laplace transform of order κ with respect to τ variable of w dp k 1 and w dp k 2 . In Proposition 4.5, we guarantee that both functions coincide in the intersection of their domain of definition. This would entail that the function L d q;1/κ (w d k 1 (τ, m, )) can be continued along direction d p , with q-exponential growth of order k 2 , see Propoposition 4.5.

The construction of the analytic solution of (1.1), u dp (t, z, ), is obtained after the application of q-Laplace transformation of order k 2 and inverse Fourier transform, providing a holomorphic function defined in T × H β × E p , where T is some well chosen bounded sector centered at 0 and {E p } 0≤p≤ς-1 represents a good covering in C * (see Definition 5.1). This result is described in Theorem 5.3. The following diagram illustrates the procedure to follow. For the attainment of Figure 1. Scheme of the different Borel levels attained in the construction of the solution the asymptotic properties of the analytic solution we make use of a Ramis-Sibuya type theorem in two levels (see Theorem 6.4), and the properties held by the difference of two analytic solutions in the intersection of their domains, whenever it is not empty. The conclusion yields two different q-Gevrey levels of asymptotic behavior of the analytic solution with respect to the formal solution depending on the geometry of the problem. The final main result states the splitting of both, the formal and the analytic solutions to the problem under study, as a sum of three terms. More precisely, if F denotes the Banach space of holomorphic and bounded functions defined in T × H β , and û(t, z, ) stands for the formal power series solution of (1.1), then it holds that û(t, z, ) = a(t, z, ) + û1 (t, z, ) + û2 (t, z, ), where a(t, z, ) ∈ F{ } and û1 (t, z, ), û2 (t, z, ) ∈ F[[ ]] and such that for every 0 ≤ p ≤ ς -1, the function u dp can be written in the form

u dp (t, z, ) = a(t, z, ) + u dp 1 (t, z, ) + u dp 2 (t, z, ), where → u dp 1 (t, z,
) is a F-valued function that admits û1 (t, z, ) as its q-Gevrey asymptotic expansion of order 1/k 1 on E p and also → u dp 2 (t, z, ) is a F-valued function that admits û2 (t, z, ) as its q-Gevrey asymptotic expansion of order 1/k 2 on E p . This corresponds to Theorem 6.6.

The paper is organized as follows. In Section 2.1, we define a weighted Banach space of continuous functions on the domain (D(0, ρ) ∪ U ) × R with q-exponential growth on the unbounded sector U with respect to the first variable, and exponential decay on R with respect to the second one. We study the continuity properties of several operators acting of this Banach space. Section 2.2. is concerned with the study of a second family of Banach spaces of functions with q-exponential growth on an infinite sector with respect to one variable and exponential decay on R with respect to the other variable. In Section 3, we recall the definitions and main properties of formal and analytic operators involved in the solution of the main equation. Namely, formal q-Borel transformation and an analytic q-Laplace transform of certain q-Gevrey orders, and inverse Fourier transform. In Section 4.1. and Section 4.2, we study the analytic solutions of two auxiliary problems in two different Borel planes and relate them via q-Laplace transformation (see Theorem 5.3). In Section 5, we describe in detail the main problem under study, and construct its analytic solution and the rate of growth of the difference of two neighboring solutions in their common domain of definition. Finally, Section 6 deals with the existence of a formal solution of the problem, and studies the asymptotic behavior relating the analytic and the formal solutions through a multilevel q-Gevrey asymptotic expansion (Theorem 6.6). This result is attained with the application of a two-level q-version of Ramis-Sibuya theorem (Theorem 6.4).

Auxiliary Banach spaces of functions

In this section, we describe auxiliary Banach spaces of functions with certain growth and decay behavior. We also provide important properties of such spaces under certain operators.

Let U d be an open unbounded sector with vertex at the origin in C, bisecting direction d ∈ R and positive opening. We take ρ > 0 and consider D(0, ρ) = {τ ∈ C : |τ | < ρ}.

We fix real numbers β, µ, δ > 0, q > 1 and α ≥ 0 through the whole section. We assume the distance from U d ∪ D(0, ρ) to the real number -δ is strictly larger than 1. Let k > 0. We denote D(0, ρ) the closure of D(0, ρ).

The next definition of a Banach space of functions, and subsequent properties have already been studied in previous works. Analogous spaces were treated in [START_REF] Lastra | On q-Gevrey asymptotics for singularly perturbed q-difference-differential problems with an irregular singularity[END_REF][START_REF]On multiscale Gevrey and q-Gevrey asymptotics for some linear q-difference differential initial value Cauchy problems[END_REF], inspired by the functional spaces appearing in [START_REF] Pravica | q-Advanced models for tsunamis and rogue waves[END_REF]. We refer the reader to [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF][START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF] for some of the proofs of the following results, whose statements are included for the sake of completeness.

2.1. First family of Banach spaces of functions with q-exponential growth and exponential decay. Definition 2.1. Let q > 1. We denote Exp q (k,β,µ,α,ρ) the vector space of complex valued continuous functions (τ, m) → h(τ, m) on (U d ∪ D(0, ρ)) × R, holomorphic with respect to τ on U d ∪ D(0, ρ) and such that

h(τ, m) (k,β,µ,α,ρ) = sup τ ∈U d ∪D(0,ρ), m∈R (1 + |m|) µ e β|m| exp - k 2 log 2 |τ + δ| log(q) -α log |τ + δ| |h(τ, m)| is finite. The space (Exp q (k,β,µ,α,ρ) , • (k,β,µ,α,ρ)
) is a Banach space. The proof of the following lemma is straightforward.

Lemma 2.2. Let (τ, m) → a(τ, m) be a bounded continuous function on (U d ∪ D(0, ρ)) × R, holomorphic with respect to τ on U d ∪ D(0, ρ). Then, it holds that a(τ, m)h(τ, m) (k,β,µ,α,ρ) ≤ sup τ ∈U d ∪D(0,ρ),m∈R |a(τ, m)| h(τ, m) (k,β,µ,α,ρ) , for every h(τ, m) ∈ Exp q (k,β,µ,α,ρ) . Proposition 2.3. Let γ 1 , γ 2 , γ 3 ≥ 0 such that γ 1 + kγ 3 ≥ γ 2 .
Let a γ 1 (τ ) be a continuous function on U d ∪ D(0, ρ), holomorphic on U d ∪ D(0, ρ), with

|a γ 1 (τ )| ≤ 1 (1 + |τ |) γ 1 ,
for every τ ∈ (U d ∪ D(0, ρ)). Then, there exists C 1 > 0, depending on k, q, α, γ 1 , γ 2 , γ 3 , such that

a γ 1 (τ )τ γ 2 σ -γ 3 q,τ f (τ, m) (k,β,µ,α,ρ) ≤ C 1 f (τ, m) (k,β,µ,α,ρ)
for every f ∈ Exp q (k,β,µ,α,ρ) . Definition 2.4. We write E (β,µ) for the vector space of continuous functions h :

R → C such that h(m) (β,µ) = sup m∈R (1 + |m|) µ exp(β|m|)|h(m)| < ∞. It holds that (E (β,µ) , • (β,µ) ) is a Banach space.
The Banach space (E (β,µ) , • (β,µ) ) can be endowed with the structure of a Banach algebra with the following noncomutative product (see Proposition 2 in [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF] for further details).

Proposition 2.5. Let Q(X), R(X) ∈ C[X] be polynomials such that deg(R) ≥ deg(Q), R(im) = 0, for all m ∈ R. Let m → b(m) be a continuous function in R such that |b(m)| ≤ 1/|R(im)|, m ∈ R. Assume that µ > deg(Q)+1. Then, there exists a constant C 2 > 0 (depending on Q(X), R(X), µ) such that b(m) +∞ -∞ f (m -m 1 )Q(im 1 )g(m 1 )dm 1 (β,µ) ≤ C 2 f (m) (β,µ) g(m) (β,µ) ,
for every f (m), g(m) ∈ E (β,µ) . In the sequel, we adopt the notation

f (m) Q g(m) := +∞ -∞ f (m -m 1 )Q(im 1 )g(m 1 )dm 1 ,
for every m ∈ R, extending the classical convolution product for Q ≡ 1. As a result, (E (β,µ) , • (β,µ) ) becomes a Banach algebra for the product b,Q defined by

f (m) b,Q g(m) := b(m)f (m) Q g(m).
The next proposition is a slightly modified version of Proposition 3 in [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF], adapted to the appearance of two different types of growth of the functions involved, which force holding some positive distance to the origin. Proposition 2.6. Let b(m), Q(X), R(X) be chosen as in Proposition 2.5. We assume

1 ≤ k ≤ κ is an integer. Let c h (m) ∈ E (β,µ) for h ≥ 0, such that (2.1) c h (β,µ) ≤ C 1 T h q h 2 2k (1-κ k ) , h ≥ 0,
for some C > 0 and T > q 1/(2k) . Let ϕ k (τ, m) be the power series

ϕ k (τ, m) = h≥0 c h (m) τ h (q 1/k ) h(h-1)/2 ∈ E (β,µ) [[τ ]],
which defines an entire function with respect to τ with values in E (β,µ) , in view of (??).

For every f (τ, m) ∈ Exp q (κ,β,µ,α,ρ) , we define a q-analog of the convolution of order k of ϕ k (τ, m) and f (τ, m) as

ϕ k (τ, m) Q q;1/k f (τ, m) := h≥0 τ h (q 1/k ) h(h-1)/2 c h (m) Q (σ -h k q,τ f )(τ, m).
Then, the function b(m)ϕ k (τ, m) Q q;1/k f (τ, m) belongs to Exp q (κ,β,µ,α,ρ) and there exists C 3 > 0, depending on µ, q, α, k, κ, Q(X), R(X), δ, T , such that

b(m)ϕ k (τ, m) Q q;1/k f (τ, m) (κ,β,µ,α,ρ) ≤ C 3 C f (τ, m) (κ,β,µ,α,ρ) .
Proof. Let f (τ, m) ∈ Exp q (κ,β,µ,α,ρ) . From the very definition of the norm

• (κ,β,µ,α,ρ) , we know that b(m)ϕ k (τ, m) Q q;1/k f (τ, m) (κ,β,µ,α,ρ) = sup τ ∈U d ∪D(0,ρ),m∈R (1 + |m|) µ e β|m| exp - κ 2 log 2 (|τ + δ|) log(q) -α log(|τ + δ|) |b(m)| × h≥0 τ h (q 1/k ) h(h-1) 2 +∞ -∞ c h (m -m 1 )Q(im 1 )f ( τ q h/k , m 1 )dm 1 = sup τ ∈U d ∪D(0,ρ),m∈R L(τ, m).
We first give upper estimates for sup τ ∈D(0,ρ),m∈R L(τ, m).

By construction, there exist two constants Q, R > 0 such that

(2.2) |Q(im 1 )| ≤ Q(1 + |m 1 |) deg(Q) , |R(im)| ≥ R(1 + |m|) deg(R) ,
for all m ∈ R. Using (2.3) and from Lemma 4 in [START_REF] Malek | On Gevrey asymptotics for some nonlinear integro-differential equations[END_REF] (see also Lemma 2.2 from [START_REF] Costin | Short time existence and Borel summability in the Navier-Stokes equation in R 3[END_REF]), we get a constant C2 > 0 with

(2.3) (1 + |m|) µ |b(m)| +∞ -∞ |Q(im 1 )| (1 + |m -m 1 |) µ (1 + |m 1 |) µ dm 1 ≤ sup m∈R Q R (1 + |m|) µ-deg(R) × +∞ -∞ 1 (1 + |m -m 1 |) µ (1 + |m 1 |) µ-deg(Q) dm 1 ≤ C2 provided that µ > deg(Q) + 1.
From the definition of f (κ,β,µ,α,ρ) and c h (β,µ) for all h ≥ 0, (2.3), and (2.4), we get that for every τ ∈ D(0, ρ) and m ∈ R, L(τ, m) is upper bounded by

(1 + |m|) µ e β|m| exp - κ 2 log 2 (|τ + δ|) log(q) -α log(|τ + δ|) |b(m)| × h≥0 |τ | h (q 1/k ) h(h-1)/2 +∞ -∞ c h (β,µ) 1 (1 + |m -m 1 |) µ e -β|m-m 1 | |Q(im 1 )| 1 (1 + |m 1 |) µ e -β|m 1 | exp κ 2 log 2 (|τ /q h/k + δ|) log(q) + α log(|τ /q h/k + δ|) dm 1 f (τ, m) (κ,β,µ,α,ρ) ≤ C2 exp - κ 2 log 2 (|τ + δ|) log(q) -α log(|τ + δ|) h≥0 |τ | h (q 1/k ) h(h-1)/2 c h (β,µ) × exp κ 2 log 2 (|τ /q h/k + δ|) log(q) + α log(|τ /q h/k + δ|) f (τ, m) (κ,β,µ,α,ρ) ≤ Ĉ2 h≥0 ρ h (q 1/k ) h(h-1)/2 c h (β,µ) f (τ, m) (κ,β,µ,α,ρ) , with Ĉ2 = C2 exp κ 2 log 2 (ρ+δ) log(q) + α log(ρ + δ) .
The assumption (??) on c h (β,µ) allows to conclude the result when restricting the domain on the variable τ to the subset D(0, ρ).

Let U d be the complementary of D(0, ρ) in U d . From what precedes we may take the supremum over U d instead of U d ∪ D(0, ρ).

By inserting terms that correspond to the ||.|| (β,µ) norm of c h (m) and to the ||.|| (κ,β,µ,α,ρ) norm of f (τ /q h/k , m), we can give the bound estimates

sup τ ∈ Ũd ,m∈R L(τ, m) ≤ sup τ ∈ U d ,m∈R (1 + |m|) µ e β|m| exp - κ 2 log 2 (|τ + δ|) log(q) -α log(|τ + δ|) |b(m)| × h≥0 +∞ -∞ (1 + |m -m 1 |) µ e β|m-m 1 | |c h (m -m 1 )| (q 1/k ) h(h-1)/2 |τ | h × |f ( τ q h/k , m 1 )|(1 + |m 1 |) µ e β|m 1 | exp - κ 2 log 2 (|τ /q h/k + δ|) log(q) -α log(|τ /q h/k + δ|) × e -β|m-m 1 | (1 + |m -m 1 |) µ |Q(im 1 )| (1 + |m 1 |) µ e -β|m 1 | exp κ 2 log 2 (|τ /q h/k + δ|) log(q) + α log(|τ /q h/k + δ|) dm 1 .
By means of the triangular inequality |m| ≤ |m -

m 1 | + |m 1 |, we deduce that (2.4) sup τ ∈ Ũd ,m∈R L(τ, m) ≤ Ĉ||f (τ, m)|| (κ,β,µ,α,ρ)
where

(2.5) Ĉ = sup τ ∈ U d ,m∈R (1 + |m|) µ exp - κ 2 log 2 (|τ + δ|) log(q) -α log(|τ + δ|) |b(m)| × h≥0 ||c h || (β,µ) |τ | h (q 1/k ) h(h-1)/2 +∞ -∞ |Q(im 1 )| (1 + |m -m 1 |) µ (1 + |m 1 |) µ dm 1 × exp κ 2 log 2 (|τ /q h/k + δ|) log(q) + α log(|τ /q h/k + δ|) .
Again, we can also apply (2.4) at this point. On the other hand, we can provide upper estimates on the following expression

(2.6) exp κ 2 log(q) log 2 τ q h/k + δ -log 2 τ q h/k × exp κ 2 log(q) log 2 τ q h/k -log 2 (|τ |) × exp κ 2 log(q) log 2 (|τ |) -log 2 (|τ + δ|) × exp α log τ q h/k -α log(|τ |) × exp α log τ q h/k + δ -α log | τ q h/k | × exp (α log(|τ |) -α log(|τ + δ|))
as follows. The proof of Proposition 3 [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF] can be applied to the second and fourth lines of (??) which yield

(2.7) exp κ 2 log(q) log 2 τ q h/k -log 2 (|τ |) = q h 2 κ 2k 2 |τ | -hκ k , and
(2.8) exp α log τ q h/k -α log(|τ |) = q -αh k ,
respectively. It is straightforward to check that the expression in the fifth line of (??) is upper bounded by (2.9)

C 31 (q h/k ) α ,
for some positive constant C 31 . We give upper bounds for the first line in (??). In the case that |τ /q h/k | ≤ 1, this expression is upper bounded by a constant which does not depend on τ nor h. Otherwise, we have

(2.10) exp κ 2 log(q) log 2 τ q h/k + δ -log 2 τ q h/k ≤ exp κ 2 log(q) log 2 τ q h/k + δ -log 2 τ q h/k ≤ sup x>1 exp( κ 2 log(q) (log 2 (x + δ)) -log 2 (x)) ≤ C 32 ,
for some C 32 > 0. We finally provide upper bounds on the third line of (??). Taking into account that

(2.11) log 2 |τ | -log 2 |τ + δ| = -log 2 |1 + δ τ | -2 log |τ | log |1 + δ τ | ≤ C 33 , τ ∈ U d
for some C 33 > 0. From (2.2), (2.4), (??), (??), (??), (??), (??), and (??) we derive the existence of C31 , C32 > 0 such that

Ĉ ≤ C31 h≥0 ||c h || (β,µ) sup τ ∈ U d |τ | h(1-κ k ) q h 2k + h 2 2k ( κ k -1) ≤ C32 C h≥0 q 1/(2k) T h ≤ C 3 C,
which yields the result, when choosing T > q 1/(2k) .

Remark: Observe that condition (??) on the coefficients c h is always satisfied in the case that only a finite number of c h is not identically zero, i.e.

ϕ k ∈ E (β,µ) [τ ].
2.2. Second family of Banach spaces of functions with q-exponential growth and exponential decay. The second family of auxiliary Banach spaces has already been studied in previous works, such as [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF][START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF]. We refer to these references for the proofs of the related results.

Let S d be an infinite sector of bisecting direction d and let ν ∈ R.

Definition 2.7. We write Exp q (k,β,µ,ν) for the vector space of continuous functions (τ, m) → h(τ, m) on S d × R, and holomorphic with respect to τ on S d , such that

h(τ, m) (k,β,µ,ν) = sup τ ∈S d ,m∈R (1 + |m|) µ e β|m| exp - k log 2 |τ | 2 log(q) -ν log |τ | |h(τ, m)| is finite. It holds that (Exp q (k,β,µ,ν) , • (k,β,µ,ν)
) is a Banach space.

Remark 2.8. Let 0 ≤ κ 1 ≤ κ 2 . For every f ∈ Exp q (κ 1 ,β,µ,ν) , it holds that f ∈ Exp q (κ 2 ,β,µ,ν) , and

f (τ, m) (κ 2 ,β,µ,ν) ≤ f (τ, m) (κ 1 ,β,µ,ν) .
The proof of the following lemma is a straightforward consequence of the definition.

Lemma 2.9. Let a(τ, m) be a bounded continuous function on S d × R, holomorphic on S d with respect to τ . Then,

a(τ, m)f (τ, m) (k,β,µ,ν) ≤ sup τ ∈S d ,m∈R |a(τ, m)| f (τ, m) (k,β,µ,ν) for every f (τ, m) ∈ Exp q (k,β,µ,ν) . Proposition 2.10. Let γ 1 , γ 2 ≥ 0 and γ 3 ∈ R such that (2.12) γ 1 + kγ 3 ≥ γ 2 .
Let a γ 1 (τ ) be holomorphic on S d , with

|a γ 1 (τ )| ≤ 1 (1 + |τ |) γ 1 , τ ∈ S d .
Then, there exists C 4 > 0, depending on k, q, ν, γ 1 , γ 2 , γ 3 such that

a γ 1 (τ )τ γ 2 σ -γ 3 q,τ f (τ, m) (k,β,µ,ν) ≤ C 4 f (τ, m) (k,β,µ,ν)
for every f ∈ Exp q (k,β,µ,ν) .

Proof. For every f ∈ Exp q (k,β,µ,ν) we have

a γ 1 (τ )τ γ 2 σ -γ 3 q,τ f (τ, m) (k,β,µ,ν) ≤ sup τ ∈S d ,m∈R (1 + |m|) µ e β|m| exp - k log 2 |τ | 2 log(q) -ν log |τ | |τ | γ 2 (1 + |τ |) γ 1 |f (τ /q γ 3 , m)| × exp - k log 2 |τ /q γ 3 | 2 log(q) -ν log |τ /q γ 3 | exp k log 2 |τ /q γ 3 | 2 log(q) + ν log |τ /q γ 3 | ≤ sup τ ∈S d q kγ 2 3 -2νγ 3 2 |τ | -kγ 3 +γ 2 (1 + |τ |) γ 1 f (τ, m) (k,β,µ,ν) .
The result follows from the condition (2.7).

Following the same lines of arguments as in Proposition 2.6, we deduce the next proposition

Proposition 2.11. Let b(m), Q(X), R(X), c h for h ≥ 0 and ϕ k (τ, m) be chosen as in Propo- sition 2.6. For every f (τ, m) ∈ Exp q (k,β,µ,ν) , it holds that b(m)ϕ k (τ, m) Q q;1/k f (τ, m) belongs to Exp q (k,β,µ,ν) and there exists C 4 > 0, depending on µ, q, ν, k, Q(X), R(X), such that b(m)ϕ k (τ, m) Q q;1/k f (τ, m) (k,β,µ,ν) ≤ CC 4 f (τ, m) (k,β,µ,ν) .

Formal and analytic operators involved in the study of the problem

The main properties of some formal and analytic transformations are displayed for the sake of completeness. In this section, E stands for a complex Banach space.

The definition and main properties of the q-analog of Borel and Laplace transformation in several different orders can be found in [START_REF] Dreyfus | Building meromorphic solutions of q-difference equations using a Borel-Laplace summation[END_REF][START_REF] Pravica | q-Advanced models for tsunamis and rogue waves[END_REF]. The proofs of the following results can be found in [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF].

Let q > 1 be a real number, and k ≥ 1 be a rational number.

Definition 3.1. For every â(T ) = n≥0 a n T n ∈ E[[T ]] we define the formal q-Borel transform of order k of â(T ) by Bq;1/k (â(T ))(τ ) = n≥0 a n τ n (q 1/k ) n(n-1)/2 ∈ E[[τ ]]. Proposition 3.2. Let σ ∈ N and j ∈ Q. Then, it holds Bq;1/k (T σ σ j q â(T ))(τ ) = τ σ (q 1/k ) σ(σ-1)/2 σ j-σ k q Bq;1/k (â(T ))(τ ) , for every â(T ) ∈ E[[T ]].
The q-analog of Laplace transformation as it is shown was developed in [START_REF] Vizio | On q-summation and confluence[END_REF]. The associated kernel of such transformation is the Jacobi theta function of order k defined by

Θ q 1/k (x) = n∈Z q -n(n-1) 2k
x n , which turns out to be a holomorphic function in C . It turns out to be a solution of the q-difference equation

(3.1) Θ q 1/k (q m k x) = q m(m+1) 2k x m Θ q 1/k (x),
for every m ∈ Z, valid for all x ∈ C . As a matter of fact, the inverse of the Jacobi theta function of order k is a function of q-Gevrey decrease of order k at infinity in the sense that for every δ > 0 there exists C q,k > 0, not depending on δ, such that

(3.2) Θ q 1/k (x) ≥ C q,k δ exp k 2 log 2 |x| log(q) |x| 1/2 , for x ∈ C under the condition that |1 + xq m k | > δ, for every m ∈ Z.
Definition 3.3. Let ρ > 0 and U d be an unbounded sector with vertex at 0 and bisecting direction d ∈ R. Let f : D(0, ρ) ∪ U d → E be a holomorphic function, continuous on D(0, ρ), such that there exist K > 0 and α ∈ R with

f (x) E ≤ K exp k log 2 |x| 2 log(q) + α log |x| , x ∈ U d , |x| ≥ ρ, and
f (x) E ≤ K, x ∈ D(0, ρ). Set π q 1/k = log(q) k n≥0 (1 - 1 q (n+1)/k ) -1 .
We define the q-Laplace transform of order k of f along direction d by

L d q;1/k (f (x))(T ) = 1 π q 1/k L d f (u) Θ q 1/k u T du u ,
where

L d := {te id : t ∈ (0, ∞)}.
We refer the reader to Lemma 4 and Proposition 6 in [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF] for the proof of the next results. The algebraic property held by q-Laplace transformation would allow to commute some operators with respect to it. Lemma 3.4. Let δ > 0. Under the hypotheses of Definition 3.3, we have that L d q;1/k (f (x))(T ) defines a bounded and holomorphic function on R d, δ ∩ D(0, r 1 ) for every 0 < r 1 ≤ q ( 1 2 -α)/k /2, where

(3.3) R d, δ := T ∈ C : 1 + re id T > δ, for all r ≥ 0 .
A different choice for d modulo 2πZ would provide the same function due to Cauchy formula.

Proposition 3.5. Let f be a function which satisfies the properties in Definition 3.3, and let δ > 0. Then, for every σ ≥ 0 one has

T σ σ j q (L d q;1/k f (x))(T ) = L d q;1/k x σ (q 1/k ) σ(σ-1)/2 σ j-σ k q f (x) (T ), for every T ∈ R d, δ ∩ D(0, r 1 ), with 0 < r 1 ≤ q ( 1 2 -α)/k /2.
Another operator which is used through the work is the inverse Fourier transform.

Proposition 3.6. Let f ∈ E (β,µ) with β > 0, µ > 1.
The inverse Fourier transform of f is defined by

F -1 (f )(x) = 1 (2π) 1/2 ∞ -∞ f (m) exp(ixm)dm
for all x ∈ R. The function F -1 (f ) extends to an analytic function on the strip

H β = {z ∈ C/|Im(z)| < β}. Let φ(m) = imf (m) ∈ E (β,µ-1)
. Then, we have

∂ z F -1 (f )(z) = F -1 (φ)(z) for all z ∈ H β . Let g ∈ E (β,µ) and let ψ(m) = 1 (2π) 1/2 f g(m)
, the convolution product of f and g, for all m ∈ R.

From Proposition 2.5, we know that ψ ∈ E (β,µ) . Moreover, we have

F -1 (f )(z)F -1 (g)(z) = F -1 (ψ)(z)
for all z ∈ H β .

4.

Formal and analytic solutions to some auxiliary convolution initial value problems with complex parameters

Let 1 ≤ k 1 < k 2 and D, D 1 , D 2 ≥ 3 be integers and define κ -1 = k -1 1 -k -1 2 .
Observe that κ > k 1 . Let q > 1 be a real number. We also consider the positive integer numbers

d D 1 , d D 2 .
For every 1 ≤ ≤ D -1 we consider nonnegative integers d , δ ≥ 1 and ∆ ≥ 0. We assume that (4.1)

δ 1 = 1, δ < δ +1 , for 1 ≤ ≤ D -2.
We also assume that

(4.2) ∆ ≥ d , d D 1 -1 κ + d k 2 + 1 ≥ δ , d k 1 + 1 ≥ δ , d D 2 -1 k 2 ≥ δ -1,
for every 1 ≤ ≤ D -1, and also

(4.3) k 1 (d D 2 -1) > k 2 d D 1 . Let Q(X), R (X) ∈ C[X] for 1 ≤ ≤ D -1 and R D 1 , R D 2 ∈ C[X], such that (4.4) deg(R D 2 ) = deg(R D 1 ), and
(4.5) deg(Q) ≥ deg(R D j ) ≥ deg(R ), µ -1 > deg(R D j ), Q(im) = 0, R D j (im) = 0, for some µ > deg(R D j ) + 1 with j = 1, 2, for all m ∈ R, 1 ≤ ≤ D -1.
We consider sequences of functions m → F n (m, ) and m → C ,n (m, ) for n ≥ 0 belonging to the Banach space E (β,µ) for some β > 0, depending holomorphically on ∈ D(0, 0 ), for some 0 > 0. We moreover assume there exist C , C F , T 0 > 0 such that

(4.6) C ,n (β,µ) ≤ C 1 T 0 n q -n 2 κ 2k 1 k 2 , F n (β,µ) ≤ C F 1 T 0 n , n ≥ 0.
We define the formal power series in

E (β,µ) [[T ]] Ĉ (T, m, ) = n≥0 C ,n (m, )T n , F (T, m, ) = n≥0 F n (m, )T n .
We consider the following initial value problem

(4.7) Q(im)σ q,T U (T, m, ) = T d D 1 σ d D 1 k 1 +1 q,T R D 1 (im)U (T, m, ) + T d D 2 σ d D 2 k 2 +1 q,T R D 2 (im)U (T, m, ) + D-1 =1 ∆ -d T d σ δ q,T 1 (2π) 1/2 +∞ -∞ Ĉ (T, m -m 1 , )R (im 1 )U (T, m 1 , )dm 1 + σ q,T F (T, m, ).
Proposition 4.1. There exists a unique formal power series

(4.8) Û (T, m, ) = n≥0 U n (m, )T n ,
solution of (4.7), where the coefficients U n (m, ) belong to E (β,µ) , for β > 0 and µ > deg(R D j )+1, j ∈ {1, 2}, given above and depend holomorphically on ∈ D(0, 0 ).

Proof. We plug the formal power series (4.8) into equation (4.7) to obtain a recursion formula for the coefficients U n , for n ≥ 0. We have

Q(im)U n (m, )q n = R D 1 (im)U n-d D 1 (m, )q d D 1 k 1 +1 (n-d D 1 ) + R D 2 (im)U n-d D 2 (m, )q d D 2 k 2 +1 (n-d D 2 ) + D-1 =1 ∆ -d q (n-d )δ   n 1 +n 2 =n-d 1 (2π) 1/2 +∞ -∞ C ,n 1 (m -m 1 , )R (im 1 )U n 2 (m 1 , )dm 1   + F n (m, )q n for every n ≥ max{d D 1 , d D 2 , max 1≤ ≤D-1 d }. Due to C ,n , F n ∈ E (β,µ)
for every n ≥ 0 and 1 ≤ ≤ D -1, we get U n ∈ E (β,µ) by recursion. We observe that Proposition 2.5. is applied in the recursion. 4.1. Analytic solutions of a first auxiliary problem in the q-Borel plane. We proceed to multiply at both sides of equation (4.7) by T k 1 and then apply the formal q-Borel transformation of order k 1 with respect to T . Let ϕ k 1 , (τ, m, ) be the formal q-Borel transform of order k 1 of Ĉ (T, m, ) with respect to T , and Ψ k 1 (τ, m, ) the formal q-Borel transform of order k 1 of F (T, m, ) with respect to T . More precisely, we have

ϕ k 1 , (τ, m, ) = n≥0 C ,n (m, ) τ n (q 1/k 1 ) n(n-1)/2 , Ψ k 1 (τ, m, ) = n≥0 F n (m, ) τ n (q 1/k 1 ) n(n-1)/2 . (4.9)
According to Lemma 5 of [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF], the expression Ψ k 1 (τ, m, ) represents an entire function of q-exponential growth of order k 1 that belongs to the Banach space Exp q (κ,β,µ,α,ρ) since κ > k 1 , provided that α satisfies T 0 > q 1 2k 1 /q α/k 1 , for any unbounded sector U d and any disc D(0, ρ). More precisely, we have

(4.10) ||Ψ k 1 (τ, m, )|| (κ,β,µ,α,ρ) ≤ C Ψ k 1 for some constant C Ψ k 1 > 0, for all ∈ D(0, 0 ).
In view of the properties of the q-Borel transformation of order k 1 , we arrive at the equation

(4.11) Q(im) τ k 1 (q 1/k 1 ) k 1 (k 1 -1) 2 w k 1 (τ, m, ) = R D 1 (im) τ d D 1 +k 1 (q 1/k 1 ) (d D 1 +k 1 )(d D 1 +k 1 -1) 2 w k 1 (τ, m, ) + R D 2 (im) τ d D 2 +k 1 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 σ d D 2 1 k 2 -1 k 1 q,τ w k 1 (τ, m, ) + τ k 1 (q 1/k 1 ) k 1 (k 1 -1) 2 Ψ k 1 (τ, m, ) + D-1 =1 ∆ -d τ d +k 1 (q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 σ δ - d k 1 -1 q,τ 1 (2π) 1/2 ϕ k 1 , (τ, m, ) R q;1/k 1 w k 1 (τ, m, )
where w k 1 (τ, m, ) stands for the formal q-Borel transformation of order k 1 of U (T, m, ) with respect to T . Observe the appearance only of negative powers of the dilation operator in each one of the terms in the sum of the right-hand side of the equation.

We assume an unbounded sector of bisecting direction

d Q,R D 1 ∈ R exists, S Q,R D 1 = z ∈ C : |z| ≥ r Q,R D 1 , | arg(z) -d Q,R D 1 | ≤ ν Q,R D 1 , for some r Q,R D 1 , ν Q,R D 1 > 0, in such a way that Q(im) R D 1 (im) ∈ S Q,R D 1 ,
for every m ∈ R. We factorize

P m,1 (τ ) = Q(im) (q 1/k 1 ) k 1 (k 1 -1) 2 - R D 1 (im) (q 1/k 1 ) (d D 1 +k 1 )(d D 1 +k 1 -1) 2 τ d D 1
in the form

P m,1 (τ ) = - R D 1 (im) (q 1/k 1 ) (d D 1 +k 1 )(d D 1 +k 1 -1) 2 d D 1 -1 =0 (τ -q (m)), with q (m) = e 2iπ d D 1 Q(im) R D 1 (im) 1/d D 1 q d D 1 +2k 1 -1 2k 1 , for every 0 ≤ ≤ d D 1 -1.
Let U d be an unbounded sector, and ρ > 0 such that the following statements hold:

1) There exists

M 1 > 0 such that |τ -q (m)| ≥ M 1 (1 + |τ |) for every 0 ≤ ≤ d D 1 -1, m ∈ R, and τ ∈ U d ∪ D(0, ρ). An appropriate choice of r Q,R D 1 and ρ yields |q (m)| > 2ρ for every m ∈ R, and 0 ≤ ≤ d D 1 -1. In the case that ν Q,R D 1 is small enough, the set {q (m) : m ∈ R, 0 ≤ ≤ d D 1
-1} stays at a positive distance to U d , and it can be chosen with the property that q (m)/τ has positive distance to 1

∈ C for every τ ∈ U d , m ∈ R and 0 ≤ ≤ d D 1 -1. 2) There exists M 2 > 0 such that |τ -q (m)| ≥ M 2 |q (m)| for every ∈ {0, . . . , d D 1 -1}, m ∈ R and τ ∈ U d ∪ D(0, ρ)
. This is a direct consequence of 1), for some small enough M 2 > 0.

In order to prove the following upper estimates in (4.13), we make use of 1) for all ∈ {0, . . . , d D 1 -1}, except for one of them, say 0 , for which 2) is applied. The previous conditions yield the existence of C P > 0 such that

(4.12) |P m,1 (τ )| ≥ M d D 1 -1 1 M 2 |R D 1 (im)|(1 + |τ |) d D 1 -1 (q 1/k 1 ) (d D 1 +k 1 )(d D 1 +k 1 -1) 2 |Q(im)| |R D 1 (im)| 1/d D 1 q d D 1 +2k 1 -1 2k 1 ≥ C P (r Q,R D 1 ) 1/d D 1 |R D 1 (im)|(1 + |τ |) d D 1 -1 ,
for every τ ∈ U d ∪ D(0, ρ), and m ∈ R.

The next result states the existence and uniqueness of a solution of (4.12) in the space Exp q (κ,β,µ,α,ρ) , provided its norm in that space is small enough.

Proposition 4.2. Under the Assumptions (4.1), (4.2), (4.3), (4.4) and (4.5), there exist

r Q,R D 1 > 0, a constant > 0 and constants ς ϕ , ς Ψ > 0 such that if (4.13) C ≤ ς ϕ C Ψ k 1 ≤ ς Ψ ,
for all 1 ≤ ≤ D -1 (see (4.6) and (??)), then the equation (4.12) admits a unique solution

w d k 1 (τ, m, ) ∈ Exp q (κ,β,µ,α,ρ) with w d k 1 (τ, m, ) (κ,β,µ,α,ρ) ≤ , for every ∈ D(0, 0 ).
Proof. Let ∈ D(0, 0 ) and consider the operator H defined by

H 1 (w(τ, m)) := R D 2 (im) P m,1 (τ ) τ d D 2 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 σ d D 2 1 k 2 -1 k 1 q,τ w(τ, m) + D-1 =1 ∆ -d τ d P m,1 (τ )(q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 σ δ - d k 1 -1 q,τ 1 (2π) 1/2 ϕ k 1 , (τ, m, ) R q;1/k 1 w(τ, m) + 1 P m,1 (τ )(q 1/k 1 ) k 1 (k 1 -1) 2 Ψ k 1 (τ, m, ).
Note that a fixed point of H 1 (w(τ, m)) will lead to a convenient solution of (4.12). To apply the fixed point theorem, we are going to prove successively two facts.

(1) One may choose small enough ς ϕ , ς Ψ , > 0, and large enough r

Q,R D 1 > 0 such that (4.14) H 1 (B(0, )) ⊆ B(0, ),
where B(0, ) stands for the closed disc centered at 0, with radius in the Banach space Exp q (κ,β,µ,α,ρ) . (2) It holds

(4.15) H 1 (w 1 (τ, m)) -H 1 (w 2 (τ, m)) (κ,β,µ,α,ρ) ≤ 1 2 w 1 (τ, m) -w 2 (τ, m) (κ,β,µ,α,ρ) ,
for every w 1 (τ, m), w 2 (τ, m) ∈ B(0, ).

Proof of (4.15).

We first check (4.15). Let w(τ, m) ∈ Exp q (κ,β,µ,α,ρ) . With (4.2) and the definition of κ, we find that

d D 1 -1 + κ(d /k 1 + 1 -δ ) ≥ d and d /k 1 + 1 -δ ≥ 0.
τ d P m,1 (τ )(q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 σ δ - d k 1 -1 q,τ 1 (2π) 1/2 ϕ k 1 , (τ, m, ) R q;1/k 1 w(τ, m) (κ,β,µ,α,ρ) ≤ ∆ -d 0 C 1 C 3 ς ϕ (q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 (2π) 1/2 w(τ, m) (κ,β,µ,α,ρ) .
Gathering Lemma 2.2, we get

1 P m,1 (τ )(q 1/k 1 ) k 1 (k 1 -1)/2 Ψ k 1 (τ, m, ) (κ,β,µ,α,ρ) ≤ 1 (q 1/k 1 ) k 1 (k 1 -1)/2 C P (r Q,R D 1 ) 1/d D 1 sup m∈R 1 |R D 1 (im)| ς Ψ . (4.17)
Condition (4.4) and the application of Proposition 2.3 and Lemma 2.2 yields

R D 2 (im) P m,1 (τ ) τ d D 2 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 σ d D 2 1 k 2 -1 k 1 q,τ w(τ, m) (κ,β,µ,α,ρ) ≤ sup m∈R |R D 2 (im)| |R D 1 (im)| C 1 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 . (4.18) An appropriate choice of r Q,R D 1 > 0, , ς Ψ , ς ϕ > 0 gives D-1 =1 ∆ -d 0 C 3 ς ϕ C 1 (q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 (2π) 1/2 + 1 (q 1/k 1 ) k 1 (k 1 -1)/2 C P (r Q,R D 1 ) 1/d D 1 sup m∈R 1 |R D 1 (im)| ς Ψ + sup m∈R |R D 2 (im)| |R D 1 (im)| C 1 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 ≤ . (4.19)
Regarding (4.17), (4.18), (4. [START_REF] Pravica | Solutions of a class of multiplicatively advanced differential equations[END_REF]) and (4.20) one concludes (4.15).

Proof of (4.16).

We proceed to prove (4.16). Let w 1 , w 2 ∈ Exp q (κ,β,µ,α,ρ) . We assume w (τ, m) (κ,β,µ,α,ρ) ≤ , = 1, 2, for some > 0. Let E(τ, m) := w 1 (τ, m) -w 2 (τ, m). On one hand, from (4.17) one has

∆ -d τ d P m,1 (τ )(q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 σ δ - d k 1 -1 q,τ 1 (2π) 1/2 ϕ k 1 , (τ, m, ) R q;1/k 1 E(τ, m) (κ,β,µ,α,ρ) ≤ ∆ -d 0 C 3 ς ϕ C 1 (q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 (2π) 1/2 E(τ, m) (κ,β,µ,α,ρ) . (4.20)
On the other hand, (4.19) yields

R D 2 (im) P m,1 (τ ) τ d D 2 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 σ d D 2 1 k 2 -1 k 1 q,τ E(τ, m) (κ,β,µ,α,ρ) ≤ sup m∈R |R D 2 (im)| |R D 1 (im)| C 1 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 E(τ, m) (κ,β,µ,α,ρ) . (4.21) We choose r Q,R D 1 > 0, ς ϕ > 0 such that D-1 =1 ∆ -d 0 C 3 ς ϕ C 1 (q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 (2π) 1/2 + sup m∈R |R D 2 (im)| |R D 1 (im)| C 1 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 C P (r Q,R D 1 ) 1/d D 1 ≤ 1 2 . (4.22)
The statement (4.16) is a direct consequence of condition (4.23) applied to (4.21) and (4.22).

Let us finish the proof of the proposition. At this point, in view of (4.15) and (4.16), one can choose > 0 such that B(0, ) ⊆ Exp q (κ,β,µ,α,ρ) , which defines a complete metric space for the norm • (κ,β,µ,α,ρ) . The map H 1 is contractive from B(0, ) into itself. The fixed point theorem states that H 1 admits a unique fixed point w d k 1 (τ, m, ) ∈ B(0, ) ⊆ Exp q (κ,β,µ,α,ρ) , for every ∈ D(0, 0 ). The construction of w d k 1 (τ, m, ) allows us to conclude that it turns out to be a solution of (4.12). The next step consists of studying the solutions of a second auxiliary problem. This problem lies in a second q-Borel plane and its solution would guarantee the extension, with appropriate growth, of the acceleration of the solution to our first auxiliary problem, described in (4.12).

We set

(4.23) Ψ k 2 (τ, m, ) = n≥0 F n (m, ) τ n (q 1/k 2 ) n(n-1)/2
the q-Borel transform of order k 2 of F (T, m, ). According to the second condition of (4.6), the expression Ψ k 2 (τ, m, ) stands for an entire function of q-exponential growth of order k 2 which belongs to the Banach space Exp q (k 2 ,β,µ,ν) provided that ν ∈ R satisfies T 0 > q 1 2k 2 /q ν/k 2 for any unbounded sector S d . More precisely, we have

(4.24) Ψ k 2 (τ, m, ) (k 2 ,β,µ,ν) ≤ C Ψ k 2
for some constant C Ψ k 2 > 0, for all ∈ D(0, 0 ).

4.2.

Analytic solutions of a second auxiliary problem in the q-Borel plane. We multiply both sides of equation (4.7) by T k 2 and apply formal q-Borel transformation of order k 2 with respect to T . In view of the properties of q-Borel transformation, the resulting problem is determined by

(4.25) Q(im) τ k 2 (q 1/k 2 ) k 2 (k 2 -1) 2 ŵk 2 (τ, m, ) = R D 1 (im) τ d D 1 +k 2 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 σ d D 1 1 k 1 -1 k 2 q,τ ŵk 2 (τ, m, ) + R D 2 (im) τ d D 2 +k 2 (q 1/k 2 ) (d D 2 +k 2 )(d D 2 +k 2 -1) 2 ŵk 2 (τ, m, ) + τ k 2 (q 1/k 2 ) k 2 (k 2 -1) 2 Ψ k 2 (τ, m, ) + D-1 =1 ∆ -d τ d +k 2 (q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 σ δ - d k 2 -1 q,τ 1 (2π) 1/2 ϕ k 2 , (τ, m, ) R q;1/k 2 ŵk 2 (τ, m, ) .
Here ŵk 2 (τ, m, ), Ψ k 2 (τ, m, ) and ϕ k 2 ,l (τ, m, ) stand for the formal q-Borel transform of order k 2 of U (T, m, ), F (T, m, ) and Ĉl (T, m, ).

We consider our second auxiliary problem, namely

(4.26) Q(im) τ k 2 (q 1/k 2 ) k 2 (k 2 -1) 2 w k 2 (τ, m, ) = R D 1 (im) τ d D 1 +k 2 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 σ d D 1 1 k 1 -1 k 2 q,τ w k 2 (τ, m, ) + R D 2 (im) τ d D 2 +k 2 (q 1/k 2 ) (d D 2 +k 2 )(d D 2 +k 2 -1) 2 w k 2 (τ, m, ) + τ k 2 (q 1/k 2 ) k 2 (k 2 -1) 2 Ψ k 2 (τ, m, ) + D-1 =1 ∆ -d τ d +k 2 (q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 σ δ - d k 2 -1 q,τ 1 (2π) 1/2 ϕ k 2 , (τ, m, ) R q;1/k 2 w k 2 (τ, m, ) .
We assume an unbounded sector of bisecting direction

d Q,R D 2 ∈ R exists, S Q,R D 2 = z ∈ C : |z| ≥ r Q,R D 2 , | arg(z) -d Q,R D 2 | ≤ ν Q,R D 2 ,
for some ν Q,R D 2 > 0, in such a way that

Q(im) R D 2 (im) ∈ S Q,R D 2 ,
for every m ∈ R. We factorize

P m,2 (τ ) = Q(im) (q 1/k 2 ) k 2 (k 2 -1) 2 - R D 2 (im) (q 1/k 2 ) (d D 2 +k 2 )(d D 2 +k 2 -1) 2 τ d D 2
in the form

P m,2 (τ ) = - R D 2 (im) (q 1/k 2 ) (d D 2 +k 2 )(d D 2 +k 2 -1) 2 d D 2 -1 =0 (τ -q ,2 (m)).
Let S d be an unbounded sector with small enough aperture in such a way that: 1) There exists In the following estimates, we apply 1) in the previous assumption to all indices ∈ {0, . . . , d D 2 -1}, except for one of them, say 0 , that we apply 2). This yields the existence of C P,2 > 0 such that

M 12 > 0 such that |τ -q ,2 (m)| ≥ M 12 (1 + |τ |) for every 0 ≤ ≤ d D 2 -1, m ∈ R,
(4.27) |P m,2 (τ )| ≥ C P,2 (r Q,R D 2 ) 1/d D 2 |R D 2 (im)|(1 + |τ |) d D 2 -1 ,
for every τ ∈ S d , and m ∈ R. (k 2 ,β,µ,ν) for ν ∈ R and depends holomorphically with respect to ∈ D(0, 0 ). Moreover,

w d k 2 (τ, m, ) (k 2 ,β,µ,ν) ≤ .
Proof. Let ∈ D(0, 0 ). We consider the map H 2 defined by

(4.28) H 2 (w(τ, m)) := R D 1 (im) P m,2 (τ ) τ d D 1 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 σ d D 1 1 k 1 -1 k 2 q,τ w(τ, m) + D-1 =1 ∆ -d τ d P m,2 (τ )(q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 σ δ - d k 2 -1 q,τ 1 (2π) 1/2 ϕ k 2 , (τ, m, ) R q;1/k 2 w(τ, m) + 1 P m,2 (τ )(q 1/k 2 ) k 2 (k 2 -1) 2 Ψ k 2 (τ, m, ).
Note that a fixed point of H 2 (w(τ, m)) will lead to a convenient solution of (4.32). To apply the fixed point theorem, we are going to prove successively two facts.

(1) One may choose small enough ς ϕ , ς Ψ , > 0, and large enough r Q,R D 2 > 0 such that (4.29) H 2 (B(0, )) ⊆ B(0, ), where B(0, ) stands for the closed disc centered at 0, with radius in the Banach space Exp q (k 2 ,β,µ,ν) . (2) It holds

(4.30) H 2 (w 1 (τ, m)) -H 2 (w 2 (τ, m)) (k 2 ,β,µ,ν) ≤ 1 2 w 1 (τ, m) -w 2 (τ, m) (k 2 ,β,µ,ν) ,
for every w 1 (τ, m), w 2 (τ, m) ∈ B(0, ).

Proof of (4.35).

We first check (4.35). Let w(τ, m) ∈ Exp q (k 2 ,β,µ,ν) . With (4.2), we find that d D 2 -1+k 2 (d /k 2 +1-δ ) ≥ d . Taking into account assumptions (4.1), (4.4), (4.5), regarding (4.33) together with Lemma 2.9, Proposition 2.10 and Proposition 2.11 we get

∆ -d τ d P m,2 (τ )(q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 σ δ - d k 2 -1 q,τ 1 (2π) 1/2 ϕ k 2 , (τ, m, ) R q;1/k 2 w(τ, m) (k 2 ,β,µ,ν) ≤ ∆ -d 0 C 4 C 3 ς ϕ (q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2 (2π) 1/2 w(τ, m) (k 2 ,β,µ,ν) . (4.31)
Gathering Lemma 2.9 we get

1 P m,2 (τ )(q 1/k 2 ) k 2 (k 2 -1)/2 Ψ k 2 (τ, m, ) (k 2 ,β,µ,ν) ≤ 1 (q 1/k 2 ) k 2 (k 2 -1)/2 C P,2 (r Q,R D 2 ) 1/d D 2 sup m∈R 1 |R D 2 (im)| ς Ψ 2 , (4.32)
for some ς Ψ 2 . Observe that ς Ψ 2 tends to 0 when ς Ψ does.

Condition (4.4), and the application of Proposition 2.10 and Lemma 2.9 yields

R D 1 (im) P m,2 (τ ) τ d D 1 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 σ d D 1 1 k 1 -1 k 2 q,τ w(τ, m) (k 2 ,β,µ,ν) ≤ sup m∈R |R D 1 (im)| |R D 2 (im)| C 4 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2 . (4.33) An appropriate choice of r Q,R D 2 > 0, , ς ϕ , ς Ψ > 0 gives D-1 =1 ∆ -d 0 C 3 ς ϕ C 4 (q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2 (2π) 1/2 + 1 (q 1/k 2 ) k 2 (k 2 -1)/2 C P,2 (r Q,R D 2 ) 1/d D 2 sup m∈R 1 |R D 2 (im)| ς Ψ + sup m∈R |R D 1 (im)| |R D 2 (im)| C 4 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2
≤ . Proof of (4.36).

We proceed to prove (4.36). Let w 1 , w 2 ∈ Exp q (k 2 ,β,µ,ν) . We assume w (τ, m) (k 2 ,β,µ,ν) ≤ , = 1, 2, for some > 0. Let E(τ, m) = w 1 (τ, m) -w 2 (τ, m). On one hand, from (4.37) one has

∆ -d τ d P m,2 (τ )(q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 σ δ - d k 2 -1 q,τ 1 (2π) 1/2 ϕ k 2 , (τ, m, ) R q;1/k 2 E(τ, m) (k 2 ,β,µ,ν) ≤ ∆ -d 0 C 3 ς ϕ C 4 (q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2 (2π) 1/2 E(τ, m) (k 2 ,β,µ,ν) .
On the other hand, (4.39) yields

R D 1 (im) P m,2 (τ ) τ d D 1 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 σ d D 1 1 k 1 -1 k 2 q,τ E(τ, m) (k 2 ,β,µ,ν) ≤ sup m∈R |R D 1 (im)| |R D 2 (im)| C 4 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2 E(τ, m) (k 2 ,β,µ,ν) . We choose r Q,R D 2 > 0, ς ϕ > 0 such that D-1 =1 ∆ -d 0 C 3 ς ϕ C 4 (q 1/k 2 ) (d +k 2 )(d +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2 (2π) 1/2 + sup m∈R |R D 1 (im)| |R D 2 (im)| C 4 (q 1/k 2 ) (d D 1 +k 2 )(d D 1 +k 2 -1) 2 C P,2 (r Q,R D 2 ) 1/d D 2 ≤ 1 2 .
We conclude (4.36). Let us finish the proof of the proposition. At this point, in view of (4.35) and (4.36), one can choose > 0 such that B(0, ) ⊆ Exp q (k 2 ,β,µ,ν) , which defines a complete metric space for the norm • (k 2 ,β,µ,ν) . The map H 2 is contractive from B(0, ) into itself. The fixed point theorem states that H 2 admits a unique fixed point w d k 2 (τ, m, ) ∈ B(0, ) ⊆ Exp q (k 2 ,β,µ,ν) , for every ∈ D(0, 0 ). The construction of w d k 2 (τ, m, ) allow us to conclude it turns out to be a solution of (4.32).

The existing link between the acceleration of w d k 1 and w d k 2 is now provided. Both functions coincide in the intersection of their domain of definition. This fact assures the extension of the acceleration of w d k 1 along direction d, with appropriate q-exponential growth in order to apply q-Laplace transformation of that order to recover the analytic solution of the main problem under study.

τ → L d q;1/κ (w d k 1 (τ, m, )) := L d q;1/κ (h → w d k 1 (h, m, ))(τ )
defines a bounded holomorphic function in R d, δ ∩ D(0, r 1 ), for 0 < r 1 ≤ q ( 1 2 -α)/κ /2. Moreover, it holds that (4.35)

L d q;1/κ (w d k 1 (τ, m, )) = w d k 2 (τ, m, ), (τ, m, ) ∈ S b d × R × D(0, 0 ),
where S b d is a finite sector of bisecting direction d.

Proof. We recall from Proposition 4.2 that w d k 1 ∈ Exp q (κ,β,µ,α,ρ) . This guarantees appropriate bounds on τ ∈ U d in order to apply q-Laplace transformation of order κ along direction d. This yields that for every δ > 0, the function L d q;1/κ (w d k 1 (τ, m, )) defines a bounded and holomorphic function in R d, δ ∩ D(0, r 1 ) for 0 < r 1 ≤ q ( 1 2 -α)/κ /2. In order to prove that (4.41) holds, it is sufficient to prove that L d q;1/κ (w d k 1 (τ, m, )) and w d k 2

are both solutions of some problem, with unique solution in certain Banach space, so they must coincide. For that purpose, we multiply both sides of equation (4.12) by τ -k 1 and take q-Laplace transformation of order κ along direction d.

The properties of q-Laplace transformation yield (4.36)

L d q;1/κ τ d D 1 w d k 1 (τ, m, ) = (q 1/κ ) d D 1 (d D 1 -1)/2 τ d D 1 σ d D 1 κ q,τ L d q;1/κ (w d k 1 )(τ, m, ), (4.37) L d q;1/κ τ d D 2 σ d D 2 1 k 2 -1 k 1 q,τ w d k 1 (τ, m, ) = (q 1/κ ) d D 2 (d D 2 -1)/2 τ d D 2 L d q;1/κ (w d k 1 )(τ, m, ), and
(4.38) L d q;1/κ τ d σ δ + d k 1 -1 q,τ 1 (2π) 1/2 ϕ k 1 , (τ, m, ) R q;1/k 1 w d k 1 (τ, m, ) = (q 1/κ ) d (d -1)/2 τ d σ δ - d k 2 -1 q,τ L d q;1/κ 1 (2π) 1/2 ϕ k 1 , (τ, m, ) R q;1/k 1 w d k 1 (τ, m, ) .
We claim that we have

(4.39) L d q;1/κ ϕ k 1 , (τ, m, ) R q;1/k 1 w d k 1 (τ, m, ) = ϕ k 2 , (τ, m, ) R q;1/k 2 L d q;1/κ (w d k 1 (τ, m, )).
This is a consequence of the change in the order of integration in the operators involved in (4.45). This situation is different from that of (60) in the proof of Proposition 12 in [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF]. Assume the variable of integration with respect to Laplace operator is r. After the change of variable r = r/q h/k 1 , we reduce the study to that of Ξ in the proof of Proposition 12 in [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF], with r replaced by r 1-h . This last argument guarantees the availability of the change of order in the integration operators involved in (4.45). We now give a proof of (4.45) under this consideration.

We have

L d q;1/κ ϕ k 1 , (τ, m, ) R q;1/k 1 w d k 1 (τ, m, ) = 1 π q 1/κ ∞ 0 ϕ k 1 , (re id , m, ) R q;1/k 1 w d k 1 (re id , m, ) 1 Θ q 1/κ re id τ dr r = 1 π q 1/κ ∞ 0   n≥0 (re id ) n (q 1/k 1 ) n(n-1)/2 C ,n (m, ) R (σ -n k 1 q,τ w d k 1 )(re id , m, )   1 Θ q 1/κ re id τ dr r = 1 π q 1/κ ∞ 0   n≥0 (re id ) n (q 1/k 1 ) n(n-1)/2 ∞ -∞ C ,n (m -m 1 , )R (im 1 )w d k 1 (re id q -n k 1 , m 1 , )dm 1   × 1 Θ q 1/κ re id τ dr r .
We make the change of variable r = r/q n/k 1 to get that the previous expression equals

1 π q 1/κ ∞ 0   ∞ -∞ n≥0 (re id ) n q n 2 /k 1 (q 1/k 1 ) n(n-1)/2 C ,n (m -m 1 , )R (im 1 )w d k 1 (re id , m 1 , )dm 1   × 1 Θ q 1/κ re id q n/k 1 τ dr r .
In view of (3.1), k -1

1 = κ -1 + k -1
2 , the change of order of the integrals and the dominated convergence theorem, the previous equation equals

= 1 π q 1/κ ∞ 0   ∞ -∞ n≥0 (re id ) n q n 2 /k 1 (q 1/k 1 ) n(n-1)/2 C ,n (m -m 1 , )R (im 1 )w d k 1 (re id , m 1 , )dm 1   × 1 Θ q 1/κ re id q n/k 2 τ q n(n+1) 2κ re id q n/k 2 τ n dr r = 1 π q 1/κ ∞ 0   ∞ -∞ n≥0 τ n q n(n-1)/(2κ) (q 1/k 1 ) n(n-1)/2 C ,n (m -m 1 , )R (im 1 )w d k 1 (re id , m 1 , )dm 1   × 1 Θ q 1/κ re id q n/k 2 τ dr r = ∞ -∞   n≥0 τ n (q 1/k 2 ) n(n-1)/2 C ,n (m -m 1 , )   R (im 1 )   1 π q 1/κ ∞ 0 w d k 1 (re id , m 1 , ) Θ q 1/κ re id q n/k 2 τ dr r   dm 1 = ϕ k 2 , (τ, m, ) R q;1/k 2 L d q;1/κ (w d k 1 (τ, m, )),
from where we conclude (4.45).

On the other hand, we observe by direct computation that (4.40)

L d q;1/κ (Ψ k 1 (τ, m, )) = Ψ k 2 (τ, m, ) for every (τ, m, ) ∈ (R d, δ ∩ D(0, r 1 )) × R × D(0, 0 ).
In view of (4.42), (4.43), (4.44), (4.45), and the last formula above (??), we derive that

Q(im) (q 1/k 1 ) k 1 (k 1 -1) 2 L d q;1/κ (w d k 1 )(τ, m, ) = R D 1 (im) (q 1/κ ) d D 1 (d D 1 -1)/2
(q 1/k 1 )

(d D 1 +k 1 )(d D 1 +k 1 -1) 2 τ d D 1 σ d D 1 κ q,τ L d q;1/κ (w d k 1 )(τ, m, ) + R D 2 (im) (q 1/κ ) d D 2 (d D 2 -1) 2 (q 1/k 1 ) (d D 2 +k 1 )(d D 2 +k 1 -1) 2 τ d D 2 L d q;1/κ (w d k 1 )(τ, m, ) + 1 (q 1/k 1 ) k 1 (k 1 -1) 2 Ψ k 2 (τ, m, ) + D-1 =1 ∆ -d τ d (q 1/κ ) d (d -1) 2 (q 1/k 1 ) (d +k 1 )(d +k 1 -1) 2 σ δ - d k 2 -1 q,τ 1 (2π) 1/2 ϕ k 2 , (τ, m, ) R q;1/k 2 L d q;1/κ (w d k 1 )(τ, m, ) , for every (τ, m, ) ∈ (R d, δ ∩ D(0, r 1 )) × R × D(0, 0 )
. We multiply at both sides of the previous equation by (

q 1/k 1 ) k 1 (k 1 -1)/2 /(q 1/k 2 ) k 2 (k 2 -1)/2 . The fact that (q 1/κ ) D(D-1) 2 (q 1/k 1 ) k 1 (k 1 -1) 2 (q 1/k 1 ) (D+k 1 )(D+k 1 -1) 2 (q 1/k 2 ) k 2 (k 2 -1) 2 = 1 (q 1/k 2 ) (D+k 2 )(D+k 2 -1) 2 , with D ∈ {d D 1 , d D 2 , d } entails that L d q;1/κ (w d k 1 (τ, m, )
) is a solution of (4.32) in its domain of definition.

Let S b d be a bounded sector of bisecting direction d such that S b d ⊆ (R d, δ ∩ D(0, r 1 )) ∩ S d , which is a nonempty set due to the assumptions on the construction of these sets. The functions

L d q;1/κ (w d k 1 (τ, m, )
) and w d k 2 (τ, m, ) are continuous complex functions defined on S b d ×R×D(0, 0 ) and holomorphic with respect to τ (resp. ) on S b d (resp. D(0, 0 )). Let ∈ D(0, 0 ) and put Ω = min{α, ν}. It is straight to check that both functions belong to the complex Banach space H (k 2 ,β,µ,Ω) , of all continuous functions (τ, m) → h(τ, m), defined on

S b d × R, holomorphic with respect to τ in S b d such that h(τ, m) H (k 2 ,β,µ,Ω) = sup τ ∈S b d ,m∈R (1 + |m|) µ e β|m| exp - k 2 2 log 2 |τ | log(q) -Ω log |τ | |h(τ, m)| is finite. It holds that L d q;1/κ (w d k 1 (τ, m, )), w d k 2 (τ,
D(0, ς) of H (k 2 ,β,µ,Ω) , whilst L d q;1/κ (w d k 1 (τ, m, )), w d k 2 (τ, m,
) are both solutions of the same problem, in the disc D(0, ς) of H (k 2 ,β,µ,Ω) , so they do coincide in the domain S b d × R × D(0, 0 ). Identity (4.41) follows from here.

Analytic solutions to a q-difference-differential equation

This section is devoted to determine in detail the main problem under study, and provide an analytic solution to it. It is worth mentioning that, although the techniques developed in previous sections are essentially novel, once the tools have been implemented, the procedure of construction of the solution coincides with that explained in Section 5 of [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF]. For the sake of completeness and a self contained work, we describe every step of the construction in detail, whilst we have decided to pass over the proofs which can be found in [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF].

Let 1 ≤ k 1 < k 2 . We define 1/κ = 1/k 1 -1/k 2 and take integers D, D 1 , D 2 larger than 3. Let q > 1 be a real number. We also consider positive integers d D 1 , d D 2 , and for every 1 ≤ ≤ D -1 we choose non negative integers d , δ ≥ 1 and ∆ ≥ 0. We make the following assumptions on the previous constants:

Assumption (A): δ 1 = 1 and δ < δ +1 for every 1 ≤ ≤ D -2.

Assumption (B):

We have

∆ ≥ d , d D 1 -1 κ + d k 2 + 1 ≥ δ , d k 1 + 1 ≥ δ , d D 2 -1 k 2 ≥ δ -1,
for every 1 ≤ ≤ D -1, and

k 1 (d D 2 -1) > k 2 d D 1 . Let Q, R D 1 , R D 2 , and R for 1 ≤ ≤ D -1 be polynomials with complex coefficients such that Assumption (C): deg(R D 2 ) = deg(R D 1 ), deg(Q) ≥ deg(R D 1 ) ≥ deg(R ). Moreover, we assume Q(im) = 0 and R D j (im) = 0 for all m ∈ R, 1 ≤ ≤ D -1. Let S Q,R D 1 and S Q,R D 2 be unbounded sectors of bisecting directions d Q,R D 1 ∈ R and d Q,R D 2 ∈ R respectively, with S Q,R D j = z ∈ C : |z| ≥ r Q,R D j , | arg(z) -d Q,R D j | ≤ ν Q,R D j ,
for some ν Q,R D j > 0, and such that

Q(im) R D j (im) ∈ S Q,R D j ,
for every m ∈ R.

Definition 5.1. Let ς ≥ 2 be an integer. A family (E p ) 0≤p≤ς-1 is said to be a good covering in C (in the plane) if the next hypotheses hold:

• E p is an open sector of finite radius 0 > 0, and vertex at the origin for every 0

≤ p ≤ ς -1. • E j ∩ E k = ∅ for 0 ≤ j, k ≤ ς -1 if and only if |j -k| ≤ 1 (we put E ς := E 0 ). • ∪ ς-1 p=0 E p = U \ {0}
for some neighborhood of the origin U.

Definition 5.2. Let (E p ) 0≤p≤ς-1 be a good covering. Let T be an open bounded sector with vertex at the origin and radius r T > 0. Given α ∈ R and ν ∈ R we assume that

0 < 0 , r T < 1, ν + k 2 log(q) log(r T ) < 0, α + κ log(q) log( 0 r T ) < 0, 0 r T ≤ q ( 1 2 -ν)/k 2 /2.
We consider a family of unbounded sectors U dp , 0 ≤ p ≤ ς -1, with bisecting direction d p ∈ R, and a family of open domains R b dp := R dp, δ ∩ D(0, 0 r T ), with

R dp, δ := T ∈ C : 1 + re idp T > δ, for every r ≥ 0 ,
for some δ < 1. We assume d p , 0 ≤ p ≤ ς -1 is chosen to satisfy the following conditions: there exist S dp ∪ D(0, ρ) and ρ > 0 such that

• Conditions 1), 2),
R dς := R d 0 ). The family {(R dp, δ ) 0≤p≤ς-1 , D(0, ρ), T } is said to be associated to the good covering (E p ) 0≤p≤ς-1 .
Let (E p ) 0≤p≤ς-1 be a good covering, and a family {(R dp, δ ) 0≤p≤ς-1 , D(0, ρ), T } associated to it. For every 0 ≤ p ≤ ς -1 we study the following equation

(5.1) Q(∂ z )σ q,t u dp (t, z, ) = ( t) d D 1 σ d D 1 k 1 +1 q,t R D 1 (∂ z )u dp (t, z, ) + ( t) d D 2 σ d D 2 k 2 +1 q,t R D 2 (∂ z )u dp (t, z, ) + D-1 =1 ∆ t d σ δ q,t (c (t, z, )R (∂ z )u dp (t, z, )) + σ q,t f (t, z, ).
The terms c (t, z, ) are determined as follows, for every 1 ≤ ≤ D -1. Let C (T, m, ) be the entire function in T , with coefficients in E (β,µ) for some β > 0 and µ ∈ R, given by

C (T, m, ) = n≥0 C ,n (m, )T n , such that µ -1 ≥ deg(R D j ), for j ∈ {1, 2}.

Assume this function depends holomorphically on

∈ D(0, 0 ) and also the existence of C , T 0 > 0 such that the left-hand side of (4.6) holds for all n ≥ 0 and ∈ D(0, 0 ). We put

c (t, z, ) := F -1 (m → C ( t, m, )) (z),
which is a holomorphic and bounded function on T × H β × D(0, 0 ), with 0 < β < β. Indeed, one can substitute T by any bounded set in C in the previous product domain.

The function f (t, z, ) is constructed as follows. Let m → F n (m, ) be a function in E (β,µ) for every n ≥ 0, depending holomorphically on ∈ D(0, 0 ). We also assume there exist C F , T 0 such that (4.6) holds and define F (T, m, ) = n≥0 F n T n . By construction, F (T, m, ) represents a holomorphic function in T on the disc D(0, T 0 /2) with values in the Banach space E (β,µ) , for all ∈ D(0, 0 ). We define

f (t, z, ) = F -1 (m → F ( t, m, ))(z)
which stands for a holomorphic and bounded function on D(0, 0 T 0 /2) × H β × D(0, 0 ), for all 0 < β < β.

Theorem 5.3. Under the construction made at the beginning of this section of the elements involved in the problem (5.1), assume that the above conditions hold. Let (E p ) 0≤p≤ς-1 be a good covering in C , for which a family {(R dp, δ ) 0≤p≤ς-1 , D(0, ρ), T } associated to this covering is considered.

Then, there exist large enough r

Q,R D 1 , r Q,R D 2 > 0 and constants ς Ψ > 0 and ς ϕ > 0 such that if C ≤ ς ϕ , C Ψ 1 ≤ ς ψ
for all 1 ≤ ≤ D -1, then for every 0 ≤ p ≤ ς -1, one can construct a solution u dp (t, z, ) of (5.1), which defines a holomorphic function on T × H β × E p , for every 0 < β < β.

Proof. Let 0 ≤ p ≤ ς -1 and consider the equation

Q(im)σ q,T U dp (T, m, ) = T d D 1 σ d D 1 k 1 +1 q,T R D 1 (im)U dp (T, m, ) + T d D 2 σ d D 2 k 2 +1 q,T R D 2 (im)U dp (T, m, ) + D-1 =1 ∆ -d T d σ δ q,T 1 (2π) 1/2 +∞ -∞ C (T, m -m 1 , )R (im 1 )U dp (T, m 1 , )dm 1 (5.2) +σ q,T F (T, m, ).
Under an appropriate choice of the constants ς Ψ and ς ϕ one can follow the construction in Section 4.1 and apply Proposition 4.2 to obtain a solution U dp (T, m, ) of (5.5).

Regarding the properties of q-Laplace transformation, and from the results obtained in Section 4.2, U dp (T, m, ) is the q-Laplace transformation of order k 2 of a function w dp k 2 along direction d p , which depends on T . Indeed,

(5.3) U dp (T, m, ) = 1 π q 1/k 2 Ld p w dp k 2 (u, m, ) Θ q 1/k 2 u T du u ,
for some L dp ⊆ S dp ∪ {0}, and w dp k 2 (τ, m, ) defines a continuous function on S dp × R × D(0, 0 ), and holomorphic with respect to (τ, ) in S dp ×D(0, 0 ). In addition to this, there exists

C w dp k 2 > 0 such that (5.4) |w dp k 2 (τ, m, )| ≤ C w dp k 2 1 (1 + |m|) µ e -β|m| exp k 2 2 log(q) log 2 |τ | + ν log |τ | ,
for some ν ∈ R. This holds for τ ∈ S dp , m ∈ R, and ∈ D(0, 0 ). Moreover, in view of Proposition 4.5, the function w dp k 2 (τ, m, ) and the q-Laplace transformation of order κ of the function w

dp k 1 (τ, m, ) along direction d 1 p , where e id 1 p R + ⊆ S dp ∪ {0}, depending on τ , coincide in (S dp ∩ D(0, r 1 )) × R × D(0, 0 ), for 0 < r 1 ≤ q ( 1 2 -α)/κ /2, for some α ∈ R. The function w dp k 1 (τ, m, ) is such that (5.5) |w dp k 1 (τ, m, )| ≤ C w dp k 1 1 (1 + |m|) µ e -β|m| exp κ 2 log(q) log 2 |τ + δ| + α log |τ + δ| , for some C w dp k 1
, δ > 0, valid for τ ∈ (D(0, ρ) ∪ U dp ), m ∈ R and ∈ D(0, 0 ). This function is the extension of a function w k 1 (τ, m, ), common for every 0 ≤ p ≤ ς -1, continuous on D(0, ρ) × R × D(0, 0 ) and holomorphic with respect to (τ, ) in D(0, ρ) × D(0, 0 ). The bounds in (5.7) with respect to m variable are transmitted to U dp (T, m, ) as defined in (5.6). This allows to define the function

u dp (t, z, ) := F -1 (m → U dp ( t, m, ))(z) = 1 (2π) 1/2 1 π 1/k 2 q ∞ -∞ Ld p w dp k 2 (u, m, ) Θ q 1/k 2 u t du u exp(izm)dm,
which turns out to be holomorphic on T × H β × E p . The properties of inverse Fourier transform allow us to conclude that u dp (t, z, ) is a solution of equation ( 5.1) defined on T × H β × E p .

Proposition 5.4. Let 0 ≤ p ≤ ς -1. Under the hypotheses of Theorem 5.3, assume that the unbounded sectors U dp and U d p+1 are wide enough so that U dp ∩ U d p+1 contains the sector

U dp,d p+1 = {τ ∈ C : arg(τ ) ∈ [d p , d p+1 ]}. Then, there exist K 1 > 0 and K 2 ∈ R such that (5.6) |u d p+1 (t, z, ) -u dp (t, z, )| ≤ K 1 exp - k 2 2 log(q) log 2 | | | | K 2 ,
for every t ∈ T , z ∈ H β , and (u, m, )

∈ E p ∩ E p+1 . Proof. Let 0 ≤ p ≤ ς -1. Taking into account that U dp,d p+1 ⊆ U dp ∩ U d p+1 ,
Θ q 1/k 2 u t is holomorphic on R b dp ∪ R b d p+1 for all (m, ) ∈ R × (E p ∩ E p+1
) and its integral along the closed path constructed by concatenation of the segment starting at the origin and with ending point fixed at ρe idp , the arc of circle with radius ρ connecting ρe idp with ρe id p+1 ⊆ R b d p+1 , and the segment from ρe id p+1 to 0, vanishes. The difference u d p+1 -u dp can be written in the form

(5.7) u d p+1 (t, z, ) -u dp (t, z, ) = 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ L d p+1 , ρ w d p+1 k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm, - 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ L dp, ρ w dp k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm + 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ C ρ,dp,d p+1 w dp,d p+1 k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm,
where L d j , ρ = [ ρ, +∞)e id j for j ∈ {p, p + 1} and C ρ,dp,d p+1 is the arc of circle connecting ρe idp with ρe id p+1 (see Figure 2). Let us put

I 1 := 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ L d p+1 , ρ w d p+1 k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm .
In view of (5.7) and (3.2), one has

I 1 ≤ C w d p+1 k 2 C q,k 2 δ(2π) 1/2 | t| 1/2 π q 1/k 2 ∞ -∞ e -β|m|-m (z) dm (1 + |m|) µ × ∞ ρ exp k 2 log 2 |u| 2 log(q) + ν log |u| |u| -3/2 exp   - k 2 log 2 |u| | t| 2 log(q)   d|u|.
We recall that we have restricted the domain on the variable z such that | (z)| ≤ β < β. Then, the first integral in the previous expression in convergent, and one derives

I 1 ≤ Cw d p+1 k 2 (2π) 1/2 ( 0 r T ) 1/2 π q 1/k 2 ∞ ρ exp k 2 log 2 |u| 2 log(q) exp   - k 2 log 2 |u| | t| 2 log(q)   |u| ν-3/2 d|u|,
for some Cw

d p+1 k 2 > 0. We derive exp k 2 log 2 |u| 2 log(q) exp   - k 2 log 2 |u| | t| 2 log(q)   = exp k 2 2 log(q) (-log 2 | | -2 log | | log |t| -log 2 |t|) × exp k 2 log(q) (log |u| log | | + log |u| log |t|) .
From the assumption that 0 < 0 < 1 and 0 < r T < 1, we get exp -

k 2 log(q) log | | log |t| ≤ | | - k 2 log(q) log(r T ) , (5.8) exp k 2 log(q) log |u| log | | ≤ | | k 2 log(q) log( ρ) ,
for t ∈ T , ∈ E p ∩ E p+1 , |u| ≥ ρ, and also (5.9) exp

k 2 log(q) log |u| log |t| ≤ |t| k 2 log(q) log( ρ) , if ρ ≤ |u| ≤ 1 exp k 2 log(q) log |u| log |t| ≤ |u| k 2 log(q) log(r T ) , if |u| ≥ 1,
for t ∈ T . In addition to that, there exists K k 2 , ρ,q > 0 such that (5.10) sup

x>0 x k 2 log(q) log( ρ) exp - k 2 2 log(q) log 2 (x) ≤ K k 2 , ρ,q .
In view of (5.11), (5.12), (5.13), and bearing in mind that the inequalities of Definition 5.2 hold, we deduce there exist K1 ∈ R, K2 > 0 such that

exp k 2 log 2 |u| 2 log(q) exp   - k 2 log 2 |u| | t| 2 log(q)   |u| ν ≤ K2 exp - k 2 2 log(q) log 2 | | | | K1 ,
for t ∈ T , r ≥ ρ, and ∈ E p ∩ E p+1 . Provided this last inequality, we arrive at (5.11)

I 1 ≤ K2 Cw d p+1 k 2 (2π) 1/2 ( 0 r T ) 1/2 π q 1/k 2 ∞ ρ d|u| |u| 3/2 exp - k 2 2 log(q) log 2 | | | | K1 = K3 exp - k 2 2 log(q) log 2 | | | | K1 ,
for some K3 > 0, for all t ∈ T , z ∈ H β , and ∈ E p ∩ E p+1 . We can estimate in the same manner the expression

I 2 := 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ L dp, ρ w dp k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm
to arrive at the existence of K4 > 0 such that (5.12)

I 2 ≤ K4 exp - k 2 2 log(q) log 2 | | | | K1 ,
for all t ∈ T , z ∈ H β , and ∈ E p ∩ E p+1 . We now provide upper bounds for the quantity

I 3 := 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ C ρ,dp,d p+1 w dp,d p+1 k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm .
From the construction of w > 0 such that

I 3 ≤ Cdp,d p+1 w k 2 ∞ -∞ e -β|m|-m (z) (1 + |m|) µ dm|d p+1 -d p ||t| 1/2 exp   - k 2 log 2 ρ | t| 2 log(q)   ,
for all t ∈ T , z ∈ H β , and ∈ E p ∩ E p+1 . We can follow analogous arguments as in the previous steps to provide upper estimates of the expression

|t| 1/2 exp   - k 2 log 2 ρ | t| 2 log(q)   .
Indeed,

|t| 1/2 exp   - k 2 log 2 ρ | t| 2 log(q)   = exp - k 2 log 2 ( ρ) 2 log(q) | | k 2 log( ρ) log(q) |t| k 2 log( ρ) log(q) × exp k 2 2 log(q) (-log 2 | | -2 log | | log |t| -log 2 |t|) |t| 1/2 .
From the assumption 0 ≤ 0 < 1 we check that exp -

k 2 log(q) log | | log |t| ≤ | | - k 2 log(q) log(r T ) ,
for t ∈ T , ∈ E p ∩ E p+1 . Gathering (5.13), we get the existence of K5 ∈ R, K6 > 0 such that

|t| 1/2 exp   - k 2 log 2 ρ | t| 2 log(q)   ≤ K6 exp - k 2 2 log(q) log 2 | | | | K5 ,
to conclude that (5.13)

I 3 ≤ K7 exp - k 2 2 log(q) log 2 | | | | K5 ,
for some K7 > 0, all t ∈ T , z ∈ H β , and ∈ E p ∩ E p+1 . We conclude the proof of this result in view of (5.14), (5.15), (5.16) and the decomposition (5.10).

Lemma 5.5. Let 0 ≤ p ≤ ς -1. Under the hypotheses of Theorem 5.3, assume that

U dp ∩ U d p+1 = ∅. Then, there exist K L p > 0, M L p ∈ R such that L d p+1 q;1/κ (w d p+1 k 1 )(τ, m, ) -L dp q;1/κ (w dp k 1 )(τ, m, ) ≤ K L p e -β|m| (1 + |m|) -µ exp - κ 2 log(q) log 2 |τ | |τ | M L p , for every ∈ (E p ∩ E p+1 ), τ ∈ (R b dp ∩ R b d p+1 ) and m ∈ R.
Proof. We first recall that, without loss of generality, the intersection R b dp ∩R b d p+1 can be assumed to be a nonempty set because one can vary δ in advance to be as close to 0 as desired.

Analogous arguments as in the beginning of the proof of Proposition 5.4 allow us to write (5.14) L d p+1 q;1/κ (w

d p+1 k 1 )(τ, m, ) -L dp q;1/κ (w dp k 1 )(τ, m, ) = 1 π q 1/κ L d p+1 , ρ w d p+1 k 1 (u, m, ) Θ q 1/κ u τ du u - 1 π q 1/κ L dp, ρ w dp k 1 (u, m, ) Θ q 1/κ u τ du u + 1 π q 1/κ C ρ,dp,d p+1 w k 1 (u, m, ) Θ q 1/κ u τ du u ,
where ρ, L dp, ρ , L d p+1 , ρ and C ρ,dp,d p+1 are constructed in Proposition 5.4. In view of (5.8) and (3.2), one has

I L 1 := 1 π q 1/κ L dp, ρ w dp k 1 (u, m, ) Θ q 1/κ u τ du u ≤ C w dp k 1 C q,κ δ |τ | 1/2 (1 + |m|) µ e -β|m| ∞ ρ exp κ log 2 |re idp +δ| 2 log(q) + α log |re idp + δ| exp κ 2 log 2 r |τ | log(q) dr r 3/2 ≤ K L p,1 |τ | 1/2 (1 + |m|) -µ e -β|m| ∞ ρ exp κ log 2 r 2 log(q) + α log r exp κ 2 log 2 r |τ | log(q) dr r 3/2
for some K L p,1 > 0. Usual calculations, and taking into account the choice of α in Definition 5.2, one derives the previous expression equals

K L p,1 |τ | 1/2 (1 + |m|) -µ e -β|m| exp - κ 2 log(q) log 2 |τ | ∞ ρ r κ log |τ | log(q) +α-3/2 dr.
Besides, we observe from direct computations that

∞ ρ r κ log |τ | log(q) +α-3 2 dr
is upper bounded by a constant times |τ | ∆ , for every ∆ < κ log(q) log(ρ). This yields (5.15)

I L 1 ≤ K L p,2 (1 + |m|) -µ e -β|m| exp - κ 2 log(q) log 2 |τ | |τ | ∆+ 1 2 ,
for some K L p,2 > 0. Analogous arguments allow us to obtain the existence of K L p,3 > 0 such that (5.16)

I L 2 := 1 π q 1/κ L d p+1 , ρ w d p+1 k 1 (u, m, ) Θ q 1/κ u τ du u ≤ K L p,3 (1 + |m|) -µ e -β|m| exp - κ 2 log(q) log 2 |τ | |τ | ∆+ 1 2 ,
for ∆ as above.

We write

I L 3 := 1 π q 1/κ C ρ,dp,d p+1 w k 1 (u, m, ) Θ q 1/κ u τ du u .
Regarding (5.8) and (3.2), one derives that

I L 3 ≤ C w dp k 1 π q 1/κ e -β|m| (1 + |m|) µ |τ | 1/2 ρ 1/2 C q,κ δ d p+1 dp exp κ log 2 | ρe iθ +δ| 2 log(q) + α log | ρe iθ + δ| exp κ 2 log 2 ρ |τ | log(q) dθ ≤ K L p,4 |τ | 1/2 e -β|m| (1 + |m|) µ exp   - κ 2 log 2 ρ |τ | log(q)   with K L p,4 = |d p+1 -d p | C w dp k 1 π q 1/κ 1 ρ 1/2 C q,κ δ exp κ log 2 ( ρ + δ) 2 log(q) + α log( ρ + δ) . Let K L p,5 = K L p,4 exp(-κ 2 log(q) log 2 ( ρ)).
It is straightforward check that (5.17)

I L 3 ≤ K L p,5 |τ | 1/2+ κ log( ρ) log(q) e -β|m| (1 + |m|) µ exp - κ 2 log 2 |τ | log(q) .
From (5.18), (5.19) and (5.20), put into (5.17), we conclude the result.

Proposition 5.6. Let 0 ≤ p ≤ ς -1. Under the hypotheses of Theorem 5.3, assume that U dp ∩ U d p+1 = ∅. Then, there exist K 3 > 0 and K 4 ∈ R such that

(5.18) |u d p+1 (t, z, ) -u dp (t, z, )| ≤ K 3 exp - k 1 2 log(q) log 2 | | | | K 4 ,
for every t ∈ T , z ∈ H β , and ∈ E p ∩ E p+1 .

Proof. Let 0 ≤ p ≤ ς -1. Under the assumptions of the statement, we observe that one can not proceed as in the proof of Proposition 5.4 for there does not exist a common function for both indices p and p + 1, defined in R b dp ∪ R b d p+1 in the variable of integration, when applying q-Laplace transform. However, one can use the analytic continuation property and write the difference u d p+1 -u dp as follows. Let ρ > 0 be such that ρe idp ∈ R b dp and ρe id p+1 ∈ R b d p+1 , and let θ p,p+1 ∈ R be such that ρe θ p,p+1 lies in both R b dp and R b d p+1 . We write u d p+1 (t, z, ) -u dp (t, z, ) as follows

J 2 := 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ L dp, ρ w dp k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm ≤ K2 exp - k 2 2 log(q) log 2 | | | | K6 , J 3 := 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ C ρ,θ p,p+1 ,d p+1 w dp,d p+1 k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm ≤ K3 exp - k 2 2 log(q) log 2 | | | | K7 , J 4 := 1 (2π) 1/2 1 π q 1/k 2 ∞ -∞ C ρ,θ p,p+1 ,dp w dp,d p+1 k 2 (u, m, ) Θ q 1/k 2 u t exp(izm) du u dm ≤ K4 exp - k 2 2 log(q) log 2 | | | | K8
We now give estimates for

J 5 := 1 (2π) 1/2 1 π q 1/k 2 × ∞ -∞ L 0, ρ,θ p,p+1 L d p+1 q;1/κ (w d p+1 k 1 )(u, m, ) -L dp q;1/κ (w dp k 1 )(u, m, ) Θ q 1/k 2 u t exp(izm) du u dm .
In view of Lemma 5.5 and (3.2), one has

J 5 ≤ K L p (2π) 1/2 1 π q 1/k 2 ∞ -∞ e -β|m|-(z)m dm (1 + |m|) µ ρ 0 exp -κ 2 log(q) log 2 |u| |u| M L p C q,k 2 δ exp k 2 2 log 2 | u t | log(q) u t 1/2 d|u| |u| .
We recall that z ∈ H β for some β < β. Then, there exists K 31 > 0 such that

J 5 ≤ K L p K 31 (2π) 1/2 | | 1/2 r 1/2 T π q 1/k 2 C q,k 2 δ ρ 0 exp -κ 2 log(q) log 2 |u| |u| M L p exp k 2 2 log 2 | u t | log(q) d|u| |u| 3/2 .
We now proceed to prove the expression

ρ 0 exp -κ 2 log(q) log 2 |u| exp k 2 2 log 2 | u t | log(q) exp k 1 2 log(q) log 2 | | d|u| |u| 3/2-M L p
is upper bounded by a positive constant times a certain power of | | for every ∈ (E p ∩ E p+1 ) and t ∈ T . This concludes the existence of K 32 > 0 such that (5.20)

J 5 ≤ K 32 | | 1/2 exp - k 1 2 log(q) log 2 | | , for every ∈ (E p ∩ E p+1 ), t ∈ T and z ∈ H β . Indeed, we have ρ 0 exp -κ 2 log(q) log 2 |u| exp k 2 2 log 2 |u| | t| log(q) exp k 1 2 log(q) log 2 | | d|u| |u| 3/2-M L p equals (5.21) exp k 1 2 log(q) log 2 | | - k 2 2 log(q) log 2 | t| × ρ 0 exp - (κ + k 2 ) 2 log(q) log 2 |u| |u| k 2 log | t| log(q) -3 2 +M L p d|u|. Given m 1 ∈ R and m 2 > 0, the function [0, ∞) x → H(x) = x m 1 exp(-m 2 log 2 (x)) attains its maximum value at x 0 = exp( m 1 2m 2 ) with H(x 0 ) = exp( m 2 1 4m
2 ). This yields an upper bound for the integrand in (5.24); the expression in (5.24) is estimated from above by (5.22)

ρ exp (M L p -3/2) 2 log(q) 2(κ + k 2 ) exp 1 2 log(q) ( k 2 2 κ + k 2 -k 2 + k 1 ) log 2 | | × exp 1 2 log(q) ( k 2 2 κ + k 2 -k 2 ) log 2 |t| |t| k 2 (M L p -3/2) κ+k 2 × exp 1 log(q) ( k 2 2 κ + k 2 -k 2 ) log | | log |t| | | k 2 (M L p -3/2) κ+k 2 .
The second line in (5.25) is upper bounded for every t because

K 33 | | k 2 (M L p -3/2) κ+k 2
for some K 33 > 0. The conclusion is achieved. The result follows from (5.22), the estimates J 1 to J 4 , and (5.23).

Existence of formal series solutions in the complex parameter and asymptotic expansion in two levels

In the first part of this section, we remind two q-analogs of Ramis-Sibuya theorem from [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF][START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF]. This result provides the tool to guarantee the existence of a formal power series in the perturbation parameter which formally solves the main problem and such that it asymptotically represents the analytic solution of that equation.

This asymptotic representation is held in the sense of q-asymptotic expansions of certain positive order. Definition 6.1. Let V be a bounded open sector with vertex at 0 in C. Let (F, • F ) be a complex Banach space. Let q ∈ R with q > 1 and let k be a positive integer. We say that a holomorphic function f : V → F admits the formal power series f ( ) = n≥0 f n n ∈ F[[ ]] as its q-Gevrey asymptotic expansion of order 1/k if for every open subsector U with (U \ {0}) ⊆ V , there exist A, C > 0 such that f ( )

- N n=0 f n n F ≤ CA N +1 q N (N +1) 2k | | N +1 ,
for every ∈ U , and N ≥ 0.

The set of functions which admit null q-Gevrey asymptotic expansion of certain positive order are characterized as follows. The proof of this result, already stated in [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF], provides the q-analog of Theorem XI-3-2 in [START_REF] Hsieh | Basic theory of ordinary differential equations[END_REF]. Lemma 6.2. A holomorphic function f : V → F admits the null formal power series 0 ∈ F[[ ]] as its q-Gevrey asymptotic expansion of order 1/k if and only if for every open subsector U with (U \ {0}) ⊆ V there exist constants K 1 ∈ R and K 2 > 0 with

f ( ) F ≤ K 2 exp - k 2 log(q) log 2 | | | | K 1 ,
for all ∈ U .

The next result leans on the one level version of the q-analog of Ramis Sibuya theorem, stated in [START_REF] Malek | Parametric Gevrey asymptotics for a q-analog of some linear initial value problem[END_REF], provides a two level result in this framework. See [START_REF]On parametric multilevel q-Gevrey asymptotics for some linear q-difference-differential equations[END_REF] for a proof. Theorem 6.3. Let (F, • F ) be a Banach space and (E p ) 0≤p≤ς-1 be a good covering in C . Let 0 < k 1 < k 2 , consider a holomorphic function G p : E i → F for every 0 ≤ p ≤ ς -1 and put ∆ p ( ) = G p+1 ( ) -G p ( ) for every ∈ Z p := E p ∩ E p+1 . Moreover, we assume:

1) The functions G p ( ) are bounded as tends to 0 on E p for every 0 ≤ p ≤ ς -1.

2) There exist nonempty sets I 1 , I 2 ⊆ {0, 1, . . . , ς -1} such that I 1 ∪ I 2 = {0, 1, . . . , ς -1} and I 1 ∩ I 2 = ∅. Also, -for every p ∈ I 1 there exist constants K 1 > 0, M 1 ∈ R such that

∆ p ( ) F ≤ K 1 | | M 1 exp - k 1 2 log(q) log 2 | | , ∈ Z p ,
-and, for every p ∈ I 2 there exist constants K 2 > 0, M 2 ∈ R such that

∆ p ( ) F ≤ K 2 | | M 2 exp - k 2 2 log(q) log 2 | | , ∈ Z p .
Then, there exists a convergent power series a( ) ∈ F{ } defined on some neighborhood of the origin and Ĝ1 ( ), Ĝ2 ( ) ∈ F[[ ]] such that G p can be written in the form

G p ( ) = a( ) + G 1 p ( ) + G 2 p ( ).
G 1 p ( ) is holomorphic on E p and admits Ĝ1 ( ) as its q-Gevrey asymptotic expansion of order 1/k 1 on E p , for every p ∈ I 1 ; whilst G 2 p ( ) is holomorphic on E p and admits Ĝ2 ( ) as its q-Gevrey asymptotic expansion of order 1/k 2 on E p , for every p ∈ I 2 .

We conclude this section with the main result in the work in which we guarantee the existence of a formal solution of the main problem (5.1), written as a formal power series in the perturbation parameter, with coefficients in an appropriate Banach space, say û(t, z, ). Moreover, it represents, in some sense to be precised, each solution u dp (t, z, ) of the problem (5.1).

From now on, F stands for the Banach space of bounded holomorphic functions defined on T × H β , with the supremum norm, where β < β, as above. formal solution of the equation (6.2) Q(∂ z )σ q,t û(t, z, )

= ( t) d D 1 σ d D 1 k 1 +1 q,t R D 1 (∂ z )û(t, z, ) + ( t) d D 2 σ d D 2 k 2 +1 q,t R D 2 (∂ z )û(t, z, ) + D-1 =1 ∆ t d σ δ
q,t (c (t, z, )R (∂ z )û(t, z, )) + σ q,t f (t, z, ).

Moreover, û(t, z, ) turns out to be the common q-Gevrey asymptotic expansion of order 1/k 1 on E p of the function u dp , seen as holomorphic function from E p into F, for 0 ≤ p ≤ ς -1. In addition to that, û is of the form û(t, z, ) = a(t, z, ) + û1 (t, z, ) + û2 (t, z, ),

where a(t, z, ) ∈ F{ } and û1 (t, z, ), û2 (t, z, ) ∈ F[[ ]] and such that for every 0 ≤ p ≤ ς -1, the function u dp can be written in the form u dp (t, z, ) = a(t, z, ) + u where → u dp 1 (t, z, ) is a F-valued function that admits û1 (t, z, ) as its q-Gevrey asymptotic expansion of order 1/k 1 on E p and also → u dp 2 (t, z, ) is a F-valued function that admits û2 (t, z, ) as its q-Gevrey asymptotic expansion of order 1/k 2 on E p .

Proof. For every 0 ≤ p ≤ ς -1, one can consider the function u dp (t, z, ) constructed in Theorem 5.3. We define G p ( ) := (t, z) → u dp (t, z, ), which is a holomorphic and bounded function from E p into F. In view of Proposition 5.4 and Proposition 5.6, one can split the set {0, 1, . . . , ς -1} in two nonempty subsets of indices, I 1 and I 2 with {0, 1, . . . , ς -1} = I 1 ∪ I 2 and such that I 1 (resp. I 2 ) consists of all the elements in {0, 1, . . . , ς -1} such that U dp ∩U d p+1 contains the sector U dp,d p+1 , as defined in Proposition 5.4 (resp. U dp ∩ U d p+1 = ∅). From (5.9) and (??) one can apply Theorem 6.4 and deduce the existence of formal power series Ĝ1 ( ), Ĝ2 ( ) ∈ F[[ ]], a convergent power series a( ) ∈ F{ } and holomorphic functions G 1 p ( ), G 2 p ( ) defined on E p and with values in F such that G p ( ) = a( ) + G 1 p ( ) + G 2 p ( ), and for j = 1, 2, one has G j p ( ) admits Ĝj ( ) as its q-Gevrey asymptotic expansion or order 1/k j on E p . We put û(t, z, ) = It only rests to prove that û(t, z, ) is the solution of (6.2). Indeed, since u dp admits û(t, z, ) as its q-Gevrey asymptotic expansion of order 1/k 1 on E p , we have that lim →0, ∈Ep sup t∈T ,z∈H β |∂ m u dp (t, z, ) -h m (t, z)| = 0, for every 0 ≤ p ≤ ς -1 and m ≥ 0. Let p ∈ {0, 1, . . . , ς -1}. By construction, the function u dp (t, z, ) solves equation (6.2). We take derivatives of order m ≥ 0 with respect to at both sides of equation (5.1) and deduce that (6.3) Q(∂ z )σ q,t (∂ m u dp )(t, z, )

= m 1 +m 2 =m m! m 1 !m 2 ! ∂ m 1 ( d D 1 )t d D 1 σ d D 1 k 1 +1 q,t R D 1 (∂ z )(∂ m 2 u dp ) + m 1 +m 2 =m m! m 1 !m 2 ! ∂ m 1 ( d D 2 )t d D 2 σ d D 2 k 2 +1 q,t R D 2 (∂ z )(∂ m 2 u dp ) + D-1 =1 m 1 +m 2 +m 3 =m m! m 1 !m 2 !m 3 ! (∂ m 1 ∆ )t d σ δ q,t (∂ m 2 c (t, z, )R (∂ z )∂ m 3 u dp (t, z, ))
+ σ q,t (∂ m f )(t, z, 0), for every (t, z, ) ∈ T × H β × E p . We let → 0 in (6.3) and obtain the recursion formula (6.4) Q(∂ z )σ q,t h m (t, z) By plugging (6.1) into (6.2) and bearing in mind (6.4) and (6.5) one concludes that the formal power series û(t, z, ) = m≥0 h m (t, z) m /m! is a solution of equation (6.2).

= m! (m -d D 1 )! t d D 1 σ d D 1 k 1 +1 q,t R D 1 (∂ z )(h m-d D 1 (t, z))+ m! (m -d D 2 )! t d D 2 σ d D 2 k 2 +1 q,t R D 2 (∂ z )(h m-d D 2 (t, z)) + D-1 =1 m 2 +m 3 =m-∆ m! m 2 !m 3 ! t d σ δ q,t (∂ m 2 c (t, z, 0)R (∂ z )h m 3 (t, z)) + σ q,t (∂ m f )(t,

  and τ ∈ S d . 2) There exists M 22 > 0 such that we have |τ -q ,2 (m)| ≥ M 22 |q ,2 (m)| for all ∈ {0, . . . , d D 2 -1}, m ∈ R and τ ∈ S d .

Proposition 4 . 3 .

 43 Under the hypotheses (4.1), (4.2), (4.3), (4.4), (4.5) and those on the geometry of the problem, there exist r Q,R D 2 , > 0 and ς ϕ , ς Ψ > 0 under (4.14) such that for every ∈ D(0, 0 ), the equation (4.32) admits a unique solution w d k 2 (τ, m, ) in the space Exp q

  37), (4.38), (4.39) and (4.40) one concludes (4.35).

Proposition 4 . 4 .

 44 We consider w d k 1 (τ, m, ) constructed in Proposition 4.2. The function

  Page 14 in Section 4.1 hold. Observe that, under this assumption, Conditions 1), 2), Page 20 in Section 4.2 hold for S dp . • For every 0 ≤ p ≤ ς -1 we have R b dp ∩ R b d p+1 = ∅, and for every t ∈ T and ∈ E p we have t ∈ R b dp (where

q; 1 /κ (w d p+1 k 1 )dp,d p+1 k 2 (

 112 we observe from the construction of the functions U dp and U d p+1 that L dp q;1/κ (w dp k 1 )(τ, m, ) and L d p+1 (τ, m, ) coincide in the domain (R b dp ∩R b d p+1 )×R×D(0, 0 ). This entails the existence of w τ, m, ), holomorphic with respect to τ on R b dp ∪R b d p+1 , continuous with respect to m ∈ R and holomorphic with respect to in D(0, 0 ) which coincides with L dp q;1/κ (w dp k 1 )(τ, m, ) on R b dp × R × D(0, 0 ) and also with L d p+1 q;1/κ (w d p+1 k 1 )(τ, m, ) on R b d p+1 × R × D(0, 0 ). Let ρ > 0 be such that ρe idp ∈ R b dp and ρe id p+1 ∈ R b d p+1 . The function u → w dp,d p+1 k 2

Figure 2 .

 2 Figure 2. Deformation of the path of integration, first case.

dp,d p+1 k 2 (τ, m, ), we have |w dp,d p+1 k 2 ( 1 1 ( 1 + 1 >

 22111 u, m, )| ≤ Cw dp k |m|) µ e -β|m| , for some Cw dp k 0, valid for u ∈ C ρ,dp,d p+1 , m ∈ R and ∈ D(0, 0 ). The estimates with (3.2) allow us to obtain the existence of Cdp,d p+1 w k 2

k 2 2 κ+k 2 2 2κ+k 2 - 2 2κ + k 2 -k 2 = -k 1 ,

 22222221 < k 2 and also, one has an upper bound for exp 1 log(q) ( k k 2 ) log | | log |t| is 1. Regarding Definition 5.2, and taking into account that k the expression (5.25) is upper bounded by

Theorem 6 . 4 .

 64 Under the hypotheses of Theorem 5.3, there exists a formal power series

dp 1 (

 1 t, z, ) + u dp 2 (t, z, ),

  z, 0), for every m ≥ max{d D 1 , d D 2 , max 1≤ ≤D-1 ∆ }, and all (t, z) ∈ T × H β .Bearing in mind that both c l and f are holomorphic w.r.t in a neighborhood of the origin, in such neighborhood one has (6.5)c (t, z, ) = m≥0 (∂ m c m )(t, z, 0) m! m , f (t, z, ) = m≥0 (∂ m f )(t, z, 0) m! mfor every 1 ≤ ≤ D -1.

  Thus, taking into account assumptions (4.1), (4.5), regarding (4.13) together with Proposition 2.3 and Proposition 2.6 we get (4.16)

	∆ -d

  m, ) and Ψ k 2 (τ, m, ) belong to H (k 2 ,β,µ,Ω) due to Proposition 4.2, Proposition 4.4. As we can see in the proof of Proposition 4.4, the operator H 2 defined in (4.34) has a unique fixed point in H (k 2 ,β,µ,Ω) provided small enough constants ς Ψ , ς ϕ > 0, for 1 ≤ ≤ D -1. Indeed, this fixed point is a solution of the auxiliary problem (4.32) in the disc

Acknowledgements: We want to express our gratitude to the anonimous referee for the valuable comments and suggestions made which helped to improve the work. Funding: This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under the Grant Agreement No 648132. The second and third authors are partially supported by the research project MTM2016-77642-C2-1-P of Ministerio de Economía, Industria y Competitividad, Spain. Competing interests: The authors declare that they have no competing interests. Availability of data and material: Data sharing is not applicable to this article as no datasets were generated or analysed during the current study. Authors' contributions: All authors contributed equally and significantly in writing this paper and typed, read, and approved the final manuscript.

 (5.19)Here, we have denoted L d j , ρ = [ ρ, +∞)e id j for j ∈ {p, p + 1}, C ρ,θ p,p+1 ,d p+1 is the arc of circle connecting ρe id p+1 with ρe iθ p,p+1 , C ρ,θ p,p+1 ,dp is the arc of circle connecting ρe idp with ρe iθ p,p+1 , L 0, ρ,θ p,p+1 = [0, ρ]e iθ p,p+1 , as it is shown in the following figure. Following the same line of arguments as those in the proof of Proposition 5.4, we can guarantee the existence of Kj > 0 and Kk ∈ R for 1 ≤ j ≤ 4 and 5 ≤ k ≤ 8 such that