Complete set of stationary solutions

[MHD93]

In the infinite square well, the stationary states solutions form a complete set, and therefore I can write a function such as ( f(x) = -x^2 + x ) as an infinite sum of them,

But this function is, clearly, not a solution of S.E., although it's written as a sum of solutions.

Why is it not a solution? is it because the sum is infinite or what?

 

[Chopin]

The eigenfunctions of the infinite square well form a complete basis only for functions which have compact support within the range of the well. In other words, you can only make functions which are zero outside of the well (since all of the eigenfunctions will be zero outside of it, there's no way you can add them up to get a nonzero number). Since the function you gave is nonzero outside of the well, it can't be represented as the sum of square well eigenfunctions.

 

[rubi]

The eigenfunctions of the infinite square well form a complete basis only for functions which have compact support within the range of the well. In other words, you can only make functions which are zero outside of the well (since all of the eigenfunctions will be zero outside of it, there's no way you can add them up to get a nonzero number). Since the function you gave is nonzero outside of the well, it can't be represented as the sum of square well eigenfunctions.

The function f:[0,1]->\mathbb R, f(x) = -x^2+x can definitely be expanded in square well eigenfunctions, but it still isn't a solution to the Schrödinger equation.

The answer to the question is that only because two functions are eigenfunctions of some operator, their sum doesn't need to be an eigenfunction anymore (unless they lie in the same eigenspace). I suggest you try it yourself with a simple matrix like M = \mathrm{diag} (1,2). However, you can use the function f as initial condition at t=0 for the time-dependent Schrödinger equation and get a completely valid solution.

 

[Chopin]

Whoops, you're right, I misread what the question was asking.

Sounds like the issue here is confusion between the time-independent and time-dependent Schrödinger equations. The former is not really a differential equation so much as a family of differential equations: (-\frac{\hbar^2\nabla^2}{2m} + V(x))\psi(x) = E\psi(x), for various values of E. The square well solutions form a complete basis for describing functions over that range, but each function has a different eigenvalue, so adding together functions with different eigenvalues does not produce another solution to the TISE.

On the other hand, the time-DEPENDENT equation is i\hbar\frac{\partial}{\partial t}\psi(x,t)=(-\frac{\hbar^2\nabla^2}{2m} + V(x))\psi(x,t). This is a single linear differential equation, not a family of equations, so you can add together any number of solutions to this, and the result will still be a solution. So you can take your parabola and decompose it into solutions of the TISE (this is basically a Fourier transform), and then add those all up and use them in the TDSE to determine the system's time evolution.

 

[tom.stoer]

Why is it not a solution??

It's a simple linear algebra problem. Suppose you have an eigenvalue equation with two solutions (eigenvectors and -values; i=1,2)

(H - E_i) u_i = 0

What about

v = a₁u₁ + a₂u₂

We apply H

Hv = H (a₁u₁ + ...) = a₁E₁u₁ + ...

but for different E₁ ≠ E₂ this is no longer an eigenvector of H

 

[stevendaryl]

In the infinite square well, the stationary states solutions form a complete set, and therefore I can write a function such as ( f(x) = -x^2 + x ) as an infinite sum of them,

But this function is, clearly, not a solution of S.E., although it's written as a sum of solutions.

Why is it not a solution? is it because the sum is infinite or what?

The way it works is this: If \psi_n(x,t) is a stationary state, that doesn't really mean that it is time-independent, it means that the time-dependency is very simple:

\psi_n(x,t) = \psi_n(x,0) e^{- \omega_n t}

where \omega_n is the characteristic frequency of state \psi_n.

Now, if you express an arbitrary function as a combination of stationary states, you can do it at a single moment (say t=0):

f(x) = \sum C_n \psi_n(x,0)

But this f(x) is not a stationary state, because if you initially choose a wave function \Psi with \Psi(x,0) = f(x), and you wait a while, the wave function will change to:

\Psi(x,t) = \sum C_n \psi_n(x,0)e^{- \omega_n t}

That's not stationary, except in the special case where all the \omega_n have the same value.