Spectral localization estimates for abstract linear Schrödinger equations
Spectral localization estimates for abstract linear Schrödinger equations
Abstract.
We study the propagation properties of abstract linear Schrödinger equations of the form , where is a self-adjoint operator and a time-dependent potential. We present explicit sufficient conditions ensuring that if the initial state has spectral support in with respect to a reference self-adjoint operator , then, for some independent of and all , the solution remains spectrally supported in with respect to , up to an remainder in norm. The main condition is that the multiple commutators of and are uniformly bounded in operator norm up to the -th order. We then apply the abstract theory to a class of nonlocal Schrödinger equations on , proving that any solution with compactly supported initial state remains approximately supported, up to a polynomially suppressed tail in -norm, inside a linearly spreading region around the initial support for all .
Key words and phrases:
Schrödinger equations; A priori estimates2020 Mathematics Subject Classification:
35Q41 (primary); 35B40 , 35B45 ,81U90 , 37K06 (secondary)1. Introduction
We consider the following abstract linear Schrödinger equation on a Hilbert space :
(1.1) |
Here is a densely defined, self-adjoint operator on and is a time-dependent potential such that admits bounded propagator on for all times.
Our aim in this paper is to control the spectral localization properties of states evolving according to (1.1). Specifically, fix a reference self-adjoint operator on and let be the spectral projection of onto the half-line . Under suitable conditions on the commutator between and the system Hamiltonian, we prove that any solution to (1.1) with initial state satisfies, for some independent of and all ,
(1.2) |
Estimate (1.2) ensures that if the initial state has spectral support in w.r.t. the reference operator , then the corresponding solution remains spectrally localized in w.r.t. up to a polynomially decaying remainder in norm. See Theorem 2.1 for precise statement and the subsequent remarks for generalizations. The physical significance is discussed in Section 2.2.
We then apply the general theory (1.2) to study spacetime localization properties for some concrete models of nonlocal dispersive equations. We work in and, for , , denote by the characteristic function of ,
We choose the reference operator , the multiplication operator by the distance functions , so that amounts to the multiplication operator by . The Hamiltonians are given by a class of nonlocal non-autonomous operators on , with satisfying suitable finite moment bounds and potential given by arbitrary uniformly bounded multiplication operators. Under condition (4.4) on the kernel of , we prove that if the wave function evolving according to such a nonlocal Schrödinger equation is supported in at , then there holds the dispersive estimate
(1.3) |
From here we conclude the Strichartz-type estimates
(1.4) |
for all and with . Furthermore, by Markov’s inequality and (1.3), we have
(1.5) |
Localization estimates (1.3)–(1.5) impose direct constraints on the size of the probability tails. See Section 4.1 for the precise results.
The novelty of our results in this paper is twofold. Firstly, we identify explicit sufficient conditions ensuring spectral localization estimate (1.2) for abstract Schrödinger equations (see conditions (C1)–(C2) below). In particular, the system Hamiltonian does not enter the conditions directly, but only through its commutators with the reference operator appearing in (1.2). Secondly, comparing to recent results [MR4254070, MR4604685] where similar propagation estimates as (1.2) are obtained for the standard Schrödinger equation and the Hartree equation, our results hold in an abstract Hilbert space setting and complement existent results in the case where the Hamiltonian involves nonlocal operators and the reference operator is not given by specific forms.
1.1. Organization of the paper
In Section 2, we present our main results, illustrate the key technical steps, and discuss the methodology in the context of relevant literature. In Section 3, we furnish the proofs of the main results, Theorems 2.2–2.4, in the general setting.
In Section 4, we illustrate applications to a large class of nonlocal Schrödinger equations on . The models under consideration have the favourable property that they satisfy the main condition laid out in Section 2 with any multiplication operator by functions in the homogeneous Sobolev space (i.e., weakly differentiable with -gradient). These could be viewed as typical models for nonlocal Hamiltonians, as the fractional Laplacian does enjoy the same property, owning to the theory of Calderón commutators.
Notation
We denote by the identity operator, the domain of an operator , and the norm of operators on and sometimes that of vectors in . denotes the space of bounded operators on . We make no distinction in our notation between a function and the associated multiplication operator on .
The commutator of two operators and is first defined as a quadratic form on (always assumed to be dense in ) and then extended to an operator. Similarly, the multiple commutators of and are defined recursively by and for .
2. Setup and main result
Let be a complex Hilbert space equipped with inner product and induced norm . We consider a dynamical system described by (1.1), with defined on a common dense domain for all times. Furthermore we assume (1.1) is globally well-posed on . By standard perturbation theory, a simple sufficient condition is that is a densely defined self-adjoint operator on and is uniformly bounded for all .
We will mainly work in the Heisenberg picture and study the Heisenberg evolution, , for differentiable families , characterized by the duality relation
(2.1) |
where is the unique global solution to the Schrödinger equation (1.1) with .
2.1. Main result
Throughout the paper, we fix a reference self-adjoint operator on , such that is dense in . Our main assumptions are then stated in terms of commutators between the reference operator and the system Hamiltonian in (1.1).
First, for some integer , we assume that ‘almost commutes’ with the potential, in the sense that extends to bounded operators for all times and satisfies, for some fixed ,
(C1) |
Second, we assume that the multiple commutators , all extend to bounded operators on ; namely, there exist such that
(C2) |
To state our main result, let
(2.2) |
so that, by (C1), amounts to the norm of the generalized group velocity operator . Recall that is the spectral projection of onto , explicitly,
(2.3) |
The main result of this paper is the following:
Theorem 2.1 (Spectral localization).
This theorem is proved at the end of Section 2.3
Remark 1.
Remark 2.
Remark 3.
Remark 4.
In [MR1720738], Hunziker, Sigal, and Soffer proved a corresponding ‘minimal escape velocity’ bound which complements estimate (2.4): Suppose in (C1) and condition (C2) holds for some . Assume further for some . Then, by [MR1720738, eq. (1.11)], for any and , there exists s.th. for all ,
(2.7) |
Combining (2.7) and (2.6) with , we obtain a refined two-sided control on the propagation of spectral supports of evolving states: Fix and assume the initial state . Then, for any , the solution generated by remains localized in the linearly expanding ‘annular regions’ , up to an remainder in norm.
2.2. Literature review and discussion
The question of locality versus nonlocality has always been a fundamental issue in the study of physical phenomena and the corresponding mathematical models. On , our results yield, through (1.3)–(1.5), spacetime localization estimates for the propagation of information in a range of non-relativistic quantum mechanical models. Establishing such estimates is a delicate matter, because the dispersive structure of the governing evolution equations generally leads to an apparent lack of locality in the restrictive sense. For instance, suppose one defines, as usual, the maximal propagation speed as the infimum of all ’s such that states initially supported in at time remains supported, for all times , in the light cone . Then, even for the simplest example of a free particle evolving according to the Schrödinger equation , infinite speed of propagation can be observed by examining the Fourier transform of the solution and using the superlinear growth of the dispersion relation. This general idea also leads to infinite speed of propagation, with the usual definition above, for typical -body quantum evolutions described by a large class of dispersive equations [MR1264524, MR1443322].
Historically, for quantum evolutions described by the standard Schrödinger equations with Hamiltonians , where the potential is sufficiently regular, one approach to recover an appropriate sense of locality is to introduce an energy cutoff adapted to the spectrum of on the initial state, and then show that the probability of finding the (microlocalized) state in the classically forbidden region vanishes asymptotically in time. Thus, localization properties of evolving states are reformulated in terms of propagation estimates for certain time-dependent observables identifying the spacetime support of states at time , and finite speed of propagation is established in terms of the resulting propagation estimates. This idea is also the starting point of the present paper.
More precisely, V. Enss proved in his seminal works [En1, En2] that if a particle with unit mass is initially localized in a ball at and has energy below , then the probability, , of finding the particle at time in the classically forbidden region vanishes as a function, i.e., . This way one obtains effective light cones (ELCs), viz., regions outside of which the probability of finding the particle vanishes asymptotically in time. Notice that the ELCs obtained this way are energy-dependent, in agreement with the physical intuition that a particle should move at a speed proportional to the square root of its energy.
The result of Enss was subsequently improved in [SigSof, Skib] and, more recently, [MR4254070, MR4604685], to Schrödinger equations with time-dependent Hamiltonians and the Hartree equation. Historically, such propagation properties have played crucial roles in scattering theory of Schrödinger operators, leading to important breakthroughs in the study of asymptotic completeness of -body problems by Enss [En3, En4], Skibsted [Skib2, Skib3], and Sigal-Soffer [SigSof1, SigSof2, SigSof3, SigSof4], among many others. For reviews of the development in scattering theory along this line, see [HunSig2, HunSig3].
Our approach to derive the main propagation estimate (1.2) is a generalization of the monotonicity method originated from the classical works of Sigal and Soffer’s [SigSof, SigSof1, SigSof2, SigSof3, SigSof4] in scattering theory, in which the authors proved that for general time-dependent Schrödinger equations on , evolving states admit ELCs that spread out in space at a finite rate. The results of the seminal works [SigSof, SigSof2] were improved in [Skib, MR1335378, MR4254070, MR4604685] and extended to non-relativistic QED in [BoFauSig], open quantum systems in [MR4529865, MVBvNL], nonlinear equations in [MR4604685] and condensed matter physics in [FLS1, FLS2, LRSZ]. Our result is motivated and built upon the works mentioned above and more.
2.3. Key steps for proving Theorem 2.1
In this section, we illustrate the key steps in the proof of our main result, the operator inequality (2.4).
Our approach to study the localization properties of quantum evolutions, pioneered in [SigSof, Skib, MR1335378, MR4254070, BoFauSig, FLS1, FLS2], is based on approximate monotonicity formulae for certain propagation-identifying observables. As the starting point, we define, for any differentiable family , the Heisenberg derivative
(2.8) |
It is readily verified that the evolution , defined by relation (2.1), solves the Heisenberg equation
(2.9) |
By the fundamental theorem of calculus and the fact that , it follows that
(2.10) |
Our goal is to establish monotonicity estimates, based on (2.10), for the evolution of suitable time-dependent observables . Roughly speaking, we design so that it satisfies the following properties:
-
•
is constrained by a time-decaying envelope, with .
-
•
is comparable with , with for while .
Combining these, we arrive at the desired estimate (2.4).
Below we elaborate on the key steps in establishing these properties.
2.3.1. Construction of ASTLOs.
We fix a test light-cone slope and let be a large adiabatic parameter. For the reference operator entering (C1)–(C2), times , and smooth cutoff functions in a suitable class (see Figure 1 below), we define the adiabatic spectral localization observables, or ASTLOs (adopting the terminology of [FLS1, FLS2]), as
(ASTLO) |
The precise construction of the class is rather flexible and not so relevant to the present exposition; we defer it to (3.2) where the detailed proof is presented. Note that is an operator-valued function defined by functional calculus of the self-adjoint operator .
Remark 5.
Regarding the nomenclature of ‘ASTLO’, we say that is adiabatic since, for a test light-cone slope and a large adiabatic parameter , the velocity varies at a slow scale. To see that identifies the spectral localization property of states, one can view as a smoothed spectral cutoff function associated to and various half-lines (see Figure 1). Indeed, for functions in appropriate classes and suitably chosen , we can arrange to have
(2.11) | ||||
(2.12) |
The detailed relations are formulated and proved in Proposition 3.4.
2.3.2. Differential inequality for ASTLOs.
By relation (2.9) and the almost-commuting assumption (C1), it is easy to verify that for some absolute constant ,
(H) |
See Lemma 5.5 and (2.3.4) for details. Further exploiting relation (H), we obtain
Theorem 2.2 (Recursive monotonicity estimate).
This theorem is proved in Section 3.1.
Remark 6.
Remark 7.
The differential inequality (RME) is ‘recursive monotone’ because the second term on the r.h.s. is of the same form as the leading term, and the latter is non-positive. Notice that (RME) (more precisely, the proof of it) is the only place where information of the evolution is used. The only property we require of the underlying evolution is the differential identity (H), which asserts that the velocity is approximately determined by the Heisenberg derivative corresponding to the autonomous part in the system Hamiltonian, up to a small time-decaying remainder due to the potential.
2.3.3. Monotonicity estimate for ASTLOs.
Since the evolution is positivity-preserving, control over along yields the same over through relations (2.11)–(2.12). Thus, our goal now is to derive monotonicity estimates for . Indeed, through (RME), we can control the growth of as follows:
Theorem 2.3 (monotonicity estimate).
This theorem is proved in Section 3.2.
Remark 8.
Here we sketch the proof of (ME) for . Integrating (RME), dropping certain non-negative terms, and using the bound on from (C1), we find
(2.14) | ||||
(2.15) |
Since the l.h.s. term in (2.14) is of the same form as the second term on the r.h.s., we can apply the same estimate to the latter, while introducing another cutoff function, say , in the same class as and . Iterating this procedure times on (2.14) until no integral is present in the r.h.s., we obtain
(2.16) |
for cutoff functions . Lastly, we apply (2.16) to bound the integral in the r.h.s. of (2.15) (which carries a prefactor of ). This, together with some additional algebraic properties of the ASTLOs, yields (ME). See Section 3.2 for detailed derivations.
2.3.4. Concluding estimate (2.4).
Theorem 2.3 establishes the approximate monotonicity for ASTLOs as alluded to in Section 2.3. It remains now to conclude the desired propagation estimate (2.4).
2.4. Discussion
We are interested in the spectral localization properties of abstract dispersive evolutions. We have proposed a modular paradigm that accomplishes this goal, as long as the underlying evolution satisfies (H) with a reference operator satisfying (C1)–(C2). Our method is geometric in nature and traces back to the line of works by Enss, Hunziker, Sigal, Skibsted, Soffer, and others, who have laid out the foundation of the geometric method for scattering theory of the Schrödinger operators (see [HunSig2, HunSig3] for reviews).
Indeed, the asymptotic localization theory of general dispersive evolutions (1.1), which we consider here, bears an intrinsic similarity to the scattering theory of the standard Schödinger operators. Both problems concern with the semiclassical behaviour of particles for large time. Notice however that the parameter in (ASTLO), essentially a semiclassical parameter, does not come with the model (1.1), but is determined by the problem directly. The precise choice of is given in (3.48).
A main technical advantage of our localization theory based on the analysis of ASTLOs lies in its flexibility. Whereas strictly monotone quantities along a given evolution equation are rare to find, the approximately monotone ASTLOs are rather easy to engineer. One reason is that the underlying evolution does not enter directly into our analysis, but only through assumptions (C1)–(C2) on the commutators between the reference operator and the system Hamiltonian.
For example, suppose (1.1) is posed on and the reference operator is a multiplication operator by a function with for some . Then, for generic pseudo-differential operators satisfying relevant domain conditions and , we have , and so (C1) is satisfied with . Consequently, the ASTLOs satisfy the differential identity (H) with equality. Since the evolution only enters our analysis through (RME) and the latter depends only on (C2) and (H), we conclude that (RME) and all subsequent modular theorems hold. See Section 4 for more details.
Moreover, since the system Hamiltonian does not enter directly into the main technical assumption (C2), but only through its commutators with in (ASTLO), we can derive conditional localization properties when the commutator assumption (C2) fails for the obvious choice of . Consider the case acting on . Let be a smoothed distance function to a smooth bounded domain . The obvious choice does not satisfies (C2), since is unbounded. However, with an energy cutoff , where and is a smooth cutoff function supported in , one can check that, with the microlocalized position operator , the (microlocal) group velocity (together with higher commutators) is bounded. Using this microlocalized version of in (ASTLO) and running the paradigm above, we obtain energy-dependent spacetime localization estimates for . See [MVBvNL] for concrete results of this nature, with applications to von-Neuman-Linblad equations in Markovian open quantum dynamics. Related propagation bounds involving microlocal cutoff are obtained in [MR4254070, MR4604685] for linear and nonlinear quantum dynamics involving standard Schrödinger operators.
Lastly, since our method is based on monotonicity estimate in the form of operator inequalities, we can reduce localization theory for quantum many-body problems to the corresponding -body problems. Consider an abstract second quantization map, , mapping -body observables acting on to many-body observables acting on a Fock space over . We assume the map is positive-preserving, i.e., for any self-adjoint -body operators , ,
(2.18) |
and, with denoting the many-body evolution of observables on ,
(2.19) |
Then, applying on both sides of (ME) yields the many-body approximate monotonicity estimate
(2.20) |
where is the number operator. See [FLS1, FLS2, LRSZ, LRZ] for related results based on this technique for quantum many-body systems arising from condensed matter physics.
3. Proofs of Theorems 2.2–2.4
In this section, we proved the main results presented in Section 2.3.
We begin with the precise definition of (ASTLO). Fix , together with a densely defined self-adjoint operator . For each , we define a class of observables by functional calculus:
(3.1) |
For a parameter , we define a class as follows:
(3.2) |
Then, for any , the operator is bounded on and non-negative definite, with . Typical examples of functions in are suitably smoothed characteristic functions of as in Figure 1.
In what follows, we will use two properties of the space , which can be readily verified:
-
(X1)
If for some with , then .
-
(X2)
For any and , there exists with and .
In principle, the class could be replaced by suitable classes of functions satisfying the abstract properties (X1)–(X2).
In view of relation (H), to prove Theorems 2.2–2.4, it suffices to derive an estimate for the Heisenberg derivative associated with the free Hamiltonian . Thus, in Sections 3.1 and 3.2, we only work with the free evolution and write and
(3.3) |
so that, with denoting the unitary evolution generated by , the Heisenberg equation (2.9) reads
(3.4) |
3.1. Proof of Theorem 2.2
Let as in (C2) and set . Recall in this subsection denotes the free evolution and .
The main result of this section is the following differential operator inequality:
Theorem 3.1.
Suppose condition (C2) holds for some . Then, for all and , there exists a constant and functions (dropped if ) depending only on , and , such that for all , the following operator inequality holds on :
(3.5) |
(The sum in the r.h.s. is dropped if .)
This theorem is proved at the end of this section. Estimate (3.5), together with property (X2) and relation (H), implies Theorem 2.2.
Remark 9.
Identity (3.4) plays a crucial role in our analysis, and it is precisely in (3.4) that the Hamiltonian structure of (4.1) is used. Indeed, for a heat-type equation with self-adjoint , we have formally, instead of (3.4),
where the brace denotes the anti-commutator. The change renders key expansion formulae below unavailable, and thus new machinery is needed to handle heat type equations. We will not seek to pursue this problem presently.
We begin with the following lemma:
Lemma 3.2.
Suppose the assumption of Theorem 3.1 holds. Then there exist (dropped if ), together with a constant such that the following operator inequality holds on :
(3.6) |
(The sum in the r.h.s. is dropped if .)
Proof.
Within this proof, we fix and write . Also, we set for . (The sign is irrelevant for our argument.)
1. By condition (C2), there exists such that
(3.7) |
This, together with the definition of (see (3.2)), implies that the hypotheses of Lemma 5.4 are satisfied for , and so there hold the commutator expansion
(3.8) |
with some such that (c.f. (5.20)–(5.21))
(3.9) |
Adding commutator expansion (3.8) to its adjoint and dividing the result by two, we obtain
(3.10) | ||||
(3.11) | ||||
(3.12) | ||||
(3.13) |
where the term is dropped for .
2. We first bound the term in line (3.11). Let , which is well defined and lies in by (3.2). Then, by (3.7) and Lemma 5.4, expansion (3.8) also holds for . This expansion, together with the fact that , implies
(3.14) | |||
(3.15) |
where line (3.14) is dropped for and, for some ,
(3.16) |
We will bound the terms in (3.14)–(3.15) using the operator estimate
(3.17) |
For the terms in line (3.14), we use (3.17) with
(3.18) |
yielding
(3.19) |
For the remainder terms in (3.15), we apply (3.17) with
(3.20) |
to obtain
(3.21) |
Combining (3.1) and (3.21) in (3.11) yields
(3.22) | ||||
This bound the term (3.11).
3. For , the term in line (3.12) is bounded similarly as in Step 2. For , we take with
(3.23) |
We claim that for some bounded operator ,
(3.24) |
For this, it suffices to show that
(3.25) |
Using relation (3.23), commutator expansion (3.8), and the fact that , we have
(3.26) |
where the -sum is dropped for and
(3.27) |
for some .
Since on , we have for all and so in line (3.26),
Estimate (3.25) follows from here. Thus we conclude claim (3.1).
Now, we apply estimate (3.17) on the first term on the r.h.s. of (3.1) with
(3.28) |
and then sum over to obtain
(3.29) |
This bounds the term in line (3.12).
4. Plugging (3.22), (3.29) back to (3.10) and using bounds (3.7), (3.9), (3.16), and (3.27), we find that for some ,
(3.30) | |||
Now, for , we choose, with , , from (3.30),
(3.31) |
which is possible since the r.h.s. of (3.31) is supported in by construction. Then the function
(3.32) |
lies in by identity (X1). Thus, by (3.30), the desired estimate (3.6) holds with the choice of from (3.32). This completes the proof of Lemma 3.2. ∎
Proof of Theorem 3.1.
To prove estimate (3.5), we first apply the differential identity (3.4) with for each . This yields
(3.33) |
By definition (3.1), we find
(3.34) |
By estimate (3.6), we find
(3.35) |
where and the second term in the r.h.s. is dropped for . Plugging (3.34) and (3.35) back to (3.33) and using the positive-preserving property of evolution yields (3.5). ∎
3.2. Proof of Theorem 2.3
Recall in this subsection denotes the free evolution and . Our main result is the following:
Theorem 3.3.
Suppose (3.5) holds. Then there exist and a function (dropped if ) depending only on , , and such that for all , the following operator inequality holds on :
(3.36) |
(The second term on the r.h.s. is dropped for .)
Using Estimate (3.36), together with relation (H) and the bound on from assumption (C1), we arrive at Theorem 2.3.
Proof of Theorem 3.3.
Within this proof, we fix and all constants depend only on , , , and . For simplicity, below we consider the case for . For negative times the argument is similar.
1. For ease of notation, for any function , we write
(3.37) |
Note in particular that .
To begin with, we claim the following holds: There exist , (dropped for ), depending only on , and , such that for all ,
(3.38) |
where the sum is dropped if .
To prove (3.38), we bootstrap the recursive monotonicity estimate (3.5). For each fixed , integrating formula (3.4) with in gives
(3.39) |
We apply inequality (3.5) to the second term on the l.h.s. of (3.39) to obtain, after transposing the -term,
(3.40) |
where , , and the second term in the r.h.s. is dropped for .
Since , estimate (3.40) implies, after dropping on the l.h.s., which is non-negative-definite due to the positive-preserving property of the evolution , and multiplying both sides by , that
(3.41) |
where the second term in the r.h.s. is dropped for .
If , then (3.41) gives (3.38). If , we proceed to apply (3.41) to the term up to -th order to get
(3.42) |
where the sum is dropped for and, for ,
Plugging (3.42) back to (3.41), we find
(3.43) |
where the third term in the r.h.s. is dropped for and the functions , for (see (X2)). Bootstrapping this procedure, we arrive at (3.38) for .
Dropping the second term in the l.h.s. of (3.40), which is non-negative since and for all , we obtain
(3.44) |
where the second term is dropped for (in which case we are done). If , then for each , we apply estimate (3.38) to the -th summand in the second term in the r.h.s. of (3.44), with remainder expanded to -th order. This way we obtain
(3.45) |
where the -sum is dropped for , the -sum is dropped if , and , are chosen according to (3.38).
3.3. Proof of Theorem 2.4
Recall that is a densely defined self-adjoint operator on and denotes the spectral cutoff operator defined in (2.3). Our goal now is to choose a function s.th. the geometric inequalities (2.11)–(2.12) hold, whereby eliminating the adiabatic parameter and the ASTLOs from (ME) so as to conclude the desired estimate (2.4).
Our main result is the following proposition:
Proposition 3.4.
Let . For functions and with
(3.47) |
let
(3.48) |
Then the following estimates hold:
(3.49) | ||||
(3.50) |
Proof.
First, by (3.47), we have for . This implies
(3.51) |
where is the characteristic function of the half-line (see Figure 3). Thus (3.49) follows.
Next, again by (3.47), we have for and so, by definition (3.48),
(3.52) |
on the subspace . Since , estimate (3.52) implies
(3.53) |
This completes the proof of Proposition 3.4. ∎
First, for as in the statement of Theorem 2.4, we set
(3.54) |
Fix any . We apply Theorem 2.3 with to get a constant and a function such that
(3.55) |
Next, we apply Proposition 3.4 with
(3.56) |
where the inequalities are ensured by the choice (3.54). The function clearly satisfies condition (3.47). If the function in (3.55), then also satisfy (3.47). (If then we drop the second term in the r.h.s. of (3.55)). Hence, applying (3.49)–(3.50) with and as in (3.48), we conclude the desired estimate, (2.4), from estimate (3.55).
This completes the proof of Theorem 2.4.∎
4. Applications to nonlocal dispersive equations
In this section, we apply the general localization theory laid out in Section 2 to study a large class of nonlocal dispersive evolution models.
We consider the following nonlocal non-autonomous Schrödinger equation:
(4.1) |
Here is a differentiable path of vectors in the Hilbert space . The Hamiltonian consists of a nonlocal part
(4.2) |
for some symmetric (and possibly singular) integral kernel with , together with a time-dependent potential
As a standing assumption, we assume that is self-adjoint on a dense domain and is uniformly bounded for all . Consequently, is self-adjoint on and so, by standard perturbation theory, admits bounded propagator with (see e.g. [GS, Theorem 25.32]).
Our main technical assumption is the following: For some integer and function (i.e., weakly differentiable with the operators
(4.3) |
satisfy, for and some ,
(4.4) |
We show in Appendix A that condition (4.4) amounts to the main technical condition (C2) in the general theory (see Section 2.1) and that a sufficient condition for (4.4) is
(4.5) |
Typical examples of the form (4.2) satisfying (4.5) include the nonlocal diffusion operators
(4.6) |
where is a radial function with profile satisfying
(4.7) |
e.g., with . By interpolation, mild singularity is allowed at , e.g., with .
Condition (4.4) are also verified by certain fractional differential operators. In his seminal work [MR0177312], Calderón proved that (4.4) holds for and , or equivalently,
(4.8) |
The boundedness of commutators of more general singular integral operators and fractional elliptic operators are subsequently established in [MR0412721, MR0358205, MR0763911, MR3286493, MR3547014, MR4443495, MR3555319], among many others, under various conditions on , and for various classes of functions (typically belonging to or BMO). As the scheme below indicates, boundedness of singular integral operators of the form (4.3) would lead to similar propagation estimates in the corresponding dynamical models.
Evolution equations involving nonlocal operators of the form (4.2) have received much research attention in recent years. For recent results concerning evolution equations involving (4.6) subject to similar conditions as (4.7), see e.g. [MR2257732, MR2542582, MR3285829, MR3289358, MR3440113, MR4182983] and, for applications to natural sciences, [MR3469920, MR4409816], as well as the references therein. For regularity theory of nonlocal evolution equations, see [MR3626038, MR3771838, MR3959442]. For an excellent recent review on nonlocal diffusion operators with integrable kernels, see [MR4187861].
Note however that all of the cited works above are concerned with, instead of Hamiltonian evolution equation as in (4.1), gradient flows of the form with of the form (4.2). This distinction should be made clear since the Hamiltonian structure of (4.1) is used crucially in proving the recursive monotonicity estimate (RME) for (wherefore in all us results in Section 2 as well). See Remark 9 below for a discussion.
Eq. (4.1) arises, among others, from the study of nonlinear nonlocal Schrödinger (NLS) equations of the form
(4.9) |
where is a bounded external potential (possibly time-dependent). Eq. (4.9) has a Hamiltonian structure inherited from the nonlocal generalization of the Ginzburg-Landau free-energy functional in the presence of external potential:
Indeed, if solves (4.9), then satisfies (4.1) with bounded for all . This convolution-type model for phase transitions was proposed in [MR1463804] and the associated -gradient flow (the nonlocal Allen-Cahn equation) has been studied in [MR1463804, MR1712445, MR1933014, MR2257732, MR2542582]. See [MR1712445, Sect. 1] for a discussion on the connection between above and the classical Ginzburg-Landau energy functional.
4.1. Results
Under the standing assumption, the evolution of a state from time according to (4.1) is given by
(4.10) |
where is the propagator for in (4.1). The evolution of an observable , dual to the evolution of states w.r.t. the coupling , is given by
(4.11) |
where is the backward propagator.
For and , we denote by and the evolution of states and observables, respectively. Let and be as in (4.4). Our main result in this section is the following:
Theorem 4.1.
Suppose (4.4) holds for and . Then, for every , there exists such that for any function and ,
(4.12) |
Proof.
We derive estimate (4.14) as a consequence of Thms. 2.2–2.3 and Proposition 3.4. Fix and (see (3.2)) with for . Below, all estimates are independent of these parameters.
First, we verify the assumptions of Theorem 2.2. Since in (4.1) with , the evolution condition (H) is satisfied with given by (4.2). By Lemma A.1, the Hamiltonian from (4.2) and verify the commutator condition (C2), with depending on . We have shown that the assumptions of Theorem 2.2 hold. Thus, by Thms. 2.2–2.3, estimate (ME) holds.
Next, define and denote by the ASTLOs from (3.1) with this choice of . Then, by estimate (ME), there exists a constant and a function such that
(4.13) |
Lastly, we use Proposition 3.4. The function clearly satisfies condition (3.47). If the function in (4.1), then also satisfy (3.47). (If then we drop the second term in the r.h.s. of (4.1)). Hence, applying (3.49)–(3.50) with in (4.1) and using that for all (see (2.3)), we conclude the desired estimate, (4.14), from estimate (4.1). This completes the proof of Theorem 4.1. ∎
Theorem 4.1 grants control over the localization of states w.r.t. to a fixed reference geometry, , and a height function . The growth rate of in turn determines the decay estimate of the probability leakage as in (4.1).
Specifically, let and . Taking and using the facts that , , we conclude from Theorem 4.1 that
Corollary 4.1 (Localization of scattering states).
Suppose (4.4) holds for and . Then, for every , there exists such that for all subset , functions , and ,
(4.14) |
To see that (4.14) controls the localization of evolving states according to (4.1), fix and define . Assuming the initial condition at is localized in in the sense that , we conclude from (4.14) that for all .
As a consequence of the localization estimate (4.14), we have the following a priori estimate on the propagation speed of traveling wave solutions to the nonlinear nonlocal Schrödinger equation (4.9):
Corollary 4.2.
Proof.
Since solves (4.1), by freezing coefficients, satisfies (4.14) and therefore we have
(4.15) |
for all . Suppose now . Then, on the one hand, we can choose such that (4.15) holds. On the other hand, since , there is a large depending only on and such that
(4.16) |
The proof of Corollary 4.2 is complete. ∎
5. Technical lemmas
5.1. Remainder estimates
In this section and the next one, we present some estimates and commutator expansions, first derived in [SigSof] and then improved in [HunSig3, Skib] etc. Below, we adapt some of the arguments from [HunSig3] and results from [MVBvNL].
Throughout this section we fix an integer . For integers and smooth functions , we define a weighted norm
(5.1) |
Note that
(5.2) |
and we have the following property:
Lemma 5.1.
Let be an integer. Suppose and there exist such that for ,
(5.3) |
Then there exists depending only on such that
(5.4) |
Proof.
We have
and the integral converges for . ∎
Write and . In what follows, as in [HunSig3, eq.(B.5)], for , we take to be an almost analytic extension of defined by
(5.5) |
where is a cutoff function with for , for , and for all . This induces a measure on as
(5.6) |
In the remainder of this section, we derive integral estimate for various functions against the measure (5.6).
The next result is obtained by adapting the argument in [HunSig3, Lem. B.1]:
Lemma 5.2 (Remainder estimate).
Proof.
For each fixed , we define
(5.10) |
by integrating (5.8) against . Using that for and for or , we find
(5.11) | ||||
(5.12) |
Since and , the integral in line (5.11) converges and can be bounded as
(5.13) |
5.2. Commutator expansions
We frequently use the following result, taken from [HunSig3, Lems. B.2]:
Lemma 5.3.
Let satisfy (5.4) for some . Then for any self-adjoint operator on ,
(5.17) |
where the integral converges absolutely in operator norm and is uniformly bounded in .
Remark 10.
We call equation (5.17) the Helffer-Sjöstrand (HS) representation. The HS representation (5.17), together with the remainder estimate (5.7), implies the following commutator expansion:
Lemma 5.4.
Let . Let satisfy (5.4) with . Let be an operator on . Let be a densely defined self-adjoint operator on . Let for some fixed and all .
Suppose
(5.18) |
Then , and we have the expansion
(5.19) | ||||
(5.20) |
where the remainders are defined by these relations and given explicitly by (5.28)–(5.29).
Moreover, there exists depending only on and , such that
(5.21) |
Remark 11.
Note that needs not to be bounded. By (5.3), it suffices for to have strictly sublinear growth.
Proof of Lemma 5.4.
Within this proof we write with .
Since is bounded, it follows that
(5.22) |
holds in the sense of quadratic forms on . Since is bounded by assumption, the r.h.s. of (5.22) is bounded and so extends to an bounded operator on . Using (5.22), we proceed by commuting successively the commutators to left and right, respectively. Iteratively, we obtain
(5.23) | ||||
(5.24) |
which hold on all of since ’s are bounded operators by assumption (5.18).
Let , be cutoff functions with for , for , and for all . Set . Since , it satisfies (5.4) for all . (Note that itself, a priori, does not satisfy (5.4) with .) Thus the HS representation 5.17 holds with and so
(5.25) |
which holds a priori on . Plugging expansions (5.23)–(5.24) into (5.25) yields
(5.26) | |||
(5.27) |
where
(5.28) | ||||
(5.29) |
Since the operator is bounded independent of , and , we have
(5.30) |
Similarly we could bound the sums in (5.26)–(5.27). Thus we see extends to a bounded operator on for each .
The following lemma is a direct consequence of the estimates proved in Lemma 5.4:
Lemma 5.5.
Let be defined as in (C1). Then for any , , the commutator extends to bounded operators for all times and satisfies, for some absolute constant ,
(5.31) |
Acknowledgments
The Author is supported by National Key R & D Program of China Grant 2022YFA100740, China Postdoctoral Science Foundation Grant 2024T170453, and the Shuimu Scholar program of Tsinghua University. He was also supported by DNRF Grant CPH-GEOTOP-DNRF151, DAHES Fellowship Grant 2076-00006B, DFF Grant 7027-00110B, and the Carlsberg Foundation Grant CF21-0680 during the completion of this paper. He thanks G. Grubb for introduction to the subject matter in Section 4, R. Frank for helpful discussion and pointing out reference [MR0177312], and C. Jiang for useful comments on earlier draft of the paper. He thanks S. Breteaux, J. Faupin, M. Lemm, C. Rubiliani, M. Sigal, and D. Ouyang for fruitful collaborations, and especially M. Lemm for pointing out (1.5) as well as other helpful comments.
Parts of this work were done while the Author was visiting MIT. Earlier version of parts of this work has appeared as a chapter in the Author’s PhD thesis at the University of Copenhagen.
Declarations
-
•
Conflict of interest: The Author has no conflicts of interest to declare that are relevant to the content of this article.
-
•
Data availability: Data sharing is not applicable to this article as no datasets were generated or analysed during the current study.
Appendix A Commutator estimates
In this appendix, we prove that condition (4.5) implies uniform estimates on multiple commutators with (multiplication operator by) Lipschitz functions. In particular, this implies that (4.4) holds with from (4.2) and .
Lemma A.1.
Let . Suppose is an operator acting on as
(A.1) |
for and integral kernel satisfying
(A.2) |
Then for every Lipschtiz function on such that for some ,
(A.3) |
there holds
(A.4) |
Proof.
We first prove that for each fixed and all , we have
(A.5) |
We prove this by a simple induction. Clearly, the term in (A.1) does not contribute to the commutators , since . Hence below we take in (A.1).
For the base case , we compute, for fixed and every ,
Taking the difference yields (A.5) with . Now assume (A.5) holds for . Then we have
Since , taking the difference of the last two expressions yields (A.5) for . This completes the induction.
Formula (A.5), together with the Schur test for integral operators, implies
(A.6) |
Now we compute, using assumptions (A.2) and (A.3), that
This bounds the first term in the r.h.s. of (A.6). Similarly we can derive the same bound for the second term in the r.h.s. of (A.6). Plugging the results back to (A.6) yields estimate (A.4). ∎