Feynmans Herleitung der Schrödingergleichung

Text auszugsweise aus

 Feynman and Hibbs Quantum Mechanics and Path Integrals,
McGraw-Hill Publishing Company 1965.  4. Kapitel  

Formeln (blau): Die "Rechnung" wurde mit Maple symbolisch durchgeführt.

 

 

4 The Schrödinger Description of Quantum Mechanics

 

The path integrals which we have discussed so far have integrands which are exponentials of actions with the property (4-1) 

S[2,1] = S[2,3]+S[3,1]

Such path integrals can be analyzed in terms of the properties of integral equations which can be deduced from them. [...]

A still more convenient method is to reduce the path integrals to differential equations if possible. This possibility exists in quantum mechanics and is, in fact, the most convenient way to present that theory. It is in almost every case easier to solve the differential equation than it is to evaluate the path integral. The conventional presentation of quantum mechanics is based on this equation, called the Schrödinger equation. Here we shall derive this equation from our formulation. [...]

In this chapter our purpose is twofold: (1) For the reader primarily interested in quantum mechanics our aim is to connect the path integral formulation with other formulations which are found in the standard literature [...]. (2) For the reader primarily interested in path integrals this chapter will show a technique which is available for a certain class of path integrals to reduce these path integrals to differential equations. This technique is best shown by the particular example of quantum mechanics which we shall develop here.

 

THE SCHRÖDINGER EQUATION

 

The Differential Equation Form. The reason that we can develop a differential equation is that the relationship of Eq. (4-1) is correct for any values of the points 1, 2 and 3. For example, the time can be only an infinitesimal time epsilon greater than the time t[3]. This will permit us to relate the value of a path integral for one time to its value a short time later. In this manner we can obtain a differential equation for the path integral.

We have already found that as a consequence of Eq. (4-1) we can define a wave function. Furthermore, we know that the equation (4-2)   

psi(x[2],t[2]) = Int(K(x[2],t[2],x[1],t[1])*psi(x[1],t[1]),x[1] = -infinity .. infinity)

gives the wave function at a time t[2] in terms of the wave function at a time t[1]. In order to obtain the differential equation we seek, we apply this relationship in the special case that the time t[2] differs only by an infinitesimal interval epsilon from t[1]. The kernel K(2,1) is proportional to the exponential of i/hq  (hq = h/(2*Pi)) times the action for the interval  t[1] to t[2]. For a short interval epsilon the action is approximately epsilon times the lagragian for this interval. That is, using the same approximation as that of Eq. (2-34)

K(i+1,i) = `1/A`*exp(i*epsilon/hq*L((x[i+1]-x[i])/epsilon,1/2*x[i+1]+1/2*x[i],1/2*t[i+1]+1/2*t[i]))

we have (4-3)       

psi(x,t+epsilon) = Int(`1/A`*exp(epsilon*i/hq*L((x-y)/epsilon,1/2*x+1/2*y))*psi(y,t),y = -infinity .. infinity)

We shall now apply this  to the special case of a particle moving in a potential V(x,t) in one dimension, i.e., that for which

L = 1/2*m*diff(x(t),t)^2-V(x,t)

In this case Eq. (4-3) becomes (4-4)

psi(x,t+epsilon) = Int(`1/A`*exp(1/2*i*m*(x-y)^2/epsilon/hq)*exp(-i/hq*epsilon*V(1/2*x+1/2*y,t))*psi(y,t),y = -infinity .. infinity)

The quantity     

(x-y)^2/epsilon

appears in the exponent of the first factor. It is clear that if y is appreciably differentfrom x, this quantity is very large and the exponential consequently oscillates very rapidly as y varies. When this factor oscillates rapidly, the integral over y gives a very small value (because of the smooth behavior of the other factors). Only if y is near x (where the exponential changes more slowly) do we get important contributions. For this reason we make the substitution

y = x+eta

with the expectation that appreciable contributions to the integral will occur only for small   eta. We obtain  (4-5)

psi(x,t+epsilon) = Int(`1/A`*exp(1/2*i*m*eta^2/epsilon/hq)*exp(-i/hq*epsilon*V(x+1/2*eta,t))*psi(x+eta,t),eta = -infinity .. infinity)

The phase of the first exponential changes by the order of 1 radian when eta is of the order  

(epsilon*hq/m)^(1/2)

so that most of the integral is contributed by values of eta in this order.

 

We may expand Psi in a power series. We need only keep terms of order epsilon. This implies keeping second-order terms in eta. The term         

epsilon*V(x+1/2*eta,t)

may be replaced by

epsilon*V(x,t)

because the error is of higher order than epsilon. Expanding the left-hand side to first order in epsilon  and the right-hand side to  first order in epsilon and second order in eta, we obtain (4-6)     

psi(x,t)+epsilon*diff(psi(x,t),t) = Int(`1/A`*exp(1/2*i*m*eta^2/epsilon/hq)*(1-i*epsilon/hq*V(x,t))*(psi(x,t)+eta*diff(psi(x,t),x)+1/2*eta^2*diff(psi(x,t),`$`(x,2))),eta = -infinity .. infinity)

If we take the leading term on the right-hand side, we have the quantity

psi(x,t)

multiplied by the integral (4-7)

`1/A`*Int(exp(1/2*i*m*eta^2/h/epsilon),eta = -infinity .. infinity) = `1/A`*2^(1/2)*(Pi*i*hq*epsilon/m)^(1/2)

On the left-hand side we just have psi(x,t). In order  that both sides agree in the limit epsilon approaches 0, it is necessary that A  be so chosen that the expression of Eq. (4-7) equals 1. That is (4-8),

A = 2^(1/2)*(1/m*h*epsilon*Pi*i)^(1/2)

[...]. This is a way of obtaining the quantity A in more complicated problems also. The A must be so chosen that the equation is correct to zero order in epsilon. Otherwise, no limit will exist as epsilon approaches 0 in the original path integral.

In order to evaluate the right-hand side of Eq. (4-6), we shall have to use two integrals  (4-9)  

Int(exp(1/2*i*m*eta^2/h/epsilon)*eta,eta = -infinity .. infinity) = 0

and (4-10)

1/A*Int(exp(1/2*i*m*eta^2/h/epsilon)*eta^2,eta = -infinity .. infinity) = i*hq*epsilon/m

Writing out the right-hand side of Eq. (4-6) gives (4-11)

psi(x,t)+epsilon*diff(psi(x,t),t) = psi(x,t)-i*epsilon/hq*V(x,t)*psi(x,t)-1/2*hq*epsilon/i/m*diff(psi(x,t),`$`(x,2))

This will be true to order epsilon if psi satisfies the differential equation (4-12)   

-hq/i*diff(psi(x,t),t) = -1/2*hq^2/m*diff(psi(x,t),`$`(x,2))+V(x,t)*psi(x,t)

This is the Schrödinger equation for our problem of a particle moving in one dimension. Corresponding equations in more complicated situations can be worked out in the same way.

 

HOME | Fächer | Physik | Elektrizität | Optik | Atomphysik | Quantenphysik | Top