Observability of a 1D Schr\"odinger equation with time-varying boundaries
Duc-Trung Hoang

TL;DR
This paper investigates the observability and controllability of a 1D Schrödinger equation with moving boundaries, providing exact conditions for observability in both linear and general moving cases.
Contribution
It offers new precise observability results for Schrödinger equations with time-dependent domains, including boundary and internal observability conditions.
Findings
Exact boundary and pointwise internal observability for linear moving boundaries.
Boundary observability conditions for general moving boundaries under certain curve conditions.
Controllability of the adjoint system established via duality theory.
Abstract
We discuss the observability of a one-dimensional Schr\"odinger equation on certain time dependent domain. In linear moving case, we give the exact boundary and pointwise internal observability for arbitrary time. For the general moving, we provide exact boundary observability when the curve satisfies some certain conditions . By duality theory, we establish the controllability of adjoint system.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsStability and Controllability of Differential Equations · Advanced Mathematical Physics Problems · Numerical methods in inverse problems
Observability of a 1D Schrödinger equation with time-varying boundaries
Duc-Trung Hoang Univ. Bordeaux, Institut de Mathématiques (IMB). CNRS UMR 5251. 351, Cours de la Libération 33405 Talence, France. Email: [email protected]
Abstract
We discuss the observability of a one-dimensional Schrödinger equation on certain time dependent domain. In linear moving case, we give the exact boundary and pointwise internal observability for arbitrary time. For the general moving, we provide exact boundary observability when the curve satisfies some certain conditions . By duality theory, we establish the controllability of adjoint system.
Keywords: Observability, controllability, Schrödinger equation, moving domain, non-autonomous evolution equation.
Mathematics Subject Classification (2010): 93B07, 93C05, 35R37.
1 Introduction
Let , and a strictly positive –function satisfying and . We consider the following system as a initial boundary value problem in a time dependent domain.
[TABLE]
For Neumann boundary observations we obtain estimates like
[TABLE]
see Theorems 2.1, 2.2 and 2.3. We refer to the first estimate as observability estimate and to the second as admissibility estimate. The two first mentioned results rely on a transformation of () to a non-autonomous equation on the fixed domain : the change of variables and new function gives an equivalent differential equation for , namely
[TABLE]
which can easily obtained by the chain rule.
To obtain Theorems 2.1 and 2.2 we apply the ’multiplier technique’: This powerful method has been developped by Morawetz [26] and was later extended by Ho [11] and Lions [17]. We extend a version of Machtyngier [18] to time-dependend multipliers. The observability estimate relies then on the “uniqueness-compacity” lemma 3.5. The pitfall of this proof strategy is that it only proves existence of some positive constant, without explicit estimates. This is in contrast with Theorem 2.3 which is as specific result for the boundary curve . In this linear moving wall case, we mimic a successful approach for a one-dimensional wave-equation obtained by Haak and the author in [9] and develop the solution of () into a series of eigenfunctions. This allows to use results from Fourier analysis; the obtained admissibility estimates are sharper than those obtained in the previous results, and the observation estimate is provided with explicit constants. Moreover, we obtain in this case admissibility and exact observability of internal point observations:
[TABLE]
see Theorem 2.5. It is remakable that the lower estimte cannot be true when on any rational point ; the fact that the considered domains extend however, seem to ’middle out’ this obstacle. Closely related to this observation are works of Castro and Khapalov [6, 14, 13] where on a fixed domain a moving point observer is considered, with similar conclusions. We also mention results from Moyano [28, 29] where in a two-dimensional circle the radius is used as a control parameter.
An additional result on -admissibility and observability of point observations are presented as well, see Theorem 2.7.
It is well-known that exact observability for an (autonomous) wave equation implies observability for the associates Schrödinger equation, see e.g. [41, Chapter 6.7 ff.]. An inspection of the proof gives several obstacles when one passes to non-autonomous problems, and we were not able to use this approach to directly infer our results from those for the wave equation in [9]. We mention that some results on the so-called Hautus-test will be subject of an independent publication [10].
2 Main Results
Before giving precise formulations of the aforementioned results, let us start by proving that the Schrödinger equation () admits a solution: to this end, we reformulate it as an abstract non-autonomous Cauchy problem in the following way: let and the family of operators be defined as
[TABLE]
wich natural domain . Moreover, by assumption, the map is continuously differentiable for all . Let . Then integration by parts gives
[TABLE]
Taking real parts and observing that
[TABLE]
we obtain
[TABLE]
For \omega>\bigl{\|}\tfrac{\ell^{\prime}}{2\ell}\bigr{\|}_{L^{\infty}}, the left hand side of (2.3) becomes positive, and the Lumer-Philips theorem asserts that generates a contraction semigroup, i.e.
[TABLE]
This ensures in particular that the family satisfies the Kato stability condition. We apply [30, Theorem V.4.8 pp.145] to conclude that generates a unique evolution family on satisfying . From this we infer a solution to () as well, by transforming the fixed domain back to the time-dependent domain.
Suppose that we are given observation operators where is another Hilbert space. Define the output function . The operator is called -admissible if there exist such that:
[TABLE]
We say that the system () is exactly -observable in time if there exist such that:
[TABLE]
If the spaces are fixed, we simply speak of admissibility and exact observability. Exact observation in time means that the knowledge of allows to recover the initial value . It is well known that exact observability is equivalent to exact controllability of the retrograde adjoint system:
[TABLE]
Moreover, it is easy to see that admissibility or observability of () is equivalent to those of ().
Results on Neumann observations
Theorem 2.1**.**
Let and be a strictly positive, twice continuously differentiable function satisfying and . Then there exists a constants such that the following admissibility inequalities hold:
[TABLE]
An explicit estimate of constant is given in the proof, see (3.15).
Concerning observability, we will have the following result. Let and be a strictly positive, twice continuously differentiable function satisfying:
[TABLE]
Integrating for [math] to of the second condition, we have . From the condition , is an increasing function, and then . It follows that , and so the condition guaranteeing admissibility is satisfied.
Theorem 2.2**.**
For all satisfying (2.4), the following observability inequality holds:
[TABLE]
Here is some positive constant depending on .
A direct application of theorem 2.2 can be used for periodic moving boundary where and . For all \tau\in\Bigl{(}0,\frac{\pi}{2\omega}\Bigr{)}, we have
[TABLE]
[TABLE]
Hence, satisfies the condition (2.4), so the curve is admissible. The problem of particles moving inside one dimensional square-well of oscillating width was proposed by Fermi and Ulam [19] in order to explain the mechanism of particles containing high energies. This model that plays an important role on theory of quantum chaos and it seems difficult to give an exact solution formula. Glasser [8] investigated the behavior of wave functions and energy in a given instantaneous eigenstate by assumptions on the smoothness of boundary. As far as we know, there are no results in the literature concerning observability and controllability with periodic boundary functions.
In the case that , the condition is ensured when and 0<t<\tfrac{1}{\varepsilon}\bigl{(}\tfrac{2}{\varepsilon\pi}-1\bigr{)}. We have the following exact analytic solution for , due to Doescher and Rice [7]
[TABLE]
where the coefficients are defined by the sine-series development of the initial value . A similar exact solution in the case of two-variable moving wall can be found in [42] where the author uses the fundamental transformation to change the moving boundary problem into a solvable one side fixed boundary problem.
Based on formula (2.5) we obtain a first result on Neumann observability at the boundary . Compared to Theorem 2.2 the admissibility constant is sharper. In contrast with Theorem 2.2, where we can only prove existence of some positive constant , we obtain now an explicit estimate for the observability constant. The proof is presented in section 3.
Theorem 2.3**.**
For every there exist explicit constants such that:
[TABLE]
In particular, the Neumann observation at the boundary of the system () is exact observable in any time . Moreover, the observability coefficient decays \thicksim\exp\Bigr{(}\frac{-2k\pi^{2}}{\varepsilon\tau}\Bigr{)} where .
Remark 2.4**.**
By Dirichlet condition for all . Differentiating yields , and so . As a result, observing or is, up to a constant, the same.
Point observations
We now focus on point observations in the case of a linearly moving wall . Observe that in the “degenerate” case that is, , the (then) autonomous Schrödinger equation has the well-known solution
[TABLE]
Clearly, there is no reasonable observability possible at rationals points since infinitely many terms in the sum vanish, independently of the leading coefficient . This changes when : from (2.5) we obtain
[TABLE]
and so
[TABLE]
Based on a remarkable result of Tenenbaum and Tucsnak we obtain the following result in section 3.
Theorem 2.5**.**
Assume . Then, for every , we have:
[TABLE]
More precisely, where and are some positive constants that appear in to proof.
Corollary 2.6**.**
For all the point observation for the system () is exactly observable in arbitrary short time.
-estimates of point observations
Finally we have to following admissibility and observability estimates.
Theorem 2.7**.**
Let . We assume that . For and , we have
[TABLE]
where are constants depending on and .
The upper estimate is a direct consequence of (2.8). Indeed, by the continuity of the embeddings and the boundedness of to obtain:
[TABLE]
Hence, it serves only to show that the lower estimate is of the right order.
3 Proof of the main results
3.1 The multiplier Lemma
We follow E. Machtyngier [18, Lemma 2.2] by using multiplier method for (): Let be a solution to () and be a real valued function. Then, due to the differential equation (),
[TABLE]
We separate the left hand side of (3.1) into three parts and simplify each of them.
Lemma 3.1**.**
The following identities hold.
[TABLE]
Proof.
To prove (3.4), we use integration by parts. Using , we have:
[TABLE]
Therefore, the left hand side of (3.4) equals
[TABLE]
Here, we already use the fact that
[TABLE]
To prove (3.8) we have
[TABLE]
since we use . Again, integration by parts shows
[TABLE]
Therefore we have:
[TABLE]
Hence, part (3.8) is proved. The last part is obvious. ∎
Now summing up the three parts and using (3.1) yields
Proposition 3.2**.**
For any real valued function and a solution to () we have
[TABLE]
3.2 Energy estimates
For a solution to () we define the first and second energy as
[TABLE]
respectively.
Lemma 3.3**.**
We have .
Proof.
Taking the derivative respected to and using , we have
[TABLE]
Now integration by parts gives
[TABLE]
whereas
[TABLE]
Therefore,
[TABLE]
so that
[TABLE]
Using , this implies easily . ∎
Lemma 3.4**.**
For all and \tau\in\Bigl{(}0,\tfrac{\pi}{2\omega}\Bigr{)}, we have:
[TABLE]
Proof.
Concerning we have
[TABLE]
The first term on the right hand side simplifies as
[TABLE]
whereas the second term simplifies as follows.
[TABLE]
We add both parts to obtain
[TABLE]
By Variation of constants, we get an explicit solution:
[TABLE]
One easily obtains an upper bound, namely . For the lower bound, we use the Poincaré (or Wirtinger) inequality on to obtain,
[TABLE]
∎
3.3 Admissibility of Neumann observations at the boundary
Proof of Theorem 2.2.
We take the function on satisfying and . By Proposition 3.2, we have
[TABLE]
Therefore, we have
[TABLE]
where we estimate all five terms separately. Concerning , we separate the products in the real part by , then use Lemmata 3.3 and 3.4 to obtain
[TABLE]
The second term is easily estimated by Lemma 3.3:
[TABLE]
Part is decoupled by Cauchy-Schwarz and then estimated using Lemma 3.4 as follows:
[TABLE]
For the forth part, we use Lemma 3.4 to obtain
[TABLE]
Finally, part is treated like part :
[TABLE]
Summing up all three estimates, we obtain
[TABLE]
where the constant is given by
[TABLE]
Replacing in yields the admissibility inequality:
[TABLE]
The second admissibility estimate follows the same lines, using on with and . ∎
3.3.1 Neumann Observability at the Boundary
Recall the following lemma
Lemma 3.5**.**
Let and be the Hilbert spaces. We consider the continuous linear operators , and such that is compact, is bounded below and:
[TABLE]
Then the kernel of has finite dimension and
Proof.
A similar proof can be found in [38, Lemma 1 pp.1] where we just replace by . ∎
Proof of Theorem 2.2.
For all satisfying , we choose two positive constants and such that:
[TABLE]
We choose where . Proposition 3.2 is then equivalent to:
[TABLE]
Taking the three last formula of the right hand side to the left, then taking the absolute to get:
[TABLE]
The sum of third and fourth terms in the right hand side of above formula can be estimated as:
[TABLE]
Due to the energy estimate in lemma 3.3 and 3.4, we have the upper bound for the second term:
[TABLE]
As a result, we combine these estimation and use (3.12) to obtain:
[TABLE]
where the last inequality come from (3.17). Therefore, there exist the constants and such that:
[TABLE]
It is sufficient to prove that there exist a constant such that
[TABLE]
Let us denote the operator from to and the operator from to that maps:
[TABLE]
[TABLE]
From admissibility and (3.19), we have:
[TABLE]
It is easy to see that is compact operator due to Rellich’s embedding lemma. In order to use the unique-compactness lemma 3.5 for , we need to check that is injective. Observe that means that satisfies () with Dirichlet conditions and zero Neumann derivative. It is well known that vanishes in this case, see for example [39, Theorem 3] or [12, Corollary 6.1]. As a consequence,
[TABLE]
for some constants .
∎
3.4 Results for linear moving walls
Recall the Doescher-Rice representation formula (2.5) that yields for
[TABLE]
and denote by
[TABLE]
For all fixed , the functions form an orthonormal basis in , since the change of variable reduces to the standard trigonometric system on .
Lemma 3.6**.**
For all finitely supported sequences we have the following relation between and the norms of the initial data .
[TABLE]
Proof.
Observe that
[TABLE]
Since is a finite sequence we may interchange differentiation and summation and obtain
[TABLE]
so that, squaring real and imaginary parts, we find
[TABLE]
Lemma 3.7**.**
Let and , then the functions for form an orthonormal system in .
Proof.
Note that \bigl{(}\tfrac{t}{\ell(t)}\bigr{)}^{\prime}=\tfrac{\ell(t)-t\ell^{\prime}(t)}{\ell(t)^{2}}=\tfrac{1}{\ell(t)^{2}}. Therefore, the obvious change of variable reduces to a standard trigonometric function on . Observe that . Now orthonormality easily follows. ∎
Observe that the above sequence is not an orthonormal basis. Indeed, with , we have for all .
3.4.1 Neumann observation at the Boundary
Proof of Theorem 2.3.
We start considering only the first term at . As in the proof of Lemma 3.6 we consider for a moment only initial data associated with finitely supported sequences . Differentiating the representation formula (2.5) term by term yields
[TABLE]
and therefore
[TABLE]
Using the monotonicity of in , we have where
[TABLE]
This allows to focus only on the integral , where we abbreviate and make a change of variable . Letting , the above double inequality rewrites as
[TABLE]
The sequence satisfies the hypotheses of [40, Theorem 3.1 and Corollary 3.3] so that, for all and
[TABLE]
On the other hand side, if , we have by periodicity and Parseval’s identity
[TABLE]
We conclude by Lemma 3.6 that
[TABLE]
This inequality being true for all leading to finitely supported sequences , it is true for any by density.
For second term at , we see for finitely supported sequences that
[TABLE]
Taking the -norm, one get the equivalent between and
[TABLE]
Clearly, the rest proof follows the lines above. ∎
3.4.2 Internal Point Observability
Proof of Theorem 2.5.
Since for all ,
[TABLE]
By definition, \sin\bigl{(}\tfrac{n\pi a}{\ell(t)}\bigr{)}=\frac{1}{2i}\bigl{(}\exp({i\tfrac{n\pi a}{\ell(t)}})-\exp({-i\tfrac{n\pi a}{\ell(t)}})\bigr{)}. Therefore,
[TABLE]
For , we extend the series by , and . The sequence is regular and satisfies the hypotheses of [40, Theorem 3.1] with and . We follow the lines of the proof of Theorem 2.3: changing the variable gives with the notation ,
[TABLE]
we write and use [40, Corollary 3.3] with :
[TABLE]
For the upper estimate, we use similar method as in theorem (2.3). More precisely,
[TABLE]
where be the integer number such that with .
3.4.3 -admissibility and observability
Proof of Theorem 2.7.
The upper estimate yielding is obtained by interpolation of the two upper estimates in Theorem 2.3. We are left with the lower estimate. Since , , and so by the Cauchy-Schwarz inequality. Let and let which is chosen to satisfy . By Hölder’s inequality we then have
[TABLE]
From trivial argument on boundedness of and :
[TABLE]
Combining with the estimate (3.25), one get:
[TABLE]
From inequalities (3.25) and (3.26) and Theorem (2.5) we deduce now
[TABLE]
Since , the result follows. ∎
4 Boundary controllability of dual problem
Since we have already stated several theorems that can be interpreted as exact observation we will briefly sketch the duality theory that allows to rephrase these assertions in terms of exact control, then the solution to adjoint problem
[TABLE]
satisfies by injection of the respective differential equations of and . Hence exact observability implies that the Gramian satisfies to the effect that has closed image. Moreover, if , taking scalar product with reveals , so is injective and hence has dense range. By the open mapping theorem, is therefore an isomorphism on . This means that the adjoint problem (4.1) can be steered to any state by an appropriate choice of the initial value . Indeed, for we have
[TABLE]
It turns out that in our case . So exact observation of the Schrödinger equation () can be reformulated as exact control for the Schrödinger equation with zero final time. We turn back to these ideas after stating our first theorem. In the case of linear moving , let be given by where . The (lower) estimate in theorems 2.3 and 2.2 then reformulates as exact observability of for the non-autonomous Cauchy problem (2.1). Some care has to be taken since is unbounded on . Indeed, is given by , then we obtain exact controllability of (4.1) in a distributional sense:
[TABLE]
Multiplying with a test function , and integrating on we obtain by partial integration
[TABLE]
This is possible for any test function only if the point evaluation vanishes. The dual statement of the lower estimate in theorems 2.3 and 2.2 is thus exact controllability of a Schrödinger equation with Dirichlet control on the right boundary,
[TABLE]
We reverse back to the moving boundary problem by taking and . Then the problem can be written as:
[TABLE]
or
[TABLE]
In general situation of satisfying condition , one take be given by . Therefore, the dual operator is given by . Using similar argument, we obtain exact controllability of a Schrödinger equation with Dirichlet control applied on both of boundaries
[TABLE]
Acknowledgements
This work was undertaken as part of the authors PhD thesis at the University of Bordeaux. The author wants to thank his advisers Bernhard Haak and El-Maati Ouhabaz for their patiently supports, great motivations, and immense knowledges.
Further, we kindly acknowledge valuable discussions with Marius Tucsnak.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] S.Agmon. On the eigenfunctions and on the eigenvalues of general elliptic boundary value problems, Comm. Pure. Appl. Math , 15, 119-147 (1962).
- 2[2] K.Beauchard. Controllability of a quantum particle in a 1D variable domain, ESAIM , Vol 14, No.1, 105-147 (2008).
- 3[3] K.Beauchard, J.M.Coron. Controllability of a quantum particle in a moving potential well, J. Functional Analysis , 232, 328-389 (2006).
- 4[4] N. Burq. Controle de l’equation des plaques en presence d’obstacles strictement convexes. Mem. Soc. Math. France , (N.S.) No. 55 (1993).
- 5[5] P.Cannarsa, G.D.Prato, J-P. Zolesio. Heat equation in non-cylindrical domains, Analisi matematica , 08, 73-77, (1989).
- 6[6] Carlos Castro, Exact controllability of the 1-D wave equation from a moving interior point , ESAIM Control Optim. Calc. Var. 19 (2013), no. 1, 301–316.
- 7[7] S.W.Doescher, M. H. Rice. Infinite Square-Well Potential with a Moving Wall, Am. J. Phys , 37, 1246, (1969).
- 8[8] L. Glasser, J. Nergo, L. Nieto. Quantum infinite square well with an oscillating wall , Chaos, Solitons, and Fractals , Volume 41, Issue 4 (2009).
