On a $q-$analog of a singularly perturbed problem of irregular type with two complex time variables

TL;DR

This paper develops analytic solutions and formal asymptotic expansions for a family of $q$-difference-differential equations with two complex time variables, extending previous work on singularly perturbed PDEs and highlighting the role of Lambert W function branches.

Contribution

It introduces a $q$-analog of singularly perturbed PDEs with two complex time variables, constructing solutions with different asymptotic behaviors and emphasizing the significance of Lambert W function branches.

Findings

Construction of outer and inner analytic solutions

Solutions exhibit different asymptotic expansions

Lambert W function's $-1$-branch is crucial

Abstract

Analytic solutions and their formal asymptotic expansions for a family of the singularly perturbed $q-$difference-differential equations in the complex domain are constructed. They stand for a $q-$analog of the singularly perturbed partial differential equations considered in our recent work [A. Lastra, S. Malek, Boundary layer expansions for initial value problems with two complex time variables, submitted 2019]. In the present work, we construct outer and inner analytic solutions of the main equation, each of them showing asymptotic expansions of essentially different nature with respect to the perturbation parameter. The appearance of the $-1$-branch of Lambert $W$ function will be crucial in this respect.

Equations

$$Q({\partial }_z)u(\mathbf{t},z,ϵ)=P(\mathbf{t},z,ϵ,{\partial }_{t_2},{\sigma }_{q;t_1},{\sigma }_{q;t_2})u(\mathbf{t},z,ϵ)+f(\mathbf{t},z,ϵ),$$

$${\sigma }_{q;t}(f(t)):=f(qt).$$

$$Q({\partial }_z)u(\mathbf{t},z,ϵ)=P(t_1^{k_1+1}{\partial }_{t_1},t_2^{k_2+1}{\partial }_{t_2},{\partial }_z,z,ϵ)u(\mathbf{t},z,ϵ)+f(\mathbf{t},z,ϵ),$$

$$u(\mathbf{t},z,ϵ)=\frac{1}{(2\pi {)}^{1/2}{\pi }_{q^{1/k_1}}}{\int }_{-\infty }^{\infty }{\int }_{L_{\gamma }}\omega (u,m,ϵ)\frac{1}{{\Theta }_{q^{1/k_1}}\left(\frac{u}{{ϵ}^{{\lambda }_1}{t}_1}\right)}{e}^{-{\left(\frac{1}{{ϵ}^{{\lambda }_2}{t}_2}\right)}^{k_2}u}{e}^{izm}dm\frac{du}{u},$$

$$ϵ\mapsto u_{h_1}(t_1,\frac{x_2}{{ϵ}^{{\mu }_2}}{e}^{\sqrt{-1}{\theta }_{h_1}},z,ϵ)$$

$$ϵ\mapsto u_{h_2}(\mathbf{t},z,ϵ),0\le h_2\le {\iota }_2-1,$$

$${\hat{\mathcal{B}}}_{q;1/k}(\hat{f}(T))(\tau )=\sum _{n\ge 0}{a}_n\frac{{\tau }^n}{(q^{1/k}{)}^{n(n-1)/2}}\in \mathbb{E}[[\tau ]].$$

$${\hat{\mathcal{B}}}_{q;1/k}(T^m{\sigma }_{q;T}^j\hat{f}(T))(\tau )=\frac{{\tau }^m}{(q^{1/k}{)}^{m(m-1)/2}}{\sigma }_{q;\tau }^{j-\frac{m}{k}}\left({\hat{\mathcal{B}}}_{q;1/k}(\hat{f}(T))(\tau )\right).$$

$${\Theta }_{q^{1/k}}(x)=\sum _{n\in \mathbb{Z}}{q}^{-\frac{n(n-1)}{2k}}{x}^n,x\in {\mathbb{C}}^{\star }.$$

$${\Theta }_{q^{1/k}}\left(q^{\frac{m}{k}}x\right)=q^{\frac{m(m+1)}{2k}}{x}^m{\Theta }_{q^{1/k}}(x),m\in \mathbb{Z},x\in {\mathbb{C}}^{\star }.$$

$$\left|{\Theta }_{q^{1/k}}(x)\right|\ge C_{q,k}\tilde{\delta }\exp \left(\frac{k}{2}\frac{{\log }^2|x|}{\log (q)}\right)|x{|}^{1/2},$$

$${\Vert f(\tau )\Vert }_{\mathbb{E}}\le K\exp \left(\frac{k}{2}\frac{{\log }^2|\tau |}{\log (q)}+\alpha \log |\tau |\right),\tau \in S_d,|\tau |\ge \rho ,$$

$${\Vert f(\tau )\Vert }_{\mathbb{E}}\le K,\tau \in \overline{D}(0,\rho ).$$

$${\mathcal{L}}_{q;1/k}^{\gamma }(f(\tau ))(T)=\frac{1}{{\pi }_{q^{1/k}}}{\int }_{L_{\gamma }}\frac{f(u)}{{\Theta }_{q^{1/k}}\left(\frac{u}{T}\right)}\frac{du}{u},$$

$${\mathcal{R}}_{\gamma ,\tilde{\delta }}=\left\{T\in {\mathbb{C}}^{\star }:\left|1+\frac{r{e}^{\gamma \sqrt{-1}}}{T}\right|>\tilde{\delta },for\; allr\ge 0\right\}.$$

$$T^{\sigma }{\sigma }_q^j({\mathcal{L}}_{q;1/k}^{\gamma }f(\tau ))(T)={\mathcal{L}}_{q;1/k}^{\gamma }\left(\frac{{\tau }^{\sigma }}{(q^{1/k}{)}^{\sigma (\sigma -1)/2}}{\sigma }_q^{j-\frac{\sigma }{k}}f(\tau )\right)(T),$$

$${\Vert h(m)\Vert }_{(\beta ,\mu )}=\underset{m\in \mathbb{R}}{\sup }(1+|m|{)}^{\mu }\exp (\beta |m|)|h(m)|<\infty .$$

$${\mathcal{F}}^{-1}(f)(x)=\frac{1}{(2\pi {)}^{1/2}}{\int }_{-\infty }^{\infty }f(m)\exp (ixm)dm,x\in \mathbb{R}.$$

$${\partial }_z{\mathcal{F}}^{-1}(f)(z)={\mathcal{F}}^{-1}(\varphi )(z),z\in H_{\beta }.$$

$$\psi (m)=\frac{1}{(2\pi {)}^{1/2}}{\int }_{-\infty }^{\infty }f(m-m_1)g(m_1)d{m}_1$$

$${\Vert f(\tau ,m)\Vert }_{(k,\beta ,\mu ,\alpha )}:=\underset{\begin{array}{c}\tau \in S_d\cup \overline{D}(0,\rho )\\ m\in \mathbb{R}\end{array}}{\sup }(1+|m|{)}^{\mu }{e}^{\beta |m|}\exp \left(\frac{-k{\log }^2(|\tau |+\delta )}{2\log (q)}-\alpha \log (|\tau |+\delta )\right)|f(\tau ,m)|$$

$${\Vert a(\tau ,m)f(\tau ,m)\Vert }_{(k,\beta ,\mu ,\alpha )}\le C_1{\Vert f(\tau ,m)\Vert }_{(k,\beta ,\mu ,\alpha )},$$

$${\gamma }_2\le k{\gamma }_3+{\gamma }_1.$$

$${\Vert a_{{\gamma }_1}(\tau ){\tau }^{{\gamma }_2}{\sigma }_{q;\tau }^{-{\gamma }_3}f(\tau ,m)\Vert }_{(k,\beta ,\mu ,\alpha )}\le C_2{\Vert f(\tau ,m)\Vert }_{(k,\beta ,\mu ,\alpha )},$$

$$\begin{gathered} {\Vert a_{{\gamma }_1}(\tau ){\tau }^{{\gamma }_2}{\sigma }_{q;\tau }^{-{\gamma }_3}f(\tau ,m)\Vert }_{(k,\beta ,\mu ,\alpha )} \\ =\underset{\begin{array}{c}\tau \in S_d\cup \overline{D}(0,\rho )\\ m\in \mathbb{R}\end{array}}{\sup }{C}_{21}(|\tau |)(1+|m|{)}^{\mu }{e}^{\beta |m|}\exp \left(\frac{-k{\log }^2(\frac{|\tau |}{q^{{\gamma }_3}}+\delta )}{2\log (q)}-\alpha \log \left(\frac{|\tau |}{q^{{\gamma }_3}}+\delta \right)\right)\left|f\left(\frac{\tau }{q^{{\gamma }_3}},m\right)\right| \\ \le {\Vert f(\tau ,m)\Vert }_{(k,\beta ,\mu ,\alpha )}\underset{\tau \in S_d\cup \overline{D}(0,\rho )}{\sup }{C}_{21}(|\tau |) \end{gathered}$$

$$\begin{gathered} C_{21}(|\tau |)=\frac{|\tau {|}^{{\gamma }_2}}{(1+|\tau |{)}^{{\gamma }_1}}\exp \left(\frac{-k\left({\log }^2(|\tau |+\delta )-{\log }^2(\frac{|\tau |}{q^{{\gamma }_3}}+\delta )\right)}{2\log (q)}\right) \\ \times \exp \left(-\alpha \log (|\tau |+\delta )+\alpha \log \left(\frac{|\tau |}{q^{{\gamma }_3}}+\delta \right)\right). \end{gathered}$$

$${\log }^2$$

$$={\log }^2(|\tau |+\delta q^{{\gamma }_3})-{\log }^2(|\tau |+\delta )+(\log (q){\gamma }_3{)}^2-2{\gamma }_3\log (q)\log (|\tau |+\delta q^{{\gamma }_3}).$$

$${\log }^2(|\tau |+\delta q^{{\gamma }_3})-{\log }^2(|\tau |+\delta )\le 2\delta (q^{{\gamma }_3}-1)\underset{x\in [|\tau |+\delta ,|\tau |+\delta q^{{\gamma }_3}]}{\max }\frac{\log (x)}{x}\le C_{22},$$

$$C_{21}(|\tau |)\le \exp \left(\frac{k{C}_{22}}{2\log (q)}+(\log (q){\gamma }_3{)}^2\right)\frac{|\tau {|}^{{\gamma }_2}}{(1+|\tau |{)}^{{\gamma }_1}}{\left(\frac{|\tau |/q^{{\gamma }_3}+\delta }{|\tau |+\delta }\right)}^{\alpha }\frac{1}{(|\tau |+\delta q^{{\gamma }_3}{)}^{k{\gamma }_3}},$$

$$\Phi (\tau ,m):=\frac{1}{R_1(im)}{\int }_{-\infty }^{\infty }f(m-m_1)R_2(i{m}_1)g(\tau ,m_1)d{m}_1$$

Full text

On a $q-$analog of a singularly perturbed problem of irregular type with two complex time variables

**A. Lastra111The author is partially supported by the project MTM2016-77642-C2-1-P of Ministerio de Economía y Competitividad, Spain, S. Malek222The author is partially supported by the project MTM2016-77642-C2-1-P of Ministerio de Economía y Competitividad, Spain.

**University of Alcalá, Departamento de Física y Matemáticas,

Ap. de Correos 20, E-28871 Alcalá de Henares (Madrid), Spain,

University of Lille 1, Laboratoire Paul Painlevé,

59655 Villeneuve d’Ascq cedex, France,

alberto.lastra@uah.es

Stephane.Malek@math.univ-lille1.fr

Abstract

Analytic solutions and their formal asymptotic expansions for a family of the singularly perturbed $q-$difference-differential equations in the complex domain are constructed. They stand for a $q-$analog of the singularly perturbed partial differential equations considered in [15]. In the present work, we construct outer and inner analytic solutions of the main equation, each of them showing asymptotic expansions of essentially different nature with respect to the perturbation parameter. The appearance of the $-1$-branch of Lambert $W$ function will be crucial in this respect.

Key words: asymptotic expansion, Borel-Laplace transform, Fourier transform, initial value problem, formal power series, $q-$difference equation, boundary layer, singular perturbation. 2010 MSC: 35C10, 35C20,35R10, 35C15

1 Introduction

The intent of this study is to provide analytic solutions and their parametric asymptotic expansions to a family of singularly perturbed $q-$difference-differential equations in the complex domain, in which two time variables act, for some fixed $q>1$.

More precisely, we consider equations of the form

[TABLE]

under null initial data $u(t_{1},0,z,\epsilon)\equiv u(0,t_{2},z,\epsilon)\equiv 0$. Here $Q(X)\in\mathbb{C}[X]$, and $P$ stands for a polynomial with complex coefficients with respect to $\boldsymbol{t}:=(t_{1},t_{2}),\partial_{t_{2}}$, a polynomial with rational powers with respect to $\sigma_{q;t_{1}}$ and $\sigma_{q;t_{2}}$; holomorphic with respect to $z$ on a horizontal strip $H_{\beta}:=\{z\in\mathbb{C}:|\hbox{Im}(z)|<\beta\}$ for some $\beta>0$, and holomorphic with respect to the perturbation parameter $\epsilon$ on a small disc centered at the origin, say $D(0,\epsilon_{0})$ for some $\epsilon_{0}>0$. Throughout the present work, $\sigma_{q;t}$ stands for the dilation operator on $t$ variable for some fixed $q>1$, i.e.

[TABLE]

We adopt the following notation $\sigma_{q;t}^{\delta}(f(t))=f(q^{\delta}t)$, for any $\delta\in\mathbb{Q}$.

The forcing term $f(\boldsymbol{t},z,\epsilon)$ is constructed under certain growth conditions, and turns out to be a holomorphic function in $\mathbb{C}^{\star}\times\mathbb{C}^{\star}\times H_{\beta^{\prime}}\times(D(0,\epsilon_{0})\setminus\{0\})$, for $0<\beta^{\prime}<\beta$.

The precise assumptions on the elements involved in the main equation under study are detailed in Section 3.

The problem under study (1) turns out to be a $q-$analog of the main problem studied in [15],

[TABLE]

under null initial data $u(t_{1},0,z,\epsilon)\equiv u(0,t_{2},z,\epsilon)\equiv 0$, and where $Q(X)\in\mathbb{C}[X]$, and $P$ is a polynomial with respect to its first three variables, with holomorphic coefficients on $H_{\beta}\times D(0,\epsilon_{0})$, and the forcing term is holomorphic on $\mathbb{C}^{\star}\times\mathbb{C}\times H_{\beta^{\prime}}\times(D(0,\epsilon_{0})\setminus\{0\})$. The analytic solutions and their asymptotic expansions of those singularly perturbed partial differential equations are obtained in [15]. More precisely, the so-called inner solutions are holomorphic solutions of (2), holomorphic on domains in time which depend on the perturbation parameter and approach infinity, admits Gevrey asymptotic expansion of certain positive order, with respect to $\epsilon$, whereas the so-called outer solutions are holomorphic solutions of (2), holomorphic on a product of finite sectors with vertex at the origin with respect to the time variables, admit Gevrey asymptotic expansion of a different positive order, with respect to $\epsilon$.

The previous phenomena of existence of different asymptotic expansions regarding different domains of the actual solution of the main problem is enhanced in the present work, in the sense that different nature on the asymptotic expansions are observed: Gevrey and $q-$Gevrey asymptotic expansions.

In this work, we search for the analytic solutions of the main problem as the inverse Fourier transform and $q-$Laplace transform of a positive order in the form

[TABLE]

(see (21) and (23)), for some appropriate $\gamma\in\mathbb{R}$ and $k_{1},k_{2},\lambda_{1},\lambda_{2}>0$ (see Section 3 for their definition). The function $\omega(\tau,m,\epsilon)$ is obtained from a fixed point argument (see Proposition 6), and belongs to $q\hbox{Exp}^{d}_{(k^{\prime}_{1},\beta,\mu,\alpha)}$, a Banach space of holomorphic functions with $q-$exponential growth and exponential decay with respect to $\tau$ and $m$, respectively (see Definition 4).

The form (3) of the analytic solutions is motivated on the shape of those of the main problem in [15], mixing both time variables in a common Laplace operator.

A first family of analytic solutions of (1) is constructed on domains of the form $\mathcal{T}_{1}\times\mathcal{T}_{2,\epsilon}\times H_{\beta^{\prime}}\times\mathcal{E}_{h_{1}}^{\infty}$, and a second on domains of the form $\mathcal{T}_{1}\times(\mathcal{T}_{2}\cap D(0,\rho_{2}))\times H_{\beta^{\prime}}\times\mathcal{E}^{0}_{h_{2}}$, where $\mathcal{T}_{1}$ is a finite sector, $\mathcal{T}_{2}$ is an unbounded sector and where $\mathcal{T}_{2,\epsilon}\subseteq\mathcal{T}_{2}$ is a bounded sector which depends on $\epsilon\in\mathcal{E}_{h_{1}}^{\infty}$, and tends to infinity with $\epsilon$ approaching the origin. The sets $\underline{\mathcal{E}}^{0}=(\mathcal{E}_{h_{2}}^{0})_{0\leq h_{2}\leq\iota_{2}-1}$ and $\overline{\mathcal{E}}^{\infty}=(\mathcal{E}_{h_{1}}^{\infty})_{0\leq h_{1}\leq\iota_{1}-1}$ represent good coverings (see Definition 5).

Different path deformations performed on the analytic solutions give rise to Theorem 2 and Theorem 3, where upper bounds on the difference of two consecutive solutions are attained (consecutive solutions in the sense that they are related to consecutive sectors in a good covering). Such bounds are related to null Gevrey and $q-$Gevrey asymptotic expansions of some positive order. As a matter of fact, the previous differences allow to apply a novel $((q,k);s)-$version of the cohomological criteria known as Ramis-Sibuya theorem. Such result is related to functions admitting $q-$Gevrey asymptotic expansions of order $k$ and a Gevrey sub-level of order $s$, see Theorem 4. We also apply a $q-$analog of Ramis-Sibuya Theorem, see Theorem 5.

The main two results of the present work are Theorem 6 and Theorem 7 relating the analytic solutions of (1) to their formal power series expansions obtaining asymptotic results of different nature. Such solutions are known as inner and outer solutions (see Definition 8 and Definition 9, resp.). Such asymptotic solutions have also been observed in the previous study [15], in the framework of singularly perturbed PDEs. However, the different nature of the asymptotic expansions regarding the outer and inner solutions is a novel phenomena which has firstly been observed in the present study.

The inner and outer expansions appear in the study of matched asymptotic expansions (see [22, 28], among others for the classical theory). In the work [6] by A. Fruchard and R. Schäfke, the method of matching is developed, studying the nature of such asymptotic expansions, under Gevrey settings.

We fix a good covering $(\mathcal{E}^{\infty}_{h_{1}})_{0\leq h_{1}\leq\iota_{1}-1}$, and consider the holomorphic solutions of the problem (1), $u_{h_{1}}(\boldsymbol{t},z,\epsilon)$ defined on $\mathcal{E}^{\infty}_{h_{1}}$ w.r.t. $\epsilon$ for all $0\leq h_{1}\leq\iota_{1}-1$. In Theorem 6, we prove that for some $\mu_{2}>\lambda_{2}$, $\theta_{h_{1}}$ and some adequate domain $\chi_{2}^{\infty}\subset\mathcal{T}_{2}$, the function

[TABLE]

with values in the Banach space of holomorphic and bounded functions on $\mathcal{T}_{1}\times\chi_{2}^{\infty}\times H_{\beta^{\prime}}$, say $\mathbb{F}_{1}$, admits a formal power series $\hat{u}^{\infty}(\epsilon)\in\mathbb{F}_{1}[[\epsilon]]$ as $((q,K);S)-$Gevrey asymptotic expansion on $\mathcal{E}_{h_{1}}^{\infty}$ for some $K,S>0$, for every $0\leq h_{1}\leq\iota_{1}-1$.

The proof of this result leans on the application of accurate estimates related to the $-1$-branch of Lambert $W$ function (see Lemma 3).

Concerning the outer solutions of the main problem under study, we consider a good covering $(\mathcal{E}^{0}_{h_{2}})_{0\leq h_{2}\leq\iota_{2}-1}$, and the holomorphic solutions of (1), $u_{h_{2}}(\boldsymbol{t},z,\epsilon)$ defined on $\mathcal{E}^{0}_{h_{2}}$ w.r.t. $\epsilon$ for all $0\leq h_{2}\leq\iota_{2}-1$. In Theorem 7, we prove that

[TABLE]

is an outer solution of (1) with values in the Banach space $\mathbb{F}_{2}$ of holomorphic and bounded functions on $\mathcal{T}_{1}\times(\mathcal{T}_{2}\cap D(0,\rho_{2}))\times H_{\beta^{\prime}}$. Moreover, there exists a formal power series $\hat{u}^{0}(\epsilon)\in\mathbb{F}_{2}[[\epsilon]]$ which is the common $q-$Gevrey asymptotic expansion of some positive order of each solution (5) on $\mathcal{E}_{h_{2}}^{0}$, for $0\leq h_{2}\leq\iota_{2}-1$.

In recent years, an increasing interest on the study of the asymptotic behavior of solutions to $q-$difference-differential equations in the complex domain has been observed. New theories giving rise to $q-$analogs of the classical theory of Borel-Laplace summability have been discussed and studied, as in the case of the work [29], by H. Tahara, where the author also provides information about $q-$analogs of Borel and Laplace transforms and related properties on convolution or Watson-type results. The use of procedures based on the Newton polygon is also exploited in recent studies, such as the work [30] by H. Tahara and H. Yamazawa. Also, it is worth mentioning the study of $q-$analogs of Briot-Bouquet type partial differential equations by H. Yamazawa, in [32]. Integral transforms involving special functions have also been considered in the study of $q-$difference-differential equations in [7, 23]. Other references in this context by the authors and collaborators are listed in the references.

Different kinds of Advanced/delayed partial differential equations are the cornerstone of mathematical models which have been recently studied. Examples of such studies have been applied to tsunamis and rogue waves which can be found in [26]. We also refer to other studies such as [24, 25], and the references therein.

The outline of the work is as follows. In Section 2, we recall some known facts about formal $q-$Borel transform, analytic $q-$Laplace transform and inverse Fourier transform together with some properties which are applied to transform the main equation under study into auxiliary problems. Afterwards, we provide the definition and related properties of some Banach spaces involved in the construction of the solution. In Section 3 we state the main problem under study (17) and two auxiliary equations. The elements involved in them in addition to the domains of existence and upper bounds of the solutions of such equations are detailed. Section 4 is devoted to the existence and description of the domain of existence for the auxiliary equation (22), and associated estimates. In the following section, Section 5, we provide analytic solutions of (17) (see Theorem 1) and estimates on the difference of two of them (see Theorem 2 and Theorem 3). After a brief summary on $q-$asymptotic expansions in the first part of Section 6, we provide formal power series expansions in the perturbation parameter of the analytic solutions and relate them asymptotically in adequate domains. These results are attained in Theorem 6 and Theorem 7. The work concludes with two technical sections, Section 7 and Section 8, left to the end of the work in order not to interfere with our reasonings.

2 Review on certain integral operators. Study of Banach spaces involved in the problem

We briefly describe the foundations of the analytic and formal operators that will allow the transformation of the main problem under study in terms of auxiliary equations. The solutions of such equations belong to certain Banach spaces which are constructed subsequently.

2.1 Review of some formal and analytic operators

This section is devoted to recall some of the basic facts on formal and analytic transformations corresponding to the $q-$analogs of those appearing in the classical Borel-Laplace summability theory. These tools were developed in [4, 27].

Through the whole section, $q>1$ stands for a real number, and $k\geq 1$ is a positive integer. $\mathbb{E}$ stands for a complex Banach space.

Definition 1

Given $\hat{f}(T)=\sum_{n\geq 0}a_{n}T^{n}\in\mathbb{E}[[T]]$, the formal $q-$Borel transform of order $k$ of $\hat{f}(T)$ is defined by

[TABLE]

The next result, whose proof can be found in Proposition 5 [19], is crucial in order to transform the main equation into an auxiliary one lying in the $q-$Borel plane.

Proposition 1

Let $m\in\mathbb{N}$ and $j\in\mathbb{Q}$. For every $\hat{f}(T)\in\mathbb{E}[[T]]$, it holds that

[TABLE]

The $q-$analog of Laplace transform used in this work was introduced in [31], and makes use of a kernel given by the Jacobi Theta function of order $k$, given by

[TABLE]

We recall that Jacobi Theta function satisfies

[TABLE]

From Lemma 4.1 [11], given any $\tilde{\delta}>0$, there exists $C_{k,q}>0$ (independent of $\tilde{\delta}$) such that

[TABLE]

for all $x\in\mathbb{C}^{\star}$ such that $|1+xq^{\frac{m}{k}}|>\tilde{\delta}$, for every $m\in\mathbb{Z}$. This last property allows to define a $q-$analog of Laplace transform with appropriate properties.

Definition 2

Let $\rho>0$ and $S_{d}$ be an unbounded sector with vertex at 0, and bisecting direction $d\in\mathbb{R}$. Let $f:S_{d}\cup D(0,\rho)\to\mathbb{E}$ be a holomorphic function, continuous up to the boundary, such that there exist $K>0$, and $\alpha\in\mathbb{R}$ with

[TABLE]

[TABLE]

We choose an argument $\gamma\in\mathbb{R}$ within the set of arguments in $S_{d}$ and define the $q-$Laplace transform of order $k$ of $f$ in direction $\gamma$ by

[TABLE]

where $L_{\gamma}=\{re^{\gamma\sqrt{-1}}:r\in(0,\infty)\}$, and $\pi_{q^{1/k}}:=\frac{\log(q)}{k}\prod_{n\geq 0}(1-q^{-\frac{n+1}{k}})^{-1}$.

The proof of the next results can be found in detail in Lemma 4 and Proposition 6 [19].

Lemma 1

Let $\tilde{\delta}>0$. In the situation of Definition 2, the integral transform $\mathcal{L}^{\gamma}_{q;1/k}(f(\tau))(T)$ defines a bounded holomorphic function on the domain $\mathcal{R}_{\gamma,\tilde{\delta}}\cap D(0,r_{1})$ for every $0<r_{1}\leq q^{(1/2-\alpha)/k}/2$, where

[TABLE]

The value of $\mathcal{L}^{\gamma}_{q;1/k}(f(\tau))(T)$ does not depend on the choice of $\gamma$ under $e^{\gamma\sqrt{-1}}\in S_{d}$.

Proposition 2

Let $f$ be as in Definition 2, and let $\tilde{\delta}>0$. Then, for all $\sigma\geq 0$ one has

[TABLE]

for all $T\in\mathcal{R}_{\gamma,\tilde{\delta}}\cap D(0,r_{1})$, where $0<r_{1}<q^{(\frac{1}{2}-\alpha)/k}/2$.

We conclude with the definition and properties regarding the inverse Fourier transform.

Definition 3

Let $\mu,\beta\in\mathbb{R}$. We write $E_{(\beta,\mu)}$ for the vector space of continuous functions $h:\mathbb{R}\to\mathbb{C}$ such that

[TABLE]

The pair $(E_{(\beta,\mu)},\left\|\cdot\right\|_{(\beta,\mu)})$ turns out to be a Banach space.

Proposition 3

Let $\mu>1,\beta>0$ and $f\in E_{(\beta,\mu)}$. The inverse Fourier transform of $f$ is defined by

[TABLE]

This function can be extended to an analytic function on the horizontal strip $H_{\beta}=\{z\in\mathbb{C}:|\hbox{Im}(z)|<\beta\}$. Moreover, it holds that the function $\phi(m)=imf(m)$ belongs to $E_{(\beta,\mu-1)}$ and

[TABLE]

In addition to this, it holds that the convolution product of $f\in E_{(\beta,\mu)}$ and $g\in E_{(\beta,\mu)}$, defined by

[TABLE]

is such that $\psi\in E_{(\beta,m)}$, and $\mathcal{F}^{-1}(f)(z)\mathcal{F}^{-1}(g)(z)=\mathcal{F}^{-1}(\psi)(z),$ for every $z\in H_{\beta}$.

2.2 Banach spaces of functions of $q-$exponential growth and exponential decay

In this section, we state the definition of the complex Banach space $qExp_{(k,\beta,\mu,\alpha)}^{d}$. The analytic solution of the main equation under study is built departing from one element in such Banach space via $q-$Laplace transform and inverse Fourier transform. Similar versions of this Banach space have already appeared in previous works by the authors such as [11, 12], and the contribution [19] of the second author.

In the whole section, we fix real numbers $\beta,\mu,k>0$, $q,\delta>1$ and $\alpha$. We also set $d\in\mathbb{R}$ and choose an infinite sector $S_{d}$ of bisecting direction $d$, with vertex at the origin; and the closed disc $\overline{D}(0,\rho)$, for some $\rho>0$.

Definition 4

Let $qExp_{(k,\beta,\mu,\alpha)}^{d}$ be the vector space of all continuous complex valued functions $(\tau,m)\mapsto f(\tau,m)$ on $(S_{d}\cup\overline{D}(0,\rho))\times\mathbb{R}$ which are holomorphic with respect to the first variable on $S_{d}\cup D(0,\rho)$ and such that

[TABLE]

is finite. The pair $(qExp_{(k,\beta,\mu,\alpha)}^{d},\left\|\cdot\right\|_{(k,\beta,\mu,\alpha)})$ is a Banach space.

In the forthcoming results, we describe some properties on the elements of the previous Banach space when combined with certain operators. The first result is a direct consequence of the previous definition.

Lemma 2

Let $(\tau,m)\mapsto a(\tau,m)$ be a bounded continuous function defined on $(S_{d}\cup\overline{D}(0,\rho))\times\mathbb{R}$, holomorphic with respect to $\tau$ on the set $S_{d}\cup D(0,\rho)$. For every $f(\tau,m)\in qExp_{(k,\beta,\mu,\alpha)}^{d}$ one has that $a(\tau,m)f(\tau,m)$ belongs to $qExp_{(k,\beta,\mu,\alpha)}^{d}$, and it holds that

[TABLE]

with $C_{1}:=\sup_{\begin{subarray}{c}\tau\in S_{d}\cup\overline{D}(0,\rho)\\ m\in\mathbb{R}\end{subarray}}|a(\tau,m)|$.

Proposition 4

Let $\gamma_{1},\gamma_{2},\gamma_{3}\geq 0$ such that

[TABLE]

Let $a_{\gamma_{1}}(\tau)$ be a bounded holomorphic function defined on $S_{d}\cup D(0,\rho)$ with $(1+|\tau|)^{\gamma_{1}}|a_{\gamma_{1}}(\tau)|\leq 1$ for every $\tau\in S_{d}\cup D(0,\rho)$. Then, for all $f\in qExp_{(k,\beta,\mu,\alpha)}^{d}$ it holds that

[TABLE]

for some $C_{2}>0$.

Proof Let $f\in qExp_{(k,\beta,\mu,\alpha)}^{d}$. Then, it holds that

[TABLE]

with

[TABLE]

Observe that $\sigma_{q;\tau}^{-\gamma_{3}}(S_{d}\cup\overline{D}(0,\rho))\subseteq S_{d}\cup\overline{D}(0,\rho)$. In order to give upper bounds for $C_{21}(|\tau|)$, we have

[TABLE]

The mean value theorem guarantees that

[TABLE]

for some positive constant $C_{22}$ which does not depend on $\tau\in S_{d}\cup\overline{D}(0,\rho)$. This entails that

[TABLE]

for every $\tau\in S_{d}\cup\overline{D}(0,\rho)$. We conclude the result by taking into account (7). $\Box$

The next result is stated without proof, heavily rests on Proposition 2 [16], and it is based on bounds stated in Lemma 2.2 [3].

Proposition 5

Let $R_{1},R_{2}\in\mathbb{C}[X]$ with $\hbox{deg}(R_{1})\geq\hbox{deg}(R_{2})$, and such that $R_{1}(im)\neq 0$ for all $m\in\mathbb{R}$. Let $\mu>\hbox{deg}(R_{2})+1$. Given $f\in E_{(\beta,\mu)}$ and $g\in q\hbox{Exp}^{d}_{(k,\beta,\mu,\alpha)}$, then the function

[TABLE]

belongs to $q\hbox{Exp}^{d}_{(k,\beta,\mu,\alpha)}$, and there exists $C^{\prime}_{2}>0$ such that

[TABLE]

3 Statement of the main problem and auxiliary equations

Let $\beta>0$ and $k^{\prime}_{1},k_{1},k_{2}\geq 1$ with $k^{\prime}_{1}>k_{1}$, and $D_{1},D_{2}\geq 2$ be integer numbers. We also consider a real number $q>1$, and choose real numbers $\Delta_{D_{1}D_{2}},d_{D_{1}},\tilde{\delta}_{D_{2}}\geq 0$, and $\Delta_{\ell_{1}\ell_{2}},\delta_{\ell_{1}},d_{\ell_{1}},\tilde{\delta}_{\ell_{2}}\geq 0$ for every $1\leq\ell_{1}\leq D_{1}-1$ and $1\leq\ell_{2}\leq D_{2}-1$. We assume that

[TABLE]

In addition to that, there exist natural numbers $\lambda_{1},\lambda_{2}$ with

[TABLE]

[TABLE]

Let us also fix polynomials $Q,R_{D_{1}D_{2}}$, and $R_{\ell_{1}\ell_{2}}$ for all $1\leq\ell_{1}\leq D_{1}-1$ and $1\leq\ell_{2}\leq D_{2}-1$, with complex coefficients such that

[TABLE]

where

[TABLE]

for some $r_{Q,R_{D_{1}D_{2}}},\eta_{Q,R_{D_{1}D_{2}}}>0$ and $d_{Q,R_{D_{1}D_{2}}}\in\mathbb{R}$. Observe this condition implies $\hbox{deg}(R_{D_{1}D_{2}})\leq\hbox{deg}(Q)$. These polynomials are chosen in such a way that $R_{D_{1}D_{2}}(im)\neq 0$ for $m\in\mathbb{R}$, and

[TABLE]

Let $\mu>1$ such that

[TABLE]

The main problem under study in this work is

[TABLE]

under initial data $u(t_{1},0,z,\epsilon)\equiv u(0,t_{2},z,\epsilon)\equiv 0$.

The function $f$ is a holomorphic function in $\mathbb{C}^{\star}\times\mathbb{C}^{\star}\times H_{\beta^{\prime}}\times(D(0,\epsilon_{0})\setminus\{0\})$ for every $0<\beta^{\prime}<\beta$, where $H_{\beta^{\prime}}$ stands for the horizontal strip $H_{\beta^{\prime}}=\{z\in\mathbb{C}:|\hbox{Im}(z)|<\beta^{\prime}\}$. It is constructed as follows. Let $\psi(\tau,m,\epsilon)$ be a continuous function, continuous in $\mathbb{C}\times\mathbb{R}\times D(0,\epsilon_{0})$, for some $\epsilon_{0}>0$, entire with respect to the first variable, and holomorphic with respect to the third one in $D(0,\epsilon_{0})$. We moreover assume there exists $C_{\psi}>0$ such that

[TABLE]

for every $(\tau,m,\epsilon)\in\mathbb{C}\times\mathbb{R}\times D(0,\epsilon_{0})$, some $\alpha\in\mathbb{R}$, $\delta>0$; and where $\mu$ satisfies (16).

In view of the definition of $q-$Laplace transform and the results described in Section 2.1, one can define

[TABLE]

for some fixed $\gamma\in\mathbb{R}$. Here, $\boldsymbol{T}:=(T_{1},T_{2})$. The function

[TABLE]

turns out to be holomorphic on the set $\mathbb{C}^{\star}\times\mathbb{C}^{\star}\times H_{\beta^{\prime}}\times(D(0,\epsilon_{0})\setminus\{0\})$. Observe that in the previous construction, given $T_{1},T_{2}\in\mathbb{C}^{\star}$, one can choose $\gamma\in\mathbb{R}$ with $\cos(\gamma-k_{2}\hbox{arg}(T_{2}))>0$ and $|1+(uq^{m/k})/T_{1}|>0$ for all $m\in\mathbb{Z}$ and with $u\in L_{\gamma}$.

For every $0\leq\ell_{1}\leq D_{1}-1$ and $0\leq\ell_{2}\leq D_{2}-1$, the function $c_{\ell_{1}\ell_{2}}(z,\epsilon)$ is constructed in the following way:

[TABLE]

$c_{\ell_{1}\ell_{2}}(z,\epsilon)$ turns out to be a holomorphic function on $H_{\beta}\times D(0,\epsilon_{0})$ whenever $m\mapsto C_{\ell_{1}\ell_{2}}(m,\epsilon)\in E_{(\beta,\mu)}$. In addition to that, we assume that uniform bounds with respect to the perturbation parameter are satisfied, i.e. there exist $\mathcal{C}_{\ell_{1}\ell_{2}}>0$ with

[TABLE]

3.1 Study of auxiliary equations

In this section we preserve the statements and constructions concerning the main problem under study (17), and the geometric and algebraic conditions held on the elements involved in the main equation.

We search for solutions of (17) of the form

[TABLE]

Assuming the solution is of the form (21), the expression $U(T_{1},T_{2},m,\epsilon)$ solves

[TABLE]

We reduce the study of solutions of (17) to those of (22), which are linked through (21). In order to solve (22), we adapt a recent approach developed in  [15] to a new situation involving both partial differential and $q-$difference operators. We seek for solutions of (22) of the special form

[TABLE]

for some appropriate function $\omega(\tau,m,\epsilon)$ and $\gamma\in\mathbb{R}$. We refer to Section 2.1 for the definitions of the elements involved in the previous expression.

Let us consider a second auxiliary equation:

[TABLE]

We define the polynomial $P_{m}(\tau)$ by

[TABLE]

whose factorization is given by

[TABLE]

with

[TABLE]

for all $0\leq\ell\leq d_{D_{1}}+\tilde{\delta}_{D_{2}}-1$, where $\Delta_{m}=\left(\frac{|Q(im)|}{k_{2}^{\tilde{\delta}_{D_{2}}}|R_{D_{1}D_{2}}(im)|}(q^{1/k_{1}})^{\frac{d_{D_{1}}(d_{D_{1}}-1)}{2}}\right)^{1/(d_{D_{1}}+\tilde{\delta}_{D_{2}})}$.

Let $d\in\mathbb{R}$ be such that the infinite sector $S_{d}$ of bisecting direction $d$ satisfies the following geometric construction: there exists $M_{1}>0$ such that $|\tau-q_{\ell}(m)|\geq M_{1}(1+|\tau|)$, for all $0\leq\ell\leq d_{D_{1}}+\tilde{\delta}_{D_{2}}-1$, $m\in\mathbb{R}$ and $\tau\in S_{d}\cup\overline{D}(0,\rho)$. The previous is a feasible condition for an appropriate choice of small enough $\eta_{Q,R_{D_{1}D_{2}}}>0$ and large enough $r_{Q,R_{D_{1}D_{2}}}>0$, in view of (14) and the definition of $q_{\ell}(m)$. More precisely, one chooses $S_{d}$ such that $q_{\ell}(m)/\tau$ has positive distance to 1, for all $0\leq\ell\leq d_{D_{1}}+\tilde{\delta}_{D_{2}}-1$, $m\in\mathbb{R}$, and $\tau\in S_{d}\cup\overline{D}(0,\rho)$. The previous choice of $d$ yields

[TABLE]

for some $C_{P}>0$, valid for all $\tau\in S_{d}\cup\overline{D}(0,\rho)$ and all $m\in\mathbb{R}$.

Proposition 6

Let $S_{d}$ be an infinite sector with vertex at the origin satisfying the previous geometric conditions. Then, if $C_{\psi}>0$ (see (18) for its definition) and $\epsilon_{0}>0$ are small enough, there exists $\varpi>0$ such that the equation (24) admits a unique solution $\omega^{d}(\tau,m,\epsilon)$ which belongs to $qExp_{(k^{\prime}_{1},\beta,\mu,\alpha)}^{d}$ with $\left\|\omega^{d}(\tau,m,\epsilon)\right\|_{(k^{\prime}_{1},\beta,\mu,\alpha)}\leq\varpi$, for all $\epsilon\in D(0,\epsilon_{0})$.

Proof Let $\epsilon\in D(0,\epsilon_{0})$. We consider the map $\mathcal{H}_{\epsilon}$ defined by

[TABLE]

Given $\varpi>0$, let $\omega(\tau,m)\in qExp_{(k^{\prime}_{1},\beta,\mu,\alpha)}^{d}$ with $\left\|\omega^{d}(\tau,m,\epsilon)\right\|_{(k^{\prime}_{1},\beta,\mu,\alpha)}\leq\varpi$. In view of (11), (12), (13), (26) and from Lemma 2, Proposition 4 and Proposition 5 one has

[TABLE]

for every $0\leq\ell_{1}\leq D_{1}-1$ and $0\leq\ell_{2}\leq D_{2}-1$. In addition to this, Lemma 2 and (26) yield

[TABLE]

Let $\epsilon_{0},C_{\psi}>0$ be small enough satisfying

[TABLE]

Then, the estimates (27) and (28) yield to $\left\|\mathcal{H}(\omega(\tau,m)\right\|_{(k^{\prime}_{1},\beta,\mu,\alpha)}\leq\varpi$. In other words, the operator $\mathcal{H}_{\epsilon}$ restricted to $\overline{D}(0,\varpi)\subseteq qExp_{(k^{\prime}_{1},\beta,\mu,\alpha)}^{d}$ is such that $\mathcal{H}_{\epsilon}(\overline{D}(0,\varpi))\subseteq\overline{D}(0,\varpi)$.

On the other hand, let $\omega_{1}(\tau,m),\omega_{2}(\tau,m)\in\overline{D}(0,\rho)\subseteq qExp_{(k^{\prime}_{1},\beta,\mu,\alpha)}^{d}$. Analogously to (27), one arrives at

[TABLE]

by choosing $\epsilon_{0}>0$ such that

[TABLE]

We conclude that the map $\mathcal{H}_{\epsilon}:\overline{D}(0,\varpi)\subseteq qExp_{(k^{\prime}_{1},\beta,\mu,\alpha)}^{d}\to qExp_{(k^{\prime}_{1},\beta,\mu,\alpha)}^{d}$ is contractive. The classical fixed point theory in complete metric spaces states the existence of a unique fixed point for $\mathcal{H}_{\epsilon}$, say $\omega^{d}(\tau,m,\epsilon)$, in $qExp_{(k^{\prime}_{1},\beta,\mu,\alpha)}^{d}$, with $\left\|\omega^{d}(\tau,m,\epsilon)\right\|_{(k^{\prime}_{1},\beta,\mu,\alpha)}\leq\varpi$. For every $\epsilon\in D(0,\epsilon_{0})$, the function $\omega^{d}(\tau,m,\epsilon)$ is a solution of (24) in view of the definition of the operator $\mathcal{H}_{\epsilon}$. Holomorphy of the map $D(0,\epsilon_{0})\ni\epsilon\mapsto\omega^{d}(\tau,m,\epsilon)$ is derived from the construction of the fixed point. $\Box$

In order to prove that the solutions of (22) and (24) are related via (23), we need to clarify how operators involved in (22) are transformed into the corresponding ones in (24). This is left to the end of the work, in Section 8 not to interfere our line of reasoning.

4 Domains of existence for the solutions of (22), and associated estimates

In this section, we describe appropriate domains on the time variables in which the solution of the main problem under consideration is well defined, within appropriate geometric conditions. Let $\tilde{\mathcal{T}}_{1}$ be a bounded sector with vertex at the origin such that there exists $\delta_{1}>0$ with

[TABLE]

for all $r\geq 0$ and all $T_{1}\in\tilde{\mathcal{T}}_{1}$, and $\gamma$ being an argument in $S_{d}$. We also fix an unbounded sector $\tilde{\mathcal{T}}_{2}$, with vertex at the origin, such that

[TABLE]

for some $\delta_{2}>0$ and well chosen $\gamma\in\hbox{arg}(S_{d})$. Observe that, in particular, there exists $\delta_{3}>0$ with $\cos(\gamma-k_{2}\hbox{arg}(T_{2}))>\delta_{3}$.

The next technical result describes accurate bounds for the solutions of (22) in different domains. The proof is left to Section 7.

Proposition 7

Let $U_{\gamma}(\boldsymbol{T},m,\epsilon)$ be defined in (23), with $\omega(\tau,m,\epsilon)=\omega^{d}(\tau,m,\epsilon)$ being the function obtained in Proposition 6. The following statements hold:

There exist small enough $\rho_{1}>0$ and large enough $\rho_{2}^{\infty}>0$ such that

[TABLE]

for every $T_{1}\in\tilde{\mathcal{T}}_{1}$ with $|T_{1}|<\rho_{1}$ and $T_{2}\in\tilde{\mathcal{T}}_{2}$ with $|T_{2}|>\rho^{\infty}_{2}$, and some $\tilde{C}_{1}>0$. 2. 2.

There exist small enough $\rho_{1},\rho_{2}>0$ such that

[TABLE]

for every $T_{1}\in\tilde{\mathcal{T}}_{1}$ with $|T_{1}|<\rho_{1}$ and $T_{2}\in\tilde{\mathcal{T}}_{2}$ with $|T_{2}|<\rho_{2}$, and some $\tilde{C}_{2}>0$.

5 Analytic solutions of the main problem: inner and outer solutions

In this section, we preserve the values of the elements involved in the main problem (17) stated in Section 3. More precisely, we assume (11)-(16), and also the hypotheses on the forcing term (19) in (18) and the coefficients in (20). Let $d\in\mathbb{R}$ and $S_{d}$ an infinite sector with vertex at $0\in\mathbb{C}$ under the geometric condition imposed in Proposition 6. Our main aim is to construct analytic solutions of (17) and their asymptotic behavior in different domains. For this purpose, we consider the analytic solutions as stated in Section 3.1.

Such solutions are defined in families of sectors with respect to the perturbation parameter, conforming good coverings of $\mathbb{C}^{\star}$ (see Definition 5). We also provide information about the difference of two solutions in consecutive sectors of the good covering, which will be crucial to determine the asymptotic behavior of the analytic solutions. We refer to consecutive solutions to solutions which are associated to consecutive elements in a fixed good covering of $\mathbb{C}^{\star}$. Let us first recall the notion of good covering in $\mathbb{C}^{\star}$.

Definition 5

Let $\iota\geq 2$ be an integer. For every $0\leq h\leq\iota-1$, we choose a finite sector with vertex at the origin $\mathcal{E}_{h}$ such that:

• •

$\mathcal{E}_{h}\subseteq D(0,\epsilon_{0})$, and $\mathcal{E}_{j_{1}}\cap\mathcal{E}_{j_{2}}=\emptyset$ if and only if $0\leq j_{1},j_{2}\leq\iota-1$ with $|j_{1}-j_{2}|\geq 2$ (under the convention that $\mathcal{E}_{\iota}:=\mathcal{E}_{0}$).

• •

$\cup_{h=0}^{\iota-1}\mathcal{E}_{h}=\mathcal{U}\setminus\{0\}$, for some neighborhood of the origin $\mathcal{U}$.

A family of sectors $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$ under these assumptions is known as a good covering in $\mathbb{C}^{\star}$.

Let $\tilde{\mathcal{T}}_{1}$ and $\tilde{\mathcal{T}}_{2}$ be sectors following the construction in Section 4.

Definition 6

Let $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$ be a good covering in $\mathbb{C}^{\star}$. We also fix a bounded sector, $\mathcal{T}_{1}$, and an unbounded sector $\mathcal{T}_{2}$ , both with vertex at the origin. For all $0\leq h\leq\iota-1$, let $S_{d_{h}}$ be an infinite sector of bisecting direction $d_{h}\in\mathbb{R}$. We say that the set $\{\mathcal{T}_{1},\mathcal{T}_{2},(\mathcal{E}_{h})_{0\leq h\leq\iota-1},(S_{d_{h}})_{0\leq h\leq\iota-1}\}$ is admissible if the following conditions hold:

• •

For every $0\leq h\leq\iota-1$, $\epsilon\in\mathcal{E}_{h}$ and $t_{1}\in\mathcal{T}_{1}$ we have $\epsilon^{\lambda_{1}}t_{1}\in\tilde{\mathcal{T}}_{1}$.

• •

For every $0\leq h\leq\iota-1$, $\epsilon\in\mathcal{E}_{h}$ and $t_{2}\in\mathcal{T}_{2}$ we have $\epsilon^{\lambda_{2}}t_{2}\in\tilde{\mathcal{T}}_{2}$.

Observe that given an admissible set $\{\mathcal{T}_{1},\mathcal{T}_{2},(\mathcal{E}_{h})_{0\leq h\leq\iota-1},(S_{d_{h}})_{0\leq h\leq\iota-1}\}$, the sectors $\tilde{\mathcal{T}}_{1}$, the choice of the sectors $\tilde{\mathcal{T}}_{1}$, $\tilde{\mathcal{T}}_{2}$ and $\gamma:=\gamma_{d}\in S_{d_{h}}$, fixed in Section 2.1, entail the existence of $\delta_{1}>0$ with

[TABLE]

In addition to this, it holds that

[TABLE]

for some $\delta_{2}>0$ and $\cos(\gamma_{d_{h}}-k_{2}\hbox{arg}(\epsilon^{\lambda_{2}}t_{2}))>\delta_{3}$, for some $\delta_{3}>0$, $\epsilon\in\mathcal{E}_{h}$, $0\leq h\leq\iota-1$.

Definition 7

Let $\{\mathcal{T}_{1},\mathcal{T}_{2},(\mathcal{E}_{h})_{0\leq h\leq\iota-1},(S_{d_{h}})_{0\leq h\leq\iota-1}\}$ be an admissible set.

Then, $\{(S_{d_{h}})_{0\leq h\leq\iota-1},\mathcal{T}_{1}\times\mathcal{T}_{2}\}$ is known as a family of sectors associated to the good covering $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$.

Theorem 1

Let $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$ be a good covering in $\mathbb{C}^{\star}$. For every $0\leq h\leq\iota-1$ we choose $d_{h}\in\mathbb{R}$ such that $S_{d_{h}}$ satisfies the geometric conditions of Proposition 6. Let $\{(S_{d_{h}})_{0\leq h\leq\iota-1},\mathcal{T}_{1}\times\mathcal{T}_{2}\}$ be a family of sectors associated to the good covering $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$. If there exist small enough $\epsilon_{0},C_{\psi}>0$ such that if

[TABLE]

then for every $0\leq h\leq\iota-1$ the problem (17) admits a solution $u_{h}(\boldsymbol{t},z,\epsilon)$, which defines a bounded and holomorphic function in $\mathcal{T}_{1}\times\mathcal{T}_{2}\times H_{\beta^{\prime}}\times\mathcal{E}_{h}$, for any fixed $0<\beta^{\prime}<\beta$.

Proof Let $0\leq h\leq\iota-1$. Proposition 6 guarantees the existence of $\epsilon_{0},C_{\psi}>0$ such that the equation (24) admits a unique solution $\omega^{d_{h}}(\tau,m,\epsilon)$ which belongs to $qExp^{d_{h}}_{(k^{\prime}_{1},\beta,\mu,\alpha)}$ for every $\epsilon\in D(0,\epsilon_{0})$, and the map $\epsilon\mapsto\omega^{d_{h}}(\tau,m,\epsilon)$ is holomorphic in $D(0,\epsilon_{0})$. Taking into account that $\{(S_{d_{h}})_{0\leq h\leq\iota-1},\mathcal{T}_{1}\times\mathcal{T}_{2}\}$ is associated to the good covering $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$ and the properties of Laplace transform stated in Section 2.1, one can construct $U_{\gamma_{h}}(\boldsymbol{T},m,\epsilon)$ in the form (23), which is well defined, bounded and continuous function on $\tilde{\mathcal{T}}_{1}\times\tilde{\mathcal{T}}_{2}\times\mathbb{R}\times D(0,\epsilon_{0})$, where $\tilde{\mathcal{T}}_{1}\subseteq\mathcal{R}_{\gamma_{h},\tilde{\delta}}\cap D(0,r_{1})$ for $0<r_{1}\leq q^{(1/2-\alpha)/k_{1}}/2$, and $\tilde{\mathcal{T}}_{2}$ is an infinite sector. The function $U_{\gamma_{h}}(\boldsymbol{T},m,\epsilon)$ is holomorphic w.r.t. $(\boldsymbol{T},\epsilon)$ on $\tilde{\mathcal{T}}_{1}\times\tilde{\mathcal{T}}_{2}\times D(0,\epsilon_{0})$. As a matter of fact, $U_{\gamma_{h}}$ is a solution of (22) in view of the properties relating both equalities in Section 8. We finally define $u_{h}(\boldsymbol{t},z,\epsilon)$ following (21):

[TABLE]

which turns out to be a holomorphic solution of (17), defined on $\mathcal{T}_{1}\times\mathcal{T}_{2}\times H_{\beta^{\prime}}\times\mathcal{E}_{h}$, for any fixed $0<\beta^{\prime}<\beta$. $\Box$

In order to provide the asymptotic behavior of the analytic solutions of (17) in different domains, with respect to the perturbation parameter $\epsilon$, we state the definition of inner and outer solutions of the problem (17).

5.1 Inner solutions of the main problem

Definition 8

Let $\iota_{1}\geq 2$. Let $(\mathcal{E}^{\infty}_{h_{1}})_{0\leq h_{1}\leq\iota_{1}-1}$ be a good covering of $\mathbb{C}^{\star}$. We also consider the admissible set $\{\mathcal{T}_{1},\mathcal{T}_{2},(\mathcal{E}^{\infty}_{h_{1}})_{0\leq h_{1}\leq\iota_{1}-1},(S^{\infty}_{d_{h_{1}}})_{0\leq h_{1}\leq\iota_{1}-1}\}$. Let $\mu_{2}>0$ be a natural number satisfying

[TABLE]

Let $\chi_{2}^{\infty}$ be a bounded domain, such that the good covering $(\mathcal{E}^{\infty}_{h_{1}})_{0\leq h_{1}\leq\iota_{1}-1}$ satisfies the following condition: for all $0\leq h_{1}\leq\iota_{1}-1$ we can select $\theta_{h_{1}}\in\mathbb{R}$ (which depends on $\mathcal{E}_{h_{1}}^{\infty}$) such that for every $x_{2}\in\chi_{2}^{\infty}$ and $\epsilon\in\mathcal{E}_{h_{1}}^{\infty}$, the complex number $t_{2}=\frac{x_{2}}{\epsilon^{\mu_{2}}}e^{\theta_{h_{1}}\sqrt{-1}}$ belongs to $\mathcal{T}_{2}$. We define the set $\mathcal{T}_{2,\epsilon,\mu_{2}}:=\{\frac{x_{2}}{\epsilon^{\mu_{2}}}e^{\theta_{h_{1}}\sqrt{-1}}:x_{2}\in\chi_{2}^{\infty}\}$.

In case that $\epsilon\in\mathcal{E}^{\infty}_{h_{1}}$, $t_{2}\in\mathcal{T}_{2,\epsilon,\mu_{2}}$, $t_{1}\in\mathcal{T}_{1}$, $z\in H_{\beta^{\prime}}$ for $0<\beta^{\prime}<\beta$, then we say that $u_{h_{1}}(\boldsymbol{t},z,\epsilon)$ represents an inner solution of (17).

Theorem 2

Under the assumptions of Theorem 1 and the constraints on the inner solutions of the main problem of Definition 8, let $(\mathcal{E}^{\infty}_{h_{1}})_{0\leq h_{1}\leq\iota_{1}}$ be a good covering of $\mathbb{C}^{\star}$. Then, there exists $C_{inn}>0$ such that for every $0\leq h_{1}\leq\iota_{1}-1$ , $\epsilon\in\mathcal{E}^{\infty}_{h_{1}}\cap\mathcal{E}^{\infty}_{h_{1}+1}$, $t_{2}\in\mathcal{T}_{2,\epsilon,\mu_{2}}$, $z\in H_{\beta^{\prime}}$ for any fixed $0<\beta^{\prime}<\beta$ and $t_{1}\in\mathcal{T}_{1}$, one has

[TABLE]

for some $\hat{C}_{1},\hat{C}_{2},\hat{C}_{3}>0$, $\Delta_{1}\in\mathbb{R}$ and a real number $k^{\prime\prime}_{1}\in(0,k_{1})$ chosen small enough.

Proof Let $0\leq h_{1}\leq\iota_{1}-1$ and consider consecutive solutions $u_{h_{1}},u_{h_{1}+1}$ of (17), constructed in Theorem 1. We recall that the function $\omega^{d_{h_{1}}}(\tau,m,\epsilon)$ stands for the analytic continuation of a function $\omega(\tau,m,\epsilon)$, holomorphic w.r.t. $\tau$ in a neighborhood of the origin $D(0,\rho)$, to the infinite sector $S_{d_{h_{1}}}$. This entails that the difference $u_{h_{1}+1}(\boldsymbol{t},z,\epsilon)-u_{h_{1}}(\boldsymbol{t},z,\epsilon)$ can be written in the following form, after an appropriate path deformation in $\tau$, avoiding the roots of $P_{m}(\tau)$ (see Figure 1).

We write $u_{h_{1}+1}(\boldsymbol{t},z,\epsilon)-u_{h_{1}}(\boldsymbol{t},z,\epsilon)=J_{1}-J_{2}+J_{3}$, where

[TABLE]

where $L_{\gamma_{j},\rho/2}=[\rho/2e^{\gamma_{j}\sqrt{-1}},\infty)$ for $j=h_{1},h_{1}+1$, and $C_{h_{1},h_{1}+1,\rho/2}$ stands for the arc of circle from $\rho/2e^{\gamma_{h_{1}}\sqrt{-1}}$ to $\rho/2e^{\gamma_{h_{1}+1}\sqrt{-1}}$.

Assume that $\epsilon\in\mathcal{E}^{\infty}_{h_{1}}\cap\mathcal{E}^{\infty}_{h_{1}+1}$ and $t_{2}\in\mathcal{T}_{2,\epsilon,\mu_{2}}$, $z\in H_{\beta^{\prime}}$ for some $0<\beta^{\prime}<\beta$ and $t_{1}\in\mathcal{T}_{1}$. Owing to the estimates leading to (29), displayed in (49), we derive that $\hbox{dist}(\chi_{2}^{\infty}\epsilon^{\lambda_{2}-\mu_{2}},0)\geq\rho_{2}^{\infty}$ for some large enough $\rho_{2}^{\infty}$. Moreover, if $\rho_{0}>0$ is small enough, then there exists $\tilde{C}_{8}>0$ such that

[TABLE]

for every $t_{1}\in\mathcal{T}_{1}$, $z\in H_{\beta^{\prime}}$, $\epsilon\in\mathcal{E}^{\infty}_{h_{1}}\cap\mathcal{E}^{\infty}_{h_{1}+1}$, and $x_{2}\in\chi_{2}^{\infty}$. An analogous upper bound is attained for $|J_{2}|$. Concerning $|J_{3}|$, one can apply the bounds stated in Proposition 6, and (6) to get the existence of $\tilde{C}_{9}>0$ such that

[TABLE]

The result follows taking into account (34) and (35), under the condition (31). $\Box$

5.2 Outer solutions of the main problem

Definition 9

Let $\iota_{2}\geq 2$ and $(\mathcal{E}_{h_{2}}^{0})_{0\leq h_{2}\leq\iota_{2}}$ be a good covering of $\mathbb{C}^{\star}$, and consider an admissible set $\{\mathcal{T}_{1},\mathcal{T}_{2},(\mathcal{E}^{0}_{h_{2}})_{0\leq h_{2}\leq\iota_{2}-1},(S^{0}_{d_{h_{2}}})_{0\leq h_{2}\leq\iota_{2}-1}\}$. Assume that $t_{2}\in\mathcal{T}_{2}$ is such that $|t_{2}|<\rho_{2}$ for some fixed $\rho_{2}>0$ which is independent of $\epsilon$, then we say that $u_{h_{2}}(\boldsymbol{t},z,\epsilon)$ represents an outer solution of (17).

Theorem 3

Under the assumptions of Theorem 1 and the constraints on the outer solutions of the main problem of Definition 9, let $(\mathcal{E}^{0}_{h_{2}})_{0\leq h_{2}\leq\iota_{2}}$ be a good covering of $\mathbb{C}^{\star}$. Then, there exists $C_{out}>0$ such that for every $0\leq h_{2}\leq\iota_{2}-1$ , $\epsilon\in\mathcal{E}^{0}_{h_{2}}\cap\mathcal{E}^{0}_{h_{2}+1}$, $t_{2}\in\mathcal{T}_{2}$ with $|t_{2}|<\rho_{2}$, $z\in H_{\beta^{\prime}}$ for any fixed $0<\beta^{\prime}<\beta$ and $t_{1}\in\mathcal{T}_{1}$, one has

[TABLE]

for some $\hat{C}_{4},\hat{C}_{5}>0$ and $\Delta_{2}\in\mathbb{R}$.

Proof Let $0\leq h_{2}\leq\iota_{2}-1$ and consider consecutive solutions $u_{h_{2}},u_{h_{2}+1}$ of (17), constructed in Theorem 1. We proceed to write the difference of two consecutive solutions in the form

[TABLE]

for all $0\leq h_{2}\leq\iota_{2}-1$, with $h_{1}$ substituted by $h_{2}$ in the expressions of $J_{1},J_{2},J_{3}$ in the proof of Theorem 2. Let $\epsilon\in\mathcal{E}^{0}_{h_{2}}\cap\mathcal{E}^{0}_{h_{2}+1}$ and $t_{2}\in\mathcal{T}_{2}$, with $|t_{2}|<\rho_{2}$, $z\in H_{\beta^{\prime}}$ for some $0<\beta^{\prime}<\beta$ and $t_{1}\in\mathcal{T}_{1}$. Owing to analogous bounds as those leading to (50) and (51), we arrive at

[TABLE]

for some $\tilde{C}_{10}>0$. Similar estimates hold for $|J_{2}|$. On the other hand, direct computations yield

[TABLE]

for some $\tilde{C}_{11}>0$. This concludes the proof. $\Box$

6 Asymptotic expansions of mixed order

This section is divided in two parts. The first part recalls some facts about $q-$asymptotic expansions and also describes $q-$asymptotic expansions which show a sub-Gevrey growth in their estimates, and related results.

In the second part of this section, we provide the existence of a formal solution of the main problem under study, written as a formal power series in the perturbation parameter, and explain the asymptotic relationship between this and the analytic solutions.

6.1 Review on $q$-asymptotic expansions

In the whole subsection, $(\mathbb{F},\left\|\cdot\right\|)$ stands for a complex Banach space.

We first recall the notion of $q-$Gevrey asymptotic expansions, which can be also found in [12] in more detail.

Definition 10

Let $V$ be a bounded open sector with vertex at 0 in $\mathbb{C}$, $q\in\mathbb{R}$ with $q>1$. We also fix a positive integer $k$. We say that a holomorphic function $f:V\to\mathbb{F}$ admits the formal power series $\hat{f}(\epsilon)=\sum_{n\geq 0}f_{n}\epsilon^{n}\in\mathbb{F}[[\epsilon]]$ as its $q-$Gevrey asymptotic expansion of order $1/k$ if for every open subsector $U$ of $V$, i.e. $\overline{U}\setminus\{0\}\subseteq V$, there exist $A,C>0$ such that

[TABLE]

for every $\epsilon\in U$ and $N\geq 0$.

Such pure $q-$asymptotic expansions have been recently studied when dealing with the asymptotic behavior of the solutions of $q-$difference-differential equations in the complex domain [5, 11, 12, 19].

We recall the definition of asymptotic expansion of mixed order as introduced in the work [20].

Definition 11

Let $V$ be a bounded open sector of the complex plane with vertex at the origin. Let $q>1$ and $s\geq 0$ be real numbers, and let $k$ be a positive integer. The function $f:V\to\mathbb{F}$ is said to admit the formal series $\hat{F}(\epsilon)=\sum_{n\geq 0}f_{n}\epsilon^{n}\in\mathbb{F}[[\epsilon]]$ as its $((q,k);s)-$Gevrey asymptotic expansion if for every open subsector $U$ of $V$ (i.e. $(\overline{U}\setminus\{0\})\subseteq V$) there exist $A,C>0$ such that

[TABLE]

for all $\epsilon\in U$ and $N\geq 0$.

Observe that the set of functions admitting $q-$Gevrey asymptotic expansions of order $1/k$ coincide with the functions admitting $((q,k);0)-$Gevrey asymptotic expansion.

We now proceed to state a $((q,k);s)-$version of the Ramis-Sibuya Theorem. The classical statement of this result involves Gevrey asymptotic expansions, and guarantee $s$-summability of some power seires. This cohomological criterion can be found in Proposition 2 [1] and Lemma XI-2-6 [9].

This version has already been stated in the work [20] where a complete proof can be found therein.

Theorem 4** ($((q,k);s)-$Ramis-Sibuya Theorem)**

Let $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$ be a good covering in $\mathbb{C}^{\star}$. For every $0\leq h\leq\iota-1$ we consider a holomorphic function $G_{h}:\mathcal{E}_{h}\to\mathbb{F}$ and define $\Delta_{h}:=G_{h+1}-G_{h}$, which turns out to be a holomorphic function in $Z_{h}:=\mathcal{E}_{h}\cap\mathcal{E}_{h+1}$ (here, $\mathcal{E}_{\iota}$ and $G_{\iota}$ stand for $\mathcal{E}_{0}$ and $G_{0}$, respectively). Assume that the following statements hold:

• •

$G_{h}$* is a bounded function in a vicinity of 0, for every $0\leq h\leq\iota-1$.*

• •

*The function $\Delta_{h}(\epsilon)$ admits null *$((q,k);s)-$Gevrey asymptotic expansion in $Z_{h}$ for every $0\leq h\leq\iota-1$, i.e. there exist $C_{1},C_{2}>0$ such that

[TABLE]

for every $\epsilon\in Z_{h}$ and all $N\geq 0$.

Then, there exists $\hat{G}\in\mathbb{F}[[\epsilon]]$ which is the common $((q,k);s)-$Gevrey asymptotic expansion of $G_{h}(\epsilon)$ on $\mathcal{E}_{h}$, for all $0\leq h\leq\iota-1$.

In [19], a $q-$Gevrey version of Ramis-Sibuya is obtained. That version is related to $q-$Gevrey asymptotic expansions of some positive order $k$, which coincide with $(q,k);0)-$Gevrey asymptotic expansions in our framework. Ramis-Sibuya theorem in that framework reads as follows:

Theorem 5** ((q-RS))**

Let $(\mathbb{F},\left\|\cdot\right\|_{\mathbb{F}})$ be a Banach space and $(\mathcal{E}_{h})_{0\leq h\leq\iota-1}$ be a good covering in $\mathbb{C}^{\star}$. For every $0\leq h\leq\iota-1$, let $G_{h}(\epsilon)$ be a holomorphic function from $\mathcal{E}_{h}$ into $\mathbb{F}$ and let the cocycle $\Delta_{h}(\epsilon)=G_{h+1}(\epsilon)-G_{h}(\epsilon)$ be a holomorphic function from $Z_{h}=\mathcal{E}_{h+1}-\mathcal{E}_{h}$ into $\mathbb{F}$ (with the convention that $\mathcal{E}_{\iota}=\mathcal{E}_{0}$ and $G_{\iota}=G_{0}$). We make further assumptions:

The functions $G_{h}(\epsilon)$ are bounded as $\epsilon$ tends to 0 on $\mathcal{E}_{h}$, for all $0\leq h\leq\iota-1$. 2. 2.

*The function $\Delta_{h}(\epsilon)$ is *$q-$exponentially flat of order $k$ on $Z_{h}$ for all $0\leq h\leq\iota-1$, meaning that there exist two constants $C_{h}^{1}\in\mathbb{R}$ and $C^{2}_{h}>0$ with

[TABLE]

for all $\epsilon\in Z_{h}$, and all $0\leq h\leq\iota-1$.

Then, there exists a formal power series $\hat{G}(\epsilon)\in\mathbb{F}[[\epsilon]]$ which is the common $q-$Gevrey asymptotic expansion of order $1/k$ of the function $G_{h}(\epsilon)$ on $\mathcal{E}_{h}$, for all $0\leq h\leq\iota-1$.

6.2 Asymptotic expansions for the analytic solutions of the main problem

In this section, we preserve all the assumptions made on the elements involved in the main problem (17), detailed in Section 3. Moreover, we depart from the geometric construction of the elements used to construct the analytic solutions of (17), described in Sections 4 and 5, and under Assumption (31).

The next result shows that the difference of two consecutive solutions of the main problem allow the application of the $((q,k);s)-$version of Ramis-Sibuya Theorem, obtained in Theorem 4, in adequate domains.

Lemma 3

Let $k^{\prime}_{1},k^{\prime\prime}_{1},k_{2},\lambda_{1},\lambda_{2},\mu_{2}$ be positive constants with $\mu_{2}>\lambda_{2}$. For every $n\in\mathbb{N}$, we consider the function

[TABLE]

Then, it holds that

[TABLE]

for $n\geq n_{0}$, where $n_{0}\geq 1$ is a large enough integer depending on $k^{\prime}_{1},k^{\prime\prime}_{1},k_{2},\lambda_{1},\lambda_{2},\mu_{2}$ with

[TABLE]

Proof We make the change of variable $\tau=\log(x)$ and consider the auxiliary function $\overline{\Psi}(\tau):=\Psi(e^{\tau})$. We search for the maximum of $\overline{\Psi}$ for $\tau\in\mathbb{R}$. It holds that

[TABLE]

It is known that the solution of the equation $\log(t)+\log(b)t=\log(a)$, for some fixed $a,b>0$, is given by $t=W(a\log(b))/\log(b)$, where $W$ is the Lambert $W$ function. The maximum of $\overline{\Psi}(\tau)$ is then attained at

[TABLE]

where $W_{-1}$ stands for the $-1$-branch of Lambert $W$ function. Let

[TABLE]

Then $\tau_{n}$ can be written in the form $\tau_{n}=\frac{-1}{C}W_{-1}(-\exp(-\frac{n\log(q)}{A}-1)C)$. Regarding Theorem 1 [2], we have

[TABLE]

which can be applied to estimate $\tau_{n}$:

[TABLE]

leading to

[TABLE]

We conclude that $\overline{\Psi}(\tau)\leq\overline{\Psi}(\tau_{n})$. Taking into account (41) we arrive at

[TABLE]

for every $0<\xi<1$, and some $\hat{C}_{1},\hat{A}_{1}>0$. The definition of $A,C$ in (40) yields

[TABLE]

Statement (38) follows directly from the last inequality. $\Box$

Theorem 6

Let the assumptions of Theorem 2 hold. Let $\mathbb{F}_{1}$ be the Banach space of holomorphic and bounded functions on $\mathcal{T}_{1}\times\chi_{2}^{\infty}\times H_{\beta^{\prime}}$. For all $0\leq h_{1}\leq\iota_{1}-1$, we consider the function

[TABLE]

which defines a holomorphic and bounded function on $\mathcal{E}^{\infty}_{h_{1}}$, with values in $\mathbb{F}_{1}$. Then, there exist a formal series $\hat{u}^{\infty}(\epsilon)\in\mathbb{F}_{1}[[\epsilon]]$ such that for all $0\leq h_{1}\leq\iota_{1}-1$, the function (43) admits $\hat{u}^{\infty}(\epsilon)$ as its $((q,K);S)-$Gevrey asymptotic expansion on $\mathcal{E}^{\infty}_{h_{1}}$, for $K$ and $S$ stated in (39).

Proof Let $0\leq h_{1}\leq\iota_{1}-1$. Regarding (32) in Theorem 2, Lemma 3 and from usual estimates, we guarantee that for every $\epsilon\in\mathcal{E}^{\infty}_{h_{1}}\cap\mathcal{E}^{\infty}_{h_{1}+1}$,

[TABLE]

Taking into account Lemma 3, we get that

[TABLE]

Theorem 4 states the existence of a formal power series $\hat{u}^{\infty}(\epsilon)\in\mathbb{F}_{1}[[\epsilon]]$ which is the common $((q,K);S)$-Gevrey asymptotic expansion of the function (43), as a function on $\mathcal{E}^{\infty}_{h_{1}}$ with values in $\mathbb{F}_{1}$, for every $0\leq h_{1}\leq\iota_{1}-1$. $\Box$

Theorem 7

Let the assumptions of Theorem 3 hold. Let $\mathbb{F}_{2}$ be the Banach space of holomorphic and bounded functions on $\mathcal{T}_{1}\times(\mathcal{T}_{2}\cap D(0,\rho_{2}))\times H_{\beta^{\prime}}$. For every $0\leq h_{2}\leq\iota_{2}-1$, we consider the function

[TABLE]

which is an outer solution of (17), holomorphic and bounded on $\mathcal{E}_{h_{2}}^{0}$, with values in $\mathbb{F}_{2}$. Then, there exists a formal power series $\hat{u}^{0}(\epsilon)\in\mathbb{F}_{2}[[\epsilon]]$ such that the function (44) admits $\hat{u}^{0}(\epsilon)$ as its $q-$Gevrey asymptotic expansion of order $1/(k_{1}\lambda_{1}^{2})$ on $\mathcal{E}_{h_{2}}^{0}$, for all $0\leq h_{2}\leq\iota_{2}-1$.

Proof Taking into account (36) in Theorem 3, and Theorem 5, we get the existence of a formal power series $\hat{u}^{0}(\epsilon)\in\mathbb{F}_{2}[[\epsilon]]$ which is the common $q$-Gevrey asymptotic expansion of the function (44), as a function on $\mathcal{E}^{0}_{h_{2}}$ with values in $\mathbb{F}_{2}$. This holds for all $0\leq h_{2}\leq\iota_{2}-1$.

$\Box$

7 Proof of Proposition 7

In this section, we give proof of the technical Proposition 7. We consider $U_{\gamma}(\boldsymbol{T},m,\epsilon)$ is constructed in the form (23).

For the first part of the proof, we take into account the property (6) on Jacobi Theta function, and Proposition 6, to arrive at

[TABLE]

where

[TABLE]

We split $I(|T_{1}|,|T_{2}|)$ into the sum of $I_{1}(|T_{1}|,|T_{2}|)$ and $I_{2}(|T_{1}|,|T_{2}|)$, where the first element is associated to the integration in $(0,\rho)$ and the second is concern with the integration restricted to $(\rho,\infty)$, for some $\rho>0$. We study each part of the splitting:

We have

[TABLE]

for some $C_{3}>0$, which after the change of variable $r=|T_{1}|r^{\prime}$ equals

[TABLE]

which is bounded for every $|T_{1}|\in\mathbb{R}_{+}$.

On the other hand, we assume $\rho,\rho_{1}>0$ are such that $\rho\geq\rho_{1}$. Then, the positive function $\phi:[\rho,\infty)\to\mathbb{R}$ defined by $\phi(r)=\log^{2}(r/|T_{1}|)$ is monotone increasing on $[\rho,\infty)$ for any choice of $|T_{1}|<\rho_{1}$. Therefore,

[TABLE]

where

[TABLE]

Let $x=|T_{2}|^{k_{2}}/\delta_{3}$. We make the change of variable $r=x\tilde{r}$ in the last integral and arrive at

[TABLE]

Taking into account that

[TABLE]

[TABLE]

we get that $I_{2.1}(x)=I_{2.2}(x)+I_{2.3}(x)$, where the splitting is done on the integral by cutting the integration path into $(0,1-\delta/x)$ and $(1-\delta/x,\infty)$ for $I_{2.2}(x)$ and $I_{2.3}(x)$, respectively.

[TABLE]

for some $C_{4}>0$ which does not depend on $x$. Concerning $I_{2.3}(x)$, we proceed analogously to arrive at

[TABLE]

for some $C_{5}>0$. In the last sequence of inequalities, we have made the change of variable $\tilde{r}+\frac{\delta}{x}=2h$. The application of Stirling formula $\Gamma(x+1)\sim_{x\to\infty}\sqrt{2\pi}x^{1/2}x^{x}e^{-x}$ leads us to

[TABLE]

for some $C_{6}>0$. From (46), (47) and (48) we conclude

[TABLE]

for some $C_{7}>0$. The first statement of Proposition 7 holds. We give proof for the second statement.

The first arguments in the proof of the first statement can be followed word by word up to the splitting of $I(|T_{1}|,|T_{2}|)$ into $I_{1}(|T_{1}|,|T_{2}|)+I_{2}(|T_{1}|,|T_{2}|)$. The quantity $I_{1}(|T_{1}|,|T_{2}|)$ is upper bounded by a constant for every $|T_{1}|>0$ and $|T_{2}|>0$. We now proceed to give upper estimates on $I_{2}(|T_{1}|,|T_{2}|)$. Let us choose $\rho_{1},\rho>0$ such that $\phi(r)$ is monotone increasing on $[\rho,\infty)$. It holds that

[TABLE]

where

[TABLE]

for some $C_{8}>0$. The conclusion follows from this last upper bound.

8 Connection of the solutions of (22) and (24)

Let $\epsilon\in D(0,\epsilon_{0})$. We consider $\omega(\tau,m,\epsilon)\in q\hbox{Exp}^{d}_{(k^{\prime}_{1},\beta,\mu,\alpha)}$. Let $\gamma\in\mathbb{R}$ be an argument of $S_{d}$.

We depart from

[TABLE]

Let $\tilde{\mathcal{T}}_{1}$ be a bounded sector and $\tilde{\mathcal{T}}_{2}$ be an unbounded sector, both with vertex at the origin, satisfying the assumptions in Section 5. $U_{\gamma}(\boldsymbol{T},m,\epsilon)$ turns out to be a continuous function on $\tilde{\mathcal{T}}_{1}\times\tilde{\mathcal{T}}_{2}\times\mathbb{R}$, holomorphic with respect to $\boldsymbol{T}$ on $\tilde{\mathcal{T}}_{1}\times\tilde{\mathcal{T}}_{2}$.

Proposition 8

In the previous framework, the following identities hold for all $T_{1}\in\tilde{\mathcal{T}}_{1},T_{2}\in\tilde{\mathcal{T}}_{2},m\in\mathbb{R}$:

• $(i)$

For every integer $\delta\geq 0$ one has

[TABLE]

• $(ii)$

For every $1\leq\ell_{1}\leq D_{1}-1$ one has

[TABLE]

• $(iii)$

It holds that

[TABLE]

• $(iv)$

For every $1\leq\ell_{1}\leq D_{1}-1$ and $1\leq\ell_{2}\leq D_{2}-1$ we have

[TABLE]

Proof Statement $(i)$ is a consequence of the derivation under the integral sign. Statement $(ii)$ is a direct consequence of Proposition 6 in [19] with $k:=k_{1}$, $\sigma:=d_{\ell_{1}}$, $j:=\delta_{\ell_{1}}$ and $f(x)=\omega(x,m,\epsilon)\exp(-\left(\frac{1}{T_{2}}\right)^{k_{2}}x)$. The third statement is a consequence of Proposition 6 in [19]. The last statement is derived from the application of $(ii)$ and the fact that

[TABLE]

$\Box$