On the local well-posedness of the nonlinear heat equation associated to the fractional Hermite operator in modulation spaces

In this note we consider the nonlinear heat equation associated to the fractional Hermite operator Hβ=(-Δ+|x|2)β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H^\beta =(-\Delta +|x|^2)^\beta $$\end{document}, 0<β≤1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0<\beta \le 1$$\end{document}. We show the local solvability of the related Cauchy problem in the framework of modulation spaces. The result is obtained by combining tools from microlocal and time-frequency analysis. As a byproduct, we compute the Gabor matrix of pseudodifferential operators with symbols in the Hörmander class S0,0m\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S^m_{0,0}$$\end{document}, m∈R\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m\in \mathbb {R}$$\end{document}.


Introduction and results
In this note we study the Cauchy problem for the nonlinear heat equation associated to the fractional Hermite operator F is a scalar function on C, with F(0) = 0. The solution u(t, x) is a complex valued function of (t, x) ∈ R × R d . We will consider the case in which F is a real analytic function with an entire extension.
Well-posedness of the heat equation has been studied by many authors, see e.g. [10,21] and the many contributions by Wong, for instance [31,32], see also [6]. In particular, heat equations associated to fractional Hermite operators were recently studied in [3], for related results see also [5]. Hermite multipliers are considered in [4], see also the textbook [25]. Recently the study of Cauchy problems in modulation spaces have been pursued by many authors, see the pioneering works [1,2]. Many deep results in this framework for nonlinear evolution equations have been obtained by B. Wang et al. in [27,28] and are also available in the textbook [30].
Following the spirit of [8,12], we shall prove the local existence and uniqueness of the solutions in modulation spaces to the Cauchy problem (1). The key arguments come from both microlocal and time-frequency analysis. In fact, we shall rely on the results related to spectral theory of globally elliptic operators developed by Helffer [20] to understand the properties of the fractional Hermite operators H β . Namely, they are pseudodifferential operators with Weyl symbols a β positive globally elliptic and in the Shubin classes 2β 1 (see Definition 2.2 and the estimate (13) below). The spectral decomposition of the Hermite operator H = − + |x| 2 is given by where P k is the orthogonal projection of L 2 (R d ) onto the eigenspace corresponding to the eigenvalue (2k +d). Namely, the range of the operator P k is the space spanned by the Hermite functions α in R d , with α multi-index in N d , such that |α| = k. The solution to the homogeneous Cauchy problem (1) (i.e., F = 0) can be formally written in terms of the heat semigroup related to H β We shall prove that the propagator K β (t) = e −t H β can be represented as a pseudodifferential operator with Weyl symbol in the Shubin class 0 1 , with related semi-norms uniformly bounded with respect to the time variable t ∈ [0, T ], for any fixed T > 0.
After that we shall leave the microlocal techniques to come to time-frequency analysis. We perform a general study concerning the boundedness of Shubin τpseudodifferential operators with symbols in the Hörmander classes S m 0,0 (for τ = 1/2 we recapture the Weyl case). The outcomes are contained in Theorem 2.4 below (see also the subsequent corollary and remark).
The main tool here is to study the decay of their related Gabor matrix representations, which we shall also control by the semi-norms in S m 0,0 . We think that such result is valuable in and of itself.
We then use the special case of Weyl operators to study (1). The integral version of the problem (1) has the form where To show that the Cauchy Problem (1) has a unique solution, we use a variant of the contraction mapping theorem (see Proposition 3.3 below). As already mentioned, the function spaces used for our results are weighted modulation spaces M p,q w , 1 ≤ p, q ≤ ∞, introduced by H. Feichtinger in 1983 [14] (then extended to 0 < p, q ≤ ∞ in [15]). We refer the reader to Sect. 2 for their definitions and main properties.
The local well-posedness results for modulation spaces read as follows: We actually do not know whether it is possible to obtain better results concerning the nonlinearity F(u) = λ|u| 2k u, k ∈ N, we refer to the work [5] for a discussion on the topic.
The tools employed follow the pattern of similar Cauchy problems studied for other equations such as the Schrödinger, wave and Klein-Gordon equations [1,2,8,12].
To compare with other results in the literature, we observe that also in [31,32], and in [6] the authors use Wigner distributions and pseudodifferential operators as tools for their main results. In the latter paper the author gives a formula for the oneparameter strongly continuous semigroup e −t H β in terms of the Weyl transforms of a L 2 -orthonormal basis made of generalized Hermite eigenfunctions. This is then used to obtain L 2 -estimates for the solution of the related initial value problem with data in L p spaces, 1 ≤ p ≤ ∞. Here the approach uses similar ideas of joining microlocal and time-frequency analysis tools, but the spaces employed are different: we use modulation spaces, which are the most common ones in time-frequency analysis.

Function spaces and preliminaries
We denote by v a continuous, positive, submultiplicative weight function on for all z 1 , z 2 ∈ R d (or for all z 1 , z 2 ∈ Z d ). We will mainly work with polynomial weights of the type Observe that, for s < 0, v s is v |s| -moderate. Moreover, we limit to weights w with at most polynomial growth, that is there exists C > 0, s > 0 such that The main features of time-frequency analysis are T x and M ξ , the so-called translation and modulation operators, defined by T x g(y) = g(y − x) and M ξ g(y) = e 2πiξ y g(y). Let g ∈ S(R d ) be a non-zero window function in the Schwartz class and consider the short-time Fourier transform (STFT) V g f of a function/tempered distribution f in S (R d ) with respect to the the window g: i.e., the Fourier transform F applied to f T x g.
(with obvious changes with p = ∞ or q = ∞) is finite.
For 1 ≤ p, q ≤ ∞ they are Banach spaces, whose norm does not depend on the window g, in the sense that different window functions in S(R d ) yield equivalent norms. Moreover, the window class S(R d ) can be extended to the modulation space , we shall simply write, using the standard notation [14], In our study, we will apply Minkowski's integral inequality to study the operator B in (3). Such inequalities do not hold whenever the indices p < 1 or q < 1, hence we shall limit ourselves to the cases 1 ≤ p, q ≤ ∞.
We do not know whether the local well-posedness is still valid in the quasi-Banach setting.
Recall that for [16,Thm. 11.3.7]). In other words, given any For this work we will use the inversion formula for the STFT (see [16,Proposition 11 and the equality holds in M p,q w (R d ). We also recall their inclusion relations: Other properties and more general definitions of modulation spaces can now be found in textbooks [9,16]. consists of all complex functions a ∈ C ∞ (R 2d ) such that for every α ∈ N 2d there exists a constant C α ≥ 0 with

Shubin classes and symbols of the operators
It immediately follows from this definition that if a ∈ m 1 (R 2n ) and α ∈ N 2n then ∂ α z a ∈ m−|α| 1 (R 2n ). Obviously m 1 (R 2n ) is a complex vector space for the usual operations of addition and multiplication by complex numbers, and we have The notion of asymptotic expansion of a symbol a ∈ m 1 (R 2d ) (cf. [24], Definition 23.2) reads as follows.
where m r = max j≥r m j we will write a ∼ ∞ j=0 a j and call this relation an asymptotic expansion of the symbol a.
The interest of the asymptotic expansion comes from the fact that every sequence of symbols (a j ) j with a j ∈ m j 1 (R 2d ), the degrees m j being strictly decreasing and such that m j → −∞ determines a symbol in some m 1 (R 2d ), that symbol being unique up to an element of S(R 2d ).
The symbol of the Hermite operator (or harmonic oscillator) H (z) = (|x| 2 + 4π 2 |ξ | 2 ) obviously belongs to 2 1 (R 2d ). From [20, Theorem 1.11.1] we infer that the fractional power H β , 0 < β < 1, can be written as a Weyl pseudodifferential operator having real symbol a β in the Shubin class 2β 1 and positive globally elliptic. Recall that a symbol a β is positive globally elliptic if there exist C > 0 and R > 0 such that Thanks to the properties of H β above, we can exploit a result by Nicola and Rodino in [22,Theorem 4.5.1] to prove that the operator e −t H β is a pseudodifferential operator with Weyl symbol in the Shubin class 0 1 , with uniform estimates with respect to t ∈ [0, T ], for any fixed T > 0.
For this purpose, we use the above theorem in the following setting: moreover we choose the parameters l = N = 0 and J = 1. If we consider the asymptotic expansion a β ∼ ∞ 0 a β, j as in Definition 2.3, together with the ellipticity condition a β,0 (z) z 2β , |z| ≥ R, then Theorem 4.5.1. guarantees that the operator e −t H β is a pseudodifferential operator with Weyl symbol b(t, z) satisfying, for every k ∈ N, T > 0, the estimate where b 0 (t, z) = e −ta β,0 (z) , and the semi-norms | · | k , k ∈ N, are defined by Defining the remainder R 1 (t, z) := b(t, z)−b 0 (t, z), we infer from (14) that R 1 (t, ·) ∈ −2 Moreover, using the ellipticity condition a β,0 (z) z , for |z| ≥ R, and the property a β,0 ∈ 2β 1 (R 2d ), one easily shows by induction that there exists a constant C > 0: that is to say, the symbol b 0 (t, z) is in the Shubin class 0 For applications of Shubin classes in the framework of Born-Jordan quantization we refer to the work [7].

Gabor analysis of -pseudodifferential operators
It can be used to define the τ -pseudodifferential operator with symbol σ via the formula For τ = 1/2 we recapture the Weyl operator. We want to consider τ -pseudodifferential operators with symbols σ in the Hörmander class S m 0,0 , m ∈ R, consisting of functions σ ∈ C ∞ (R 2d ) such that, for every α ∈ N 2d , The related semi-norms are denote by Fix g ∈ S(R d )\{0}. We define the Gabor matrix of a linear continuous operator This is a slightly abuse of notation, since originally Gabor matrices we defined for time-frequency shifts π(λ), with λ varying in a lattice ⊂ R 2d . We observe that the almost diagonalization of Gabor matrices of pseudodifferential operators with symbols in the modulation space M ∞,1 (R 2d ) treated in [17] (and in many subsequent papers on the topic) are valid in both the continuous and discrete case. So we adopt this terminology in the continuous framework.
For m = 0 we are reduced to the Hörmander class S 0 0,0 , whose Gabor matrix characterization for Weyl operators was shown in [18, Theorem 6.1], see also [23]. Even though m = 0 is our case of interest, for our goal we need to control such matrix by the semi-norms of S 0 0,0 . Moreover, for further references, we shall formulate our result in the case of τ -pseudodifferential operators having symbols in the more general class S m 0,0 , m ∈ R. We are going to use the following result for τ -pseudodifferential operators [11,Lemma 4.1]. Lemma 2.1 Fix a non-zero window g ∈ S(R d ) and set τ = W τ (g, g) for τ ∈ [0, 1]. Then, for σ ∈ S R 2d , where z = (z 1 , z 2 ), y = (y 1 , y 2 ) and T τ and J are defined as follows: The Gabor matrix for a τ -pseudodifferential operator Op τ (σ ) with symbol σ ∈ S m 0,0 enjoys the following decay.
Proof Using the representation in (21) and we can write (observe that the above integral is absolutely convergent since τ ∈ S(R 2d )). Now we estimate T τ (z, y)))| ≤ C N |σ | 2N ,m λ m λ − T τ (z, y) −s for every s ≥ 0 since τ ∈ S(R 2d ). Choose s = |m| + 2d + 1. Then the submultiplicativity of · |m| allows us to control from above the right-hand side of the last inequality by Hence, for every N ∈ N we can find C(N ) > 0 such that (22) is satisfied. This concludes the proof.
For fixed τ ∈ (0, 1) we observe that T τ (z, y) z + y , hence the matrix decay can be controlled by a function which does not depend on the τ -quantization. Namely, ∈ (0, 1). Consider a τ -pseudodifferential operator Op τ (σ ) with symbol σ ∈ S m 0,0 . Then, for every N ∈ N, there exists C = C(τ, N ) > 0 such that Remark 2. 2 We conjecture that pseudodifferential operators in the Hörmander class S m 0,0 , m ∈ R, can be characterized via the Gabor matrix in (22), extending the case m = 0 already shown in [18]. To be precise, we allow to write Studying the Gabor matrix decay for M ∞ · m ⊗ · s and following the pattern of the proofs as in the paper [18] one should get the result easily. Since this subject is outside the scope of the paper, we will write the details in a separate work.

Local well-posedness in modulation spaces
For m = 0 the semi-norms on 0 1 are exactly the ones in (15). Observe that (the inclusion is continuous).
The results in the previous yields the boundedness of the Weyl operator e −t H β on modulation spaces.
Then for every T > 0 there exists C = C(T ) > 0 such that Proof Consider u 0 ∈ M p,q w (R d ). Fix g ∈ S(R d )\{0} such that g 2 = 1. Then using the inversion formula for u 0 in (8) we can write In the previous section we showed that e −t H β is a Weyl operator with symbol b(t, ·) in 0 1 with semi-norms uniformly bounded w.r.t. t ∈ [0, T ]. The continuous embedding in (24) and Theorem 2.4 let us write Since w(y) v(y − z)w(z), and v(z) z s for some s > 0, we can write Choosing N such that 2N −s > 2d+1 and using the convolution relations L 1 * L p,q → L p,q we obtain the claim.
Corollary 3.2 Consider 0 < β ≤ 1, s ∈ R. For any fixed T > 0 there exists C = C(T ) > 0 such that The same result holds by replacing the Sobolev space H s with the Shubin-Sobolev space Q s .
As already done in [12], in order to show the local existence of the solution we will make use of the following variant of the contraction mapping theorem (cf., e.g., [26, Proposition 1.38]).

Proposition 3.3 Let N and T be two Banach spaces. Consider a linear operator
for some C 0 > 0, and suppose to have a nonlinear operator F : T → N with F(0) = 0 and Lipschitz bounds for all u, v in the ball B μ := {u ∈ T : u T ≤ μ}, for some μ > 0. Then, for all u lin ∈ B μ/2 there exists a unique solution u ∈ B μ to the equation u = u lin + BF(u), with the map u lin → u Lipschitz continuous with constant at most 2.
Proof of Theorem 1. 1 We apply Theorem 3.1 with T = 1, q = 1 and w(x, ξ) = ξ s . For every 1 ≤ p < ∞, the operator K β (t) in (3)  Observe that similar results were obtained in [21, Theorem 1.1]. We conclude this note by addressing the reader to open problems in this field. First, it is still not clear whether better results can be obtained when considering the nonlinearity F(u) = F k (u) = λ|u| 2k u = λu k+1ūk , λ ∈ C, k ∈ N.
In fact, this was the case for the wave and vibrating place equation, cf. [8,12], where more general modulation spaces were considered.
Moreover, another open question is the well-posedness of the Cauchy problem (1) with initial datum u 0 ∈ M p,q s (R d ), 0 < p ≤ ∞, 0 < q ≤ 1. We conjecture that the result holds true as well, but the techniques employed so far do not apply in this case.