Quantum harmonic analysis on lattices and Gabor multipliers

We develop a theory of quantum harmonic analysis on lattices in R. Convolutions of a sequence with an operator and of two operators are defined over a lattice, and using corresponding Fourier transforms of sequences and operators we develop a version of harmonic analysis for these objects. We prove analogues of results from classical harmonic analysis and the quantum harmonic analysis of Werner, including Tauberian theorems and a Wiener division lemma. Gabor multipliers from time-frequency analysis are described as convolutions in this setting. The quantum harmonic analysis is thus a conceptual framework for the study of Gabor multipliers, and several of the results include results on Gabor multipliers as special cases.


Introduction
In time-frequency analysis, one studies a signal ψ ∈ L 2 (R d ) by considering various time-frequency representations of ψ. An important class of time-frequency representations is obtained by fixing ϕ ∈ L 2 (R d ) and considering the short-time Fourier transform V ϕ ψ of ψ with window ϕ, which is the function on the time-frequency plane R 2d given by V ϕ ψ(z) = ψ, π(z)ϕ L 2 for z ∈ R 2d , where π(z) : L 2 (R d ) → L 2 (R d ) is the time-frequency shift given by π(z)ϕ(t) = e 2πiω·t ϕ(t − x) for z = (x, ω). The intuition is that V ϕ ψ(z) carries information about the components of the signal ψ with frequency ω at time x.
Compared to (1) we see that G ϕ,ξ m modifies the time-frequency content of ψ in a simple way, namely by multiplying the samples of its timefrequency representation with a mask m. Gabor multipliers have been studied by [4,8,13,15,19,20,28,32] among others.
Gabor multipliers are the discrete analogues of the much-studied localization operators [1,8,9,31]. In [44] we showed that the quantum harmonic analysis developed by Werner and coauthors [38,54] provides a conceptual framework for localization operators, leading to new results and interesting reinterpretations of older results on localization operators. The goal of this paper is therefore to develop a version of quantum harmonic analysis for lattices to provide a similar conceptual framework for Gabor multipliers.
With this aim we introduce two convolutions of operators and sequences in Section 4. Following [17,39,54] we first define the translation of an operator S on L 2 (R d ) by λ ∈ Λ to be the operator for S ∈ B and trace class T , where · T is the trace class norm. While there are other classes of operators that would ensure that S ⋆ Λ T ∈ ℓ 1 (Λ), see for instance the Schwartz operators [37], B has the advantage of being a Banach space, hence allowing the use of tools such as Banach space adjoints. The space B has previously been studied by [13,16,17] among others.
Equipped with both convolutions and Fourier transforms, we naturally ask whether the Fourier transforms turn convolutions into products. We show for z ∈ R 2d that where Λ • is the adjoint lattice of Λ defined in Section 5, and (3) F W (c ⋆ Λ S)(z) = F Λ σ (c)(z)F W (S)(z). These results include as special cases the so-called fundamental identity of Gabor analysis [18,35,48,52] and results on the spreading function of Gabor multipliers due to [13]. We show versions of (2) and (3) for general classes of operators and sequences, taking care to give a precise interpretation of the objects and equalities in all cases.
A fruitful approach to Gabor multipliers due to Feichtinger [15] is to consider the so-called Kohn-Nirenberg symbol of operators. The Kohn-Nirenberg symbol of an operator S on L 2 (R d ) is a function on R 2d , and Feichtinger used this to reduce questions about Gabor multipliers in the Hilbert Schmidt operators to questions about functions in L 2 (R 2d ). This approach has later been used in other papers on Gabor multipliers [4,13,19]. As Gabor multipliers are examples of convolutions, we show in Section 6 that this approach can be generalized and phrased in terms of our quantum harmonic analysis, and that one of the main results of [15] finds a natural interpretation as a Wiener's lemma in our setting.
In Section 7 we show the extension of some deeper results of harmonic analysis on R d to our setting. We obtain an analogue of Wiener's classical Tauberian theorem, similar to the results of Werner and coauthors [38,54] in the continuous setting. As an example we have the following equivalent statements for S ∈ B : (i) The set of zeros of F Λ σ (S⋆ ΛŠ * ) contains no open subsets in R 2d /Λ • . (ii) If c ⋆ Λ S = 0 for c ∈ ℓ 1 (Λ), then c = 0. (iii) B ′ ⋆ Λ S is weak*-dense in ℓ ∞ (Λ). These results are related to earlier investigations of Gabor multipliers by Feichtinger [15]. In particular, he showed that if S = ξ ⊗ϕ is a rankone operator and F Λ σ (S ⋆ ΛŠ * ) has no zeros, then any m ∈ ℓ ∞ (Λ) can be recovered from the Gabor multiplier G ϕ,ξ m . Since Gabor multipliers are given by convolutions, the equivalence (i) ⇐⇒ (ii) shows that we can recover m ∈ ℓ 1 (Λ) from G ϕ,ξ m under the weaker condition (i) -this holds in particular for finite sequences m.
Finally, we apply our techniques to prove a version of Wiener's division lemma. At the level of Weyl symbols this turns out to reproduce a result by Gröchenig and Pauwels [30], but in our context it has the following interpretation: If F W (S) has compact support for some operator S, and the support is sufficiently small compared to the density of Λ, then there exists a sequence m ∈ ℓ ∞ (Λ) such that S = m ⋆ Λ A for some A ∈ B. If S belongs to the Schatten-p class of compact operators, then m ∈ ℓ p (Λ). The above result fits well into the common intuition that operators S with compactly supported F W (S) (so-called underspread operators) can be approximated by Gabor multipliers [13] -i.e. by operators c⋆ Λ T where T is a rank-one operator. The result shows that if we allow T to be any operator in B, then any underspread operator S is precisely of the form S = c ⋆ Λ T for a sufficiently dense lattice Λ.
We end this introduction by emphasizing the hybrid nature of our setting. In [54], Werner introduced quantum harmonic analysis of functions on R 2d and operators on the Hilbert space L 2 (R d ). We are considering the discrete setting of sequences on a lattice instead of functions on R 2d . If we had modified the Hilbert space L 2 (R d ) accordingly, many of our results would follow by the arguments of [54], as already outlined in [38]. However, we keep the same Hilbert space L 2 (R d ) as in the continuous setting. We are therefore mixing the discrete (lattices) and the continuous (L 2 (R d )), which leads to some extra intricacies.

Conventions
By a lattice Λ we mean a full-rank lattice in R 2d , i.e. Λ = AZ 2d for A ∈ GL(2d, R). The volume of Λ = AZ 2d is |Λ| := det(A). For a lattice Λ, the Haar measure on R 2d /Λ will always be normalized so that R 2d /Λ has total measure 1.
If X is a Banach space and X ′ its dual space, the action of y ∈ X ′ on x ∈ X is denoted by the bracket y, x X ′ ,X , where the bracket is antilinear in the second coordinate to be compatible with the notation for inner products in Hilbert spaces. This means that we are identifying the dual space X ′ with anti linear functionals on X. For two Banach spaces X, Y we use L(X, Y ) to denote the Banach space of continuous linear operators from X to Y , and if X = Y we simply write L(X). The notation P Q means that there is some C > 0 such that P ≤ C · Q.

Spaces of operators and functions
3.1. Time-frequency shifts and the short-time Fourier transform. For z = (x, ω) ∈ R 2d we define the time-frequency shift operator π(z) by (π(z)ψ)(t) = e 2πiω·t ψ(t − x) for ψ ∈ L 2 (R d ). Hence π(z) can be written as the composition M ω T x of a translation operator (T x ψ)(t) = ψ(t − x) and a modulation operator (M ω ψ)(t) = e 2πiω·t ψ(t). The time-frequency shifts π(z) are unitary operators on L 2 (R d ). For ψ, ϕ ∈ L 2 (R d ) we can use the time-frequency shifts to define the short-time Fourier transform V ϕ ψ of ψ with window ϕ by The short-time Fourier transform satisfies an orthogonality condition, sometimes called Moyal's identity [21,26].
, and the relation holds, where the leftmost inner product is in L 2 (R 2d ) and those on the right are in L 2 (R d ).
By replacing the inner product in the definition of V ϕ ψ by a duality bracket, one can define the short-time Fourier transform for other classes of ψ, ϕ. The most general case we need is that of a Schwartz function ϕ ∈ S(R d ) and a tempered distribution ψ ∈ S ′ (R d ); we define

3.2.
Feichtinger's algebra. An appropriate space of functions for our purposes will be Feichtinger's algebra S 0 (R d ), first introduced by Feichtinger in [14].
With the norm above, S 0 (R d ) is a Banach space of continuous functions and an algebra under multiplication and convolution [14]. By [26,Thm. 11.3.6], the dual space We get the following chain of continuous inclusions: . One important reason for using Feichtinger's algebra is that it consists of continuous functions, and that sampling them over a lattice produces a summable sequence [14, Thm. 7C)].
where the implicit constant depends only on the lattice Λ.

3.4.
Banach spaces of operators on L 2 (R d ). The results of this paper concern operators on various function spaces, and we will pick operators from two kinds of spaces: the Schatten-p classes T p for 1 ≤ p ≤ ∞ and a space B of operators defined using the Feichtinger algebra.
3.4.1. The Schatten classes. Starting with the Schatten classes, we recall that any compact operator S on L 2 (R d ) has a singular value decomposition [6, Remark 3.1], i.e. there exist two orthonormal sets {ψ n } n∈N and {φ n } n∈N in L 2 (R d ) and a bounded sequence of positive numbers {s n (S)} n∈N such that S may be expressed as with convergence of the sum in the operator norm. Here ψ ⊗ φ for ψ, φ ∈ L 2 (R d ) denotes the rank-one operator ψ ⊗ φ(ξ) = ξ, φ L 2 ψ.
For 1 ≤ p < ∞ we define the Schatten-p class T p of operators on L 2 (R d ) by T p = {T compact : {s n (T )} n∈N ∈ ℓ p }. To simplify the statement of some results, we also define T ∞ = L(L 2 ) with · T ∞ given by the operator norm. The Schatten-p class T p is a

3.4.2.
The Weyl transform and operators with symbol in S 0 (R 2d ). The other class of operators we will use will be defined in terms of the Weyl transform. We first need the cross-Wigner distribution W (ξ, η) of two functions ξ, η ∈ L 2 (R d ), defined by For f ∈ S ′ 0 (R 2d ), we define the Weyl transform L f of f to be the operator L f : f is called the Weyl symbol of the operator L f . By the kernel theorem for modulation spaces [26,Thm. 14.4.1], the Weyl transform is a bijection from S ′ 0 (R 2d ) to L(S 0 (R d ), S ′ 0 (R d )). Notation. In particular, any S ∈ L(S 0 (R d ), S ′ 0 (R d )) has a Weyl symbol, and we will denote the Weyl symbol of S by a S . By definition, this means that L a S = S.
It is also well-known that the Weyl transform is a unitary mapping from L 2 (R 2d ) to HS [46]. This means in particular that S, T HS = a S , a T L 2 for S, T ∈ HS, which often allows us to reduce statements about Hilbert Schmidt operators to statements about L 2 (R 2d ). We then define B to be the Banach space of continuous operators B consists of trace class operators L 2 (R d ) and we have a norm-continuous inclusion ι : B ֒→ T [24,29].
The dual space B ′ can also be identified with a Banach space of operators. By definition, τ : B → S 0 (R 2d ) given by τ (S) = a S is an isometric isomorphism. Hence the Banach space adjoint τ * : . In this paper we will always consider elements of B ′ as operators Remark. For more results on B and B ′ we refer to [16,17]. In particular we mention that we could have defined B using other pseudodifferential calculi, such as the Kohn Nirenberg calculus, and still get the same space B with an equivalent norm. We would also like to point out that the statements of this section may naturally be rephrased using the notion of Gelfand triples, see [17].
3.5. Translation of operators. The idea of translating an operator S ∈ L(L 2 ) by z ∈ R 2d using conjugation with π(z) has been utilized both in physics [54] and time-frequency analysis [17,39]. More precisely, we define for z ∈ R 2d and S ∈ B ′ the translation of S by z to be the operator α z (S) = π(z)Sπ(z) * . We will also need the operation S →Š = P SP , where P is the parity The main properties of these operations are listed below, note in particular that part (i) supports the intuition that α z is a translation of operators. See Lemmas 3.1 and 3.2 in [44] for the proofs.
By the last part we can unambiguously writeŠ * .

Convolutions of sequences and operators
In [54], the convolution of a function f ∈ L 1 (R 2d ) and an operator S ∈ T was defined by the operator-valued integral and the convolution of two operators S, T ∈ T was defined to be the function These definitions, along with a Fourier transform defined for operators, have been shown to produce a theory of quantum harmonic analysis with non-trivial consequences for topics such as quantum measurement theory [38] and time-frequency analysis [44]. The setting where R 2d is replaced by some lattice Λ ⊂ R 2d is frequently studied in timefrequency analysis, and our goal is therefore to develop a theory of convolutions and Fourier transforms of operators in that setting. For a sequence c ∈ ℓ 1 (Λ) and S ∈ T , we define the operator and for operators S ∈ B and T ∈ T we define the sequence Hence S ⋆ Λ T is the sequence obtained by restricting the function S ⋆ T to Λ.
Remark. We use the same notation ⋆ Λ for the convolution of an operator and a sequence and for the convolution of two operators. The correct interpretation of ⋆ Λ will always be clear from the context.
The fact that S⋆ Λ T is a well-defined and summable sequence on Λ is less straightforward.

4.1.
Gabor multipliers and sampled spectrograms. If we consider rank-one operators, these convolutions reproduce well-known objects from time-frequency analysis. First consider the rank-one operator ξ 1 ⊗ξ 2 for ξ 1 , ξ 2 ∈ L 2 (R d ). The operators c⋆ Λ (ξ 1 ⊗ξ 2 ) are well-known in time-frequency analysis as Gabor multipliers [4,13,15,19]: it is simple to show that which is the definition of the Gabor multiplier G ξ 2 ,ξ 1 c used in timefrequency analysis [19] Remark. In this sense, operators of the form c ⋆ Λ S are a generalization of Gabor multipliers. We mention that this is a different generalization from the multiple Gabor multipliers introduced in [13].
Finally, if S ∈ T is any operator, then one may calculate that often called the lower symbol of S with respect to ϕ 1 , ϕ 2 and Λ [15].
Remark. In particular, Proposition 4.1 does not hold for all S ∈ T . By Remark 4.6 in [4], there exists a function ψ ∈ L 2 (R) such that Since ψ ⊗ ψ,ψ ⊗ψ ∈ T , this shows that the assumption S ∈ B in Proposition 4.1 is necessary.

4.2.
Associativity and commutativity of convolutions. Since the convolution S ⋆ T of two operators S, T ∈ T is commutative in the continuous setting [54,Prop. 3.2], it follows from the definitions that the convolutions (4) and (5) are commutative. It is also a straightforward consequence of the definitions that the convolutions are bilinear.
In the original theory of Werner [54], the associativity of the convolution operations is of fundamental importance. Associativity still holds in some cases when moving from R 2d to Λ, but we will later see in Corollary 7.2.2 that the convolution of three operators over a lattice is not associative in general. In what follows, c * Λ d denotes the usual convolution of sequences Proof. For the proof of (i), we write out the definitions of the convolutions and use the commutativity (4) and (5) = (c ⋆ Λ S) ⋆ Λ T by commutativity.
We have used the easily checked relation α −λ ′ (Š) = P α λ ′ (S)P . For the second part, we find that To pass to the last line we have used the relation , which is easily verified.
Remark. Part (ii) of this result along with the trivial estimate c ⋆ Λ T T ≤ c ℓ 1 T T shows that T is a Banach module (see [23]) over ℓ 1 (Λ) if we define the action of c ∈ ℓ 1 (Λ) on T ∈ T by c ⋆ Λ T . The same proofs also show that this is true when T is replaced by B or any Schatten class T p for 1 ≤ p ≤ ∞.
Example 4.1. Let ϕ, ξ ∈ L 2 (R d ) and c ∈ ℓ 1 (Λ), and define S = ξ ⊗ ξ and T =φ ⊗φ. If we use (8) to simplify S ⋆ Λ T and (9) to simplify (c ⋆ Λ S) ⋆ Λ T , the first part of the result above becomes In words, the convolution of a sequence c with samples of a spectrogram |V ϕ ξ| 2 can be described using the action of a Gabor multiplier c⋆(ξ ⊗ξ).
In applications of convolutional neural networks to audio processing, one often considers the spectrogram of an audio signal as the input to the network. Convolutions of sequences with samples of spectrograms therefore appear naturally in such networks, and the connection (10) has been exploited in this context -see the proof of [12, Thm. 1].

4.3.
Young's inequality. The convolutions in (4) and (5) [54] using the L p -norms of functions and Schatten-p-norms of operators. In the discrete case, it is not always possible to use the Schatten-p-norms, since Proposition 4.1 requires S ∈ B. We will therefore always require that one of the operators belongs to B. A Young's inequality for Schatten classes can then be established by first extending the domains of the convolutions by duality. If S ∈ B and c ∈ ℓ ∞ (Λ), we define c ⋆ Λ S ∈ L(L 2 ) by (11) c ⋆ Λ S, R L(L 2 ),T := c, R ⋆ ΛŠ * ℓ ∞ ,ℓ 1 for any R ∈ T . and if S ∈ B and T ∈ L(L 2 ) = T ∞ we define T ⋆ Λ S ∈ ℓ ∞ (Λ) by (12) T ⋆ Λ S, c ℓ ∞ (Λ),ℓ 1 (Λ) := T, c ⋆ ΛŠ * L(L 2 ),T for any c ∈ ℓ 1 (Λ). It is a simple exercise to show that these definitions define elements of and that they agree with (4) and (5) when c ∈ ℓ 1 (Λ) or T ∈ T . A standard (complex) interpolation argument then gives the following result, since (  ( Remark. If 1 ∈ ℓ ∞ (Λ) is given by 1(λ) = 1 for any λ, then Feichtinger observed in [15,Thm. 5.15] that φ ∈ S 0 (R d ) generates a so-called tight Gabor frame if and only if the Gabor multiplier 1 ⋆ Λ (φ ⊗ φ) is the identity operator I in L(L 2 ). A similar result holds in the more general case: if S ∈ B, then 1 ⋆ Λ S * S = I if and only if S generates a tight Gabor g-frame, recently introduced in [50].
We may also use duality to define the convolution We end this section by showing that the space c 0 (Λ) of sequences vanishing at infinity corresponds to compact operators under convolutions with S ∈ B. The second part of this statement is due to Feichtinger [15,Thm. 5.15] for the special case of Gabor multipliers.
Proof. By [44,Prop. 4.6], the function T ⋆ S belongs to the space C 0 (R 2d ) of continuous functions vanishing at infinity. Since T ⋆ Λ S is simply the restriction of T ⋆ S to Λ, it follows that T ⋆ Λ S ∈ c 0 (Λ). For the second part, let c N be the sequence Then c N ⋆ Λ S = |λ|<N c(λ)α λ (S) is a compact operator for each N ∈ N, and by Proposition 4.3 and the bilinearity of convolutions Hence c ⋆ Λ S is the limit in the operator topology of compact operators, and is therefore itself compact.

Fourier transforms
In [54], Werner observed that if one defines a Fourier transform of an operator S ∈ T to be the function then the formulas hold for f ∈ L 1 (R 2d ) and S, T ∈ T . The transform F W , called the Fourier-Wigner transform (or the Fourier-Weyl transform [54]) is an isomorphism F W : B → S 0 (R 2d ), can be extended to a unitary map F W : HS → L 2 (R 2d ), and to an isomorphism F W : B ′ → S ′ 0 (R 2d ) by defining F W (S) for S ∈ B ′ by duality [17, Cor. 7.6.3]: for any f ∈ S 0 (R 2d ).
Here ρ : S 0 (R 2d ) → B is the inverse of F W . In fact, F W and the Weyl transform are related by a symplectic Fourier transform: for any S ∈ B ′ we have F W (S) = F σ (a S ), where a S is the Weyl symbol of S. As an important special case, the Fourier-Wigner transform of a rank-one operator φ ⊗ ψ is Since we have defined convolutions of operators and sequences, it is natural to ask whether a version of (14) holds in our setting. We start by defining a suitable Fourier transform of sequences.
Symplectic Fourier series. For the purposes of this paper, we identify the dual group R 2d with R 2d by the bijection R 2d ∋ z → χ z ∈ R 2d , where χ z is the symplectic character 1 χ z (z ′ ) = e 2πiσ(z,z ′ ) . Given a lattice Λ ⊂ R 2d , it follows that the dual group of Λ is identified with R 2d /Λ • (see [11,Prop. 3.6.1]), where Λ • is the annihilator group The group Λ • is itself a lattice, namely the so-called adjoint lattice of Λ from [17,48]. Given this identification of the dual group of Λ, the Fourier transform of c ∈ ℓ 1 (Λ) is the symplectic Fourier series Hereż denotes the image of z ∈ R 2d under the natural quotient map is a function on R 2d /Λ • . If we denote by A(R 2d /Λ • ) the Banach space of functions on R 2d /Λ • with symplectic Fourier coefficients in ℓ 1 (Λ), the Feichtinger algebra has the following property [14, Thm. 7 B)].
Lemma 5.1. If Λ is a lattice, the periodization operator P Λ : is continuous and surjective.
1 Phase space, which in this paper is R 2d , is more properly described by (the isomorphic) space R d × R d . The symplectic characters appear because they are the natural way of identifying the group R d × R d with its dual group.
(ii) One may define Feichtinger's algebra S 0 (G) for any locally compact abelian group G [14]. In fact, all our function spaces besides When we identify the dual group of Λ with R 2d /Λ • , the Poisson summation formula for functions in S 0 (R 2d ) takes the following form.
Since F Λ σ is a Fourier transform it extends to a unitary mapping 5.1. The Fourier transform of S ⋆ Λ T . We now consider a version of (14) for sequences. The formula for F Λ σ (S ⋆ Λ T ) is a simple consequence of the Poisson summation formula.
for any z ∈ R 2d .
Proof. From [44, Thm. 8.2], we know that S ⋆ T ∈ S 0 (R 2d ). Hence is an isomorphism. By applying Poisson's summation formula from where we used that F σ is its own inverse to conclude that Since F W (S)F W (T ) ∈ S 0 (R 2d ), Theorem 5.2 says that the equation holds for any z ∈ R 2d .
Remark. Theorem 5.3 has also been proved and used in [43,Cor. A.3] in noncommutative geometry, with stronger assumptions on S, T .
Theorem 5.3 has many interesting special cases. We will frequently refer to the following version, which follows since a short calculation using the definition of the Fourier-Wigner transform shows that (18) F W (Š * )(z) = F W (S)(z).

5.2.
The Fourier transform of c ⋆ Λ S. When c ∈ ℓ 1 (Λ), we obtain the expected formula for F W (c ⋆ Λ S).
Proposition 5.4. If c ∈ ℓ 1 (Λ) and S ∈ T , then Proof. One easily verifies the formula F W (α λ (S))(z) = e 2πiσ(λ,z) F W (S)(z), showing that the Fourier transform of a translation is a modulation. Hence To move F W inside the sum, we use that the sum λ∈Λ c(λ)α λ (S) converges absolutely in T , and F W is continuous from T to L ∞ (R 2d ) by the Riemann-Lebesgue lemma for F W [44, Prop. 6.6].

Technical intermezzo. Let
, where q : R 2d → R 2d /Λ • is the natural quotient map. To extend this map to distributions f ∈ A ′ (R 2d /Λ • ), one can use Weil's formula [25, (6.2.11)] to show that for f ∈ A(R 2d /Λ • ) and g ∈ S 0 (R 2d ) one has . This shows that the map f → f •q agrees with the Banach space adjoint P * Λ • : The natural way to extend f ∈ A ′ (R 2d /Λ • ) is therefore to consider P * Λ • f ∈ S ′ 0 (R 2d ), and by an abuse of notation we will use f to also denote the extension P * Λ • f -by definition this means that when f ∈ . We also remind the reader that for c ∈ ℓ ∞ (Λ) one defines F Λ σ (c) as an element of where (F Λ σ ) −1 (g) are the symplectic Fourier coefficients of g. This is [34, Example 6.8] for the group G = R 2d /Λ • . Finally, recall that we can multiply f ∈ S ′ 0 (R 2d ) with g ∈ S 0 (R 2d ) to obtain an element f g ∈ S ′ 0 (R 2d ) given by (21) f g, h S ′ 0 ,S 0 := f, gh S ′ 0 ,S 0 for h ∈ S 0 (R 2d ).

5.2.2.
The case c ∈ ℓ ∞ (Λ). The technical intermezzo allows us to make sense of the following generalization of Proposition 5.4. Recall in particular that F Λ σ (c) is shorthand for the distribution P * Λ • (F Λ σ (c)) ∈ S ′ 0 (R 2d ). Proposition 5.5. If c ∈ ℓ ∞ (Λ) and S ∈ B, then For h ∈ S 0 (R 2d ), we get from (15), (11) and (20) (in that order) By Theorem 5.3 we find using (18) that where we also used that ρ is the inverse of F W . On the other hand we find using (21) and (19) that Remark. For Gabor multipliers c ⋆ Λ (ψ ⊗ ψ), Propositions 5.4 and 5.5 were proved in [13,Lem. 14], and have been used in the theory of convolutional neural networks [12].

Riesz sequences of translated operators in HS
Two of the useful properties of the Weyl transform f → L f are that it is a unitary transformation from L 2 (R 2d ) to the Hilbert-Schmidt operators HS, and that it respects translations in the sense that As a consequence, statements concerning translates of functions in L 2 (R 2d ) can be lifted to statements about translates of operators and convolutions ⋆ Λ in HS. This approach was first used for Gabor multipliers in [15,19], and has later been explored in other works [4,13] we include these results for completeness, and because the proofs and results find natural formulations and generalizations in the framework of this paper.
For fixed S ∈ HS and lattice Λ, we will be interested in whether {α λ (S)} λ∈Λ is a Riesz sequence in HS, i.e. whether there exist A, B > 0 such that for all finite sequences c ∈ ℓ 2 (Λ) Since the Weyl transform is unitary and preserves translations, if we let a S be the Weyl symbol of S, then (22) is clearly equivalent to the fact that {T λ (a S )} λ∈Λ is a Riesz sequence in L 2 (R 2d ), meaning that for finite c ∈ ℓ 2 (Λ). Following [4,13,15,19] we can use a result from [3] to give a characterization of when (22) holds in terms of an expression familiar from Corollary 5.3.1.
Theorem 6.1. Let Λ be a lattice and S ∈ B. Then the following are equivalent.
Proof. The equality in (i) is Corollary 5.3.1. By the preceding discussion, {α λ (S)} λ∈Λ is a Riesz sequence in HS if and only if {T λ (a S )} λ∈Λ is a Riesz sequence in L 2 (R 2d ). The result from [3] (see [4] for a statement for general lattices and symplectic Fourier transform) says that {T λ (a S )} λ∈Λ is a Riesz sequence if and only if there exist A, B > 0 such that Since the Weyl transform and Fourier-Wigner transform are related by F σ (a S ) = F W (S), we we may restate this condition as Note that the middle term is P Λ • (|F W (S)| 2 )(ż), and since S ∈ B we know that |F W (S)| 2 ∈ S 0 (R 2d ). Therefore P Λ • (|F W (S)| 2 ) ∈ A(R 2d /Λ • ) by Lemma 5.1, which in particular means that P Λ • (|F W (S)| 2 ) is a continuous function on the compact space R 2d /Λ • . For a continuous function on a compact space, condition (23) is equivalent to having no zeros. This completes the proof.
Remark. (i) Since we assume S ∈ B, the first condition above is in fact equivalent to {α λ (S)} λ∈Λ generating a frame sequence in HS, which is a weaker statement than (2) above. The proof of this in [4] for Gabor multipliers works in our more general setting. (ii) As mentioned in the introduction, Feichtinger [15] used the Kohn-Nirenberg symbol rather than the Weyl symbol. This makes no difference for our purposes -we have opted for the Weyl symbol as it is related to F W by a symplectic Fourier transform.
6.1. The biorthogonal system and best approximation. Any Riesz sequence has a so-called biorthogonal sequence and, by the theory of frames of translates [7,Prop. 9.4.2], if the Riesz sequence is of the form {α λ (S)} λ∈Λ for some S ∈ B, then the biorthogonal system has the same form. This means that there exists S ′ ∈ HS such that the biorthogonal system is and biorthogonality means that where δ λ,λ ′ is the Kronecker delta. Now note that for T ∈ HS the definition (5) of T ⋆ Λ S ′ implies that With this observation we can formulate the standard properties of the biorthogonal sequence using convolutions with R.
(iii) For any T ∈ V 2 , Proof. This is simply a restatement of the properties of the biorthogonal sequence of a Riesz sequence using the relation T, α λ (S ′ ) HS = T ⋆ Λ R(λ) -with this observation, parts (i), (ii) and (iii) follow from [7, Thm. 3.6.2].
Remark. (i) If the convolution of three operators were associative, we could find for any T ∈ HS (not just T ∈ V 2 as above) that However, we will soon see that the convolution of three operators is not associative. (ii) For T, R ∈ HS, we have strictly speaking not defined T ⋆ Λ R (since (5) has stronger assumptions than simply HS). However, it is clear by the Cauchy Schwarz inequality for HS that (5) also in this case.
We will now answer two natural questions. First, to what extent does R inherit the nice properties of S -is it true that R ∈ B? Then, how is R related to S? The answer is provided by the following theorem, first proved by Feichtinger [15,Thm. 5.17] for Gabor multipliers, and the proof finds a natural formulation using our tools. Theorem 6.3. Assume that S ∈ B and that {α λ (S)} λ∈Λ is a Riesz sequence in HS. If R is defined as above, then R ∈ B and R = b ⋆ ΛŠ * where b ∈ ℓ 1 (Λ) are the symplectic Fourier coefficients of .
Proof. By [7, Thm. 3.6.2], the generator S ′ of the biorthogonal system belongs to V 2 , hence there exists some b ′ ∈ ℓ 2 (Λ) such that S ′ = b ′ ⋆ Λ S. Since R =Š ′ * , one easily checks by the definitions ofˇand * that . By part (i) of Lemma 6.2 and the associativity of convolutions, we have Taking the symplectic Fourier series of this equation using (17) and Corollary 5.3.1, we find for a.e.ż ∈ R 2d /Λ • , and by assumption on S (see Theorem 6.1 and its proof) the denominator is bounded from below by a positive constant. Since S ∈ B, we know that |F W (S)| 2 ∈ S 0 (R 2d ), and therefore Lemma 5.1 implies that To prepare for the next result, fix S ∈ B and let hence V ∞ is the set of operators given as a convolution c ⋆ Λ S for c ∈ ℓ ∞ (Λ). The first part of the next result says that when {α λ (S)} λ∈Λ is a Riesz sequence, then the Schatten-p class properties of c ⋆ Λ S are precisely captured by the ℓ p properties of c. This result appears to be a new result even for Gabor multipliers. We also determine for any T ∈ HS the best approximation (in the norm · HS ) of T by an operator of the form c ⋆ Λ S. See [15,Thm. 5.17] and [13,Thm. 19] for the statement for Gabor multipliers. Corollary 6.3.1. Assume that S ∈ B and that {α λ (S)} λ∈Λ is a Riesz sequence in HS, and let R be as above.
Hence V ∞ ∩ T p = ℓ p (Λ) ⋆ Λ S and c ℓ p c ⋆ Λ S T p c ℓ p . (ii) For any T ∈ HS, the best approximation in · HS of T by an operator c ⋆ Λ S with c ∈ ℓ 2 (Λ) is given by Equivalently, the symplectic Fourier series of c is given by Proof. (i) By Proposition 4.3 part (i) we get C R (T ) ℓ p ≤ T T p R B , and by part (ii) of the same proposition we get D S (c) T p c ℓ p S B . Hence both maps in the statement are continuous. It remains to show that the two maps are inverses of each other, which will follow from the associativity of convolutions. First assume that c ∈ ℓ p (Λ). Then where we have used associativity and part (i) of Lemma 6.2. Then assume T ∈ V ∞ ∩ T p , so that T = c ⋆ Λ S for c ∈ ℓ ∞ (Λ). We find Hence D S and C R are inverses. In particular V ∞ ∩T p = ℓ p (Λ)⋆ Λ S as D S is onto V ∞ ∩ T p , and V ∞ ∩ T p is closed in T p (hence a Banach space) since D S : ℓ p (Λ) → T p has a left inverse C R and therefore has a closed range in T p . (ii) We claim that the map T → (T ⋆ Λ R) ⋆ Λ S is the orthogonal projection from HS onto ℓ 2 (Λ) ⋆ Λ S, which is a closed subset of HS = T 2 by part (i) (or (24)). If T = c ⋆ Λ S for some c ∈ ℓ 2 (Λ), then c = T ⋆ Λ R by part (i) -therefore T = (T ⋆ Λ R) ⋆ Λ S. Then assume that T ∈ (ℓ 2 (Λ) ⋆ Λ S) ⊥ . As we saw in (25), we can write From the proof of Theorem 6.3, One easily checks that (26) shows that (T ⋆ Λ R) ⋆ Λ S = 0. Finally, to obtain the equivalent expression recall from Theorem 6.3 that R = b ⋆ ΛŠ * for b ∈ ℓ 1 (Λ). Hence by associativity and commutativity of convolutions, It follows from (17) that we get We have a known expression for F Λ σ (b) from Theorem 6.3, and a known expression for F Λ σ (T ⋆ ΛŠ * ) from Theorem 5.3 -inserting these expressions into the equation above yields the desired result.
The key to the results of this section is Wiener's lemma, used in the proof of Theorem 6.3. In fact, we may interpret these results as a variation of Wiener's lemma. To see this, recall that V 2 = span{α λ (S) : λ ∈ Λ} = ℓ 2 (Λ) ⋆ Λ S ⊂ HS. Then {α λ (S)} λ∈Λ is a Riesz sequence if and only if the convolution map D S : ℓ 2 (Λ) → V 2 given by has a bounded inverse [7, Thm. 3.6.6]. Corollary 6.3.1 therefore says the following: if S ∈ B and the convolution map D S : ℓ 2 (Λ) → V 2 has a bounded inverse, then the inverse is given by the convolution The similarities with Wiener's lemma are evident when we compare this to the following formulation of Wiener's lemma [27,Thm. 5.18]: If b ∈ ℓ 1 (Z) and the convolution map ℓ 2 (Z) → ℓ 2 (Z) defined by c → c * Z b has a bounded inverse on ℓ 2 (Z), then the inverse is given by the convolution map

Tauberian theorems
In the continuous setting, where one considers functions on R 2d and the convolutions briefly introduced at the beginning of Section 4, a version of Wiener's Tauberian theorem for operators was obtained by Kiukas et al. [38], building on earlier work by Werner [54]. This theorem consists of a long list of equivalent statements for T p and L p (R 2d ) for p = 1, 2, ∞, and as a starting point for our discussion we state a shortened version for p = 2 below. We wish to obtain versions of this theorem when R 2d is replaced by a lattice Λ, functions on R 2d are replaced by sequences on Λ and we still consider operators on L 2 (R d ). In this discrete setting, statements (3) and (4) in Theorem 7.1 are still equivalent, mutatis mutandis, while the analogues of (1) and (5) can never be true. First we show that the discrete version of statement (1) can never hold.
Proposition 7.2. Let Λ be any lattice in R 2d and let S ∈ HS. Then the linear span of {α λ (S)} λ∈Λ is not dense in HS.
(1) The following are equivalent: If any of the statements above holds, c ∈ ℓ ∞ (Λ) is recovered from c ⋆ Λ S by c = (c ⋆ Λ S) ⋆ Λ R for some R ∈ B. In particular, the map c ⋆ Λ S → c is continuous L(L 2 ) → ℓ ∞ (Λ).
Let us rewrite the statements of the theorem in the case that S is a rank-one operator S = ϕ ⊗ ϕ for ϕ ∈ S 0 (R d ). By (8) we find that and by (6) c ⋆ Λ S is the Gabor multiplier Hence the equivalences (i) ⇐⇒ (ii) provides a characterization using the symplectic Fourier series of V ϕ ϕ| Λ of when the symbol c of a Gabor multiplier is uniquely determined.
for z ∈ supp(f ). The next result is a version of this statement for the convolutions and Fourier transforms of operators and sequences. At the level of Weyl symbols, this result is due to Gröchenig and Pauwels [30] (see also the thesis of Pauwels [45]) using different techniques. We choose to include a proof using the techniques of this paper to show how the the statement fits our formalism. Note that apart from the function g -introduced to ensure A ∈ B -Theorem 7.4 is obtained by replacing the convolutions and Fourier transforms in Wiener's division lemma by the convolutions and Fourier transforms of sequences and operators.
To show that S = (S ⋆ Λ T ) ⋆ Λ A, we will show that their Fourier-Wigner transforms are equal. Using Proposition 5.4 and Theorem 5.3 we find that To show that this equals F W (S), we consider three cases.
A similar argument using duality brackets shows that essentially the same result even holds for S ∈ B ′ . Theorem 7.5. Assume that S ∈ B ′ satisfies supp(F W (S)) ⊂ (1 − 2ǫ) Λ • for some 0 < ǫ < 1/2. Pick g ∈ C ∞ c (R 2d ) such that g| (1−ǫ) Λ • ≡ 1 and supp(g) ⊂ Λ • . If T ∈ B satisfies F W (T )(z) = 0 for z ∈ supp(g), then S = (S ⋆ Λ T ) ⋆ Λ A, where A ∈ B is given by F W (A) = g F W (T ) . Proof. We have already seen that A ∈ B. Let f ∈ S 0 (R 2d ). Then In the last line we multiplied the right hand side by a bump function b ∈ C ∞ c (R 2d ) such that b| (1−2ǫ) Λ • ≡ 1 and supp(b) ⊂ (1 − ǫ) Λ •this does not change anything by the assumptions on the supports of F W (S) and b. We find using Theorem 5.3 and Proposition 5.4 that (18).
Operators S such that supp(F W (S)) ⊂ [− a 2 , a 2 ] d × [− b 2 , b 2 ] d where ab ≤ 1 are called underspread, and provide realistic models of communication channels [13,30,40,41,51]. We immediately obtain the following consequence. Corollary 7.5.1. Any underspread operator S ∈ B ′ can be expressed as a convolution T = c⋆ Λ A with c ∈ ℓ ∞ (Λ) and A ∈ B for a sufficiently dense lattice Λ. In particular, S is bounded on L 2 (R d ).
It is known (see [13]) that for an operator S to be well-approximated by Gabor multipliers -i.e. operators c ⋆ Λ (ψ ⊗ ψ) for ψ ∈ L 2 (R d ) -S should be underspread. The result above shows that any underspread operator S is given precisely by a convolution S = c ⋆ Λ A if we allow A to be any operator in B, not just a rank-one operator. In fact, A as constructed in the theorem will never be a rank-one operator, since F W (A) has compact support -this is not possible for rank-one operators [36]. If S satisfies S ∈ T p in addition to the assumptions of Theorem 7.5, then c = S ⋆ T ∈ ℓ p (Λ) by Proposition 4.3. Hence the p-summability of c in S = c ⋆ Λ A reflects the fact that S ∈ T p . Theorem 7.5 also implies that underspread operators S are determined by the sequence S ⋆ Λ T when T ∈ B is chosen appropriately. This was a major motivation for [30], since when T is a rank-one operator T = ϕ ⊗ ϕ, the sequence S ⋆ ΛŤ is the diagonal of the so-called channel matrix of S with respect to ϕ -see [30,45] for a thorough discussion and motivation of these concepts. Finally, note that Theorem 7.5 gives a (partial) discrete analogue of part (5) of Theorem 7.1.