Physics 582, Problem Set 4 Solutions TAs: Hart Goldman and Ramanjit Sohal Fall 2018 1. PATH INTEGRAL FOR A PARTICLE IN A DOUBLE POTENTIAL WELL 1. In real time, q 0, T/2 q 0, T/2 = q 0 e i HT q 0 = where Dq e q(t f )=q 0 q(t 0 )= q 0 i t f t dtl(q, q) i (1.1) L = 1 2 m q2 λ(q 2 q 2 0) 2. (1.2) On Wick rotating t iτ, we obtain the imaginary time expression q 0, T/2 q 0, T/2 = Dq e 1 t f t dτl i E (q, q) (1.3) where L E = 1 2 m q2 + λ(q 2 q0) 2 2. (1.4) So, in imaginary time, the potential is V (q) = λ(q 2 q0) 2 2 which is an inverted double potential well. 2. In real time, the equations of motion are readily found to be m q = 4λq(q 2 q0) 2 (real time) (1.5) while, in imaginary time, m q = 4λq(q 2 q0) 2 (imaginary time). (1.6) Again, we see that the EOM in imaginary time describe a point particle moving in an inverted potential well. In order to find the trajectory in imaginary time, we make use of the conservation of energy, E = T + V = 1 2 m q2 λ(q 2 q0) 2 2 (1.7)
2 with the initial and final conditions q( T/2) = q 0 and q(t/2) = q 0, respectively. This expression is obtained by multiplying Eq. (1.6) by q and integrating with respect to τ. Now, we are interested in the limit T. The initial/final conditions q(±t/2 ± ) = ±q 0 can only be enforced if we also require that q(±t/2 ± ) = 0 (otherwise the particle would travel past ±q 0 and not be able to return). So, E i = E f = 1 2 m(0)2 λ((±q 0 ) 2 q0) 2 2 = 0. (1.8) Thus, in general, 0 = 1 2λ 2 m q2 λ(q 2 q0) 2 2 = q = ± m (q2 q0) 2 2. (1.9) This ODE has solutions [ ] 2λ q c (τ) = q 0 tanh q 0 m (τ τ 0). (1.10) In order to satisfy q(± ) = ±q 0, we must choose the positive root in the last equation. This solution is unique up to the choice of τ 0 as it is the solution to a second order ODE with initial/final conditions for q and q specified. In real time, this trajectory would be classically forbidden. In particular, the amplitude (in real time) we calculated corresponds to the amplitude for the particle to tunnel from the well centred at q 0 to the well at q 0. 3. Without loss of generality, we can set τ 0 = 0 in the imaginary time trajectory. We compute, S[q c ] = dτ [ ] 1 2 q2 c + λ(qc 2 q0) 2 2 = 4 3 0 q3 2mλ. (1.11) 4. We write q(τ) = q c (τ)+ξ(τ) where ξ(τ) satisfies the initial/final conditions ξ(± ) = 0 ] 2λ and q c (τ) = q 0 tanh [q 0 (τ τ m 0) is the classical solution to the EOM obtained above. Now, expanding the action to quadratic order in ξ(τ), we find 1 1 S 2 m q2 c V (q c )dτ + 2 m ξ 2 1 2 mv (q c )ξ 2 dτ (m q c + V (q c ))ξdτ. (1.12) The last integral vanishes as q c satisfies the EOM. So, we can write S[q, q] S[q c, q c ] + S eff [ξ, ξ]. (1.13)
3 Since the effective action is quadratic in ξ, we can absorb the factor of 1/ in the path integral into ξ by defining ξ = ξ. Dropping the tilde on xi to lighten notation, we have that lim q 0, T/2 q 0, T/2 = e Sc/ e T ξ(± )=0 S eff[ξ, ξ]. (1.14) On integrating by parts, we can write S eff [ξ, ξ] = 1 ] ξ [ m d2 2 dτ V (q 2 c ) ξdτ 1 2 ξâξdτ. (1.15) Using the explicit form of q c (τ), we can express the differential operator  as ]  = m d2 2λm dτ + 2 12λq2 0 tanh [q 2 0 τ 4λq0. 2 (1.16) In order to evaluate the path integral, we follow the standard procedure and expand ξ in orthonormal eigenstates, ψ n, of  with eigenvalues A n: ξ(τ) = n c nψ n (τ). Thus, Defining the integration measure as Dξ = N n ξâξdτ = 1 c 2 2 na n. (1.17) n 2π dc n we find that dc n 2π e 1 2 Anc2 n = [An ] 1/2 (1.18) and so ξ(± )=0 e S eff[ξ, ξ] = Det[Â] 1/2 (1.19) where the explicit form of  is given in Eq. (1.16). 2. PATH INTEGRAL FOR A CHARGED PARTICLE MOVING ON A PLANE IN THE PRESENCE OF A PERPENDICULAR MAGNETIC FIELD 1. The derivation of the path integral follows that in the free particle case very closely and so we will be a little terse. We have that Ĥ = 1 2m (p + e c A)2 = p2 2m + e2 2mc 2 A2 + e (p A + A p). (2.1) 2mc
4 Now, the transition amplitude is given by F (q f, t f q i, t i ) = d 2 q q f, t f q, t q, t q i, t i (2.2) = d 2 q 1... d 2 q N q f, t f q N, t N q N, t N q N 1, t N 1... q 1, t 1 q i, t i (2.3) where we have inserted resolutions of identity at intermediate times. Since the Hamiltonian is time-independent, we can approximate q j, t j q j 1, t j 1 = q f e ih(t j t j 1 )/ q N (2.4) = δ(q j q j 1 ) i t q j H q j 1 + O( t 2 ). (2.5) Now, inserting a resolution of identity in terms of momentum eigenstates, we find q j H q j 1 = d 2 p j (2π ) 2 eip j (q j q j 1 )/ (2.6) { 1 i t [p 2j + e2 2m c A(r j) 2 + e 2 c p j A(r j 1 ) + e ] c p j A(r j ) (2.7) We now adopt the mid-point rule and replace r j, r j 1 (r j + r j 1 )/2 (note that in the derivation of the path integral for a free particle, terms depending on r j 1 did not appear). Plugging this back into Eq. (2.5), re-exponentiating, and then plugging the result back into Eq. (2.3), we find F (q f, t f q i, t i ) (2.8) N N+1 { = lim d 2 d 2 N+1 p j q j N (2π ) exp i t [ p 2 j (q ( j q j 1 ) H p j, q )] } j + q j 1. t 2 j=1 j=1 j=1 (2.9) } + O( t 2 ). In the limit N, we find F (q f, t f q i, t i ) = DpDqe i t f t i dt[p q H(p,q)] (2.10) where DpDq = lim N N j=1 d 2 q j d 2 p j (2π ) 2. (2.11)
5 In order to integrate out p (which will amount to effecting a Legendre transform), we must discretize the path integral again. We have that d 2 p ( e p2 q A(q)) ei[p ( mc 2m e2 2mc (2π ) 2 2 A2 ] t m ) [ = e i t 1 2 m( q e mc A)2 e2 2mc ]. 2 A2 (2.12) 2πi t Thus, where F (q f, t f q i, t i ) = q(t f )=q f q(t i )=q i i t f e t L(q, q)dt i (2.13) as expected. L(q, q) = 1 2 ( q e mc A)2 e2 2mc 2 A2 = 1 2 m q2 e q A (2.14) c 2. As derived in the previous part, S = tf In the ultra-quantum limit (m 0), this becomes S tf t i 1 t i 2 mṙ2 e A ṙdt. (2.15) c e c A dr dt dt = e c Γ A dr (2.16) where Γ is the trajectory of the particle which begins and ends at r 0. Using Stokes theorem, we obtain S = e c Ω A ds = e c Ω B ds = e c Φ (2.17) where Ω is the region enclosed by Γ and Φ is the magnetic flux through this region. 3. There is no ambiguity involved in evaluating the above formula in the path integral for a particle moving in a uniform magnetic field on the plane of size L L. However, if we take L, then the plane plus the point at infinity is topologically equivalent to the two-sphere. Now consider a closed trajectory on the sphere. Using the result of the previous part, the amplitude of this trajectory is given by the Aharanov-Bohm flux it encloses. However, on the sphere, a closed trajectory bounds two regions and so there is an ambiguity as to what is considered to be inside. In particular, if we denote the surface area of the sphere as A and that of the smaller of the two regions bounded by the trajectory as a, we have that A ds = Ba OR B(A a). (2.18) Ω
6 The extra minus for the second possibility comes from the orientation of the trajectory. If we want the path integral to be well defined, we must require that the actions corresponding to these two regions differ by a multiple of 2π so that e is 1/ = e is 2/ +2πim = e is 2/. So, we must have e ( c B(A a) e ) Ba = 2πm = BA = hc c e m Φ 0m. (2.19) Hence the total flux must be quantized to an integer multiple of the fundamental flux quantum Φ 0. 3. PATH INTEGRALS FOR A SCALAR FIELD 1. The vacuum persistence amplitude in Minkowski space is J 0 0 J = N DφDφ e i d 4 x[ µφ µ φ m 2 φ φ+j φ+jφ ]. (3.1) We obtain the Euclidean space expression by Wick rotating the action (sending t ix D ) so that is S E = µ φ µ φ + m 2 φφ d 4 x. (3.2) So, J 0 0 J = N DφDφ e 1 d 4 x[ µφ µφ+m 2 φφ J φ Jφ ]. (3.3) 2. We will first do the calculation in Euclidean space and then Wick rotate back to Minkowski space. We first define a shifted field ξ as φ(x) = φ c (x) + ξ(x) = φ (x) = φ c(x) + ξ (x) (3.4) where φ c is a fixed field configuration (which we fix below). On integrating by parts, we can write the Lagrangian in terms of these fields: L E = φ [ 2 + m 2 ]φ J φ Jφ (3.5) = φ c[ 2 + m 2 ]φ c J φ c Jφ c + ξ [ 2 + m 2 ]ξ + ξ [ 2 + m 2 ]φ c (3.6) + ξ[ 2 + m 2 ]φ c J ξ Jξ. (3.7)
7 In order to make the cross terms vanish, we fix φ c to satisfy [ 2 + m 2 ]φ c = J = [ 2 + m 2 ]φ c = J (3.8) and so φ c (x) = d 4 x G E (x x )J(x ), φ c(x) = d 4 x G E (x x )J (x ) (3.9) where G E (x x ) is the Euclidean Green function. Hence the Euclidean action becomes S E = d 4 x(φ c[ 2 + m 2 ]φ c J φ c Jφ c) + d 4 x(ξ [ 2 + m 2 ]ξ) (3.10) = d 4 xd 4 x [J (x)g E (x x )J(x )] + d 4 x(ξ [ 2 + m 2 ]ξ) (3.11) where we have used the explicit forms of φ c and φ c in terms of the Green function and we have used the fact that the Green function is symmetric: G E (x x ) = G E (x x). Hence, ( Z E [J] = Z E [0] exp ) d 4 xd 4 x [J (x)g E (x x )J(x )] (3.12) where Z E [0] DξDξ e d 4 xξ [ 2 +m 2 ]ξ. (3.13) In order to evaluate Z E [0], we expand ξ in orthonormal eigenfunctions, ψ n, of  2 + m 2 with eigenvalues A n : ξ(x) = n c n ψ n (x) = ξ (x) = n c nψ n(x). (3.14) We define the measure as DξDξ = N n dc n dc n 2π 2π = N n dre(c n )dim(c n ) 2π 2πi (3.15) where Re(x) and Im(x) denote the real and imaginary parts of x and N is a normalization constant which we will use to absorb unimportant constant factors. Evaluating the integrals in the usual way, we find Z E [J] = N [ Det( 2 + m 2 ) ] 1 exp ( ) d 4 xd 4 x [J (x)g E (x x )J(x )]. (3.16)
8 In order to find the corresponding expression in Minkowski space, we Wick rotate back: 2 + m 2 2 + m 2 (3.17) G E (x x ) ig(x x ). (3.18) Hence, Z[J] = N [ Det( 2 + m 2 ) ] 1 exp (i ) d 4 xd 4 x [J (x)g(x x )J(x )]. (3.19) Note that the functional determinant is raised to the power 1 in contrast to the single real scalar field case in which the determinant is raised to the power 1/2. This is a reflection of the fact that a single complex scalar is equivalent to two real scalars. 3. By inspection of the Minkowski space path integral (before performing the field shift and completing the square), we see that G 2(x x ) = 0 T φ(x)φ(x ) 0 = 1 1 δz[j]. (3.20) i 2 Z[J] δj (x)δj (x ) Using the form of the path integral in Eq. (3.19), we find J=J =0 G 2(x x ) = 0 T φ(x)φ(x ) 0 (3.21) = 1 [ ( ) δ i d 4 yg(x y)j(y) Z[J]] (3.22) Z[J] δj (x ) J=0 = 1 ( 2 d 4 yg(x y)j(y)) Z[J] (3.23) Z[J] J=0 = 0. (3.24) By an identical calculation, we also find G 2 (x x ) = 0 T φ (x)φ (x ) 0 = 1 1 i 2 Z[J] δz[j] δj(x)δj(x ) = 0. (3.25) J=J =0
9 On the other hand, G 2 (x x ) = 0 T φ (x)φ(x ) 0 (3.26) = 1 1 δz[j] i 2 Z[J] δj (x)δj(x ) (3.27) J=J =0 = δ [ ( ) i d 4 yg(x y)j(y) Z[J]] (3.28) δj(x ) J=0 = 1 [ ( ) ( ) ig(x x )Z[J] + i d 4 yg(x y) i d 4 yg(y x )J (y) Z[J]] Z[J] So, (3.29) = ig(x x ). (3.30) By a similar calculation, one finds G 2 (x x ) = 0 T φ (x)φ(x ) 0 = ig(x x ). (3.31) G 2(x x ) = 0 T φ(x)φ (x ) 0 = ig(x x). (3.32) 4. The propagator for the complex scalar field is identical to that of a real scalar field and so all the results in the notes carry apply to the current situation. Here we will simply restate these results as they are derived in the notes (but the reader is expected to work though the derivations in detail). Euclidean space: The Euclidean Green function satisfies [ 2 x + m 2 ]G E (x x ) = δ (4) (x x ) = G E (x x d 4 p e ipµ(x x )µ ) = (2π) 4 p 2 + m. (3.33) 2 Following the derivations in the notes, we obtain G E (x x ) = 1 ( ) m K (2π) 2 x x 1 (m x x ) (3.34) where K ν (z) is the modified Bessel function. For large distances (m x x 1), we can use the asymptotic behaviour of the Bessel function to find that π/2m G E (x x 2 e m x x [ ) 1 + O((m x x ) 1 ) ]. (3.35) (2π) 2 (m x x ) 3/2 Hence the Green function decays exponentially. On the other hand, for short distances ((m x x 1), we find that G E (x x ) 1 4π 2 1 x x 2 (3.36) J=0
10 and so the Green function decays as a power law. Minkowski space: The Minkowski and Euclidean space propagators are related as follows: G(x x ) = ig E (x xi ) x4 ix 0. (3.37) Under a Wick rotation, the Euclidean and Minkowski metrics are related as So, in our case, G(x x ) = x x s 2. (3.38) i m s 4π K 1(m s 2 ). (3.39) 2 2 In order to understand the asymptotic behaviour, we must consider two cases corresponding to when s represents two space-like or time-like separated points. For the space-like separated case, s 2 < 0, we have that s 2 is a real, positive number and so the argument of K ν (z) in the Green function is also a positive number. Hence the Minkowski and Euclidean Green functions are identical (up to an overall factor of i) under the replacement x x s 2. Hence the long-distance behaviour of the Green function is G(x x ) i π/2m 2 e m s 2 [ (2π) 2 (m 1 + O((m ] s 2 ) 1 ). (3.40) s 2 ) 3/2 whereas the short-distance behaviour is G(x x ) i 1 4π 2 s. (3.41) 2 In the time-like separated case (s 2 > 0), we have that s 2 = i s 2 is imaginary. As in the notes, we use the relation K 1 (iz) = π 2 H(1) 1 ( z) where H 1 (z) is a Hankel function. Using the asymptotic behaviour of the Hankel function, we find that in the long-distance regime, G(x x ) π/2 m 2 4π 2 (m s 2 s 2 eim. (3.42) ) 3/2 Hence the Green function exhibits oscillatory behaviour and power-law decay. Conversely, in the short-distance regime, and so we again obtain power-law decay. G(x x ) 1 4π 2 s 2, (3.43)
11 5. We can express the four-point functions in terms of the two-point functions using Wick s theorem. For instance, 1 G a 4(x 1, x 2, x 3, x 4 ) = 0 T φ(x 1 ) φ(x 2 ) φ(x 3 )φ(x 4 ) 0 (3.44) = 0 T φ 1φ 2 0 0 T φ 3 φ 4 0 + 0 T φ 1φ 3 0 0 T φ 2φ 4 0 (3.45) + 0 T φ 1φ 4 0 0 T φ 2φ 3 0 (3.46) = G 2(x 1 x 3 )G 2(x 2 x 4 ) + G 2(x 1 x 4 )G 2(x 2 x 3 ) (3.47) = G 2 (x 1 x 3 )G 2 (x 2 x 4 ) + G 2 (x 1 x 4 )G 2 (x 2 x 3 ), (3.48) where we have used the fact that G 2(x x ) = G 2 (x x ). procedure, we can compute Following the same G 4(x 1, x 2, x 3, x 4 ) = 0 (3.49) G 4 (x 1, x 2, x 3, x 4 ) = 0 (3.50) and G b 4(x 1, x 2, x 3, x 4 ) = G 2 (x 1 x 2 )G 2 (x 3 x 4 ) + G 2 (x 1 x 4 )G 2 (x 2 x 3 ) (3.51) G c 4(x 1, x 2, x 3, x 4 ) = G 2 (x 1 x 2 )G 2 (x 3 x 4 ) + G 2 (x 1 x 3 )G 2 (x 2 x 4 ). (3.52) Comparing the above expressions, we find that G a 4(x 3, x 2, x 1, x 4 ) = G c 4(x 1, x 2, x 3, x 4 ) (3.53) G a 4(x 1, x 4, x 3, x 2 ) = G c 4(x 1, x 2, x 3, x 4 ) (3.54) G a 4(x 1, x 3, x 2, x 4 ) = G b 4(x 1, x 2, x 3, x 4 ) (3.55) G a 4(x 4, x 2, x 3, x 1 ) = G b 4(x 1, x 2, x 3, x 4 ) (3.56) G b 4(x 1, x 2, x 4, x 3 ) = G c 4(x 1, x 2, x 3, x 4 ) (3.57) G b 4(x 2, x 1, x 3, x 4 ) = G c 4(x 1, x 2, x 3, x 4 ). (3.58) These relations follow from using the symmetry of the two-point function. We can also note that the four-point functions are symmetric under permutations of the spacetime positions of the φ s or φ s. For instance, G a 4(x 1, x 2, x 3, x 4 ) = G a 4(x 2, x 1, x 3, x 4 ). Using this fact, the above set of relations are equivalent to the single relation G a 4(x 1, x 4, x 3, x 2 ) = G c 4(x 1, x 2, x 3, x 4 ) = G b 4(x 1, x 2, x 4, x 3 ). (3.59) 1 For brevity, we will write φ i in place of φ(x i ).