Propagator and delta functions

jostpuur
Messages
2,112
Reaction score
19
http://en.wikipedia.org/wiki/Propagator

What does this equation mean:

<br /> \Big(H_x - i\hbar\partial_t\Big) K(x,t,x&#039;,t&#039;) = -i\hbar\delta(x-x&#039;)\delta(t-t&#039;)<br />

Wouldn't it be more relevant to emphasize these equations:

<br /> \Big(H_x - i\hbar\partial_t\Big) K(x,t,x&#039;,t&#039;) = 0,\quad\quad t\neq t&#039;<br />
<br /> K(x,t,x&#039;,t&#039;) = \delta(x-x&#039;),\quad\quad t=t&#039;<br />
 
Physics news on Phys.org
jostpuur said:
http://en.wikipedia.org/wiki/Propagator

What does this equation mean:

<br /> \Big(H_x - i\hbar\partial_t\Big) K(x,t,x&#039;,t&#039;) = -i\hbar\delta(x-x&#039;)\delta(t-t&#039;)<br />

Wouldn't it be more relevant to emphasize these equations:

<br /> \Big(H_x - i\hbar\partial_t\Big) K(x,t,x&#039;,t&#039;) = 0,\quad\quad t\neq t&#039;<br />
<br /> K(x,t,x&#039;,t&#039;) = \delta(x-x&#039;),\quad\quad t=t&#039;<br />
<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<,,

The second formulation, apart from the missing h bar, and the first are really the same. But, the latter form was used prior to Dirac to define Green's functions, and the techniques for solution are somewhat tedious -- see any E&M text written before 1925, and "elementary" discussions of Green's functions.

Four dimensional Green's functions, that is a source at x' and t', are basic to solutions of Maxwell's equations. Once QM arrived, physicists were comfortable with inverse operators, and your first fomulation became the approach of choice.

If I recall correctly, Feynman discussed this issue early on in his first QED papers; Schwinger, likewise. Schwinger, however is very difficult to read.

Regards,
Reilly
 
Once the propagator and the time evolution are defined:

<br /> K(x,t,x&#039;,t&#039;) = \frac{1}{\sqrt{2\pi\hbar |t-t&#039;|}}\exp\Big(\frac{i(x-x&#039;)^2}{2\hbar(t-t&#039;)}\Big)<br />

<br /> \psi(x,t) = \int\limits_{-\infty}^{\infty}K(x,t,x&#039;,t&#039;)\psi(x&#039;,t&#039;)dx&#039;<br />

we are usually interested in two things. One is

<br /> i\hbar\partial_t\psi(x,t) = -\frac{\hbar^2}{2}\partial_x^2\psi(x,t)<br />

so that the Schrödinger equation is satisfied, and the other one is

<br /> \lim_{t\to t&#039;} \psi(x,t) = \psi(x,t&#039;)<br />

so that the correct initial value is satisfied. In order to prove that the Schrödinger equation is satisfied, we need this:

<br /> \Big(H_x - i\hbar\partial_t\Big) K(x,t,x&#039;,t&#039;) = 0,\quad\quad t\neq t&#039;<br />

And in order to prove that the initial value is satisfied, we need this:

<br /> K(x,t,x&#039;,t&#039;) = \delta(x-x&#039;),\quad\quad t=t&#039;<br />

How could this equation:

<br /> \Big(H_x - i\hbar\partial_t\Big) K(x,t,x&#039;,t&#039;) = -i\hbar\delta(x-x&#039;)\delta(t-t&#039;)<br />

replace the previous two? I see that if you set \delta(t-t&#039;)=0 for t\neq t&#039;, then you get the Schrödinger equation satisfied, but how about the initial value? What is the delta function with respect to time supposed to do there?

Not only do I not understand how to use this delta function with respect to time, I also don't understand where this equation is supposed to come from.

When t\neq t&#039;, the following can be calculated:

<br /> i\hbar\partial_t K(x,t,x&#039;,t&#039;) \;=\; -\frac{1}{2}\sqrt{\frac{\hbar}{2\pi}} \frac{1}{|t-t&#039;|^{5/2}}<br /> \Big(i(t-t&#039;) \;-\; \frac{(x-x&#039;)^2}{\hbar}\Big) \exp\Big(\frac{i(x-x&#039;)^2}{2\hbar (t-t&#039;)}\Big)<br /> \;=\; -\frac{\hbar^2}{2}\partial_x^2 K(x,t,x&#039;,t&#039;)<br />

You cannot get the delta function with respect to time with a limit t\to t&#039; at least, because there is no divergent quantity there. It's pure zero:

<br /> \lim_{t\to t&#039;} \Big(-\frac{\hbar^2}{2}\partial_x^2 - i\hbar\partial_t\Big) K(x,t,x&#039;,t&#039;) = \lim_{t\to t&#039;} 0 = 0<br />
 
I see (roughly) what the identity means now. A following calculation can be carried out:

<br /> \int\limits_{-\infty}^{\infty} dx \int\limits_{-\infty}^{\infty} dt \Big(\big(-\frac{\hbar^2}{2}\partial_x^2<br /> \;+\;i\hbar\partial\big)\psi(x,t)\Big) K(x,t,x&#039;,t&#039;)<br /> \;=\; \lim_{\epsilon\to 0^+} \int\limits_{-\infty}^{\infty} dx<br /> \int\limits_{]-\infty,t&#039;-\epsilon]\cup [t&#039;+\epsilon,\infty[} dt<br /> \Big(\psi(x,t)\big(-\frac{\hbar^2}{2}\partial_x^2\big)K(x,t,x&#039;,t&#039;)<br /> \;+\; \big(i\hbar\partial_t\psi(x,t)\big) K(x,t,x&#039;,t&#039;)\Big)<br />

<br /> = \lim_{\epsilon\to 0^+}<br /> i\hbar \int\limits_{-\infty}^{\infty} dx\Big(\psi(x,t&#039;-\epsilon)K(x,t&#039;-\epsilon,x&#039;,t&#039;) \;-\; \psi(x,t&#039;+\epsilon)<br /> K(x,t&#039;+\epsilon,x&#039;,t&#039;)\Big)<br />

<br /> = \lim_{\epsilon\to 0^+}<br /> \frac{i\hbar}{\sqrt{2\pi\hbar\epsilon}} \int\limits_{-\infty}^{\infty} dx\big(\psi(x,t&#039;) + o(\epsilon)\big)<br /> (-2i)\sin\Big(\frac{(x-x&#039;)^2}{2\hbar\epsilon}\Big)<br /> \;=\; \frac{(i\hbar)(-2i)}{\sqrt{\pi}} \psi(x&#039;,t&#039;) \int\limits_{-\infty}^{\infty} \sin(y^2)dy \;=\; \hbar\psi(x&#039;,t&#039;)<br />

Which would be the same thing as

<br /> \Big(-\frac{\hbar^2}{2}\partial_x^2 - i\hbar\partial_t\Big)K(x,t,x&#039;,t&#039;) = \hbar\delta(x-x&#039;)\delta(t-t&#039;)<br />

I see the coefficient -i vanished somewhere. Maybe there's a mistake somewhere...

Anyway, I still don't see what's the point of this. In what kind of situation are we interested in integrating with respect to time like this? In order to get the time evolution, there's no need to integrate anything with respect to time.
 
jostpuur said:
<br /> \psi(x,t) = \int\limits_{-\infty}^{\infty}K(x,t,x&#039;,t&#039;)\psi(x&#039;,t&#039;)dx&#039;<br />
You've got a free t' parameter on the rhs, but not on the lhs.
The point of this sort of thing is that, given a \psi(x&#039;,t&#039;) you can
find \psi(x,t). To do that, a t' integration is needed.

K(x,t,x',t') should be considered as an inverse operator to
<br /> \left(H_x - i \hbar \partial_t \right)<br />
(in the distributional sense).
 
strangerep said:
jostpuur said:
<br /> \psi(x,t) = \int\limits_{-\infty}^{\infty}K(x,t,x&#039;,t&#039;)\psi(x&#039;,t&#039;)dx&#039;<br />
You've got a free t' parameter on the rhs, but not on the lhs.

The initial time t&#039; can be chosen arbitrarily. Or if you know the wave function only at some given instant, then you should use that given instant.

The propagator has the following transitivity property

<br /> \int\limits_{-\infty}^{\infty} K(x,t,x&#039;,t&#039;) K(x&#039;,t&#039;,x&#039;&#039;,t&#039;&#039;)dx&#039; = K(x,t,x&#039;&#039;,t&#039;&#039;)<br />

and it doesn't matter what t&#039; you use in the middle. It follows from this property, that the time-evolution equation gives the same final wave function \psi(x,t) for all initial times t&#039;. The time-evolution equation, in which integration is carried out with respect to x&#039;, and not with respect to a fixed variable t&#039;, is correct after all, so no integration with respect to t&#039; is needed.
 
Take the FT of your first eq, and get K(E,p) = 1/(E-p*p) -- I'm dropping constants, except those that are really important.

Now get K(x,t) by inverting. So, K(x,t) == Integral dE dp{exp(-iEt +ipx)/(E-p*p)}. Do the E integration first, to handle the pole, push it below the E axis -- add -i(epsilon) to E, and close the contour in the lower half plane, Then

K = Integral dp exp{ -ip*pt +ipx } So, as t->0, K does turn into a delta function.
(This type of contour integration is common in scattering theory, and in EE as well.)

Regards,
Reilly Atkinson
 
My final question is that what is this identity good for.

It seems like an interesting result, but I cannot see any use for it.
 
jostpuur said:
The propagator has the following transitivity property

<br /> \int\limits_{-\infty}^{\infty} K(x,t,x&#039;,t&#039;) K(x&#039;,t&#039;,x&#039;&#039;,t&#039;&#039;)dx&#039; = K(x,t,x&#039;&#039;,t&#039;&#039;)<br />
And what ever happened to integral over t'?

If you are talking about classical limit, it should actually look like this:

<br /> \int\limits_{-\infty}^{\infty} dx&#039; \int\limits_{t}^{t&#039;&#039;} dt&#039; K(x,t,x&#039;,t&#039;) K(x&#039;,t&#039;,x&#039;&#039;,t&#039;&#039;) = K(x,t,x&#039;&#039;,t&#039;&#039;)<br />
 
  • #10
K^2 said:
jostpuur said:
The propagator has the following transitivity property

<br /> \int\limits_{-\infty}^{\infty} K(x,t,x&#039;,t&#039;) K(x&#039;,t&#039;,x&#039;&#039;,t&#039;&#039;)dx&#039; = K(x,t,x&#039;&#039;,t&#039;&#039;)<br />

And what ever happened to integral over t'?

If you are talking about classical limit, it should actually look like this:

<br /> \int\limits_{-\infty}^{\infty} dx&#039; \int\limits_{t}^{t&#039;&#039;} dt&#039; K(x,t,x&#039;,t&#039;) K(x&#039;,t&#039;,x&#039;&#039;,t&#039;&#039;) = K(x,t,x&#039;&#039;,t&#039;&#039;)<br />

The equation that I wrote is correct, and it has fixed t&#039;. It is a property of propagators which appear in QM and also in the study of the heat equation.
 
  • #11
The propagator has the following transitivity property

<br /> \int\limits_{-\infty}^{\infty} K(x,t,x&#039;,t&#039;) K(x&#039;,t&#039;,x&#039;&#039;,t&#039;&#039;)dx&#039; = K(x,t,x&#039;&#039;,t&#039;&#039;)<br />

The meaning and importance of this identity become clear if you realize that the propagator (Green function) K(x,t : y,t_{0})can be interpreted as transition amplitude. That is the probability amplitude for the system prepared at t_{0} with position eigenvalue y to be found at a later time t at x ;

K(x,t : y,t_{0}) = \langle x | e^{-iH(t-t_{0})}| y \rangle = \langle x,t | y,t_{0}\rangle \ \ \ (1)

where \langle x,t| and |y,t_{0}\rangle are the eigenket and the eigenbra of the position operator in the Hiesenberg picture. Because at any given time these eigen vectors form a complete set, we can insert the identity operator;

\int dz |z,t_{1}\rangle \langle z,t_{1}| = \hat{1}

at any place we want. So, by dividing the evolution interval (t_{0},t_{2}) into two parts (t_{0},t_{1}) and (t_{1},t_{2}), we can decompose the transition amplitude as :

\langle x , t_{2} | y , t_{0} \rangle = \int dz \langle x ,t_{2} | z , t_{1} \rangle \langle z , t_{1}| y , t_{0}\rangle \ \ (2)

(t_{2} &gt; t_{1} &gt; t_{0})

This composition property of the transition amplitude is an important consistency requirement underlying the whole formalism of path integration. Indeed, we can use it to show that the transition amplitude \langle x , t| y , t_{0} \rangle satisfies Schrodinger equation in the variables (x,t), just as the propagater K(x,t : y,t_{0}) ;

i\partial_{t}\langle x , t | y , t_{0} \rangle = H(t_{0})\langle x ,t | y ,t_{0} \rangle \ \ \ (3)

( argue that you can write \langle x , t + \epsilon | y , t_{0} \rangle = \delta (x - y) - i \epsilon\delta(x-y)H(t_{0})
and use eq(2) for t_{1} = t and t_{2} = t + \epsilon)

When he path integration method is carried over to Brownian motion, eq(2) is known as the Chapman-Kolmogorov equation, and in diffusion theory, the Smoluchowsky equation.

regards

sam
 
Back
Top