Mean Field Derivation of DNLS from the Bose–Hubbard Model

We prove that the flow of the discrete nonlinear Schrödinger equation (DNLS) is the mean field limit of the quantum dynamics of the Bose–Hubbard model for N interacting particles. In particular, we show that the Wick symbol of the annihilation operators evolved in the Heisenberg picture converges, as N becomes large, to the solution of the DNLS. A quantitative Lp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^p$$\end{document}-estimate, for any p≥1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p \ge 1$$\end{document}, is obtained with a linear dependence on time due to a Gaussian measure on initial data coherent states.

The aim of the present paper is to show that the quantum expectation ofb k / √ N is close to u k (t, ω) in a suitable L p -measure sense when N is large, which is precisely what we mean by the mean field limit of (2). This is achieved thanks to an explicit estimate in terms of the parameters of the model, global in time.
The literature on the mean field derivation of the nonlinear Schrödinger equation (NLS), the Hartree equation, and more in general about the study of many-body quantum-mechanical systems, is quite rich (see Sect. 2). However, it seems that a direct mean field derivation of the DNLS (1) together with quantitative, explicit estimates is missing. In some works (see for example [1,33,37] and references therein) the DNLS is obtained directly from the NLS equation, in the framework of the tight-binding approximation. However, combining these two kinds of results, a growth exponential in time of the mean field estimate for DNLS follows (essentially due to the Grönwall lemma). In the present paper, a growth linear in time of the mean field estimate is provided. Our first step to deal with the mean field asymptotics is to consider the rescaled operatorsâ and the Heisenberg equation whereâ k (0) :=â k and 1 ≤ k ≤ L. Notice that Eq. (4) is clearly the operator counterpart of (1). We now rewrite Eq. (4) in terms of the Wick symbols (see Sect. 3.1) of the operatorsâ k (t) and H. Let φ ω (z) := e ω·z− 1 2 |ω| 2 be the normalized coherent states in F B (C L ) and recall thatâ k φ √ Nω = ω k φ √ Nω . Define the symbols ρ k (t,ω, ω) := φ √ Nω ,â k (t)φ √ Nω ; (5) H N (ω, ω) := φ √ Nω , Hφ √ Nω = N H(ω, ω) (6) = N 1≤j≤L E j |ω j | 2 + J(ω j+1 ω j +ω j , ω j+1 ) + U 2 |ω j | 4 . (7) Then, by the Wick bracket (see [9,15]) we get the equation with initial data ρ k (0, ω,ω) = ω k . We recall that, as an asymptotic series, that here the Wick bracket is a finite sum { ·, H } Wick = L 1 + L 2 where (10) Notice that L 1 is the Poisson bracket in the variables (ω, ω). Thus, it is easily seen that the DNLS (1) exactly reads By denoting Δ := {(ω, ω) | ω ∈ C L } ⊂ C 2L , and (Φ t , Φ t ) : Δ ⊂ C 2L → C 2L the flow ofγ = i(∂ ω H(γ), −∂ωH(γ)), it follows that The equality (12), together with Eq. (8), tells us that ρ k − u k is a kind of semiclassical perturbation term, and thus we expect ρ k − u k → 0 as N → +∞. Indeed, we will prove such a result with respect to an L p (μ N )-norm, where p ≥ 1 and μ N is a suitable Gaussian measure, invariant under the DNLS flow. With respect to this target, recall that the total number operator defined as is conserved by the quantum flow, i.e., Moreover, the well-known 2 -conservation law for the DNLS can be rewritten as Both these two important properties will be used in the proof of Theorem 1, and for this reason we define the invariant Gaussian probability measure where ω = x + iy, dω ∧ dω := π −L dxdy and c N,L := N L is the normalization constant. This measure is linked (see Proposition 1) to a weighted trace formula involving Wick operators that will be an important tool to our approach.
We are now ready to state the main result of the paper.

Theorem 1.
Let u(t, ω) be the flow of the DNLS Eq. (1) and let ρ k (t, ω) be the solution of (8) Here, Γ and S(α, β) denote the Gamma function and the Stirling numbers of second kind, respectively.
We notice that (17) can be written with the condition L/N → 0 as N → +∞, this means that the number of particles can be supposed to be large with respect to the number of L "sites" of the Bose-Hubbard model, which is the regime considered in some experiments, see for example [40]. We also stress that the L p -norm used above allows to discuss, in the measure sense, the pointwise estimate for |ρ k −u k |(t,ω, ω). Indeed, we have the following Corollary 1. Fix a parameter 0 < < 1 2 and define the set Then, for any Notice that p → A p is an increasing function, whence inequality (21) provides, when N is fixed and p is large, a measure of the region where |ρ k −u k | is large.
On the other hand, in the case of a fixed p ≥ 1 and large values of N we have a vanishing measure of the region where |ρ k − u k | is super-linear in time. We underline an important consequence of this observation. Indeed this means that, from the viewpoint of the Gaussian measure, if this superlinear (in time) mean field estimate is sharp then it is associated with a set of coherent states which is negligible as N → +∞ with the rate shown in (21). Of course, any exponential (in time) upper bound gives rise to the same conclusion. An interesting open problem is to show that the same feature holds for more general quantum dynamics than the one associated with our manybody operator (2). Our paper deals with the Bose-Hubbard model, a simpler setting with respect to that of quantum field theory. However, the explicit estimate in terms of the parameters of the model and its linear dependence on time in (17) seem to be a novel and promising result with respect to other kind of mean field estimates on the NLS equation.
Furthermore, we stress that Theorem 1 can be seen as an Egorov-type result, written to the first order and with respect to the L p -norm, for Wick symbols. With respect to this observation, we recall Proposition 5.1 in [16] where it is proved the convergence, as → 0, of the Wick symbol of an evolved quantum observable toward the Weyl symbol composed with the Hamiltonian flow. In this result, the well-known bound of the Ehrenfest time |t| < T is shown. We also recall Proposition 5.10 and Theorem 5.6 in [5] where, in the framework of evolved Wick operators on the Fock space and with a quantum dynamics much more general than our, it is proved the convergence toward the solution of the Hartree equation as 1 N → 0, but the estimate on the remainder in Theorem 5.6 is again local in time. In our main result, we avoid locality in time by making use of the L p (μ N )-norm (the meaning and the properties of the measure μ N are clarified in Proposition 1 and 2).
In Sect. 7.2 of [26], the authors discuss, thanks to the Wick quantization for a class of symbols, how the many-body quantum mechanics of bosons can be viewed as a deformation quantization of the Hartree theory. We stress that our paper also makes use of Wick operators, but deals with a different class of quantum dynamics and another way to get the derivation of the mean field dynamics, which is here a discrete NLS.
To conclude, we stress the absence in (17) of the parameters E and J involved in the quadratic part of the operator H in (2), in agreement with a well-known elementary result: any quantum expectation of the Heisenberg equations of a linear system (quadratic Hamiltonian) yields the classical equations of motion. Thus, the distance between the Wick symbol ρ k solving Eq. (8) and the k-th component u k of the flow for Eq. (11) is ruled only by nonlinearity, namely by the parameter U . As a consequence, Theorem 1 holds for Hamilton operators H with a completely general quadratic part. This is not the primary target of the present work, but we observe here that a more general setting of H ensures a larger set of invariant measures for the DNLS flow and whence an interesting open problem is to study the link between this kind of mean field estimates and the possible various invariant measures.
We also remark that the equation in (1) with general quadratic part, that usually describes particles in one-dimensional periodic lattice, can be used also to modelize two-and three-dimensional lattices with different topologies (see for example [22] and references therein).
The paper is organized as follows. In Sect. 2, we shortly comment on the most influential results in the literature on the subject, to finally stress the main innovations characterizing our work. Sect. 3 is devoted to the proofs of Theorem 1 and Corollary 1 stated in the Introduction; the proof is divided into various technical steps. Sect. 3.1 is "Appendix" consisting of three subsections.

Synopsis of the Literature and Motivations of the Work
We here provide a commented list of some papers involving NLS, Hartree equation, and more in general the study of many-body quantum-mechanical systems in mean field and semiclassical limit, relevant and somehow connected to our work.
A first reference work in the field, is the review [39], where the author discusses a variety of classical as well as quantum models for which kinetic equations can be derived rigorously, and where the probabilistic nature of the problem is emphasized.
A second reference paper, of particular interest to our work for its use of coherent states, is [29], where Hepp shows that in the many-body framework the classical limit of the expectation values of products of Weyl operators, translated in time by the quantum dynamics and taken on coherent states centered in x-space and p-space are shown to become the exponentials of coordinate functions of the classical orbit in phase space. The Hepp results have been extended in [27]. For a recent review and discussion of Hepp's method, we also address the reader to Sect. 10.4.2 of [21].
The convergence of the N -particle Schrödinger dynamics of bosons toward the Hartree dynamics is proved in the mean field limit in [24]. The authors work in the Heisenberg picture (as in our paper) with a class of bounded operators, whereas we consider annihilation operators, that are unbounded, and another type of convergence.
A rigorous derivation of the cubic NLS in dimension one is shown in [4]. An approach for deriving higher-order corrections to the mean field limit for the quantum systems is provided in [11], a simple and effective method is given in [35].
A reference role in the literature is played by those recent works dealing with the rigorous version of the Bogoliubov theory of superfluids; see, e.g., the review [38]. The Gross-Pitaevskii equation is rigorously deduced, for example, in [14], whereas the fluctuations around it are studied in [13,17]. The convergence to a limiting Hartree dynamics is instead studied in [6,36], whereas in [7,10] and the Hartree-Fock-Bogoliubov and the Bogoliubov-de Gennes equations are derived by the method of the quasi-free reduction.
Further results are, for example, a derivation of the 1D focusing cubic NLS obtained in [20], where the difficulties due to the attractive interaction are discussed and new energy estimates are shown, and a mean field derivation of the defocusing 2D cubic NLS is provided in [28]; the mean field dynamics of a mixture of bosons is treated in [32].
The literature on the subject is actually huge, and the above description represents just a short summary of it. However, with the respect to the existing framework of methods and results, our contribution here is characterized by a certain number of aspects deserving a short discussion.
1. We consider the Bose-Hubbard lattice model (2). This is certainly a context much simpler than the quantum field theory of a boson gas generally considered in the literature quoted above. However, its interest is motivated by the modern experiments on many-body effects in optical lattices, where the lattice models are the basic tool for the theoretical interpretation of the results [12,19]. 2. The privileged physical quantity considered here is the coherent expectation of the local annihilation operator, which is shown to satisfy the DNLS equation within a certain approximation limit. As specified below, we would be able to do the same with any observable, say any polynomial of the Dirac operators. The choice of the annihilation operator, besides its simplicity, is quite natural if, along an Ehrenfest-like line of thought, one wants to compare the quantum expectation of the Heisenberg equation of a certain operator with the "classical" scalar equation obtained by replacing the operator with its quantum expectation in the Heisenberg equation itself, as heuristically done in physics. 3. A clear innovation of our approach consists in distributing the initial data (i.e., coherent states) according to an invariant probability measure, which then calls naturally for the use of the L 2 norm (which is then generalized to L p ), as typical and meaningful in statistical mechanics, the right framework for this kind of problems; for the experimental relevance of coherent states see [12,19]. However, in certain applications, e.g., quantum computing, special initial conditions, and thus point-wise estimates, may play the relevant role. Indeed, we refer for example to the quantum walks in the Bose-Hubbard model, that are unitary processes describing the evolution of initially localized wave functions on a lattice potential, see [30].
Concerning the choice of the measure, we take the Gaussian one, inherited by the quantum trace measure with density e −λ N , λ being a suitable parameter. Of course, in a statistical mechanical framework, one would like to work with the Gibbs-Von Neumann density, namely e −β H , β being the inverse temperature. The latter point makes part of a work in progress. Here, we only stress that the Hamiltonian (2) reads H = E N + · · · , the dots denoting the other two terms, N commuting with both of them. As will be further discussed below, and expected from the tight-binding assumptions made to deduce the Bose-Hubbard model, the first term E N is the leading one with respect to the other two. In a sense, we are thus considering an invariant measure that is approximately connected to the Gibbs one. 4. The joint use of an invariant measure on the initial coherent states and of the Wick formalism allows us to bound distance between the symbol of an observable at time t and the classical evolution of its initial symbol by a constant growing linearly with t. The linear dependence on time is not a particular feature of the Gaussian measure, the latter being instead quite convenient in order to get an explicit estimate of the overall constant multiplying time. Within the framework of results in measure on coherent states, the linear growth in time of our bound represents an interesting news, since most of time dependencies obtained in the literature up to now are typically exponential, which is an unavoidable consequence of the Grönwall lemma.

Results
In this section, we provide the proof of the main theorem we have stated in the Introduction. To such a purpose, we will need some preliminary lemmas and propositions. Among them, Lemmas 1 and 2 are just quoted and used, their statements and proofs being reported at the end of the section.
The following result provides a weighted trace formula for Wick operators, involving the positive definite operator e −λ N with λ > 0. In order to make a link with the Gaussian measure μ N given in (16), we have to write a bijective relation between the parameters λ and N . This result will be useful for the subsequent result on expectation values of Wick operators under quantum dynamics. Proposition 1. Let μ N be as in (16). Let Op W (g) be a Wick operator on where γ λ := Tr(e −λ N ) and e λ = N + 1.
Proof. We begin by the equality and notice that the Wick symbol of e −λ N reads Equality (25) allows to write the constant Thanks to the Wick-product, (24) can be rewritten as We also remind formula (2.38) in [15] that provides a link between Wick and anti-Wick symbols where ∂ω k ∂ω k , and recall that Wick and anti-Wick symbols of e −λ N are unique. Now write explicitly the Wick-product and integrate by parts the integral in (27). This gives (formally) Recall that where c N,L := N L . Our target is thus to prove the well-defined equation namely which is solved by μ = N/(N + 1), and since μ = 1 − e −λ we recover Remark 1. We now recall thatb † kb μ = Op W (g) when g =ω k ω μ , (see Sect. 3.1). For these Wick operators, Proposition 1 reads Since {z α √ α! , α ∈ Z n + } is an orthonormal set in the Fock-Bargmann space (see [15]), an easy computation shows that Thus, equality (37) can be considered as the version, in the Fock-Bargmann space, of the Quantum Wick Theorem showed in [25] that works in the Fock space and with the related bosonic creation and annihilation operators of quantum field theory.
In the next, we provide a kind of quantum mean value formula for the time evolved G(s) := U † (s) GU (s) where U (s) = e −i Hs with H as in (2) and G = Op W (g) are Wick operators (see Sect. 3.1). This result will be applied within the proof of Theorem 1 for operators of type G = (â † kâ k + 1 N ) p with creation and annihilation operators as in (3). This tool allows to avoid, in our setting and for our estimates, the well-known problem of Ehrenfest time, as well as to avoid the application of Grönwall Lemma (and thus exponential in time upper bounds) used in many papers on mean field estimates for NLS equations.
Proof. We apply Proposition 1 and recall that the trace is invariant by unitary conjugations of operators, so that Now we recall that [ N, H] = 0 and whence [ N, U (s)] = 0, which gives and applying again Proposition 1 for G(s) := U † (s)Op W (g)U (s) and, in view of Remark 2, we conclude This observation will be useful in the application of this equality with a = √ N .
We provide two technical lemma used in the next.

Lemma 1.
Let P(ω, ω) be as in (68), then for any 1 ≤ p < ∞ there exists a positive constant C 1,p such that Moreover, for a positive constant C 2,p .
Proof. . We first notice that for any fixed ω ∈ C L , P(ω, ω) is a sum of real nonnegative numbers , so by using Hölder inequality we get Since for any integrating with respect to Gaussian measure and performing the change of variables ω j = √ Nω j we have For each j = 1, . . . , L, we factorize the integrals not containing ω j , so introducing the variable v ∈ C and its corresponding measure dv ∧ dv we have where the Euler Gamma function has been introduced. Notice that in the first line we exploited the definition of c N,L = N L , while in the second line we used the fact that we have exactly L equal integrals. Taking the fourth of root in the last expression, we get inequality (44) with To get (45), we need to compute the mean value ϕ √ Nω ,n α k ϕ √ Nω =: n α k for any positive integer α, wheren k =â † k (0)â k (0). We find that from the definition of Wick- * product there exists a recurrence relation between these quantities where 1 is the constant function 1(ω, ω) ≡ 1, so that in general where S(α, β) is the Stirling number of the second kind with integer parameters α and β (see computations below). Sincen k and N −1 commute as operators we can expand (n k + N −1 ) 2p using the binomial theorem Taking again the fourth root, we get (45) with constant We now complete the proof of this lemma, showing that the coefficients of the polynomial in (52) are the Stirling number of the second kind (see [3], Par. 24.1.4 for their definition and properties). By the recurrence relation (51), we have where we used the fact that S(α, α) = S(α, 1) = 1, as is easy verified using (51). Comparing the last expression with the general expansion of n α+1 k as in (52) with exponent α + 1, we see that which is precisely the recurrence relation defining Stirling numbers.
Proof. We begin by In particular, the second term can be rewritten and this last form equals φ √ Nv , Op W (g) In what follows, we get an estimate for |ρ k (t,ω, ω) − u k (t, ω)| for any fixed ω ∈ C L . This will be used, in the proof of Theorem 1, to have the L p (μ N ) estimate.

Proposition 3. Let
with H as in (6). Let u(t, ω) := (u 1 , . . . u L )(t, ω) be the solution of (1), and Then, ds. (69) Proof. The semigroup identity applied to our case gives where the operator L 2 reads and thus We now recall the definition where φ √ Nv denotes the complex conjugated coherent state. Thanks to Lemma 2, Applying twice this formula, we get Applying the same computations for the derivatives onv j , we get The sum in (73) can now be rewritten as which simplifies to The sum exhibits the following upper bound namely We need to get an estimate for â † and thus We now look at By using againâ and hence Inserting (90)-(93) into (88), we get Thus, Ann. Henri Poincaré We observe that and As a consequence, and recall equalities (71)-(73) which imply the statement (69).
In view of previous propositions, we can now provide the proof of the main result of the paper.
Proof of Theorem 1. Recalling (69), we define the positive function and for the sake of simplicity we avoid to write the dependence on (ω, ω).
More in details, The Hölder inequality fg L 1 ≤ f L p g L q with 1/q + 1/p = 1, allows and hence t 0 ψ(s) ds This gives We now focus our attention to The invariance of μ N under the flow Φ t−s implies where we have just defined the positive definite operator Now assume that p = 2 m with m ∈ N so that We get, thanks to the normalization of μ N , and recalling (109), The Cauchy-Schwarz inequality gives Sinceâ † k (s)â k (s) + 1 N is positive definite, we have the upper bound Observe that, since U (s)U (s) = Id, Now apply Proposition 2 and Remark 3 in order to rewrite (114) as Integrating these terms, see Lemma 1, we have two constants C 1,p and C 2,p > 0 such that and Thus, We are now in the position to conclude so that by defining we have, in the case p = 2 m with m ∈ N, Now observe that, thanks to normalization of μ N and a simple application of Hölder inequality, we have ρ for any α ≥ p. Thus, fix α := 2 p so that ensures now for all p ≥ 1 the inequality where the last inequality is guaranteed since μ N is a Gaussian type measure and |ω k | p is a polynomial term. Inequality (125) gives ρ k (t) − u k (t) L p (μN ) < +∞ and thus ρ k L p (μN ) < +∞.
An immediate consequence of Theorem 1 is the next corollary.
Proof of Corollary 1. Let 0 < < 1 2 and define the set Then, Recalling inequality (125), we get hence under for example the assumption that |U rs kj | ≤Ū . Indeed, in this more general case the dependence from the parameter L into the L p -estimate of Theorem 1 will be L 4 /N in place of L/N , namely where the constant A p will be changed by a more general formula.
In the current paper, we have considered the most simple setting involving the quantum dynamics driven by H as in (2) in order to make more clear and direct the exposition of our approach of mean field estimates for Wick operators on Bargmann space. Moreover, concerning our focussing on the dynamics of a k as the main physical quantity, we stress that along the same lines, we can treat the dynamics of polynomials in a k and a † k . For example, the quantity could be considered in the same way, and the asymptotics as N → +∞ recovers, in L p (μ N )-norm, the quantity related to the components j and k of the DNLS flow at time t.
Finally, for what concerns the problem of reduced density matrices (see for example [35,36]), we recall the setting of the operator Γ Ψ : L 2 → L 2 (see Sects. 10 where Ψ = Ψ (N ) is fixed in the N -particle sector F (N ) L 2 s (R N ; C) of the Fock space, and where b * (u j ) and b(u k ) are the bosonic creation and annihilation operators on the Fock space, and associated with a fixed orthonormal set {u j } j∈N of L 2 (R; C). A preliminary treatment of the link between Γ Ψ with the contents of our paper is reported in the PhD Thesis of one of the authors [34].

Fock-Bargmann Space and Wick Quantization
In this subsection, we provide an overview of Fock-Bargmann space and the Wick quantization, recalling the basic notions we need in the framework of this paper. Here, we mainly follow the notations of Sect. 5.2 in [15], but we also address the reader to Sects. 1.6-2.7 of [23] and Sect. 2.6 in [21]. LetĀ(C n ) be the set of the anti-analytic functions f : C n → C. The Fock-Bargmann space is defined as with a scalar product here z := x + iy and dz ∧ dz := π −n dxdy, thus the integral can be written over R 2n . In this paper, we consider n = L where L is a parameter of the Bose-Hubbard model. The creation and annihilation operators are defined as: Notice thatb † k ,b k are well defined on F B (C n ) and [b k ,b † μ ] = δ kμ Id. Coherent states are represented, with its normalization e − 1 2 |ω| 2 , as For a given operator A : whereas outside the diagonal (ω, ω) the Wick symbol reads The Wick quantization of an entire function σ : C n × C n → C is given as In view of these settings, we have A = Op W (σ( A)), and we call this a Wick operator.
To be more precise about the set of these operators, suppose that A (possibly unbounded) is defined on F B (C n ) together with its adjoint A † , and assume that for all ω ∈ C L , φ ω belongs to the domains of A and A † . Then, ω → σ( A)(ω, ω) is a smooth function on C n and moreover σ( A)(ω, ω) is the restriction on the diagonal of σ( A)(z, ω) as in (137), which is furthermore an entire function (see [23]-pg. 139). As shown in Proposition 1.69 of [23], any entire function K(z, ω) is uniquely determined by its restriction to {z =ω}.
Thanks to these observations, A = Op W (σ) is uniquely related to the symbol on the diagonal, and for this reason frequently in the literature one refers to A as the Wick quantization of (136).
We also stress that if the starting point is the definition (138) one must prove, for a given entire function σ, that Op W (σ) is well defined on F B (C n ) is the sense we have just described.
A simple computation shows that These equalities directly allow to write the Wick symbol of the Bose-Hubbard operator H in (2) φ ω , and the rescaled Wick symbol We also recall that any bounded operator on F B (C n ) is a well-defined Wick operator. Furthermore, a large class of Weyl operators (see Sect. 2.1 in [23]) can be rewritten as a Wick operator by the following link of symbols σ Wick = e Δ/2 σ Weyl , namely (for = 1) see also Proposition 2.97 in [23] for the link with standard quantization, and more detailed setting for the allowed symbols. The set of Wick operators is closed under composition, and the Wickproduct is defined as the symbol of the composition of two operators, It can be shown (see [15]) the following asymptotics (in multi-index notation) where ∂ω r := ∂ω i1 ∂ω i2 . . . ∂ω ir . About the convergence of the right-hand side, we address the reader to [9]. We stress that in the asymptotics (144) the semiclassical parameter is = 1 and the absence of the factor 2 r used in [9] is a consequence of the setting of the scalar product in (133).
The Wick bracket is defined as the symbol of the commutator For trace class Wick operators, a very useful formula gives Tr (Op W (σ)) = σ(ω, ω) dω ∧ dω.
Remark 4. Let U (t) := e −i Ht andb k (t) := U † (t)b k U (t). Recalling Remark 2, we have thatb k (t) is a Wick operator for any t ≥ 0. We denote now its symbol by σ k (t,ω, ω). Moreover, it is easy to see that ρ k (t,ω, ω) hence the operator Op W (ρ k ) is well defined in the Wick quantization.

Remarks on Phase Space Analysis
Let Ψ (t) be the solution of the quantum dynamics for some fixed initial data Ψ (0) ∈ F B (C L ), for example a normalized function in the N -sector of F B (C L ) whence representing N -particle states (see [23], pg. 48) and Π Ψ (t) the related projection operator. Let us consider the related Wick symbol σ 0 (ω, ω) := σ W ( Π Ψ (0) )(ω, ω) = | φ ω , Ψ(0) | 2 . A way to study the semiclassical localization of the operator as the parameter h = 1/N → 0 is to describe the essential support (see Theorem 8.16 in [41]). In order to do this, one has to study the semiclassical Wave Front set, also called Frequency Set (see for example [31], or [18] for the periodic setting which is the one of the current paper) where ψ α belongs to an orthonormal set, given for example by eigenfunctions labeled by the multi-index α ∈ N L of the number operator N := k b † b k . The phase space localization of this set can be described by the use of the coherent states φ √ Nω and study of the FBI transform φ √ Nω , B(t)ψ α (see Sect. 3.6point 3 of [31]) which reads (thanks to coherent states decomposition) Notice that these are Gaussian-type integrals. Our estimate in Proposition 3 can be obtained in the same way for | φ √ Nω , B(t)φ √ Nv | (i.e., also outside the diagonal ω = v), where instead of the evolved annihilation operator a k (t) we have Π Ψ (t) . Thus, our L 2 -estimates in Theorem 1 with respect to Gaussian measures can be easily adapted to get a semiclassical estimate for the integral (152) with a more general invariant measure. To conclude, we stress that higher-order corrections of (150) implies better localization estimates for the function ω → φ √ Nω , B(t)ψ α in (152). This can be done in the spirit of Egorov Theorem (see [8,16]) by iterating the semigroup identity (70) working on the time evolved Wick symbols governed by Eq. (8).