![Arxiv:1902.03413V3 [Math.FA]](http://data.docslib.org/img/004007d96d3936b443f1829b86bd8a8a-1.webp)
arXiv:1902.03413v3 [math.FA] 11 Aug 2020 rnfr,qaiBnc pcs ouainsae,Wee ama Wiener spaces, modulation spaces, quasi-Banach transform, nw ssottm ore rnfr SF) e stog int modulation though us the Let first (STFT). Recall transform Fourier short-time as known opera Anty-Wick as known therein. already references 10, the were 8, and they 7, 4, mechanics 2, quantum t [1, f instance from In mathematics: for and see, of 12] calculus, fields [11, pseudodifferential different to Daubechies in I. authors several by by papers vestigated the with popular become For S (1) f ( S ∈ e od n phrases. and words Key 2010 oaiainoeaoscnb nrdcdvatetm-rqec r time-frequency the via introduced be can operators Localization T tradition. long-standing a has operators localization of study The R z d EA N MOHESFREGNUCIN OF EIGENFUNCTIONS FOR SMOOTHNESS AND DECAY ) ′ V ( = ( { \ ahmtc ujc Classification. Subject Mathematics Abstract. dows A lcto eain o inraagmsae nteqaiBnc set quasi-Banach a the spaces in modulation spaces for amalgam relations convolution Wiener exhibited. new for tool, relations a plication As spaces. ouainspace modulation atlclzto operators localization pact M L e ao tm.Smlry o symbols for Similarly, atoms. Gabor few ϕ R 2 f a ϕ v d ∞ egnutoswt o-eoegnausaeide ihycompres highly indeed are eigenvalues non-zero with -eigenfuctions 0 nipratrl spae yqaiBnc inraagmadmodula and amalgam Wiener quasi-Banach by played is role important An 1 z s ( as ) EEIOBSINN,EEACREO N AI NICOLA FABIO AND CORDERO, ELENA BASTIANONI, FEDERICO ,ϕ } ⊗ 1 ,ω x, z , edfietesottm ore rnfr fatmee distrib tempered a of transform Fourier short-time the define We . 1 ϕ 2 M ( 1 r culyShat functions. Schwartz actually are R 2 ϕ , ) ω = ) 2 d f 2 ∈ ), nteShat class Schwartz the in ( esuydcyadsotns rprisfregnucin of eigenfunctions for properties smoothness and decay study We t > s R h = ) ,π f, 2 d edfietetm-rqec shift time-frequency the define we , sbpcsof (subspaces 0 OAIAINOPERATORS LOCALIZATION M ( e ,ω x, 2 iefeunyaayi,lclzto prtr,sottm Four short-time operators, localization analysis, Time-frequency p, πitω ∞ M ) f ϕ cnann h eegespace Lebesgue the (containing ω ( i A t = a ϕ n translation and ) 1. 1 T , ,ϕ F Introduction 2 S M Operators . ( f r nw ob opc.W hwta their that show We compact. be to known are p, 73;3S5 63;47B10. 46E35; 35S05; 47G30; T ∞ x x ( ϕ 1 f R )( 2 ( d t a ω ), = ) = ) ntewihe ouainspaces modulation weighted the in > p A T ϕ a f 1 x Z ,ϕ 2 ( R prtr fafunction a of operators d/s t 2 d − ihsymbols with f the ) ( x gmspaces. lgam y ) L ) ,x ω, , p ϕ π ), L ( ( 2 y < p z egnucin of -eigenfunctions − 1 2 3 4 5 50]. 45, 44, 43, 32, 31, = ) o inlanalysis signal rom x a ∞ oueteSTFT. the roduce ∈ ) M ntewide the in e n win- and , e noa onto sed R dmulti- nd − z e ieyin- widely hen d igare ting 2 2 epresentation πiyω . T os f [41] cf. tors, z com- 1 tion Fix . e have hey dy. f on ution ϕ R ier d ∈ : 2 FEDERICO BASTIANONI,ELENA CORDERO,ANDFABIO NICOLA
ϕ1,ϕ2 The localization operator Aa with symbol a and windows ϕ1,ϕ2 is formally defined to be
ϕ1,ϕ2 (2) Aa f(t)= a(x, ω)Vϕ1 f(x, ω)MωTxϕ2(t) dxdω. R2d Z 2 −πt ϕ1,ϕ2 If ϕ1(t) = ϕ2(t) = e , then Aa = Aa is the classical Anti-Wick operator and ϕ1,ϕ2 the mapping a 7→ Aa is a quantization rule in quantum mechanics [5, 13, 41, 50]. If one considers a symbol a in the Lebesgue space Lq(R2d) (1 ≤ q < ∞) and 1 d window functions ϕ1,ϕ2 in the Feichtinger’s algebra M (R ) (see below for its defi- ϕ1,ϕ2 nition) then the localization operator Aa is in the Schatten class Sq (see [7]). This ϕ1,ϕ2 2 Rd implies, in particular, that Aa is a bounded and compact operator on L ( ). Sharp results for localization operators on modulation spaces (Banach case) were obtained in [9], thus concluding the open issues related to this problem. The focus of this paper is the properties of eigenfunctions of compact localization operators. The study of eigenvalues and eigenfunctions of a restrict class of compact and ϕ,ϕ self-adjoint localization operators, namely of the type AχΩ , where Ω is a compact domain of the time-frequency plane and the window ϕ is in L2(Rd), was pursued in [1, 2, 3]. The focus of the previous papers, as well as new recent contributions [37, 38], extending the previous results, is the asymptotic behavior of the eigenvalues, depending on the domain Ω. Our perspective is different: properties of eigenfunctions of a compact localization operator having a general symbol a, without any requirement on the geometry of the function a. In particular, the symbol a does not need to have a compact support. ϕ1,ϕ2 Besides, the related localization operator Aa is not necessarily a self-adjoint op- ϕ1,ϕ2 erator (but Aa is compact). It is easy to check that the adjoint of a localization operator is given by ϕ1,ϕ2 ∗ ϕ2,ϕ1 (Aa ) = Aa¯ ;
hence the self-adjointness property forces the choice ϕ1 = ϕ2 and the symbol a real ϕ,ϕ valued, as for the case AχΩ mentioned above. Our framework can allow the use of two different windows ϕ1 and ϕ2 to analyze and synthesize the signal f, respectively. Moreover, the symbol a can be a complex-valued function. For a linear bounded operator T on L2(Rd) we denote by σ(T ) the spectrum of T , that is the set {λ ∈ C : T − λI is not invertible}; in particular, the set σP (T ) denotes the point spectrum of T , that is 2 d σP (T )= {λ ∈ C : ∃ a non-zero f ∈ L (R ) such that T f = λf}. ϕ1,ϕ2 2 Rd Under our assumptions, the operator Aa is a compact mapping on L ( ). Hence, the spectral theory for compact operators yields
ϕ1,ϕ2 ϕ1,ϕ2 σ(Aa ) \{0} = σP (Aa ) \{0}. DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 3
For sake of completeness, we recall that for compact operators on L2(Rd) the number ϕ1,ϕ2 ϕ1,ϕ2 zero is always in the spectrum σ(Aa ) and that the point spectrum σP (Aa )\{0} (possibly empty) is at most countably infinite. 2 Rd ϕ1,ϕ2 A function f ∈ L ( ) is called an eigenfunction of the operator Aa if there exists λ ∈ C such that ϕ1,ϕ2 Aa f = λf. ϕ1,ϕ2 We are interested in the properties of eigenfuctions of Aa related to eigenvalues ϕ1,ϕ2 ϕ1,ϕ2 λ ∈ σP (Aa ) \{0}, whenever σP (Aa ) \{0}= 6 ∅. To chase our goal, the new idea is to use properties of quasi-Banach modulation and Wiener amalgam spaces. Some of these issues are investigated for the first time in this paper. To give a flavour of these results, we recall the definition of modulation spaces (un-weighted case: for the more general cases see below). Fix a non-zero window function g ∈ S(Rd) and p, q ∈ (0, ∞). The modulation space M p,q(Rd) consists of all tempered distributions f ∈ S′(Rd) such that
1 q q p p kfkM p,q := kVgfkLp,q = |Vgf(x, ω)| dx dω < ∞ Rd Rd Z Z ! (with natural modifications when p = ∞ or q = ∞). For p, q ≥ 1 the function k·kM p,q is a norm, and different window functions g ∈ S(Rd) yield equivalent norms, thus the same space. Modulation spaces M p,q(Rd), p, q ≥ 1, are Banach spaces, invented by H. Feichtinger in [19], where many of their properties were already investigated. ϕ1,ϕ2 p,∞ We remark that a localization operator Aa with symbol a ∈ M ,1 ≤ p< ∞, 1 and windows ϕ1,ϕ2 ∈ M , belongs to the Schatten class Sp and, in particular, it is a compact operator on L2(Rd) (cf. [8, Theorem 1]). The (quasi-)Banach space M p,q, p,q > 0, were first introduced and studied by Y.V. Galperin and S. Samarah in [27], see the next section for more details. Roughly speaking, the mapping f is in M p,q(Rd) if it locally behaves like a function in FLq(Rd) and “decays” as a function in Lp(Rd) at infinity. In this paper we extend convolution relations for Banach modulation spaces ex- hibited in [7, 47] to the quasi-Banach ones, see Proposition 3.1 below. This result seems remarkable by itself. Such convolution relations will be crucial for our main result, that can be simplified as follows (cf. Theorem 3.7). Theorem 1.1. Consider a symbol a ∈ M p,∞(R2d), 0
we use Gabor frames. Let Λ = αZd × βZd be a lattice of the time-frequency plane. The set of time-frequency shifts G(g, Λ) = {π(λ)g : λ ∈ Λ}, for a non-zero g ∈ L2(Rd), is called a Gabor system. The set G(g, Λ) is a Gabor frame if there exist constants A,B > 0 such that 2 2 2 2 Rd (3) Akfk2 ≤ |hf, π(λ)gi| ≤ Bkfk2, ∀f ∈ L ( ). Xλ∈Λ If A = B = 1 we call G(g, Λ) a Parseval Gabor frame. In this case (3) reduces to 2 2 2 Rd (4) kfk2 = |hf, π(λ)gi| , ∀f ∈ L ( ). Xλ∈Λ From now on we restrict to Parseval Gabor frames. If the eigenfunction f above γ Rd satisfies f ∈ γ>0 M ( ), then f is highly compressed onto a few Gabor atoms π(λ)g. Indeed, for N ∈ N+, its N-term approximation error σN (f) presents super- polynomial decay.T We define d d (5) ΣN = {p = ck,nπ(αk,βn)g : ck,n ∈ C, F ⊂ Z × Z , card F ≤ N} k,nX∈F (the set of all linear combinations of Gabor atoms consisting of at most N terms). Note that ΣN is not a linear subspace since ΣN + ΣN = Σ2N . That is why the approximation of a signal f by elements of ΣN is often referred to as non-linear approximation. Given a function f ∈ L2(Rd), the N-term approximation error in L2(Rd) is
(6) σN (f) = inf kf − pk2. p∈ΣN
That is, σN (f) is the error produced when f is approximated optimally by a linear combination of N Gabor atoms. γ Rd We shall show in Corollary 3.9 that, if f ∈ γ>0 M ( ), then, for every r > 0 there exists a positive constant C = C(r, f) such that T −r σN (f) ≤ CN .
Since σN (f) is the error produced when f is approximated optimally by a linear combination of N Gabor atoms, the decay above shows the high compression of the eigenfunction onto such atoms, cf. [28, Subsection 12.4]. Another main result (Theorem 3.10) states, in the same spirit, that for symbols in ∞ R2d the weighted modulation space Mvs⊗1( ), s> 0 (see Definition 2.4 below), the cor- 2 ϕ1,ϕ2 Rd responding L eigenfunctions of Aa , with eigenvalues λ =6 0, are actually in S( ). ϕ1,ϕ2 ∞ R2d p,∞ R2d Also in this case Aa is a compact operator, since Mvs⊗1( ) ⊂ M ( ), when- ever p> 2d/s, cf. the subsequent Remark 2.7. These symbol classes include certain measures and the result applies, in particular, to Gabor multipliers (see e.g., [21]). We leave the precise statement to the interested reader. DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 5
In short, the paper is organized as follows. Section 2 is devoted to the function spaces involved in our study. In particular, we show new multiplication relations for Wiener amalgam spaces in the quasi-Banach setting and we prove convolution relations for quasi-Banach modulation spaces. Section 3 represents the core of the paper. We first exhibit continuity results for Weyl operators on modulation spaces (involving the quasi-Banach setting), cf. Theorem 3.3. Then we study the eigen- functions’ properties for Weyl operators (Propositions 3.5 and 3.6 below). Next we show our main result on the eigenfunctions’ regularity and smoothness of local- ization operators: Theorem 3.7 and related consequences in terms of the N-term approximation (Proposition 3.8 and Corollary 3.9). The last Section 4 is devoted to the study of eigenfunctions of localization operators with symbols in the weighted q Rd Lebesgue spaces Lm( ), 1 ≤ q < ∞. This section suggests a wider study of the topic for localization operators on groups [50]. We shall pursue this issue in a subse- quent paper. Let us conclude this introduction with another open problem, that can be worth investigating. Namely, in this paper we work with compact localization op- ∞ R2d p,∞ R2d erators having symbols a in Mvs⊗1( ), for s> 0, or in the larger space M ( ), 0 2. Preliminaries Notation. We define xy = x · y, the scalar product on Rd. The Schwartz class is denoted by S(Rd), the space of tempered distributions by S′(Rd). We use the brackets hf,gi to denote the extension to S′(Rd) × S(Rd) of the inner product hf,gi = f(t)g(t)dt on L2(Rd). The Fourier transform of a function f on Rd is normalized as R Ff(ω)= e−2πixωf(x) dx. Rd Z The involution g∗ is given by g∗(t) = g(−t). Given a measurable and positive R2d p Rd weight function m on and 1 ≤ p< ∞, we denote by Lm( ) the Banach space of measurable functions f : Rd → C satisfying 1/p p p p kfkLm = |f(x)| m(x) dx , Rd Z up to the equivalence relation f ∼ g if and only if f(x) = g(x) for a.e. x. For ∞ d R ∞ d p = ∞, we have f ∈ Lm ( ) if kfkLm = ess supx∈R |f(x)|m(x) < ∞, again up to 6 FEDERICO BASTIANONI,ELENA CORDERO,ANDFABIO NICOLA the equivalence relation above. For s ∈ R, recall the Sobolev spaces Hs(Rd)= {f : fˆ(ω)hωis ∈ L2(Rd)}, where hωis =(1+ |ω|2)s/2. We use T ∗ for the adjoint of an operator T . Observe that π(z)∗ = π(−z) and the following commutation relations hold (7) π(z)π(w)= e−2πiz1w2 π(w)π(z). For f,g ∈ L2(Rd), the cross-Wigner distribution is defined by: t t (8) W (f,g)(x, ω)= f(x + )g(x − )e−2πitω dt. Rd 2 2 Z Given two normed spaces A and B, we denote by A ֒→ B the continuous embedding of A into B. Finally we recall the definition of the Schatten classes Sp,0 2.1. Weight functions. In the sequel v will always be a continuous, positive, sub- d d multiplicative weight function on R (or on Z ), i.e., v(z1 + z2) ≤ v(z1)v(z2), for all d d d d z1, z2 ∈ R (or for all z1, z2 ∈ Z ). We say that m ∈Mv(R ) (or m ∈Mv(Z )) if m d d is a positive, continuous weight function on R (or on Z ) v-moderate: m(z1 + z2) ≤ d d Cv(z1)m(z2) for all z1, z2 ∈ R (or for all z1, z2 ∈ Z ). We will mainly work with polynomial weights of the type s 2 s/2 d d (9) vs(z)= hzi =(1+ |z| ) , s ∈ R, z ∈ R (or Z ). Observe that, for s< 0, vs is v|s|-moderate. d d Given two weight functions m1, m2 on R (or Z ), we write d d (m1 ⊗ m2)(x, ω)= m1(x)m2(ω), x,ω ∈ R (or Z ). 2.2. Spaces of sequences. Z2d p,q Z2d Definition 2.1. For 0 < p,q ≤∞, m ∈Mv( ), the space ℓm ( ) consists of all sequences a =(ak,n)k,n∈Zd for which the (quasi-)norm 1 q q p p p kak p,q = |a | m(k, n) ℓm k,n Zd Zd ! nX∈ kX∈ (with obvious modification for p = ∞ or q = ∞) is finite. DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 7 p,q Z2d p Z2d For p = q, ℓm ( ) = ℓm( ), the standard spaces of sequences. Namely, in Zd dimension d, for 0 p p p kakℓm = |ak| m(k) < ∞. Zd ! kX∈ Here there are some properties we need in the sequel [26, 27]: p1 Zd p2 Zd i) Inclusion relations: If0 < p1 ≤ p2 ≤ ∞, then ℓm ( ) ֒→ ℓm ( ), for any) positive weight function m on Zd. d (ii) Young’s convolution inequality: Consider m ∈Mv(Z ), 0 < p,q,r ≤∞ with 1 1 1 (10) + =1+ , for 1 ≤ r ≤∞ p q r and (11) p = q = r, for 0 r p q (12) ka ∗ bkℓm ≤ Ckakℓm kbkℓv , where C is independent of p,q,r, a and b. If m ≡ v ≡ 1, then C = 1. (iii) H¨older’s inequality: For any positive weight function m on Zd, 0 < p,q,r ≤ ∞, with 1/p +1/q =1/r, p Zd q Zd r Zd ,( ) ℓm( ) · ℓ1/m( ) ֒→ ℓ (13) .where the symbol ֒→ denotes that the inclusion is a continuous mapping 2.3. Wiener Amalgam Spaces [16, 17, 18, 25, 27, 34, 39, 40]. Definition 2.2. Consider p, q ∈ (0, ∞], a weight function m ∈Mv and the compact d p q Rd set Q = [0, 1] . The Wiener amalgam space W (L , Lm)( ) consists of the functions Rd C p Rd f : → such that f ∈ Lloc( ) and for which the control function: Q q Zd Zd (14) Ff (k) := kf · TkχQkLp ∈ ℓm , k ∈ . p q The quasi-norm on W (L , Lm) is given by Q kfk p q := F (k) W L ,Lm f q ( ) ℓm = kf · Tkχ Qk p q L ℓm 1 q q p p q (15) = |f (t)| χQ (t − k) dt m (k) , Rd Zd ! kX∈ Z with suitable adjustments for the cases p, q = ∞. 8 FEDERICO BASTIANONI,ELENA CORDERO,ANDFABIO NICOLA This special definition allows us to grasp the sense of the amalgam: we first view f “locally” through translations TkχQ of the sharp cutoff function χQ, and measure those local pieces in the Lp-norm, then we measure the global behavior of those local q pieces according to the ℓm-norm. The “window” through which we view f locally need not be a unit d-dimensional cube, cf. [17, 27, 34, 40]. In the sequel we shall use the following properties: (i) Inclusion relations: For 0 (iii) For m ∈Mv, 0 Proposition 2.3 (Multiplication relations). Consider m, w ∈Mv, 0 p1 q1 p2 q2 p3 q3 .( W (L , Lm) · W (L , Lw/m) ֒→ W (L , Lw (19) Proof. The result is well known for 1 ≤ pi, qi ≤ ∞, cf. [16, 34]. Here we show that the same proof works for quasi-Banach spaces. Indeed, since the standard p1 q1 H¨older inequality holds for Lebesgue exponents in (0, +∞], for f1 ∈ W (L , Lm), p2 q2 f2 ∈ W (L , Lw/m) we have kf1f2TkχQkLp3 = k(f1TkχQ)(f2TkχQ)kLp3 ≤kf1TkχQkLp1 kkf2TkχQkLp2 . Defining ak = kf1TkχQkp1 and bk = kf2TkχQkp2 and using H¨older’s inequality for q1 q2 q3 sequences ℓ ℓ ֒→ ℓ , for 1/q1 +1/q2 =1/q3 (0 < qi ≤∞, i =1, 2, 3), we obtain kakbkw(k)kℓq3 = k(akm(k))(bkw(k)/m(k))kℓq3 ≤kakm(k)kℓq1 kbkw(k)/m(k)kℓq2 . This completes the proof. 2.4. Modulation Spaces. We use the extension to quasi-Banach spaces introduced first by Y.V. Galperin and S. Samarah in [27]. d Definition 2.4. Fix a non-zero window g ∈ S(R ), a weight m ∈ Mv and 0 < p,q Rd p, q ≤ ∞. The modulation space Mm ( ) consists of all tempered distributions f ∈ S′(Rd) such that the (quasi-)norm 1 q q p p,q p,q p p (20) kfkMm = kVgfkLm = |Vgf(x, ω)| m(x, ω) dx dω Rd Rd Z Z ! DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 9 (obvious changes with p = ∞ or q = ∞) is finite. p,q Rd The most famous modulation spaces are those Mm ( ) with 1 ≤ p, q ≤ ∞, invented by H. Feichtinger in [19]. In that paper he proved they are Banach spaces, whose norm does not depend on the window g, in the sense that different window functions in S(Rd) yield equivalent norms. Moreover, the window class S(Rd) can 1,1 Rd be extended to the modulation space Mv ( ) (so-called Feichtinger algebra). p Rd p,p Rd p,q Rd For shortness, we write Mm( ) in place of Mm ( ) and M ( ) if m ≡ 1. p,q Rd The modulation spaces Mm ( ), 0 < p,q < 1, where introduced almost twenty years later by Y.V. Galperin and S. Samarah in [27] and then studied in [35, 39, 48, 49] (see also references therein). In this framework, it appears that the largest p,q Rd natural class of windows universally admissible for all spaces Mm ( ), 0 < p,q ≤∞ (with m having at most polynomial growth) is the Schwartz class S(Rd). There are thousands of papers involving modulation spaces with indices 1 ≤ p, q ≤ ∞, whereas very few works deal with the quasi-Banach case 0 < p,q ≤ 1. Indeed, many properties related to the latter case are still unexplored. In this paper their contribution is fundamental, since they are the key tool for un- derstanding the properties of eigenfunctions of localization operators having symbols with some decay at infinity, measured in the Lp-mean, 0 Proposition 2.1. Consider m ∈Mv, r ∈ (0, 1], 0 0: ∞ R2d p,∞ R2d Mvs⊗1( ) ⊂ M ( ) if p> 2d/s. 10 FEDERICO BASTIANONI,ELENA CORDERO, AND FABIO NICOLA This follows from the recent contribution [33, Theorem 1.5]. Hence localization oper- ϕ1,ϕ2 ∞ R2d Rd ators Aa with symbols a in Mvs⊗1( ) and windows ϕ1,ϕ2 in S( ) are compact operators belonging to the Schatten class Sp, with p> 2d/s. We will repeatedly use the following result, cf. [27, Theorem 3.3] (see also [28, Theorem 12.2.1] for p ≥ 1). Theorem 2.8. Assume that m ∈ Mv. For 0 ∞ p p (21) kVgfkW (L ,Lm) ≤ CkVgfkLm . We also need to recall the inversion formula for the STFT (see [28, Proposition 1 Rd p,q Rd 11.3.2]): assume g ∈ Mv ( ) \{0}, 1 ≤ p, q ≤∞, f ∈ Mm ( ), then 1 (22) f = 2 Vgf(x, ω)MωTxg dx dω, kgk R2d 2 Z p,q Rd and the equality holds in Mm ( ). 2.5. Gabor Frames. Let Λ = αZd × βZd be a lattice of the time-frequency plane. The set of time-frequency shifts G(g, Λ) = {π(λ)g : λ ∈ Λ}, for a non-zero g ∈ L2(Rd), is called a Gabor system. The set G(g, Λ) is a Gabor frame, if there exist constants A,B > 0 such that (3) holds true. If G(g, Λ) is a Gabor frame, then it can be shown that the frame operator Sf = hf, π(λ)giπ(λ)g, f ∈ L2(Rd) Xλ∈Λ is a topological isomorphism on L2(Rd). Moreover, the system G(γ, Λ), where the function γ = S−1g ∈ L2(Rd) is the canonical dual window of g, is a Gabor frame and we have the reproducing formula (23) f = hf, π(λ)giπ(λ)γ = hf, π(λ)γiπ(λ)g, f ∈ L2(Rd) Xλ∈Λ Xλ∈Λ 2 d with unconditional convergence in L (R ). In particular, if γ = g and kgk2 = 1 then A = B = 1, the frame operator is the identity S = I and the Gabor frame is called Parseval Gabor frame. Observe that formula (4) holds true. More generally, any window function γ ∈ L2(Rd), such that (23) is satisfied, is called alternative dual window for g. In general, given two functions g,γ ∈ L2(Rd), it is customary to extend the notion of Gabor frame operator Sg,γ, related to g,γ, as follows 2 d Sg,γf = hf, π(λ)giπ(λ)γ, f ∈ L (R ), Xλ∈Λ DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 11 whenever the previous operator is well-defined. With this notation the reproducing 2 d formula (23) can be rephrased as Sg,γ = I on L (R ), with I being the identity operator. Modulation spaces provide a natural setting for time-frequency analysis, thanks to discrete equivalent norms produced by means of Gabor frames. The key result is the following (see [28, Chapter 12] for 1 ≤ p, q ≤ ∞, and [27, Theorem 3.7] for 0 2d d d d Theorem 2.9. Assume m ∈ Mv(R ), Λ = αZ × βZ , g,γ ∈ S(R ) such that 2 d Sg,γ = I on L (R ). Then p,q Rd (24) f = hf, π(λ)giπ(λ)γ = hf, π(λ)γiπ(λ)g, f ∈ Mm ( ) Xλ∈Λ Xλ∈Λ p,q Rd with unconditional convergence in Mm ( ) if 0 < p,q < ∞ and with weak-* con- ∞ Rd vergence in M1/v( ) otherwise. Furthermore, there are constants 0 < A ≤ B such p,q Rd that, for all f ∈ Mm ( ), 1 q q p p p (25) Akfk p,q ≤ |hf, π(αk,βn)gi| m(αk,βn) ≤ Bkfk p,q , Mm Mm Zd Zd ! nX∈ kX∈ independently of p, q, and m. Similar inequalities hold with g replaced by γ. In other words, p,q d p,q p,q kfkMm (R ) ≍k(hf, π(λ)gi)λkℓm (Λ) = k(Vgf(λ))λkℓm (Λ). 3. Main Results We first study convolution relations for modulations spaces. Let us recall that, for the Banach cases, convolution relations were studied in [7] and [46, 47]. Our approach is general, the techniques use Gabor frames via the equivalence (25), plus H¨older’s and Young’s inequalities for sequences. Proposition 3.1. Let ν(ω) > 0 be an arbitrary weight function on Rd, 0 < p,q,r,t,u,γ ≤ ∞, with 1 1 1 (26) + = , u t γ and 1 1 1 (27) + =1+ , for 1 ≤ r ≤∞ p q r whereas (28) p = q = r, for 0 2d For m ∈ Mv(R ), m1(x) = m(x, 0) and m2(ω) = m(0,ω) are the restrictions to Rd ×{0} and {0} × Rd, and likewise for v. Then p,u d q,t d r,γ d R −1 R R ( ) Mm1⊗ν( ) ∗ Mv1⊗v2ν ( ) ֒→ Mm (29) with norm inequality kf ∗ hk r,γ . kfk p,u khk q,t . Mm Mm ⊗ν M −1 1 v1⊗v2ν Proof. We use the key idea in [7, Proposition 2.4] to measure the modulation space −πx2 −d/2 −πx2/2 norm with respect to the Gaussian windows g0(x)= e and g(x)=2 e = d (g0 ∗ g0)(x) ∈ S(R ). −2πixω ∗ A straightforward computation shows Vgf(x, ω) = e (f ∗ Mωg )(x) (recall that g∗(x)= g(−x)). Thus, using the same argument as in the proof [7, Proposition ∗ ∗ ∗ ∗ 2.4], we recall the easily verified identity Mω(g0 ∗ g0)= Mωg0 ∗ Mωg0 and write the STFT of f ∗ h as follows: −2πixω ∗ −2πixω ∗ ∗ Vg(f ∗ h)(x, ω)= e (f ∗ h) ∗ Mωg (x)= e (f ∗ Mωg0) ∗ (h ∗ Mωg0) (x) . In the sequel, we first use the norm equivalence (25), written in terms of the STFT p,q p,q Zd Zd as kfkMm ≍k(Vgf(λ))λ∈Λkℓm (Λ), where Λ = α × β . Then we majorize m by m(αk,βn) . m(αk, 0)v(0, βn)= m1(αk)v2(βn), and finally use Young’s convolution inequality for sequences in the k-variable and H¨older’s one in the n-variable. The indices p,q,r,γ,t,u fulfil the equalities in the assumptions. We show in details the case when r,γ,t,u < ∞: r,γ r,γ 2d kf ∗ hkMm ≍k((Vg(f ∗ h))(αk,βn)m(αk,βn))k,nkℓ (Z ) γ/r 1/γ . |(f ∗ M g∗) ∗ (h ∗ M g∗)(αk)|rm (αk)r v (βn)γ βn 0 βn 0 1 2 Zd Zd ! nX∈ kX∈ 1/γ ∗ ∗ γ γ = k(f ∗ Mβng0) ∗ (h ∗ Mβng0)k r Zd v2(βn) ℓm1 (α ) Zd ! nX∈ 1/γ ∗ γ ∗ γ γ . kf ∗ Mβng0k p Zd kh ∗ Mβng0k q Zd v2(βn) ℓm1 (α ) ℓv1 (α ) n∈Zd ! X 1 1 u t t ∗ u u ∗ t v2(βn) . kf ∗ Mβng0k p Zd ν(βn) kh ∗ Mβng0k q Zd ℓm1 (α ) ℓv1 (α ) ν(βn)t Zd ! Zd ! nX∈ nX∈ p,u = k((Vg0 f)(λ))λk k((Vg0 h)(λ))λk q,t ℓm ⊗ν (Λ) ℓ −1 (Λ) 1 v1⊗v2ν ≍kfk p,u khk q,t . Mm ⊗ν M −1 1 v1⊗v2ν DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 13 The cases when one among the indexes r,γ,t,u is equal to ∞ are done similarly. This concludes the proof. 3.1. Weyl Operators. Every continuous operator from S(Rd) to S′(Rd) can be represented as a pseudodifferential operator in the Weyl form Lσ, with Weyl symbol σ ∈ S′(R2d). The operator is formally defined by x + y 2πi(x−y) ω (30) Lσf(x)= σ ,ω e f(y) dydω. R2d 2 Z The crucial relation between the action of the Weyl operator Lσ on time-frequency shifts and the short-time Fourier transform of its symbol is contained in [29, Lemma 3.1]. Lemma 3.1. Consider g ∈ S(Rd), Φ= W (g,g). Then, for σ ∈ S′(R2d), z + w (31) |hL π(z)g, π(w)gi| = V σ , j(w − z) , z,w ∈ R2d, σ Φ 2 where j(z1, z2)=(z2, −z1). We first recall Schatten class results for the Weyl calculus in terms of modulation spaces, initially proved for 1 ≤ p ≤∞ in [6, Theorem 4.5], for 0 p,1 kLσkSp ≤kσkM . 2 d In particular, Lσ is a compact operator on L (R ). Theorem 3.3. (i) Assume p,q,γ ∈ (0, ∞] such that 1 1 1 (32) + = . p q γ p,min{1,γ} 2d d If σ ∈ M (R ), then the pseudodifferential operator Lσ, from S(R ) to S′(Rd), extends uniquely to a bounded operator from M q(Rd) to M γ (Rd). ∞,1 R2d (ii) If s,r ≥ 0, t ≥ r + s, and the symbol σ ∈ Mvs⊗vt ( ), then the pseudodiffe- d ′ d rential operator Lσ, from S(R ) to S (R ), extends uniquely to a bounded operator 2 Rd 2 Rd from Mvr ( ) into Mvr+s ( ). Proof. (i) Assume γ ≥ 1, then, by (32), p ≥ γ ≥ 1 and q ≥ γ ≥ 1. The claim was proved by Toft in [46, Theorem 4.3]. The case γ < 1, p, q ∈ (0, ∞] was again proved by Toft in [48, Theorem 3.1]. d (ii) Let g ∈ S(R ) with kgkL2 = 1. From the inversion formula (22), Vg(Lσf)(w)= hLσπ(z)g, π(w)gi Vgf(z) dz. R2d Z 14 FEDERICO BASTIANONI,ELENA CORDERO, AND FABIO NICOLA The desired result thus follows if we can prove that the map M(σ) defined by M(σ)G(w)= hLσπ(z)g, π(w)gi G(z) dz R2d Z 2 R2d 2 R2d is continuous from Lvr ( ) into Lvr+s ( ). Using (31), we see that it is sufficient to prove that the integral operator with integral kernel z + w V σ , j(w − z) hzi−rhwir+s Φ 2 is bounded on L2(R2d). This follows from Schur’s test (see, e.g., [28, Lemma 6.2.1 ∞,1 (b)]). Indeed, by assumption σ ∈ Mvs⊗vt , so that z + w s t sup VΦσ , j(w − z) hz + wi hw − zi dz < ∞ R2d R2d 2 w∈ Z and similarly z + w s t sup VΦσ , j(w − z) hz + wi hw − zi dw < ∞. R2d R2d 2 z∈ Z Hence it is sufficient to prove that for some positive constant C > 0 we have (33) hz + wi−shw − zi−thzi−rhwir+s ≤ C, ∀w, z ∈ R2d. Let us prove the estimate (33). Setting x = z + w, y = w − z, the inequality (33) can be rephrased as (34) hxi−shyi−thx − yi−rhx + yir+s ≤ C, ∀x, y ∈ R2d. For |x| < 2|y|, observe that |x+y| < 3|y| and since t ≥ r+s we get the estimate (34). For |x| ≥ 2|y|, we use hx + yi≍hx − yi≍hxi and (34) immediately follows. We remark that 2 Rd Rd Mvs ( )= Qs( ), the Shubin-Sobolev spaces, cf. [6, 41]. In particular, for s = 0, M 2(Rd) = L2(Rd). Thus, we recover the known result, c.f. [47, Theorem 4.3]: ∞,1 R2d Corollary 3.4. If s,r ≥ 0, t ≥ r + s, and the symbol σ ∈ Mvs⊗vt ( ), then the d ′ d pseudodifferential operator Lσ, from S(R ) to S (R ), extends uniquely to a bounded d d operator from Qr(R ) into Qr+s(R ). An application of the previous theorem concerns the study of eigenfunctions’ properties for Weyl operators. Observe that a Weyl operator having symbol σ ∈ M p,γ(R2d) with p< ∞ is a compact operator on L2(Rd) by Theorem 3.2. DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 15 Proposition 3.5. Consider a Weyl symbol σ ∈ M p,γ for some 0 (observe γ2 < γ1 since p < ∞). Continuing this way we construct a decreasing sequence of indices γn > 0 and such that limn→∞ γn = 0. This proves the claim. Recall that, for s,t > 0, by inclusion relations for modulation spaces we have ∞,1 R2d p,1 R2d Mvs⊗vt ( ) ⊂ M ( ), for any p > 2d/s, [33, Theorem 1.5 and Lemma 4.10]. Hence the operators below are compact operators by Theorem 3.2. ∞,1 R2d Proposition 3.6. Consider a Weyl symbol σ ∈ Mvs⊗vt ( ) for some s > 0 and 2 d every t> 0. If λ ∈ σP (Lσ)\{0}, then any eigenfunction f ∈ L (R ) with eigenvalue λ is in S(Rd). ∞,1 R2d Proof. By Theorem 3.3, if the symbol σ is in Mvs⊗vt ( ), for some s> 0 and every 2 Rd 2 Rd Rd t> 0, then the Weyl operator acts continuously from L ( ) into Mvs ( )= Qs( ). d Starting now with the eigenfunction f in Qs(R ) and repeating the same argument d with t ≥ s we obtain that the eigenfunction is in Q2s(R ). Proceeding this way we Rd infer that f ∈ ∩n∈N+ Qns( ). The inclusion relations for Shubin-Sobolev spaces and the property (see e.g., [7, 20, 47]) d d ∩r>0Qr(R )= S(R ), prove the claim. 3.2. Localization Operators. The study of eigenfunctions for a localization op- ϕ1,ϕ2 erator Aa uses its representation as a Weyl one: ϕ1,ϕ2 Aa = La∗W (ϕ2,ϕ1) (cf. [7] and references therein) where W (ϕ2,ϕ1) is the the cross-Wigner distribution ϕ1,ϕ2 defined in (8). Therefore, the Weyl symbol of Aa is given by (35) σ = a ∗ W (ϕ2,ϕ1). ϕ1,ϕ2 We will deduce properties for Aa using its Weyl form La∗W (ϕ2,ϕ1), as detailed be- low. Precisely, we shall focus on properties of eigenfunctions of compact localization operators. 16 FEDERICO BASTIANONI,ELENA CORDERO, AND FABIO NICOLA Theorem 3.7. Consider a symbol a ∈ M p,∞(R2d), for some 0 d 2d Proof. Since the windows ϕ1,ϕ2 ∈ S(R ), the cross-Wigner distributionT is in S(R ) ⊂ M q,γ (R2d), for every 0 < q,γ ≤ ∞. We next apply the convolution relations for modulation spaces (29). Namely, if p ≥ 1, choose q = 1, so that r = p, and p,∞ 2d 1,γ 2d p,γ 2d .σ = a ∗ W (ϕ2,ϕ1) ∈ M (R ) ∗ M (R ) ֒→ M (R ), ∀γ > 0 If 0 p,∞ 2d p,γ 2d p,γ 2d .σ = a ∗ W (ϕ2,ϕ1) ∈ M (R ) ∗ M (R ) ֒→ M (R ), ∀γ > 0 In both cases we obtain that σ ∈ M p,γ(R2d), for every γ > 0. Hence the claim immediately follows by Proposition 3.5. As a consequence, the eigenfunctions are extremely concentrated on the time- frequency space, having very few Gabor coefficients large whereas all the others are negligible. Consider a Parseval Gabor frame G(g, Λ) for L2(Rd), with Λ = αZd × βZd and d g ∈ S(R ). For N ∈ N+, we defined in (5) the space ΣN to be the set of all linear combinations of Gabor atoms consisting of at most N terms. Given a function f ∈ L2(Rd), the N-term approximation error in L2(Rd) is recalled in (6). Namely, σN (f) is the error produced when f is approximated optimally by a linear combination of N Gabor atoms. Assume f ∈ M p(Rd) for some 0 |cι(1)|≥|cι(2)|≥···≥|cι(m)|≥|cι(m+1)|≥··· The sequence (˜cm)m∈N+ = (|cι(m)|)m∈N+ is called the non-increasing rearrangement of (ck,n)k,n above. With this notations, denoting by popt the best approximation of f in ΣN , the N-term approximation error becomes N (36) σN (f) = inf kf − pk2 = kf − poptk2 ≤ f − cι(m)π(ι(m))g p∈ΣN m=1 2 1 X ∞ 2 2 ≤ C |cι(m)| . ! m=XN+1 DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 17 By abuse of notation, given a = (am)m, am ≥ 0 for every m, a non-increasing sequence (a1 ≥ a2 ≥···≥ am ≥ am+1 ≥··· ) we write 1 ∞ 2 2 σN (a)= am . m=N+1 ! X The key tool is now the following lemma, see [42] and [15] (we also refer to [28, Lemma 12.4.1]): Lemma 3.2. Let a =(am)m, am ≥ 0 for every m, be a non-increasing sequence and consider 0 kfkM p ≍k(hf, π(αk,βn)gi)k,nkℓp(Λ) and the sequence (|hf, π(αk,βn)gi|)k,n can be rearranged in a non-increasing one (am), as explained above. Applying Lemma 3.2 to such a sequence, from the right- hand side inequality in (38) and using the majorization in (36) we infer (39). Corollary 3.9. Consider a Parseval Gabor frame G(g, Λ), with g ∈ S(Rd), and Λ= αZd ×βZd, α,β > 0. Under the assumptions of Theorem 3.7, any f eigenfunction of ϕ1,ϕ2 Aa (with eigenvalue λ =06 ) is highly compressed onto a few Gabor atoms π(λ)g, in the sense that its N-term approximation error satisfies the following property: for every r> 0 there exists C = C(r, f) > 0 with −r (40) σN (f) ≤ CN . Proof. By Theorem 3.7, the eigenfunction fulfils f ∈ M p(Rd), for every p > 0. Hence the assumptions of Proposition 3.8 are satisfied for every 0 We next consider the case of compact localization operators with symbols a ∈ ∞ R2d 2 Mvs⊗1( ), s > 0. In this case L eigenfunctions reveal to be Schwartz functions, as shown below. ∞ R2d Theorem 3.10. Consider a symbol a ∈ Mvs⊗1( ), for some s> 0, and non-zero Rd ϕ1,ϕ2 ϕ1,ϕ2 windows ϕ1,ϕ2 ∈ S( ). Assume that σP (Aa )\{0}= 6 ∅. If λ ∈ σP (Aa )\{0}, any eigenfunction f ∈ L2(Rd) with eigenvalue λ belongs to the Schwartz class S(Rd). Rd R2d 1,1 R2d Proof. The assumption ϕ1,ϕ2 ∈ S( ) implies W (ϕ2,ϕ1) ∈ S( ) ⊂ Mvr⊗vt ( ), for every r,t > 0. We next apply the convolution relations for modulation spaces ϕ1,ϕ2 ∞,1 R2d (29), obtaining that Aa = Lσ with σ ∈ Mvs⊗vt ( ), for some s > 0 and every t> 0. Hence the claim immediately follows by Proposition 3.6. Remark 3.11. The nice properties of eigenfunctions for localization operators stud- ied so far seem to depend on the fact that such operators are not only compact but belong to the Schatten class Sp, 0 (41) lim sup Vg0 a(x, ω)=0. |x|→∞ |ξ|≤R It seems that for symbols satisfying condition (41) the techniques developed above do not work anymore. It would be very interesting to know whether for compact 2 operators that are not in the Schatten class Sp, 0 4. Symbols in Lp(R2d) spaces We now consider localization operators with symbols in weighted Lebesgue spaces. ϕ1,ϕ2 Rd Let us recall that any localization operator Aa with windows in S( ) and symbol a ∈ Lq(R2d), with 1 ≤ q < ∞, is a compact operator, cf. [50, Proposition 13.3]. The case of weighted Lebesgue spaces and, more generally, Potential Sobolev spaces ϕ1,ϕ2 was treated in [6]: let us stress that any localization operator Aa with Schwartz q R2d windows and symbol a in Lm( ), with 1 ≤ q < ∞ is a compact operator on L2(Rd). R2d q R2d Theorem 4.1. Let m ∈Mv, m(z) ≥ 1 for every z ∈ , a ∈ Lm( ), 1 ≤ q < ∞, Rd ϕ1,ϕ2 and non-zero windows ϕ1,ϕ2 ∈ S( ). Assume that σP (Aa ) \ {0}= 6 ∅. If ϕ1,ϕ2 2 Rd λ ∈ σP (Aa ) \ {0}, any eigenfunction f ∈ L ( ) with eigenvalue λ satisfies p Rd f ∈ p>0 Mm( ). T DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 19 q R2d Proof. By assumption and using (18), we start with a symbol a in Lm( ) = q q R2d 2 Rd W (L , Lm)( ). Consider the eigenvector f ∈ L ( ) and the window ϕ1 ∈ Rd S( ). Then by Theorem 2.8 the STFT Vϕ1 f is in the Wiener amalgam space ∞ 2 R2d q p1 R2d W (L , L )( ). Proposition 2.3 yields that aVϕ1 f ∈ W (L , Lm )( ), with 1 1 1 + = , q 2 p1 so that the index p1 satisfies p1 < min{q, 2}. Consider now a non-zero window g ∈ S(Rd). Using (7), ϕ1,ϕ2 ϕ1,ϕ2 Vg(Aa f)(w)= hAa f, π(w)gi = (aVϕ1 f)(z)hπ(z)ϕ2, π(w)gi dz R2d Z = (aVϕ1 f)(z)hϕ2, π(−z)π(w)gi dz R2d Z −2πiz1w2 = (aVϕ1 f)(z)e hϕ2, π(w − z)gi dz R2d Z so that, ϕ1,ϕ2 (42) |Vg(Aa f)(w)| ≤ |(aVϕ1 f)(z)||Vgϕ2(w − z)|dz = |aVϕ1 f|∗|Vgϕ2|(w). R2d Z We estimate ϕ1,ϕ2 q p1 R2d ∞ 1 R2d |Vg(Aa f)(w)|≤|aVϕ1 f|∗|Vgϕ2|(w) ∈ W (L , Lm )( ) ∗ W (L , Lv)( ). ∞ 1 R2d q′ 1 R2d -Observing that W (L , Lv)( ) ֒→ W (L , Lv)( ) and applying the convolu ϕ1,ϕ2 ∞ p1 p1 tion relations (17) we infer |Vg(Aa f)| ∈ W (L , Lm ) ֒→ Lm . This proves that ϕ1,ϕ2 p1 Rd Aa f ∈ Mm ( ). ϕ1,ϕ2 p1 Rd Recalling the assumption Aa f = λf, λ =6 0, we infer f ∈ Mm ( ). p1 Rd We now repeat the previous argument starting with f ∈ Mm ( ). By Theorem ∞ p1 R2d q p2 q p2 ,( the STFT Vϕ1 f ∈ W (L , Lm )( ) and aVϕ1 f ∈ W (L , Lm2 ) ֒→ W (L , Lm 2.8 (since m2 ≥ m), with 1 1 1 + = , q p1 p2 ϕ1,ϕ2 ∞ p2 p2 . so that p2 < p1. Arguing as above we infer |Vg(Aa f)(w)| ∈ W (L , Lm ) ֒→ Lm p2 Thus, the eigenfunction f belongs to the smaller space Mm . Continuing this way we construct a strictly decreasing sequence of indices pn > 0 and such that lim pn =0. n→∞ By induction and using the same argument as above one immediately obtains that pn Rd pn+1 Rd if f ∈ Mm ( ) then f ∈ Mm ( ). This concludes the proof. 20 FEDERICO BASTIANONI,ELENA CORDERO, AND FABIO NICOLA Acknowledgments The last two authors were partially supported by the Gruppo Nazionale per l’Analisi Matematica, la Probabilit`ae le loro Applicazioni (GNAMPA) of the Is- tituto Nazionale di Alta Matematica (INdAM). References [1] L. D. Abreu and M. D¨orfler. An inverse problem for localization operators. Inverse Problems, 28(11):115001, 16, 2012. [2] L. D. Abreu, K. Gr¨ochenig, and J. L. Romero. On accumulated spectrograms. Trans. Amer. Math. Soc., 368(5):3629–3649, 2016. [3] L. D. Abreu, J. a. M. Pereira, and J. L. Romero. Sharp rates of convergence for accumulated spectrograms. Inverse Problems, 33(11):115008, 12, 2017. [4] D. Bayer and K. Gr¨ochenig. Time-frequency localization operators and a Berezin transform. Integral Equations Operator Theory, 82(1):95–117, 2015. [5] F. A. Berezin. Wick and anti-Wick symbols of operators. Mat. Sb. (N.S.), 86(128):578–610, 1971. [6] P. Boggiatto, E. Cordero, and K. Gr¨ochenig. Generalized anti-Wick operators with symbols in distributional Sobolev spaces. Integral Equations Operator Theory, 48(4):427–442, 2004. [7] E. Cordero and K. Gr¨ochenig. Time-frequency analysis of localization operators. J. Funct. Anal., 205(1):107–131, 2003. [8] E. Cordero and K. Gr¨ochenig. Necessary conditions for Schatten class localization operators. Proc. Amer. Math. Soc., 133(12):3573–3579, 2005. [9] E. Cordero and F. Nicola. Sharp integral bounds for Wigner distributions. Int. Math. Res. Not. IMRN, (6):1779–1807, 2018. [10] E. Cordero, S. Pilipovi´c, L. Rodino, and N. Teofanov. Localization operators and exponential weights for modulation spaces. Mediterr. J. Math., 2(4):381–394, 2005. [11] I. Daubechies. Time-frequency localization operators: a geometric phase space approach. IEEE Trans. Inform. Theory, 34(4):605–612, 1988. [12] I. Daubechies and T. Paul. Time-frequency localization operators—a geometric phase space approach. II. The use of dilations. Inverse Problems, 4(3):661–680, 1988. [13] M. A. de Gosson. Symplectic methods in harmonic analysis and in mathematical physics, vol- ume 7 of Pseudo-Differential Operators. Theory and Applications. Birkh¨auser/Springer Basel AG, Basel, 2011. [14] M. A. de Gosson. The canonical group of transformations of a Weyl-Heisenberg frame; appli- cations to Gaussian and Hermitian frames. J. Geom. Phys., 114:375–383, 2017. [15] R. A. DeVore and V. N. Temlyakov. Some remarks on greedy algorithms. Adv. Comput. Math., 5(2-3):173–187, 1996. [16] H. G. Feichtinger. Banach spaces of distributions of Wiener’s type and interpolation. In Func- tional analysis and approximation (Oberwolfach, 1980), volume 60 of Internat. Ser. Numer. Math., pages 153–165. Birkh¨auser, Basel-Boston, Mass., 1981. [17] H. G. Feichtinger. Banach convolution algebras of Wiener type. In Functions, series, operators, Vol. I, II (Budapest, 1980), pages 509–524. North-Holland, Amsterdam, 1983. [18] H. G. Feichtinger. Generalized amalgams, with applications to Fourier transform. Canad. J. Math., 42(3):395–409, 1990. DECAY AND SMOOTHNESS FOR EIGENFUNCTIONS OF LOCALIZATION OPERATORS 21 [19] H. G. Feichtinger. Modulation spaces on locally compact abelian groups. In Technical report, University of Vienna, 1983, and also in “Wavelets and Their Applications”, pages 99–140. M. Krishna, R. Radha, S. Thangavelu, editors, Allied Publishers, 2003. [20] H. G. Feichtinger and K. Gr¨ochenig. Gabor frames and time-frequency analysis of distributi- ons. J. Funct. Anal., 146(2):464–495, 1997. [21] H. G. Feichtinger and K. Nowak. A first survey of Gabor multipliers. In Advances in Gabor analysis, Appl. Numer. Harmon. Anal., pages 99–128. Birkh¨auser Boston, Boston, MA, 2003. [22] C. Fern´andez and A. Galbis. Compactness of time-frequency localization operators on L2(Rd). J. Funct. Anal., 233(2):335–350, 2006. [23] C. Fern´andez and A. Galbis. Some remarks on compact Weyl operators. Integral Transforms Spec. Funct., 18(7-8):599–607, 2007. [24] C. Fern´andez, A. Galbis. and E. Primo. Compactness of Fourier integral operators on weighted modulation spaces. Trans. Amer. Math. Soc., 372(1):733–753, 2019. [25] J. J. F. Fournier and J. Stewart. Amalgams of Lp and lq. Bull. Amer. Math. Soc. (N.S.), 13(1):1–21, 1985. [26] Y. V. Galperin. Young’s convolution inequalities for weighted mixed (quasi-) norm spaces. J. Inequal. Spec. Funct., 5(1):1–12, 2014. p,q [27] Y. V. Galperin and S. Samarah. Time-frequency analysis on modulation spaces Mm , 0 < p, q ≤ ∞. Appl. Comput. Harmon. Anal., 16(1):1–18, 2004. [28] K. Gr¨ochenig. Foundations of time-frequency analysis. Applied and Numerical Harmonic Anal- ysis. Birkh¨auser Boston, Inc., Boston, MA, 2001. [29] K. Gr¨ochenig. Time-frequency analysis of Sj¨ostrand’s class. Rev. Mat. Iberoam., 22(2):703–724, 2006. [30] K. Gr¨ochenig and Y. Lyubarskii. Gabor (super)frames with Hermite functions. Math. Ann., 345(2):267–286, 2009. [31] K. Gr¨ochenig and J. Toft. Isomorphism properties of Toeplitz operators and pseudo-differential operators between modulation spaces. J. Anal. Math., 114:255–283, 2011. [32] K. Gr¨ochenig and J. Toft. The range of localization operators and lifting theorems for modu- lation and Bargmann-Fock spaces. Trans. Amer. Math. Soc., 365(8):4475–4496, 2013. [33] W. Guo, J. Chen, D. Fan and G. Zhao. Characterizations of Some Properties on Weighted Modulation and Wiener Amalgam Spaces Michigan Math. J., 68:451–482, 2019. [34] C. Heil. An introduction to weighted wiener amalgams. In In: Wavelets and their Applications, M. Krishna, R. Radha and S. Thangavelu, eds.,, pages 183–216. Allied Publishers, New Delhi, 2003. [35] M. Kobayashi. Modulation spaces M p,q for 0 < p,q ≤ ∞. J. Funct. Spaces Appl., 4(3):329–341, 2006. [36] M. Kobayashi. Dual of modulation spaces. J. Funct. Spaces Appl., 5(1):1–8, 2007. [37] F. Luef and E. Skrettingland. Mixed-state localization operators: Cohen’s class and trace class operators. J. Fourier Anal. Appl., 25(4):2064–2108, 2019. [38] F. Luef and E. Skrettingland. On accomulated Cohen’s class distributions and mixed-state localization operators. e-prints, arXiv:1808.06419, 2018. [39] H. Rauhut. Coorbit space theory for quasi-Banach spaces. Studia Math., 180(3):237–253, 2007. [40] H. Rauhut. Wiener amalgam spaces with respect to quasi-Banach spaces. Colloq. Math., 109(2):345–362, 2007. [41] M. A. Shubin. Pseudodifferential operators and spectral theory. Springer-Verlag, Berlin, second edition, 2001. 22 FEDERICO BASTIANONI,ELENA CORDERO, AND FABIO NICOLA [42] S. B. Stechkin. On absolute convergence of orthogonal series. Dokl. Akad. Nauk SSSR, 102:37– 40, 1955. [43] N. Teofanov. Gelfand-Shilov spaces and localization operators. Funct. Anal. Approx. Comput., 7(2):135–158, 2015. [44] N. Teofanov. Continuity and Schatten–von Neumann properties for localization operators on modulation spaces. Mediterr. J. Math., 13(2):745–758, 2016. [45] N. Teofanov. Bilinear localization operators on modulation spaces. J. Funct. Spaces, pages Art. ID 7560870, 10, 2018. [46] J. Toft. Continuity properties for modulation spaces, with applications to pseudo-differential calculus. I. J. Funct. Anal., 207(2):399–429, 2004. [47] J. Toft. Continuity properties for modulation spaces, with applications to pseudo-differential calculus. II. Ann. Global Anal. Geom., 26(1):73–106, 2004. [48] J. Toft. Continuity and compactness for pseudo-differential operators with symbols in quasi- Banach spaces or H¨ormander classes. Anal. Appl. (Singap.), 15(3):353–389, 2017. [49] B. Wang, Z. Huo, C. Hao, and Z. Guo. Harmonic analysis method for nonlinear evolution equations. I. World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2011. [50] M. W. Wong. Wavelet transforms and localization operators, volume 136 of Operator Theory: Advances and Applications. Birkh¨auser Verlag, Basel, 2002. Dipartimento di Matematica, Universita` di Torino, via Carlo Alberto 10, 10123 Torino, Italy E-mail address: [email protected] Dipartimento di Matematica, Universita` di Torino, via Carlo Alberto 10, 10123 Torino, Italy E-mail address: [email protected] Dipartimento di Scienze Matematiche, Politecnico di Torino, corso Duca degli Abruzzi 24, 10129 Torino, Italy E-mail address: [email protected]