ITP-SB-98-04 Mode regularization, time slicing, Weyl ordering and phase space path integrals for quantum mechanical nonlinear sigma models. 8 9 9 1 n Fiorenzo Bastianelli1 a J 5 Dipartimento di Fisica, Universit`a di Modena 1 via Campi 213/A, I-41100, Modena, Italy 1 v 5 Koenraad Schalm2 and Peter van Nieuwenhuizen3 0 1 1 Institute for Theoretical Physics 0 8 State University of New York at Stony Brook 9 / Stony Brook, New York, 11794-3840, USA h t - p e Abstract h : v i X A simple, often invoked, regularization scheme of quantum mechanical path in- r tegrals in curved space is mode regularization: one expands fields into a Fourier a series, performs calculations with only the first M modes, and at the end takes the limit M . This simple scheme does not manifestly preserve reparametrization → ∞ invariance of the target manifold: particular noncovariant terms of order h¯2 must be added to the action in order to maintain general coordinate invariance. Regulariza- tion by time slicing requires a different set of terms of order h¯2 which can be derived from Weyl ordering of the Hamiltonian. With these counterterms both schemes give the same answers to all orders of loops. As a check we perform the three-loop calcu- lation of the trace anomaly in four dimensions in both schemes. We also present a diagrammatic proof of Matthews’ theorem that phase space and configuration space path integrals are equal. 1Also INFN, Sezione di Bologna, Italy. E-mail: [email protected] 2E-mail: [email protected] 3Research supported by National Science Foundation Grant Phy 9722101. E-mail: [email protected] 1 I Introduction Quantum mechanical pathintegrals incurved target spaces have variousapplica- tions. Two decades ago they were used to quantize collective coordinates of solitons [1], and one decade ago to calculate chiral [2] and trace [3, 4] anomalies in quantum field theories. In the former case, the path integrals were regularized by expanding the fluctuations about the solitons into normal modes, and cutting the sum over modes off at a maximum energy or at a maximum number of modes [5]. In the trivial vacuum, on the other hand, two methods have been widely used: the time discretization method1 which has been pioneered by Dirac and Feynman, and the mode truncation method already described in [6]. In the literature, both methods have been used on a par. In this article we point out the relation and the differences between the two methods. Both the time discretization method and the mode truncation method are particular regularization schemes which define the path integral. One expects therefore that the results for physical objects, for example the transition element T = xµ,t xµ,t , can be obtained from path integrals whose actions differ at most h f f| i ii by finite local counterterms. In quantum gauge field theories, one can fix such am- biguities by requiring that Ward identities are satisfied. Similarly, for theories in curved target space which are general coordinate invariant one can require that the transition element be a “bi-scalar” (a scalar in xµ and xµ), but this fixes the ambigu- f i ities only up to covariant terms, namely a term proportional to the scalar curvature R. Only an experiment, for example the measurement of the trace anomaly, can fix the coefficient of this term. One can also view the Hamiltonian operator Hˆ as an observable, and by writing T = xµ,t xµ,t = xµ exp( βHˆ) xµ where β = t t h f f| i ii h f| −¯h | ii f− i in Euclidean space, one can fix the ambiguities in T by requiring that it satisfies the Schroedinger equation with Hˆ as Hamiltonian. If Hˆ itself contains no term proportional to R, then the action will contain a term 1R as we will show, and −8 T will have an R term in the exponent with coefficient 1 . If Hˆ contains a term −12 αR, the coefficient of R in T will be ( 1 +α). In general different regularization − 12 schemes can lead to different finite local counterterms of any order in h¯, but finite local counterterms of order h¯ and h¯2 only are generated in the path integral action if one considers different operator orderings in the Hamiltonian. A Hamiltonian which is a target space scalar2 is 1 Hˆ = g−41(xˆ)pˆµg21(xˆ)gµν(xˆ)pˆνg−41(xˆ) (1) 2 but changing the order of the operators will in general destroy the invariance under targetspacediffeomorphisms. Iftheambiguities duetousing different regularization 1Time discretization in the sector with solitons is complicated because the canonical mo- menta (of the non-zero-modes) satisfy the equal-time commutation relation [π(x,t),ϕ(y,t)] = ′ ′ i¯h[δ(x y) φ (x)φ (y)/M ] where φ (x) is the classicalsoliton solution with mass M . − 2The−gene−ratosrolfor tasorlget spascoel diffeomorspolhisms xˆµ xˆµ +ξµ(xˆ) is a sum of an orbital psaorlt Gˆorb(ξ)= 1 (pˆ ξµ(xˆ)+ξµ(xˆ)pˆ ) and a spin part Gˆspin→(ξ)=(2∂ ξλg ξλ∂ g ) ∂ . Closure 2ih¯ µ µ µ λν − λ µν ∂gµν of the commutator algebra fixes Gˆspin(ξ). An example of a target space diffeomorphism scalar is (1), namely [Gˆorb(ξ)+Gˆspin(ξ),Hˆ]=0, as one may check by an explicit (tedious) calculation. 2 schemes correspond to the ambiguities due to different operator orderings, the terms of order h¯3 and higher in the Hamiltonian should be the same in all cases. This conclusion is corroborated by the observation that in these non-linear sigma models, N-loop graphs are convergent by power counting for N 3. In the time discretization method, a very clear conn≥ection between Hˆ and the transitionelementexistsifoneusesphasespacepathintegrals(Feynman’sapproach) instead of configuration space path integrals (Dirac’s approach), and rewrites the HamiltonianinWeylorderedform. ThenonemayreplacetheWeylorderedoperator Hˆ (pˆ,xˆ) in the kernel by the corresponding function at the midpoint (Berezin’s W theorem) ǫ ǫ x +x ˆ 1 2 dp x exp( H ) p p x exp H(p, ) dp x p p x . (2) 1 W 2 1 2 h | −h¯ | ih | i → −h¯ 2 h | ih | i Z (cid:18) (cid:19)Z The Weyl ordering of Hˆ in (1) leads to the following local finite correction to the naive Hamiltonian H = 1gµνp p to be used in the path integral 2 µ ν h¯2 ∆V = (R+gµνΓ σΓ ρ), (3) W µρ νσ 8 where Γ σ are the Christoffel symbols.3 This counterterm has been extensively µν discussed and used in the literature [7]. Both in the time discretization and in the mode truncation method, one writes the paths as xµ(τ) = xµ (τ)+qµ(τ) for 1 τ 0, where xµ (τ) is a solution of the free field equations witbhgthe correct bou−nda≤ry c≤onditions (xµbg(τ) = xµ+τ(xµ xµ)) bg f f − i andoneexpandsthequantumfluctuationsqµ(τ)intoacompletesetofeigenfunctions ofthe freefield equations with vanishing boundaryvalues andintegrates over afinite number of coefficients qµ. The difference between the time discretization method m and mode regularization is that in the latter one uses the continuum action and naive continuum Feynman rules for a finite number of modes and one chooses a measure for the integrals over these modes which is usually normalized such that it reproduces the standard result for the free particle if the interaction part of the action, S , is set to zero. In the time discretization method, on the other hand, one int uses a discretized action whose discretized Feynman rules are derived (in explicit form [8]) from the Hamiltonian starting point. Up to this point we have made fairly obvious statements. Both regularization schemes are well-defined, and this is how they are used in the literature. However, it is already clear that this cannot be the whole truth, because with time slicing one has the nontrivial ∆V in (3), whereas with mode truncation such a ∆V seems at W first to be absent. The precise way in which also mode truncation leads to a ∆V will be presented in section II. We will see that one needs the following counterterm for mode regularization h¯2 1 ∆V = R gµνgαβg Γ γΓ δ . (4) MR γδ µα νβ 8 − 3 (cid:18) (cid:19) 3OurconventionsfortheRiemannandRiccitensorare: R σ =∂ Γ σ+Γ σΓ τ (ρ µ); R = R µ. Hence, at the linearized level, R = 1(∂ρµν∂ h ∂ρ hµν ∂ ρhτ +µ2νh− ) w↔here ρν ρµν µν 2 µ ν − µ ν − ν µ µν h=ηµνg and h =ηνρ∂ h . µν µ ρ µν 3 Asacheck weshallobtainthecorrecttraceanomalyinfourdimensions fromathree- loop calculation. In Riemann normal coordinates, the Christoffel symbols Γ ρ in µν ∆V vanish at the origin of the coordinate system, but at the three-loop level one finds a contribution from ∆V by expanding each Christoffel symbol into qτ∂ Γ ρ τ µν and contracting the two quantum fields. Since ∂ Γ ρ contains a part proportional τ µν totheRiemanntensor, itisclearthat∆V yields anonvanishing contribution. Hence dropping ∆V in Riemann normal coordinates is incorrect. If one eliminates the momenta by integrating over them in the time discretized path integral, one obtains N factors (detgµν)−1 at the N midpoint coordinates 2 1(qµ +qµ). Exponentiating these in the familiar Faddeev-Popov way one obtains 2 k+1 k what we have called “Lee-Yang ghosts”, namely commuting aµ and anticommuting bµ and cµ ghost fields [4]. Phase space path integrals are free from ambiguities because they are finite. Clearly, upon transition from phase space to configuration space, ambiguities are created; technically this is due to the fact that the momenta are replaced by the q˙’s and ghosts, each of which introduces divergences and hence ambiguities into the theory. In the phase space approach, the vertices are different (for example 1p gµν(x)p 2 µ ν instead of 1x˙µg (x)x˙ν), as well as the propagators ( p p is not proportional to 2 µν h µ νi g g q˙αq˙β ), but the transition element should be the same (Matthews’ theorem). µα νβ h i In section III we present a graphical proof. In section IV we draw conclusions and we show that the coefficient of the R term in ∆V is scheme independent and equals 1. 8 II Mode regularization We now describe how one can define mode regularized path integrals in curved space. Ideally, one would like to derive mode regularization from first principles, i.e., starting from the transition amplitude defined as a matrix element of the evolution operatoracting intheHilbert space ofphysical states: xµ exp( βHˆ) xµ . However, h f| −¯h | ii thisderivationlooksquitedifficult, sothatweprefertotakeamorepragmaticcourse ofaction, andattempt a direct definition of themoderegularized pathintegral. This definitionwillbesupplemented bycertainconsistency requirements whichwespecify later on. The transition amplitude can formally be written as follows x(0)=xf 1 xµ,t xµ,t = ˜x exp S ; (5) h f f| i ii D −βh¯ Zx(−1)=xi (cid:20) (cid:21) 0 1 S = dτ g (x)(x˙µx˙ν +aµaν +bµcν), µν Z−1 2 ˜x = dDx(τ)dDa(τ)dDb(τ)dDc(τ). D −1<τ<0 Y We have shifted and rescaled the time parameter t t = βτ, and since all terms in f − the action only depend on βh¯, we set h¯ = 1 from now on. Note that β then counts the number of loops. We will evaluate the path integral in a perturbative expansion 4 in β and in the coordinate displacements ξµ about the final point: ξµ = xµ xµ. i − f Formally the path integral is a scalar since the action is a scalar and the ghost fields make up a scalar measure on the space of paths.4 Mode regularization will destroy this formal covariance. However, we will see that covariance can be recovered by adding a suitable noncovariant counterterm ∆V to the action S. We start parametrizing xµ(τ) = xµ (τ)+qµ(τ), (6) bg where xµ (τ) is a background trajectory and qµ(τ) the quantum fluctuations. The bg background trajectory is taken to satisfy the free field equations of motion and is a function linear in τ connecting xµ to xµ in the chosen coordinate system, thus i f enforcing the boundary conditions xµ (τ) = xµ ξµτ; ξµ = xµ xµ. (7) bg f − i − f Then the quantum fields qµ(τ) should vanish at the time boundaries and can be expanded into sines. For the Lee-Yang ghosts we use the same Fourier expansion; this may be considered as part of our definition of mode regularization.5 Hence ∞ ∞ qµ(τ) = qµ sin(πmτ); aµ(τ) = aµ sin(πmτ); m m m=1 m=1 X X ∞ ∞ bµ(τ) = bµ sin(πmτ); cµ(τ) = cµ sin(πmτ). (8) m m m=1 m=1 X X At this point the formal measure ˜x can be defined in terms of integration over D the Fourier coefficients ∞ D ˜x = q a b c = A dqµdaµdbµdcµ, (9) D D D D D m m m m m=1µ=1 Y Y which fixes the path integral for a free particle up to the constant A 1 0 1 ˜x exp Q = A; Q = dτ δ (q˙µq˙ν +aµaν +bµcν). (10) µν Z D (cid:20)−β (cid:21) Z−1 2 The constant A will be fixed later on from a consistency requirement. Note that limiting the integration over the number of modes to a finite number M gives a nat- ural regularization of the path integral. This regularization resolves the ambiguities that show up in the continuum limit. 4The factors (detg(x ))1/2dxµ in the discretized path integral are target space scalars and i i exponentiating them leads to ghosts. 5AnotherarguQmenttojustify thattheghostsshouldbeexpandedintosinesisthatthe classical solutions of their field equations are aµ =bµ =cµ =0, and that the quantum fluctuations do not modify the boundary conditions of the classical solutions. In [8] it was shown that the results for the transition amplitude do not change if one uses cosines for the ghosts. 5 Now we expand the action about the final point xµ and obtain f S = S +S = S +S +S +... (11) 2 int 2 3 4 where 0 1 S = dτ g (ξµξν +q˙µq˙ν +aµaν +bµcν), 2 µν Z−1 2 0 1 S = dτ ∂ g (qα ξατ)(ξµξν +q˙µq˙ν +aµaν +bµcν 2q˙µξν), (12) 3 α µν Z−1 2 − − 0 1 S = dτ ∂ ∂ g (qαqβ +ξαξβτ2 2qαξβτ)(ξµξν +q˙µq˙ν +aµaν +bµcν 2q˙µξν). 4 α β µν Z−1 4 − − All geometrical quantities, like g or ∂ g , are evaluated at the final point xµ, but µν α µν f for notational simplicity we do not exhibit this dependence. S is taken as the free 2 part and defines the propagators while S gives the vertices as usual. Therefore, int the quantum perturbative expansion reads: 1 hxµf,tf|xµi,tii = (D˜x) exp −βS = Ae−21βgµνξµξνhe−β1Sinti Z (cid:20) (cid:21) 1 1 1 = Ae−21βgµνξµξν h1− βS3 − βS4 + 2β2S32i+O(β23) . (13) (cid:18) (cid:19) Aiming at a two-loop computation, we have kept only those terms contributing up to O(β), taking into account that ξµ O(√β), as follows from the exponential ∼ appearing in the last line of (13). The propagators that follow from S are given by 2 qµ(τ)qν(σ) = βgµν(x )∆(τ,σ) f h i − aµ(τ)aν(σ) = βgµν(x )••∆(τ,σ) (14) f h i bµ(τ)cν(σ) = 2βgµν(x )••∆(τ,σ) f h i − where ∆ is regulated by the mode cut-off introduced below (10): M 2 ∆(τ,σ) = sin(πmτ)sin(πmσ) , (15) −π2m2 m=1(cid:20) (cid:21) X and has as continuum value ∆(τ,σ) = τ(σ +1)θ(τ σ)+σ(τ +1)θ(σ τ). (16) − − (A dot on the left(right) of ∆(τ,σ) indicates differentiation with respect to τ(σ).) Using standard Wick contractions, we computed: 1 1 1 S = ∂ g ξαξµξν, (17) 3 α µν h−β i −β 4 1 β 1 S = ∂ ∂ g (gµνgαβ gµαgνβ) (gµνξαξβ +gαβξµξν 2gµαξνξβ) 4 α β µν h−β i 24 − − 24 − (cid:20) 1 1 ξµξνξαξβ , (18) −β 12 (cid:21) 6 1 β S2 = ∂ g ∂ g (gαβgµνgλρ 4gαρgµνgβλ 6gαβgµλgνρ h2β2 3i α µν β λρ 96 − − (cid:20) +4gαρgβµgνλ+4gαµgβλgνρ) 1 + gµλgνρξαξβ +2(gαβgµλ gαλgµβ)ξνξρ+(2gαλgβρ gαβgλρ)ξµξν 48 − − (cid:18) +(2gµβgλρ 4gµλgβρ)ξαξν − (cid:19) 1 1 + (gαβξµξνξλξρ 4gαλξµξνξβξρ +4gµλξαξνξβξρ) β96 − 1 1 + ξαξµξνξβξλξρ . (19) β232 (cid:21) This completes the calculation of the transition amplitude in the two-loop approxi- mation using the mode regularized path integral. Atthispointweshouldmakecontactwithotherschemes andtesttheconsistency ofourrules. Todothat,weuseourtransitionamplitudetoobtainthetimeevolution of an arbitrary wave function Ψ(x ,t ) = dDx g(x ) xµ,t xµ,t Ψ(x ,t ). (20) f f i i h f f| i ii i i Z q We need the factor g(x ) because the transition element is formally a biscalar as i we explained before,qbut then Ψ(x ,t ) is a scalar and hence also Ψ(x ,t ), which in f f i i turn implies that the measure must be a scalar as well. Since the transition amplitude (13) is given in terms of an expansion around the final point (x ,t ), we Taylor expand also the wave function Ψ(x ,t ) and the f f i i measure g(x ) in eq. (20) about this point i q 1 Ψ(xi,ti) = Ψ(xf,tf) β∂tΨ(xf,tf)+ξµ∂µΨ(xf,tf)+ ξµξν∂µ∂νΨ(xf,tf)+O(β23) − 2 1 g(x ) = g(x ) 1+ξµΓ α + ξµξν(∂ Γ α +Γ αΓ β)+O(ξ3) . (21) i f µα µ να µα νβ 2 q q (cid:18) (cid:19) We then perform the integration over dDxµ = dDξµ in (20) and match the various terms. The leading term fixes A D −D Ψ = A(2πβ)2Ψ A = (2πβ) 2, (22) → and the terms of order β give 1 1 1 1 β ∂ Ψ+ 2Ψ+ RΨ gµνgαβgγδ∂ g ∂ g Ψ+ gµνgαβgγδ∂ g ∂ g Ψ = 0. t µ αγ ν βδ µ αγ β νδ − 2∇ 8 −32 48 (cid:20) (cid:21) (23) This last equation means that the wave function Ψ satisfies the following Schroe- dinger equation at the final point (xµ,t ): f f 1 ∂ Ψ = (H +∆V )Ψ = 2Ψ+∆V Ψ, (24) t 0 eff eff − −2∇ 7 where the effective potential is given by 1 1 1 ∆V = R+ gµνgαβgγδ∂ g ∂ g gµνgαβgγδ∂ g ∂ g eff µ αγ ν βδ µ αγ β νδ −8 32 − 48 1 1 1 = R gµνΓ βΓ α + gµνgαβgγδ∂ g ∂ g (25) µα νβ µ αγ β νδ −8 − 8 24 1 1 = R+ gµνgαβg Γ γΓ δ. γδ µα νβ −8 24 Clearly, to obtain the “free” Hamiltonian H from a path integral, one should sub- 0 tract the potential ∆V , and thus use in (5) the following classical action eff 1 0 1 S = dτ g (x)x˙µx˙ν +β∆V (x) . (26) µν eff − β Z−1 (cid:18)−2β (cid:19) With this counterterm mode regularization gives the same results as the time dis- cretization scheme described in the introduction. Note that (25) is different from the counterterm for time discretization. The difference is the last term in the second line of (25). Finally, as a non-trivial test, which was one of our motivations to carry out the present research, we have computed the trace anomaly for a real conformal scalar field in four dimensions using mode regularization andthe newly found counterterm. Anomalies in Fujikawa’s approach are given by the regulated trace of the Jacobian J of the symmetry transformation. −βR An = limTr(J e h¯ ). (27) β→0 Atheoryofconsistentregulators exists[9]anditgivesusforconformalscalarfields the Hamiltonian (1) minus an imRprovement potential h¯2 D−2 R. The Jacobian for 8(D−1) trace anomalies is proportional to unity. In mode regularization (27) translates into the evaluation of the expectation value of the Jacobian with respect to the path integral based on the formal transition element (5) including the just calculated counterterms andthe improvement term. Taking the trace of (5) means first putting x = x which gives that x = x and then integrating over x with dDx g(x ). f i bg f f f f As the measure of (5) includes the constant A = (2πβ)−D, the β-independeqnt term 2 of the r.h.s. of (27) is given by the (D/2 + 1)-loop term in the path integral. Expanding the metric in Riemann normal coordinates we find the first nonvanishing contributions at the two-loop level: graphs with the topology of the number 8 and the background values of the various potential terms. From (18) with ξµ = 0 one obtains6 1 1 (1 (D 2)/(D 1)) + = ( βh¯) R+( βh¯) − − − R. (28) • − 6 −4 − 8 (cid:18) (cid:19) The result is the same as obtained from time slicing and yields the correct trace anomaly in D = 2 dimensions. 6In equations (28), (29), (30), (31) only the topology of the graphs is indicated, e.g., the figure 8 in (28) stands for all possible graphs of that shape, including ghost loops. 8 At the three-loop level we find three regular graphs 1 1 = ( βh¯)2 R2 , (29) 72 − −6 µν (cid:18) (cid:19) 1 1 = ( βh¯)2 R2 , (30) 72 − −6 µνλρ (cid:18) (cid:19) 1 1 1 = ( βh¯)2 2R+ R2 + R2 ,(31) − 480∇ 720 µνλρ 1080 µν (cid:18) (cid:19) as well as a graph coming from the potential terms 1 1 (D 2)/(D 1) 1 = (βh¯)2 − − − 2R R2 . (32) −6 16 ∇ − 72 µνλρ! There are only two differences with respect to time slicing: with time slicing, the factor 1 in (30) becomes 1, and the factor 1 in (32) becomes 1 . These −6 −4 −72 −48 two modifications lead to the same final expression. Adding all contributions of connected and disconnected graphs we find the correct result d4x 1 An(Weyl)(spin0,D = 4) = g(x)σ(x) R2 R2 2R . (33) (2π)2 720 µνλρ − µν −∇ Z q (cid:16) (cid:17) where σ(x) is the arbitrary function appearing in the Jacobian: J = σ(x)δD(x y). − This shows that mode regularization works after all.7 The fact that the right answer isobtainedsuggeststhatnonewcountertermsareneededinthisscheme. Acomplete three-loop calculation in arbitrary coordinates and with a non-vanishing ξµ could be used to test eq. (20) at order β2. As we noted in the introduction, however, all three- and higher-loop graphs are power-counting convergent. This means that at these orders in h¯ any consistent regularization scheme will yield the same answer. We therefore claim that with the presently found counterterm, the mode regularized path integral is consistent to all orders. This suggests that the difference is indeed due to operator-orderings. It would be interesting to justify mode regularization from first principles. Before closing this section, we comment on how loop integrations in Feynman diagrams are done. Mode cut-off allows one to disentangle ambiguities that appear in the continuum limit of certain integrals over the ∆’s. Resorting to the mode regulated expression for ∆, one can use partial integration and take boundary terms into account if they are non-vanishing. Using partialintegration repeatedly, one gets • • expressions containing only ∆’s, ∆ ’s and ∆’s which are unambiguous in the con- tinuum limit, and computes them there. Useful identities obtained at the regulated 7 In [4] the noncovariantpart of the counterterm (25) was missed. In Riemann normal coordi- natesthistermaffectsonlythecoefficientoftheR2 terminthefourdimensionaltraceanomaly, µνλρ which was erroneously calculated. 9 level from (15) are • • •• • ∆ (τ,τ)+ ∆(τ,τ) = ∂ (∆(τ,τ)), τ •• •• ∆(τ,σ) = ∆ (τ,σ). (34) The first identity was used to compute (28). Here is a list of integrals that are easily computed using partial integration8 and whose values are different from the time discretization method: • • •• I = dτ τ (∆ + ∆) = 0, 1 σ=τ | Z 1 • • • • I = dτdσ (∆) (∆ ) (∆ ) = , (35) 2 −12 ZZ 1 • • • • I = dτdσ ∆ (∆) (∆ ) (∆ ) = . 3 180 ZZ Time discretization would give I = 1, I = 1 and I = 7 .9 1 −2 2 −6 3 360 The non-vanishing value of I in time discretization is needed to compensate the 1 explicit factors of g−1/4(x ) appearing in the measure, see (36), and together with i I it is responsible for the different counterterms required in the two regularization 2 schemes. Finally, I leads to the different values for the coefficient of R2 in 3 µνλρ eq. (30) in the two schemes. The counterterms with two Christoffel symbols are also different in both schemes but the final result for the transition element (and hence the trace anomaly) is the same. This confirms our approach to mode regularization. III Phase space path integrals in curved space and Matthews’ theorem In the phase space approach to path integrals in curved space the vertices and propagators are different from those used in the configuration space approach. For example, the leading terms in the action are 1p gµν(x)p and 1x˙µg (x)x˙ν, respec- 2 µ ν 2 µν tively, and the q˙µ(τ)q˙ν(σ) propagator contains a Dirac delta function but the h i p (τ)p (σ) propagator is completely regular. Although the vertices and propa- µ ν h i gators differ, the result for the transition element should be the same (Matthews’ theorem10) [10], and the way this comes about is due to a new kind of ghosts [3, 4]. 8To obtain I2, use (∆•)(•∆•) = 12∂τ(∆•)2, partially integrate, use ••∆ = ∆••, and then use (∆••)(∆•)2 = 31∂σ(∆•)3. ToobtainI3,writetheintegrandas∆(•∆)12∂τ(∆•)2,partiallyintegrate, and use then that ∆(••∆)(∆•)2 = ∆(∆••)(∆•)2 = ∆1∂ (∆•)3. Then partially integrate once 3 τ more. 9In the time discretization scheme δ(σ τ) is a Kronecker delta function and θ(0) = 1. This 2 − leadstoafullyconsistentschemeasshownin[8]. Forexample, dσdτδ(σ τ)θ(σ τ)θ(σ τ)= 1 − 1 − 1 −1 1 dσdτδ(σ τ)θ(σ τ)θ(τ σ) = whereas mode regularization would give and = − •− − 4 • • RR •• 3 2 − 3 6 respectively. Using ∆(τ,σ) = σ+θ(τ σ), ∆ (τ,σ) = 1 δ(τ σ) and ∆ = δ(τ σ), the RR − − − − results below (35) follow. 10Matthews’ theorem holds for connected graphs and originally only applied to nonderivative interactions. 10