linear Schrodinger equation

This post is about the existence and uniqueness of the solution to the linear Schrodinger equation. That is to say, $latex i\partial_{t}u+\Delta u=0, (t,x)\in\mathbb{R}\times\mathbb{R}^{d}, u(t,x)\in\mathbb{C}$ with the initial condition u |_{t=0}=u_{0} The main tool is due to Strichartz, in particular his norms. This post, serves to be an introduction to this method in using the linear Schrodinger equation.

First of all, since this equation is linear, the first idea comes to mind is to use Fourier transforms with respect to x. In fact, there is a group structure in the evolution process. That is what we are going to do right now.

Take the Fourier transform of the above Schrodinger equation, then we obtain that

i \partial_t\widehat{u} + (-i \xi )^{2} \hat{u}=0 where we denote \hat{u} the Fourier transform of u. So, we get that \partial_{t} \hat{u} +i \xi^{2} \hat{u} =0. We solve this ordinary equation very easily, that is \hat{u}(t,\xi ) =e^{i \xi^{2} t}\widehat{u_{0}}(\xi ). Then we take the inverse Fourier transform, we obtain that u(t,x)=S_{t} \ast u_{0} where \widehat{S_{t}} =e^{i \xi^{2} t} ( t>0 ) ,\widehat{S_{0}} =\delta_{0} the Dirac distribution at the origin. The group structure of this evolution can be seen very easily using the Fourier transform. The real problem is the regularity.

This is not very hard. Note that, we have $latex |\hat{u} (t,\xi )|=|\overline{u_{0}}(\xi)|$. Have we seen something like this before which concerns the norm of the Fourier transform of a function? Yes, in the Sobolev spaces. In deed, if we suppose that u_{0} \in H^{s} ( \mathbb{R}^{d} ), then we see easily that u(t,x) \in H^{s} (\mathbb{R}^{d} ) due to the above comment for fixed time t. Note that here we do not require that the time t be non-negative. This also guarantees that the evolution process hasa group structure. Now we can have a closer look at the operator S_{t}. We have seen that it is in fact an evolution operator. This said, we have that if there are two initial conditions, u_{0} ,v_{0}, then using the norm <u,v>=\int_{\mathbb{R}^{d}} ( 1+ |\xi |^{2} )^{s} \bar{u}vd\xifor the Hilbert space H^{s}, we have that < S_{t} u,S_{t} v>=< u,v>. Using the group structure, we see that the conjugate of S_{t} is just S_{t}^{\ast}=S_{-t}. All the way around, we have not seen any explicite formulation of S_{t}. In fact, we have seen that \widehat{S_{t}} =e^{i \xi^{2} t}. In general, we can not calculate the inverse Fourier transform of e^{i \xi^{2}t} directly, since it doesn’t lie in L^{1} (\mathbb{R}^{d} ). So, we must solve this problem in the sense of distribution. We can construct a sequenceof Gaussian-like functions f_{n} (\xi ) =e^{({it} -1/n )\xi^{2}}.This f_{n} has an inverse Fourier transform. And, this sequence converges to f=e^{i \xi^{2} t} in the topology of S' ( \mathbb{R}^{d}) (that is to say, for any g \in S ( \mathbb{R}^{d}), we have (f_{n} ,g ) \rightarrow ( f,g ) due to the dominant convergence theorem).After some calculations, we get that S_{t} = \frac{1}{( 4 \pi{it})^{d/2}} e^{i \frac{x^{2}}{4t}}. What is more, we verify easily that as t\rightarrow 0, we have that S_{t} \rightarrow S_{0} = \delta_{0} in the topology of S' ( \mathbb{R}^{d} ). This, combined with the comments above about the isometricity of S_{t}, shows that the application \mathbb{R}\rightarrow H^{s} ( \mathbb{R}^{d} ) ,t \mapsto S_{t} u_{0} is continuous. This result is very natural, and in some sense it has to be so, since we aresolving a differential equation, so the solution u_{t} must be continuous with respect to the temporal variable. Now we have that \hat{u} ( t, \xi ) =e^{i\xi^{2} t} \widehat{u_{0}} ( \xi). What does this expression mean? Note that, for the initial data, u_{0}, its Fourier transform has an amplitude \widehat{u_{0}} ( \xi ) for the frequency \xi. But for the time t, the absolute value of this amplitude doesn’t change, yet it oscillates much, in fact it oscillates in a fashion proportional to \xi^{2}. We can interpret this oscillation as the propagation of the component at the frequency \xi of the original wave. In other words, the higher the frequency, the faster its propagation.

Yet the problem is to control this ‘dispersion’. But how? We should use which norms? It is a very good yet hard question. Perhaps we can use the L^p norms. Which L^p? For example, for the L^2, we know already that the norm doesn’t change with time, that is ||S_tu_0||_{L^2}=||u_0||_{L^2}. There is another norms, which we can use, it is the L^1 or L^{\infty}. Note that S_tu_0 is in fact a convolution operator, so this reminds us of the inequality of Young, that is ||S_tu_0||_{L^{\infty}}\leq ||S_t||_{L^{\infty}}||u_0||_{L^1}. Note that, ||S_t||_{L^{\infty}}=\frac{1}{(4\pi t)^{d/2}}. So, we have some estimations on the L^2 and L^{\infty}, so now we can use the interpolation theorems on exposants, that is, if 1/p=a/2+(1-a)/\infty, 1/p'=a/2+(1-a)/1, then we have that ||S_tu_0||_{L^p}\leq 1^a(\frac{1}{(4\pi t)^{d/2}})^{1-a}||u_0||_{L^{p'}}. Note that 1/p'-1/p=1-a, so we have that ||S_tu_0||_{L^p}\leq 1^a(\frac{1}{(4\pi t)^{d/2}})^{1/p'-1/p}||u_0||_{L^{p'}}.

This estimation is for every instant t. Can we say something global? That is to say, can we have some norm to measure the whole time global behavior of u(t,x)?

This is what Strichartz did. He defines a norm on the space C^{\infty}(\mathbb{R}\times\mathbb{R}^d). That is, for a smooth function f(t,x), we define ||f||_{L^q_tL^r_x}=(\int_{\mathbb{R}}||u(t, .)||_{L^r_x}^qdt)^{1/q} for q<\infty and similarly for q=\infty. Is this norm surprising? Not really. First, we want to know the L^r information of u(t,x)for each instant t. So, we get a positive valued function on t. Then we want to control this function using some norm L^q. So, this is a very natural approach. Now perhaps some dimensional analysis is helpful. Suppose that f_{a,b}(t,x)=f(at,bx), then we can verify easily that ||f_{a,b}||_{L^q_tL^r_x}=a^{-1/q}b^{-d/r}. So, if we want that there be some inequality like ||S_tu_0||_{L^q_tL^r_x}\leq C||u_0||_{L^p}, we must also have it valid for u_{0,b}(x)=u_0(bx). Note that, for the second term, we have that ||u_{0,b}||_{L^p}=b^{-d/p}||u_0||_{L^p}. While for the first term, we have that S_tu_{0,b}=S_{tb^2}u_0(bx)(this is really due to the fact that in the expression of S_t(x), in the exponent, there is x^2/t, while in the rational term, there is 1/t^{d/2}). Thus, according to the above remark, we have that ||S_tu_{0,b}||_{L^q_tL^r_x}=b^{-2/q-d/r}||S_tu_0||_{L^q_tL^r_x}. So, for the inequality to be held for any reasonable u_0, we must have that -2/q-d/r=-d/p, or 2/p+d/r=d/p. The easiest case for L^p is of course L^2. And thus, we have the inequality of Strichartz,

If u_0\in L^2(\mathbb{R}^d), then we have that ||S_tu_0||_{L^q_tL^r_x}\leq C||u_0||_{L^2} for some constant C independent of u_0 with 2/p+d/r=d/2 and (q,r,d)\neq(2,\infty,2),q>2.

Note that this result says that the linear operator, T:H=L^2(\mathbb{R}^d)\rightarrow B=L^q(\mathbb{R},L^q(\mathbb{R}^d)) is continuous, where it is easy to verify that H is a Hilbert space and B is a Banach space with the norms defined above. This important observation reminds us to use something operator theory. To show that T is bounded, we can in some way consider its adjoint, T^*. That is to say, T^*:B'\rightarrow H'=\overline(H). There is a general result from operator theory,

Suppose that T:H\rightarrow B is a linear operator from a Hilbert space(its scalar product is <u,av>=a<u,v>(a\in\mathbb{C}), and for the scalar product induced on H', we write it <u',av'>=\overline{a}<u',v'>(u',v'\in H'=\overline{H})) to a Banach space(the action of a linear functional on the elements is (f,b)(f\in B',b\in B)), and its adjoint T^*:B'\rightarrow H' defined by <T^*(b'),h'>=(b',\overline{T(h')}). Then we have that ||TT^*||_{L(B',B)}=||T||_{L(H,B)}^2=||T^*||_{L(B',H')}.

This is not hard to prove. We just have to utilize the expression like ||T(h)||=\sup_{||b'||\leq1}|(b',T(h))| and things like that. We will not prove this result here. One remark is that the composition T^*T doesn’t make sense in general.

Now we want to see what T^*,TT^* are in our context. It is not hard to see that B'=L^{q'}(\mathbb{R},L^{r'}(\mathbb{R}^d)). One one hand, <T^*(b'),h'>_{H'}=\int_{\mathbb{R}^d}T^*(b')\overline{h'}dx. One the other hand, <T^*(b'),h'>_{H'}=(b',\overline{T(h')})_{B',B}=\int_{\mathbb{R}\times\mathbb{R}^d}=\int_{\mathbb{R}}<S_th',b'>_Hdt.But we have seen above that <S_th',b'>_H=<h',S_{-t}b'>_H. So, we get that <T^*(b'),h'>_{H'}=\int_{\mathbb{R}} <h',S_{-t}b'>_H=<h',\int_{\mathbb{R}}S_{-t}b'(t,.)dt>, thus we see that T^*(b')=\int_{\mathbb{R}} S_{-t}b'(t,.)dt. Thus, TT^*(b')=S_t(\int S_{-s}b'(s,.))ds=\int S_{t-s}b'(s,.)ds, just a convolution. So to prove the continuity of T, we have to show the continuity of T^* or TT^*. But for T^*, perhaps the difficulty is the same as T. So, we first consider TT^*, the most complex one. For f\in B', we have that TT^*(f)=\int S_{t-s}f(s,.)ds. So, we first evaluate ||TT^*(f)||_{L^r_x}\leq\int_{\mathbb{R}}||S_{t-s}f(s,.)||_{L^r_x}dt. Note that for almost all s, f(s,.)\in L^{r'}(\mathbb{R}^d), so using the estimation on the dispersion above, we get that ||S(t-s)f(s,.)||_{L^r}\leq c\frac{1}{|t-s|^{d/2}})^{1/r'-1/r}||f(s,.)||_{L^{r'}}. Now we get another convolution, ||TT^*(f)||_{L^r_x}\leq c\int_{\mathbb{R}})\frac{1}{|t-s|^{d/2}})^{1/r'-1/r}||f(s,.)||_{L^{r'}_x}ds=c\int_{\mathbb{R}}(\frac{1}{|t-s|^{d/2}})^{1/r'-1/r}h(s)ds. Note that \frac{d}{2}(1/r'-1/r)=\frac{d}{2}(1-2/r)=d/2-d/r=2/q, so all in all, we have that ||TT^*(f)(t)||_{L^r_x}\leq c\int \frac{1}{|t-s|^{2/q}}h(s)ds=c(\frac{1}{|s|^{2/q}}\ast h(s))(t). Now comes the last thing, ||TT^*(f)||_B\leq c||\frac{1}{|t|^{2/q}}\ast h(t)||_{L^q_t}.

Now we have a convolution which involves a special function 1/|t|^{2/q}. This function doesn’t lie in any L^p. Perhaps we can use the inequality of Hardy-Littlewood-Sobolev,

If a\in(0,n), (p,w)\in(1,\infty)^2 such that 1/p+a/n=1+1/w then for all f\in L^p(\mathbb{R}^n), we have that ||\frac{1}{|x|^a}\ast f(x)||_{L^w(\mathbb{R}^n)}\leq C||f||_{L^p(\mathbb{R}^n)} for some constant  C independent of f.

Now apply this theorem to our context, here n=1,a=2/q, w=q,p=q'(indeed h(t)=||f(t,.)||_{L^r_x} lies in L^q(\mathbb{R})). Besides, we have that 1/p+a/n=1/q'+2/q=1+1/q=1+1/w. What a relief, it satisfies this inequality, and we conclude that ||TT^*(f)||_B\leq cC||h(t)||_{L^q_t}=cC||f||_{B'}. That is to say, the operator TT^*:B'\rightarrow B is continuous, and so is T:H\rightarrow B. And we proved the result. Note that in the above argument, we use implicitly that q\neq\infty(a condition also in the Hardy-Littlewood-Sobolev inequality). For the case, q=\infty, we must have that r=2, thus ||TT^*(f)||_{L^{\infty}_tL^2_x}=\sup_t||TT^*(f)(t,.)||_{L^2}=\sup_t||\int S_{t-s}f(s,.)ds||_{L^2}\leq \sup_t\int ||S_{t-s}f(s,.)||_{L^2}ds=\sup_{t}\int ||f(s,.)||_{L^2}ds=||f||_{L^1_tL^2_x} where we have used the isometricity of S_t for the norm L^2. So, this case is also proved. And thus the whole result is finished. What is more, reexamine the condition in the inequality of Hard-Littlewood-Sobolev, we have to ensure that a=2/q\in (0,n)=(0,1). This means that q>2. So, the above result is really valid for the case q>2, just like what the condition says in the theorem.

Now we want to generalize the above analysis for the case where the equation is not homogeneous, that is:

i\partial_tu+\Delta u=f(f\in C(\mathbb{R},S(\mathbb{R}^d))) with initial condition u|{t=0}u_0\in S(\mathbb{R}^d)

First we can write the explicit formula for this solution u(t)=S_tu_0-i\int_0^tS_{t-s}f(s)ds. A little digression here. We see that S_t serves as an integral kernel, which means that it takes into account the effects of stimulation. For example, if f=0. At time t, the effect of the stimulation u_0 is taken to be S_tu_0. If at first, we add some stimulation, f(0), then its effect at time t becomes S_tf(0), and if at time t_1 we add another stimulation f(t_1), then this stimulation will have an effect S_{t-t_1}f(t_1) at time t. So, the whole effect of all the stimulation f(s)(0\leq s\leq t) will be \int_0^t S_{t-s}f(s)ds. Then why don’t we take a similar integral for the stimulations u(s)(0\leq s\leq t)? The point is that S_{t-s}u(s)=S_{t-s}S_su_0=S_tu_0(\forall s), so the effect of these u(s) are uniform, that is why we consider only the starting point, that is enough. Now return to our work, our result for the inhomogeneous case is

Suppose that (q_2,r_2) satisfy 2/q_2+d/r_2=d/2 and q_2,r_2\in[2,\infty],q_2>2. Then the solution u to the equation i\partial_tu+\Delta u=f(f\in L^{q_2'},\mathbb{R},L^{r'_2}(\mathbb{R}^d)) with initial condition u|_{t=0}=0 satisfies that u\in L^{q_1}(\mathbb{R},L^{r_1}(\mathbb{R}^d)) for any (q_1,r_1) satisfying the same condition as (q_2,r_2) as above. What is more, for each such pair, there is a constant C independent of f such that ||u||_{L^{q_1}_tL^{r_1}_x}\leq C||f||_{L^{q_2'}_tL^{r'_2}_x}.

The main idea of the proof for this result is that, we first show that ||u||_{L^{q_2}_tL^{r_2}_x}\leq c||f||_{L^{q_2'}_tL^{r_2'}_x} for some constant c using the same argument as above, and then ||u||_{L^{\infty}_tL^2_x}\leq c'||f||_{L^{q_2'}_tL^{r_2'}_x} for another constant c' using the continuity of the operator T^*:B'\rightarrow H'. So, we can apply the interpolation theorems to obtain that for any q_1\in [q_2,\infty], ||u||_{L^{q_1}_tL^{r_1}_x}\leq c''||f||_{L^{q_2'}_tL^{r_2'}_x} for some constant c'' depending on q_1. As for the case q_1\in[2,q_2], we will prove it later.

So, all this way, we have shown the global property of the solution to the linear Schrodinger equation.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s