# Backstepping Control for the Schrödinger Equation with an Arbitrary Potential in a Confined Space

Authors: Pedro Franco; A. Schaum; Roberto Bernal Jaquez
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

In this work the control design problem for the Schrödinger equation with an arbitrary potential is addressed. In particular a controller is designed which (i) for a space-dependent potential steers the state probability density function to a prescribed solution and (ii) for a space and state-dependent potential exponentially stabilizes the zero solution. The problem is addressed using a backstepping controller that steers to zero the deviation between the initial probability wave function and the target probability wave function. The exponential convergence property is rigorously established and the convergence behavior is illustrated using numerical simulations for the Morse and the Pöschl-Teller potentials as well as the semilinear Schrödinger equation with cubic potential.


---

## Body

## 1. Introduction

The problem of quantum control is an important problem in today’s quantum technology with relevant applications in quantum information systems, molecular chemistry, and atomic physics, among others [1–3]. Even though control theory of classical mechanical systems has been extensively developed and applied in recent years, its application to quantum mechanical systems bears some additional hurdles due to particularities of the quantum world which in particular impose restrictions on the controllability of these systems [4–6]. In spite of the progress made in this field, still more research and understanding of the underlying quantum phenomena are necessary in order to fully achieve control in the quantum mechanical realm [7–9]. A particularly interesting approach to the boundary control of Partial Differential Equations (PDEs) in general and Schrödinger equations in particular is the backstepping method [10, 11], yielding exponential stabilization of the L2-norm in the case of a potential-free Schrödinger equation. This method is used in a variety of engineering problems [12, 13] and has been recently extended to the case of regulator design for a Schrödinger equation with a general purpose space-dependent potential and subject to distributed and boundary perturbations [14]. It is an interesting question, which to the authors knowledge still remains open, how these approaches perform for particular potentials of interest, like the Morse potential [15–17] or the Pöschl-Teller potential [18–21]. Furthermore, the aforementioned results are restricted to potentials which are not a function of the state, i.e., to the case of linear Schrödinger equations. Nevertheless, important application examples require the consideration of state-dependent potentials, e.g., to describe the evolution of solitons in optical fibers [22, 23] or the description of deep water waves [24]. This leads to a semilinear Schrödinger equation. As it is well-known, the application of the backstepping approach to semilinear systems is only resolved for particular cases [25–29]. Accordingly, the problem of controlling a semilinear Schrödinger equation is still an open task.In this work, we consider the problem of steering a Probability Density Function (PDF) into a target PDF. It should be noted that the use of the PDF control perspective has not been widely explored [30]. In particular, for space-dependent potentials convergence to an arbitrary desired solution is achieved, and for state-dependent potentials, i.e., the case of a semilinear Schrödinger equation, the zero solution is exponentially stabilized. This task is resolved by using boundary actuation following the backstepping method [10]. In contrast to transforming the complete Schrödinger equation into a target PDE, the system is viewed as an interconnection of a linear dynamic with a nonlinear static subsystem, just as in the framework of absolute stability analysis [31, 32], and the backstepping transformation is performed only for the linear subsystem.The paper is organized as follows. In Section2, the problem is stated from the perspective of PDF control and the initial and target systems are given. In Section 3, the boundary control is designed and the main results of this work are given. In Section 4, to illustrate our findings, we set the control for three significant physical systems in order to show the validity of the method developed in this paper. In Section 4.1, we control a Schrödinger system with a Morse potential, and a system with a Pöschl-Teller potential is examined in Section 4.2. In Section 4.3 the exponential stabilization of the solution for the semilinear Schrödinger equation is considered. Conclusions are given in Section 5.


## 2. Problem Formulation

Consider the Schrödinger equation with a nonlinear potentialV and boundary control(1a)i∂tψx,t=-∂x2ψx,t+Vx,ψx,tψx,t,ψx,0=ψ0x(1b)∂xψ0,t=0(1c)ψ1,t=utwith the restriction(1d)∫01ψx,tψx,t¯dx=1,∀t≥0.In (1a), (1b), (1c), and (1d), x∈[0,1] is the spatial coordinate describing the position of a confined particle within a bounded space or cavity, t is the time, ψ(x,t)∈L2([0,1]×R)→C is the complex-valued wave function that gives the quantum state of the system and ψ(x,t)¯ is its conjugate complex counterpart, V(x,ψ(x,t)) is the real-valued nonlinear potential function determining the movement of the particle within the cavity, and u(t)∈C∞(R+) is the control input. The function V is continuously differentiable with respect to both arguments.Equation (1d) states that the system is confined into x∈[0,1] and that the probability of finding the particle within the cavity is conserved for all times. Stated in terms of the PDF(2)ρx,t=ψx,tψx,t¯this means that(3)∫01ρx,tdx=1,∀t≥0.The problems addressed in this paper consist in the design of a control lawu=ϖ(t,ψ) such that the wave function ψ(x,t) converges exponentially in the L2-norm to a target wave function ψd(x,t), i.e., that there exist constants M≥1,κ>0 such that, for all t≥0,(4)ψ-ψd=∫01ψx,t-ψdx,t2dx1/2≤Mψ0-ψd0e-κt,where ψ0(x)=ψ(x,0), ψd0(x)=ψd(x,0) denote the initial profiles of ψ(x,t) and ψd(x,t), respectively, and the target wave function ψd(x,t) is a solution of the same Schrödinger equation (1a), (1b), (1c), and (1d) with the control input u(t)=ud(t). Two different cases are considered:(I)
For the case of a space-dependent potentialV(x) an arbitrary target system can be considered with the state ψd(x,t) associated with a quantum system with the same potential but different boundary and initial conditions ψd(1,t)=ud(t), ∂xψd(0,t)=0, and ψd(x,0)=ψd0.(II)
For the case of a space and state-dependent potentialV(x,ψ(x,t)) the target state is given by the zero solution ψd(x,t)=0∀x∈[0,1],t≥0, which corresponds to the solution of (1a), (1b), (1c), and (1d) with initial condition ψd0(x)=0,∀x∈[0,1], and the control input ud(t)=0,∀t≥0.Introducing the deviation variables(5)ψ~x,t=ψx,t-ψdx,t,u~t=ut-udtwhere in case (I) ψd(x,t)≠0 and in case (II) ψd(x,t)=0, ud(t)=0, it can easily be verified that the associated Schrödinger equations for ψ~(x,t) have the form(6a)i∂tψ~x,t=-∂x2ψ~x,t+Vx,ψ~x,tψ~x,t,ψ~x,0=ψ~0x(6b)∂xψ~0,t=0(6c)ψ~1,t=u~t.The associated error PDF is given by(7)ρ~x,t=ψ~x,tψ~x,t¯=ψx,t-ψdx,t2≥0together with the probability of a nonzero estimation error(8)Pt=∫01ρ~x,tdx.In terms of the PDFρ~(x,t) and the probability P(t) the convergence requirement (4) reads(9)Pt≤Me-κtP0,P0=∫01ρ~x,0dxand corresponds to the convergence in the L1-norm of ρ~ to zero.


## 3. Control Design

For the purpose of control design, the dynamics (6a), (6b), and (6c) are interpreted as an interconnection of a linear dynamical subsystem L given by (10a)-(10d) with distributed input ν~(x,t) and output σ~(x,t), and a nonlinear static subsystem N given by (10e) is defined as follows:(10a)i∂tψ~x,t=-∂x2ψ~x,t+ν~x,t,ψ~x,0=ψ~0x(10b)∂xψ~0,t=0(10c)ψ~1,t=u~t(10d)σ~x,t=ψ~x,t(10e)ν~x,t=Vx,σ~x,tσ~x,t.Note that this kind of interpretation of semilinear (or nonlinear) systems in terms of interconnected subsystems is standard in absolute stability theory [31, 32] and corresponds to the nonlinear feedback loop sketched in Figure 1.Figure 1
Interconnection scheme (10a), (10b), (10c), (10d), and (10e) with L given by (10a)-(10d) and N given by (10e).Accordingly, the main idea behind the approach in the present paper consists in designing a feedback control for the linear part in such a way that the semilinear system given by the two-subsystem interconnection has an exponentially stable zero solution.


### 3.1. Open-Loop Dynamics

Before addressing the control design problem, the stability properties of the open-loop dynamics are analyzed. The following result establishes the instability ofρ~(x,t)=0 and will further be very useful in the subsequent developments.Lemma 1.
Foru~(t)=0, the probability that the wave function ψ(x,t) does not converge to the desired wave function ψd(x,t) is conserved over time; i.e.,(11)∀t≥0:Pt=∫01ρ~x,tdx=K,K=constant.Accordingly, the solution ρ~(x,t)=0 for the error PDF (7) is unstable.Proof.
A direct calculation shows that(12)dPtdt=∫01∂tρ~x,tdx=∫01∂tψ~x,tψ~x,t¯+ψ~x,t∂tψ~x,t¯dx=∫01i∂x2ψ~x,t-iVx,ψ~x,tψ~x,tψ~x,t¯++ψ~x,t-i∂x2ψ~x,t¯+iVx,ψ~x,tψ~x,t¯dx. Taking into account the fact that V(x,ψ(x,t))∈R, it holds that (13)Vx,ψx,t=Vx,ψx,t¯∈R. Thus, it follows that (14)dPtdt=i∫01∂x2ψ~x,tψ~x,t¯-ψ~x,t∂x2ψ~x,t¯dx+-i∫01Vx,ψ~x,tψ~x,tψ~x,t¯-ψ~x,tVx,ψ~x,tψ~x,t¯dx. The second term is clearly zero, and integration by parts of the first term and substitution of the boundary conditions (6b), (6c) yield (15)∫01i∂x2ψ~x,tψ~x,t¯-ψ~x,t∂x2ψ~x,t¯dx=i∂xψ~x,tψ~x,t¯-ψ~x,t∂xψ~x,t¯01-i∫01∂xψ~x,t∂xψ~x,t¯-∂xψ~x,t∂xψ~x,t¯dx. The second term is again zero and from the first term it holds that(16)dPtdt=2I∂xψ~1,tu~t,and for u~(t)=0 it results that the probability P is conserved over time, implying that ρ~(x,t) does not converge to zero or, equivalently, that ρ~d(x,t)=0 is unstable.Remark 2.
From the preceding proof it follows that the conservation property (11) holds independently of the particular shape of the potential V(x,ψx,t).


### 3.2. Control of the Linear Subsystem

Note that the control of the linear subsystem withν~(x,t)=0 corresponds to controlling the linear Schrödinger equation without potential, as depicted schematically in Figure 2.Figure 2
Interconnection scheme for the control of the linear Schrödinger systemL with ν~(x,t)=0 and controller C.The problem of stabilizing the zero solution of the linear Schrödinger equation has been solved in [10, 11]. The purpose of this section is to recall the results from [10, 11] and to put them in perspective for the analysis of the linear-nonlinear subsystem interconnection (10a), (10b), (10c), (10d), and (10e).Forψ~, a solution of (10a), (10b), (10c), (10d), and (10e) (with ν~(x,t)=0) considers the state transformation T(17)ϕ~x,t=Tψ~x,t=ψ~x,t-∫0xkx,ξψ~ξ,tdξand the target dynamics for ϕ~(x,t) given by(18a)i∂tϕ~x,t=-∂x2ϕ~x,t-iγϕ~x,t,ϕ~x,0=ϕ~x,0(18b)∂xϕ~0,t=0(18c)ϕ~1,t=0.A straightforward calculation (see [10, 11] for details) shows that the dynamics of ϕ~(x,t) corresponds to (18a), (18b), and (18c) if and only if the kernel k(x,ξ) of the integral transformation (17) satisfies the PDE (19a)∂x2kx,ξ+∂ξ2kx,ξ=-iγkx,ξ(19b)∂ξkx,0=0(19c)kx,x=-γ2ix.The solution of this PDE is given by (see [10, 11, 33])(20)kx,ξ=iγxI1-iγx2-ξ2-iγx2-ξ2,where I1(x) represents the first-order modified Bessel function of the first kind. From (17) it follows that the control input u~(t) is given by(21)u~t=∫01k1,ξψ~ξ,tdξ=∫01iγI1-iγ1-ξ2-iγ1-ξ2ψ~ξ,tdξ.For the dynamics (18a), (18b), and (18c), it is known from [11] that the following result holds true.Lemma 3.
The solutionϕ~(x,t) of (18a), (18b), and (18c) converges exponentially to zero; i.e., there exist constants N≥1,κ>0 such that (22)ϕ~x,t≤Nϕ~x,0e-κt.Proof.
Consider theL2 norm of ϕ~(23)ϕ~x,t=Et,Et=∫01ϕ~x,tϕ~x,t¯dx and recall that it corresponds to the L1 norm of the probability density function (24)ρϕ~x,t=ϕ~x,tϕ~x,t¯.From Lemma 1, it follows that, for γ=0, the conservation of probability implies that dE(t)/dt=0. For γ≠0, it follows from a short calculation that (25)dEtdt=∫01∂tϕ~x,tϕ~x,t¯+ϕ~x,t∂tϕ~x,t¯dx=-2γ∫01ϕ~x,tϕ~x,t¯dx=-2γEtimplying that (26)Et=E0e-2γt or equivalently that(27)ϕ~x,t=ϕ~x,0e-γtmeaning that ϕ~(x,t) converges exponentially to zero; i.e., inequality (4) holds for ϕ~ with N=1 and rate κ=γ>0.Given the invertibility of the integral transformation (17) (see [11]), it holds that ψ~ also converges exponentially to zero in the L2-norm or, equivalently, ρ~ converges exponentially to zero in the L1-norm. This is summarized in the following corollary.Lemma 4.
Let the controlu~(t) be chosen according to (21) for some γ>0. Then, for ν~(x,t)=0, it holds that the deviation state ψ~(x,t) defined in (5) converges exponentially to zero and the probability P(t) defined in (8) satisfies(28)Pt≤Me-2γtP0for some M>0; i.e., (9) is satisfied.Proof.
From [11] it is known that the inverse transformation for (17) exists, is unique, and is given by (29)ψ~x,t=T-1ϕ~x,t=ϕ~x,t+∫0xlx,ξϕ~ξ,tdξ with a bounded kernel l(x,ξ) (see, e.g., [11] for details). Both transformations, T and its inverse T-1, are linear and bounded operators, which implies that the following inequalities hold: (30)ψ~x,t2=T-1ϕ~x,t2≤T-1op2ϕ~x,t2≤Ne-2γtT-1op2ϕ~x,02=Ne-2γtT-1op2Tψ~x,02≤Ne-2γtT-1op2Top2ψ~x,02, where ·op denotes the norm of an operator, defined by (see [34]) (31)Top=sup⁡Tψ:ψ=1, and in the second inequality, Lemma 3 was used. This implies that (32)ψ~x,t2≤Me-2γtψ~x,02 or, equivalently, that (33)Pt=ψ~x,t2≤Me-2γtP0 with M=NT-1op2Top2≥1.From the preceding lemma the following useful corollary is obtained.Corollary 5.
Foru~(t) given by (21) the differential equation (34)dPtdt=2u~I∂xψ~1,twith ψ~(x,t) being the solution of the Schrödinger equation (10a), (10b), (10c), (10d), and (10e) with ν~(x,t)=0 has a solution P(t) which satisfies (28).Proof.
It follows from (16) and Lemma 4.To illustrate the functioning of the backstepping controller, in Figure3 the real and imaginary parts of the solution ϕ~(x,t) of (18a), (18b), and (18c) are shown for the initial condition(35)ϕ~x,0=Ae-σx-x02e-ix,with A=2, σ=10, and x0=0.5. It can be seen that both real and imaginary parts of ϕ~(x,t) converge exponentially to zero, as stated in Lemma 3.Figure 3
Real and imaginary parts of the solutionϕ~ of (18a), (18b), and (10c).
(a)
Real part ofϕ~ (b)
Imaginary part ofϕ~


#### 3.2.1. Arbitrary Potential Case

In this subsection, the main result of this work is obtained by proving that, in the presence of an arbitraryV(x,ψ(x,t)), the controller (21) stabilizes the origin for the dynamics (6a), (6b), and (6c), which in turn means that ψ(x,t) converges exponentially to ψd(x,t). In terms of the system-interconnection structure presented in Figure 4, this means that the controller yields absolute stability of the linear subsystem with respect to the nonlinear feedback.Figure 4
Interconnection scheme (10a), (10b), (10c), (10d), and (10e) with L given by (10a)-(10d), N by (10e), and C by (21).The main result is stated in the next theorem.Theorem 6.
For an arbitrary potentialV(x,ψ~(x,t)), the error state ψ~(x,t) satisfying (6a), (6b), and (6c) converges exponentially to zero if the control u~(t) is chosen according to (21), i.e.,(36)ψ~x,t2=∫01ρ~x,tdx≤Me-2γtψ0~x2with M given in (28).Proof.
From (16) in the proof of Lemma, 1 it follows that (37)dPtdt=2u~I∂xψ~1,t holds for (10a), (10b), (10c), (10d), and (10e) for an arbitrary potential V(x,ψ~(x,t)). Thus, the rate of change of the probability does not depend on the potential and, as a consequence of Corollary 5, by applying the control (21) with γ>0 it follows that (38)Pt≤Me-2γtP0, with M given in (28). This implies that ψ~(x,t) converges exponentially to zero.


## 3.1. Open-Loop Dynamics

Before addressing the control design problem, the stability properties of the open-loop dynamics are analyzed. The following result establishes the instability ofρ~(x,t)=0 and will further be very useful in the subsequent developments.Lemma 1.
Foru~(t)=0, the probability that the wave function ψ(x,t) does not converge to the desired wave function ψd(x,t) is conserved over time; i.e.,(11)∀t≥0:Pt=∫01ρ~x,tdx=K,K=constant.Accordingly, the solution ρ~(x,t)=0 for the error PDF (7) is unstable.Proof.
A direct calculation shows that(12)dPtdt=∫01∂tρ~x,tdx=∫01∂tψ~x,tψ~x,t¯+ψ~x,t∂tψ~x,t¯dx=∫01i∂x2ψ~x,t-iVx,ψ~x,tψ~x,tψ~x,t¯++ψ~x,t-i∂x2ψ~x,t¯+iVx,ψ~x,tψ~x,t¯dx. Taking into account the fact that V(x,ψ(x,t))∈R, it holds that (13)Vx,ψx,t=Vx,ψx,t¯∈R. Thus, it follows that (14)dPtdt=i∫01∂x2ψ~x,tψ~x,t¯-ψ~x,t∂x2ψ~x,t¯dx+-i∫01Vx,ψ~x,tψ~x,tψ~x,t¯-ψ~x,tVx,ψ~x,tψ~x,t¯dx. The second term is clearly zero, and integration by parts of the first term and substitution of the boundary conditions (6b), (6c) yield (15)∫01i∂x2ψ~x,tψ~x,t¯-ψ~x,t∂x2ψ~x,t¯dx=i∂xψ~x,tψ~x,t¯-ψ~x,t∂xψ~x,t¯01-i∫01∂xψ~x,t∂xψ~x,t¯-∂xψ~x,t∂xψ~x,t¯dx. The second term is again zero and from the first term it holds that(16)dPtdt=2I∂xψ~1,tu~t,and for u~(t)=0 it results that the probability P is conserved over time, implying that ρ~(x,t) does not converge to zero or, equivalently, that ρ~d(x,t)=0 is unstable.Remark 2.
From the preceding proof it follows that the conservation property (11) holds independently of the particular shape of the potential V(x,ψx,t).


## 3.2. Control of the Linear Subsystem

Note that the control of the linear subsystem withν~(x,t)=0 corresponds to controlling the linear Schrödinger equation without potential, as depicted schematically in Figure 2.Figure 2
Interconnection scheme for the control of the linear Schrödinger systemL with ν~(x,t)=0 and controller C.The problem of stabilizing the zero solution of the linear Schrödinger equation has been solved in [10, 11]. The purpose of this section is to recall the results from [10, 11] and to put them in perspective for the analysis of the linear-nonlinear subsystem interconnection (10a), (10b), (10c), (10d), and (10e).Forψ~, a solution of (10a), (10b), (10c), (10d), and (10e) (with ν~(x,t)=0) considers the state transformation T(17)ϕ~x,t=Tψ~x,t=ψ~x,t-∫0xkx,ξψ~ξ,tdξand the target dynamics for ϕ~(x,t) given by(18a)i∂tϕ~x,t=-∂x2ϕ~x,t-iγϕ~x,t,ϕ~x,0=ϕ~x,0(18b)∂xϕ~0,t=0(18c)ϕ~1,t=0.A straightforward calculation (see [10, 11] for details) shows that the dynamics of ϕ~(x,t) corresponds to (18a), (18b), and (18c) if and only if the kernel k(x,ξ) of the integral transformation (17) satisfies the PDE (19a)∂x2kx,ξ+∂ξ2kx,ξ=-iγkx,ξ(19b)∂ξkx,0=0(19c)kx,x=-γ2ix.The solution of this PDE is given by (see [10, 11, 33])(20)kx,ξ=iγxI1-iγx2-ξ2-iγx2-ξ2,where I1(x) represents the first-order modified Bessel function of the first kind. From (17) it follows that the control input u~(t) is given by(21)u~t=∫01k1,ξψ~ξ,tdξ=∫01iγI1-iγ1-ξ2-iγ1-ξ2ψ~ξ,tdξ.For the dynamics (18a), (18b), and (18c), it is known from [11] that the following result holds true.Lemma 3.
The solutionϕ~(x,t) of (18a), (18b), and (18c) converges exponentially to zero; i.e., there exist constants N≥1,κ>0 such that (22)ϕ~x,t≤Nϕ~x,0e-κt.Proof.
Consider theL2 norm of ϕ~(23)ϕ~x,t=Et,Et=∫01ϕ~x,tϕ~x,t¯dx and recall that it corresponds to the L1 norm of the probability density function (24)ρϕ~x,t=ϕ~x,tϕ~x,t¯.From Lemma 1, it follows that, for γ=0, the conservation of probability implies that dE(t)/dt=0. For γ≠0, it follows from a short calculation that (25)dEtdt=∫01∂tϕ~x,tϕ~x,t¯+ϕ~x,t∂tϕ~x,t¯dx=-2γ∫01ϕ~x,tϕ~x,t¯dx=-2γEtimplying that (26)Et=E0e-2γt or equivalently that(27)ϕ~x,t=ϕ~x,0e-γtmeaning that ϕ~(x,t) converges exponentially to zero; i.e., inequality (4) holds for ϕ~ with N=1 and rate κ=γ>0.Given the invertibility of the integral transformation (17) (see [11]), it holds that ψ~ also converges exponentially to zero in the L2-norm or, equivalently, ρ~ converges exponentially to zero in the L1-norm. This is summarized in the following corollary.Lemma 4.
Let the controlu~(t) be chosen according to (21) for some γ>0. Then, for ν~(x,t)=0, it holds that the deviation state ψ~(x,t) defined in (5) converges exponentially to zero and the probability P(t) defined in (8) satisfies(28)Pt≤Me-2γtP0for some M>0; i.e., (9) is satisfied.Proof.
From [11] it is known that the inverse transformation for (17) exists, is unique, and is given by (29)ψ~x,t=T-1ϕ~x,t=ϕ~x,t+∫0xlx,ξϕ~ξ,tdξ with a bounded kernel l(x,ξ) (see, e.g., [11] for details). Both transformations, T and its inverse T-1, are linear and bounded operators, which implies that the following inequalities hold: (30)ψ~x,t2=T-1ϕ~x,t2≤T-1op2ϕ~x,t2≤Ne-2γtT-1op2ϕ~x,02=Ne-2γtT-1op2Tψ~x,02≤Ne-2γtT-1op2Top2ψ~x,02, where ·op denotes the norm of an operator, defined by (see [34]) (31)Top=sup⁡Tψ:ψ=1, and in the second inequality, Lemma 3 was used. This implies that (32)ψ~x,t2≤Me-2γtψ~x,02 or, equivalently, that (33)Pt=ψ~x,t2≤Me-2γtP0 with M=NT-1op2Top2≥1.From the preceding lemma the following useful corollary is obtained.Corollary 5.
Foru~(t) given by (21) the differential equation (34)dPtdt=2u~I∂xψ~1,twith ψ~(x,t) being the solution of the Schrödinger equation (10a), (10b), (10c), (10d), and (10e) with ν~(x,t)=0 has a solution P(t) which satisfies (28).Proof.
It follows from (16) and Lemma 4.To illustrate the functioning of the backstepping controller, in Figure3 the real and imaginary parts of the solution ϕ~(x,t) of (18a), (18b), and (18c) are shown for the initial condition(35)ϕ~x,0=Ae-σx-x02e-ix,with A=2, σ=10, and x0=0.5. It can be seen that both real and imaginary parts of ϕ~(x,t) converge exponentially to zero, as stated in Lemma 3.Figure 3
Real and imaginary parts of the solutionϕ~ of (18a), (18b), and (10c).
(a)
Real part ofϕ~ (b)
Imaginary part ofϕ~


### 3.2.1. Arbitrary Potential Case

In this subsection, the main result of this work is obtained by proving that, in the presence of an arbitraryV(x,ψ(x,t)), the controller (21) stabilizes the origin for the dynamics (6a), (6b), and (6c), which in turn means that ψ(x,t) converges exponentially to ψd(x,t). In terms of the system-interconnection structure presented in Figure 4, this means that the controller yields absolute stability of the linear subsystem with respect to the nonlinear feedback.Figure 4
Interconnection scheme (10a), (10b), (10c), (10d), and (10e) with L given by (10a)-(10d), N by (10e), and C by (21).The main result is stated in the next theorem.Theorem 6.
For an arbitrary potentialV(x,ψ~(x,t)), the error state ψ~(x,t) satisfying (6a), (6b), and (6c) converges exponentially to zero if the control u~(t) is chosen according to (21), i.e.,(36)ψ~x,t2=∫01ρ~x,tdx≤Me-2γtψ0~x2with M given in (28).Proof.
From (16) in the proof of Lemma, 1 it follows that (37)dPtdt=2u~I∂xψ~1,t holds for (10a), (10b), (10c), (10d), and (10e) for an arbitrary potential V(x,ψ~(x,t)). Thus, the rate of change of the probability does not depend on the potential and, as a consequence of Corollary 5, by applying the control (21) with γ>0 it follows that (38)Pt≤Me-2γtP0, with M given in (28). This implies that ψ~(x,t) converges exponentially to zero.


## 3.2.1. Arbitrary Potential Case

In this subsection, the main result of this work is obtained by proving that, in the presence of an arbitraryV(x,ψ(x,t)), the controller (21) stabilizes the origin for the dynamics (6a), (6b), and (6c), which in turn means that ψ(x,t) converges exponentially to ψd(x,t). In terms of the system-interconnection structure presented in Figure 4, this means that the controller yields absolute stability of the linear subsystem with respect to the nonlinear feedback.Figure 4
Interconnection scheme (10a), (10b), (10c), (10d), and (10e) with L given by (10a)-(10d), N by (10e), and C by (21).The main result is stated in the next theorem.Theorem 6.
For an arbitrary potentialV(x,ψ~(x,t)), the error state ψ~(x,t) satisfying (6a), (6b), and (6c) converges exponentially to zero if the control u~(t) is chosen according to (21), i.e.,(36)ψ~x,t2=∫01ρ~x,tdx≤Me-2γtψ0~x2with M given in (28).Proof.
From (16) in the proof of Lemma, 1 it follows that (37)dPtdt=2u~I∂xψ~1,t holds for (10a), (10b), (10c), (10d), and (10e) for an arbitrary potential V(x,ψ~(x,t)). Thus, the rate of change of the probability does not depend on the potential and, as a consequence of Corollary 5, by applying the control (21) with γ>0 it follows that (38)Pt≤Me-2γtP0, with M given in (28). This implies that ψ~(x,t) converges exponentially to zero.


## 4. Application Examples

In this section, the stabilization of the Schrödinger equation is studied for three different potentials. In the first case, the Schrödinger equation with a Morse potential [15] is considered. This potential has been extensively used to model the interatomic potential energy of the diatomic molecule in the study of anharmonic vibrations [15–17]. Coupled Morse potentials have been used to model bond interactions and describe the vibrational states of complex molecules [20]. In Section 4.2 the control of the Schrödinger equation with a Pöschl-Teller potential [18] is studied. The study of this potential is important in the description of out of plane molecular vibrations [20], diatomic vibrations [21], and atomic and neutron scattering modeling, among others [19]. Both potentials are shown in Figure 5.Figure 5
Potentials used in examples.
(a)
Morse potential withDe = 10, a = 4, and xe = 0.3 (b)
Pöschl-Teller potential withλ = 6, a = 10, and xe = 0.5In Section4.3 a cubic state-dependent potential is considered which corresponds to the study of solitons in optical fibers [22, 23] or the description of deep water waves [24].For the associated stabilization tasks, the control law (21) is applied to the initial system (1a), (1b), (1c), and (1d) by setting (39)ut=udt+u~twith(40)udt=2sin⁡2t,γ=10.The resulting Schrödinger equations are solved using the time-splitting finite difference (TSFD) method [35].


### 4.1. Morse Potential

Consider the Schrödinger equation with a Morse potential(41)Vx=De1-e-ax-xewith depth De and width a, and an initial condition as defined in (35) (with A=100/π4, σ=50 and x0=0.3), i.e., a Gaussian function modeling a wave packet.The target wave function is defined as the solution of the same equation with a translated initial conditionψ0d with respect to the initial system ψ0 (x0=0.7) and boundary condition ud(t)=2sin⁡(2t)Making the state transformation (17) on the error variable (5) with k(x,ξ) given in (20) and choosing the control u~(t) according to (21), ψ~(x,t) satisfies (36); i.e., ψ~x,t→0. This, in turn, means that ψ(x,t)-ψd(x,t)→0.Figure6 shows the initial system, the target system, and the error system, along with the controlled system for ψ which can be seen to be converging to ψd.Figure 6
Numerical simulations for the case in which the Morse potential is used. (a) Initial system. (b) Target system. (c) Error system. (d), (e), (f) Comparison between|ψ|2 (red) and |ψd|2 (dotted in blue) at times t=0.04, t=0.6, and t=1, respectively.
(a) (b) (c) (d) (e) (f)


### 4.2. Pöschl-Teller Potential

We consider now the stabilization of the Schrödinger equation (1a), (1b), (1c), and (1d) with a Pöschl-Teller potential given by(42)Vx=-λλ+12sech2⁡ax-xewith the same initial and boundary conditions used in the example of the Morse potential. Using the same target wave function as in the previous example and error variable (5), we make the state transformation (17) with k(x,ξ) given in (20). Once more, choosing the control u~(t) according to (21), ψ~(x,t) satisfies (36); i.e., ψ~x,t→0. This, in turn, means that ψx,t-ψdx,t→0 as in our last example. The stabilization for this system is shown in Figure 7.Figure 7
Numerical simulations for the case in which the Pöschl-Teller potential is used. (a) Initial system. (b) Target system. (c) Error system. (d), (e), (f) Comparison between|ψ|2 (red) and |ψd|2 (dotted in blue) at times t=0.04, t=0.6, and t=1, respectively.
(a) (b) (c) (d) (e) (f)


### 4.3. Semilinear Schrödinger Equation

In the semilinear Schrödinger equation, the potential is given by(43)Vx,ψx,t=-ψx,t2ψx,tand the target wave function is the zero solution ψd(x,t)=0. The initial condition chosen for this simulation is a soliton centered in x=0.2(44)ψ0=2e-ix-.2/2cosh⁡20x-.2,and the kernel was chosen equal to the preceding cases. Simulation results are shown in Figure 8.Figure 8
Numerical simulations for the case of the semilinear Schrödinger equation with cubic potential (43).
(a)
Open-loop system (b)
Controlled system


## 4.1. Morse Potential

Consider the Schrödinger equation with a Morse potential(41)Vx=De1-e-ax-xewith depth De and width a, and an initial condition as defined in (35) (with A=100/π4, σ=50 and x0=0.3), i.e., a Gaussian function modeling a wave packet.The target wave function is defined as the solution of the same equation with a translated initial conditionψ0d with respect to the initial system ψ0 (x0=0.7) and boundary condition ud(t)=2sin⁡(2t)Making the state transformation (17) on the error variable (5) with k(x,ξ) given in (20) and choosing the control u~(t) according to (21), ψ~(x,t) satisfies (36); i.e., ψ~x,t→0. This, in turn, means that ψ(x,t)-ψd(x,t)→0.Figure6 shows the initial system, the target system, and the error system, along with the controlled system for ψ which can be seen to be converging to ψd.Figure 6
Numerical simulations for the case in which the Morse potential is used. (a) Initial system. (b) Target system. (c) Error system. (d), (e), (f) Comparison between|ψ|2 (red) and |ψd|2 (dotted in blue) at times t=0.04, t=0.6, and t=1, respectively.
(a) (b) (c) (d) (e) (f)


## 4.2. Pöschl-Teller Potential

We consider now the stabilization of the Schrödinger equation (1a), (1b), (1c), and (1d) with a Pöschl-Teller potential given by(42)Vx=-λλ+12sech2⁡ax-xewith the same initial and boundary conditions used in the example of the Morse potential. Using the same target wave function as in the previous example and error variable (5), we make the state transformation (17) with k(x,ξ) given in (20). Once more, choosing the control u~(t) according to (21), ψ~(x,t) satisfies (36); i.e., ψ~x,t→0. This, in turn, means that ψx,t-ψdx,t→0 as in our last example. The stabilization for this system is shown in Figure 7.Figure 7
Numerical simulations for the case in which the Pöschl-Teller potential is used. (a) Initial system. (b) Target system. (c) Error system. (d), (e), (f) Comparison between|ψ|2 (red) and |ψd|2 (dotted in blue) at times t=0.04, t=0.6, and t=1, respectively.
(a) (b) (c) (d) (e) (f)


## 4.3. Semilinear Schrödinger Equation

In the semilinear Schrödinger equation, the potential is given by(43)Vx,ψx,t=-ψx,t2ψx,tand the target wave function is the zero solution ψd(x,t)=0. The initial condition chosen for this simulation is a soliton centered in x=0.2(44)ψ0=2e-ix-.2/2cosh⁡20x-.2,and the kernel was chosen equal to the preceding cases. Simulation results are shown in Figure 8.Figure 8
Numerical simulations for the case of the semilinear Schrödinger equation with cubic potential (43).
(a)
Open-loop system (b)
Controlled system


## 5. Conclusions

In this work a control law for the Schrödinger equation with an arbitrary space- and state-dependent potential has been designed to steer the state evolution to match a target wave function in a bounded domain. For the case of a state-independent potential the target wave function is an arbitrary solution of the Schrödinger equation for a given initial state and input signal, whereas for the case of a state-dependent potential the target wave function is given by the zero solution. The control design is carried out by interpreting the semilinear dynamics as a two-subsystem interconnection of a linear dynamic subsystem with a nonlinear static subsystem and following the backstepping design approach for the linear subsystem. The closed-loop stability of the target wave function is rigorously established in terms of the probability of convergence.In order to illustrate our findings, three different cases of the Schrödinger equation with physically significant potentials have been considered: (i) a Morse potential, (ii) a Pöschl-Teller potential, and (iii) a cubic state-dependent potential. In these examples, numerically solved using the time-splitting finite difference method, it is shown that the designed control makes an initial probability density function to approach and match a target probability density function.Future research will focus on the question of how to stabilize nonzero target wave functions for the case of the semilinear Schrödinger equation and to include the observer design problem for output-feedback control.



---

*Source: 1826783-2018-07-19.xml*