Short Time Quantum Propagator and Bohmian Trajectories
Total Page:16
File Type:pdf, Size:1020Kb
Short Time Quantum Propagator and Bohmian Trajectories Maurice de Gosson∗ Universit¨at Wien, Fakult¨at f¨ur Mathematik, NuHAG Wien 1090 Basil Hiley University of London Birkbeck College, Theoretical Physics Unit London WC1E 7HX September 11, 2018 Abstract We begin by giving correct expressions for the short-time action; following the work of one of us and Makri–Miller. We use these es- timates to derive a correct expression modulo ∆t2 for the quantum propagator and we show that the quantum potential is negligible mod- ulo ∆t2 for a point source. We finally prove that this implies that the quantum motion is classical for very short times. Keywords: generating function, short-time action, quantum potential, Bohmian trajectories arXiv:1304.4771v1 [quant-ph] 17 Apr 2013 1 Introduction Explicit approximate expressions for the short-time action play an essential role in various aspects of quantum mechanics (for instance the Feynman path integral, or semi-classical mechanics), and so does the associated Van Vleck determinant. Unfortunately, as already observed by Makri and Miller [15, 16], the literature seems to be dominated by formulas which are wrong ∗Financed by the Austrian Research Agency FWF (Projektnummer P20442-N13). 1 even to the first order of approximation! Using the calculations of Makri and Miller, which were independently derived by one of us in [3] using a slightly different method, we show that the correct approximations lead to precise estimates for the short-time Bohmian quantum trajectories for an initially sharply located particle. We will see that these trajectories are classical to the second order in time, due to the vanishing of the quantum potential for small time intervals. In this paper we sidestep the philosophical and ontological debate around the “reality” of Bohm’s trajectories and rather focus on the mathematical issues. 2 Bohmian trajectories Consider a time-dependent Hamiltonian function n p2 H(x,p,t)= j + U(x,t) (1) =1 2mj Xj and the corresponding quantum operator n ~2 2 ~ ∂ H(x, i x,t)= − 2 + U(x,t). (2) − ∇ =1 2mj ∂xj Xj b The associated Schr¨odinger equation is ∂Ψ i~ = H(x, i~ ,t)Ψ , Ψ(x, 0)=Ψ0(x). (3) ∂t − ∇x Let us write Ψ in polarb form ReiΦ/~; here R = R(x,t) 0 and Φ= S(x,t) ~ ≥ are real functions. On inserting ReiS/ into Schr¨odinger’s equation and separating real and imaginary parts, one sees that the functions R and S satisfy, at the points (x,t) where R(x,t) > 0, the coupled system of non- linear partial differential equations n 2 n ∂S 1 ∂S ~2 ∂2R + + U(x,t) 2 = 0 (4) ∂t =1 2mj ∂xj − =1 2mjR ∂xj Xj Xj 2 n ∂R 1 ∂ 2 ∂Φ + R = 0. (5) ∂t =1 mj ∂xj ∂xj Xj 2 The crucial step now consists in recognizing the first equation as a Hamilton– Jacobi equation, and the second as a continuity equation. In fact, introduc- ing the quantum potential n 2 2 Ψ ~ ∂ R Q = 2 (6) − =1 2mjR ∂xj Xj (Bohm and Hiley [1]) and the velocity field Ψ 1 ∂Φ 1 ∂Φ v (x,t)= , ..., (7) m1 ∂x1 m ∂x n n the equations (4) and (5) become ∂Φ Ψ + H(x, Φ,t)+ Q (x,t) = 0 (8) ∂t ∇x ∂ρ Ψ 2 + div(ρv )=0 , ρ = R (9) ∂t The main postulate of the Bohmian theory of motion is that particles follow quantum trajectories, and that these trajectories are the solutions of the differential equations Ψ ~ 1 ∂Ψ x˙ j = Im . (10) mj Ψ ∂xj The phase space interpretation is that the Bohmian trajectories are de- termined by the equations Ψ Ψ 1 Ψ Ψ ∂U Ψ ∂Q Ψ x˙ j = pj ,p ˙j = (x ,t) (x ,t). (11) mj −∂xj − ∂xj It is immediate to check that these are just Hamilton’s equations for the Hamiltonian function n 2 Ψ p Ψ H (x,p,t)= j + U(x,t)+ Q (x,t) (12) =1 2mj Xj which can be viewed as a perturbation of the original Hamiltonian H by the quantum potential QΨ (see Holland [13, 14] for a detailed study of quantum trajectories in the context of Hamiltonian mechanics). The Bohmian equations of motion are a priori only defined when R = 0 6 (that is, outside the nodes of the wavefunction); this will be the case in our constructions since for sufficiently small times this condition will be satisfied by continuity if we assume that it is case at the initial time. An important feature of the quantum trajectories is that they cannot cross; thus there will be no conjugate points like those that complicate the usual Hamiltonian dynamics. 3 3 The Short-Time Propagator The solution Ψ of Schr¨odinger’s equation (3) can be written Ψ(x,t)= K(x,x0; t)Ψ0(x0)dx0 Z where the kernel K is the “quantum propagator”: K(x,x0; t)= x exp( iHt/~) x0 . h | − | i Schr¨odinger’s equation (3) is then equivalentb to ∂K i~ = H(x, i~ ,t)K , K(x,x0;0) = δ(x x0) (13) ∂t − ∇x − where δ is the Diracb distribution. Physically this equation describes an isotropic source of point-like particle emanating from the point x0 at initial time t0 = 0. We want to find an asymptotic formula for K for short time intervals ∆t. Referring to the usual literature, such approximations are given by expressions of the type 2 1 n/ i K(x,x0;∆t)= ρ(x,x0;∆t)exp S(x,x0;∆t) 2πi~ ~ p where S(x,x0;∆t) is the action along the classical trajectory from x0 to x in time ∆t and 2 ∂ S(x,x0;∆t) ρ(x,x0;∆t) = det − ∂x ∂x j k 1≤j,k≤n is the corresponding Van Vleck determinant. It is then common practice (especially in the Feynman path integral literature) to use the following “midpoint approximation” for the generating function S: n mj 2 1 S(x,x0;∆t) (xj x0) (U(x,t0)+ U(x0,t0))∆t (14) ≈ =1 2∆t − − 2 Xj or, worse, n mj 2 1 S(x,x0; t,t0) (xj x0) (U( 2 (x + x0),t0)∆t (15) ≈ =1 2∆t − − Xj However, as already pointed out by Makri and Miller [15, 16], these “approx- imations” are wrong; they fail to be correct even to first orderin∆t! In fact, 4 Makri and Miller, and one of us [3] have shown independently, and using different methods, that the correct asymptotic expression for the generating function is given by n mj 2 2 S(x,x0;∆t)= (xj x0) U(x,x0)∆t + O(∆t ) (16) =1 2∆t − − Xj e where U(x,x0, 0) is the average value of the potential over the straight line joining x0 at time t0 to x at time t with constant velocity: e 1 U(x,x0)= U(λx + (1 λ)x0, 0)dλ. (17) 0 − Z For instance whene 1 2 2 2 2 H(x,p)= (p + m ω x ) 2m is the one-dimensional harmonic oscillator formula (16) yields the correct expansion 2 m 2 mω 2 2 S(x,x0; t)= (x x0) (x + xx0 + x0)∆t + O(∆t ); (18) 2∆t − − 6 the latter can of course be deduced directly from the exact value mω 2 S(x,x0; t,t0)= ((x + x0) cos ω∆t 2xx0) (19) 2sin ω∆t − by expanding sin ω∆t and cos ω∆t for ∆t 0. This correct expression is → of course totally different from the erroneous approximations obtained by using the “rules” (14) or (15). Introducing the following notation, n 2 (xj x0) S(x,x0;∆t)= mj − U(x,x0)∆t, (20) =1 2∆t − Xj e e leads us to the Makri and Miller approximation (formula (17c) in [15]) for the short-time propagator: n/2 1 i 2 K(x,x0;∆t)= ρ(x,x0;∆t)exp S(x,x0;∆t) + O(∆t ) 2πi~ ~ p (21) e 5 where 2 ∂ S(x,x0;∆t) ρ(x,x0;∆t) = det . ∂x ∂x − j k !1 e ≤j,k≤n It turns out that this formula can be somewhat improved. The Van Vleck determinant ρ(x,x0;∆t) is explicitly given, taking formula (20) into account, by 1 ′′ ρ(x,x0;∆t) = det M U (x,x0)∆t −∆t − x,x0 where M is the mass matrix (the diagonal matrixe with positive entries the masses mj) and 2 ∂ U(x,x0) U ′′ = . x,x0 ∂x ∂x − j k !1 e ≤j,k≤n Writing e 1 ′′ 1 −1 ′′ 2 M U (x,x0)∆t = M[I × M U (x,x0)∆t ] −∆t − x,x0 −∆t n n − x,x0 1 2 e = M[I × + O(∆te)], −∆t n n we have by taking the determinant of both sides m1 mn 2 ρ(x,x0;∆t)= · · · det I × + O(∆t ) . (∆t)n n n 2 2 Noting that det In×n + O(∆t ) =1+ O(∆t ), we thus have m1 mn 2 ρ(x,x0;∆t)= · · · (1 + O(∆t )). (22) (∆t)n Writing m1 m ρ(∆t)= · · · n (23) (∆t)n which is just the Van Vleck density for the free particle Hamiltonian. We e thus have 2 ρ(x,x0;∆t)= ρ(∆t)+ O(∆t ) (24) and hence we can rewrite formula (21) as e n/2 1 i 2 K(x,x0;∆t)= ρ(∆t)exp S(x,x0;∆t) + O(∆t ). (25) 2πi~ ~ p We will see below that this formulae allows ane easy study of the quantum potential for K.