Gaussian integrals and Feynman diagrams. February 28

Size: px
Start display at page:

Download "Gaussian integrals and Feynman diagrams. February 28"

Transcription

1 Gaussian integrals and Feynman diagrams February 28

2 Introduction A mathematician is one to whom the equality e x2 2 dx = 2π is as obvious as that twice two makes four is to you. Lord W.T. Kelvin to his students

3 Introduction A mathematician is one to whom the equality e x2 2 dx = 2π is as obvious as that twice two makes four is to you. A physicist must be one to whom the formula Lord W.T. Kelvin to his students is equally obvious.

4 Introduction Generalizing P (x)e x2 2 dx e x2 2 dx = 2π: P (x)e x,ax 2 dx Here, A is a positively-defined R n symmetric matrix and, is the standard inner product on R n. P (x)e x,ax 2 r 3 r/2 1 1 r! Br(x,...,x) dx, where is a R n parameter.

5 Introduction Generalizing P (x)e x2 2 dx e x2 2 dx = 2π: P (x)e x,ax 2 dx Here, A is a positively-defined R n symmetric matrix and, is the standard inner product on R n. P (x)e x,ax 2 r 3 r/2 1 1 r! Br(x,...,x) dx, where is a R n parameter.

6 Introduction Generalizing P (x)e x2 2 dx e x2 2 dx = 2π: P (x)e x,ax 2 dx Here, A is a positively-defined R n symmetric matrix and, is the standard inner product on R n. P (x)e x,ax 2 r 3 r/2 1 1 r! Br(x,...,x) dx, where is a R n parameter.

7 Motivation In classical mechanics, the motion of a point-particle of mass m in a potential field V = V (x) is determined by the Newton s equation mx = V (x). Alternatively, one can state the law of motion in the form of the stationary action principle. One starts by introducting the Lagrangian of the system L(x, x, t) = (kinetic energy) (potential energy) and the action functional S[x] = b a L(x, x, t) dt. Stationary action principle The trajectory of a particle x = x(t) (t [a, b]) has to be an extremum of the action functional S.

8 Motivation In classical mechanics, the motion of a point-particle of mass m in a potential field V = V (x) is determined by the Newton s equation mx = V (x). Alternatively, one can state the law of motion in the form of the stationary action principle. One starts by introducting the Lagrangian of the system L(x, x, t) = (kinetic energy) (potential energy) and the action functional S[x] = b a L(x, x, t) dt. Stationary action principle The trajectory of a particle x = x(t) (t [a, b]) has to be an extremum of the action functional S.

9 Motivation In classical mechanics, the motion of a point-particle of mass m in a potential field V = V (x) is determined by the Newton s equation mx = V (x). Alternatively, one can state the law of motion in the form of the stationary action principle. One starts by introducting the Lagrangian of the system L(x, x, t) = (kinetic energy) (potential energy) and the action functional S[x] = b a L(x, x, t) dt. Stationary action principle The trajectory of a particle x = x(t) (t [a, b]) has to be an extremum of the action functional S.

10 Motivation In classical mechanics, the motion of a point-particle of mass m in a potential field V = V (x) is determined by the Newton s equation mx = V (x). Alternatively, one can state the law of motion in the form of the stationary action principle. One starts by introducting the Lagrangian of the system L(x, x, t) = (kinetic energy) (potential energy) and the action functional S[x] = b a L(x, x, t) dt. Stationary action principle The trajectory of a particle x = x(t) (t [a, b]) has to be an extremum of the action functional S.

11 Motivation Stationary action principle The trajectory of a particle x = x(t) (t [a, b]) has to be an extremum of the action functional S. This means that such a trajectory has to be a solution of the variational problem δs = 0. The latter can be reduced to solving the Euler-Lagrange equation (or rather a system of EL equations). In the one-dimensional case it has the form L x ( ) L t x = 0, x(a) = x 0, x(b) = x 1. Example The Lagranian for a mass on a spring system is L = mx 2 2 kx2 2. The EL equation reads: kx t mx = 0. That is, mx = kx (with some boundary conditions).

12 Motivation Stationary action principle The trajectory of a particle x = x(t) (t [a, b]) has to be an extremum of the action functional S. This means that such a trajectory has to be a solution of the variational problem δs = 0. The latter can be reduced to solving the Euler-Lagrange equation (or rather a system of EL equations). In the one-dimensional case it has the form L x ( ) L t x = 0, x(a) = x 0, x(b) = x 1. Example The Lagranian for a mass on a spring system is L = mx 2 2 kx2 2. The EL equation reads: kx t mx = 0. That is, mx = kx (with some boundary conditions).

13 Motivation Stationary action principle The trajectory of a particle x = x(t) (t [a, b]) has to be an extremum of the action functional S. This means that such a trajectory has to be a solution of the variational problem δs = 0. The latter can be reduced to solving the Euler-Lagrange equation (or rather a system of EL equations). In the one-dimensional case it has the form L x ( ) L t x = 0, x(a) = x 0, x(b) = x 1. Example The Lagranian for a mass on a spring system is L = mx 2 2 kx2 2. The EL equation reads: kx t mx = 0. That is, mx = kx (with some boundary conditions).

14 Motivation Remark The stationary action principle is often referred to as the principle of the least action, but, in general, classical trajectories do not have to minimize the action. In classical field theory one studies not the motion of point-particles, but rather motion of a continuum of particles (e.g. a string, a membrane, a jet of fluid etc.) Now, the trajectory of a system is described by a classical field φ = φ(x, t). The stationary action principle still holds: one starts with the Lagrangian L = L(φ,... ) and looks for φ s that give the extremum to the action functional S[φ] = D L(φ,... ) dx dt. Example The Lagrangian[ for a string has the form L(z, z x, z t) = 1 2 m ( ) z 2 ( t λ z ) 2 ] x. The EL equation for the Gaussian integrals 2 z and Feynman λ 2 z diagrams

15 Motivation In classical field theory one studies not the motion of point-particles, but rather motion of a continuum of particles (e.g. a string, a membrane, a jet of fluid etc.) Now, the trajectory of a system is described by a classical field φ = φ(x, t). The stationary action principle still holds: one starts with the Lagrangian L = L(φ,... ) and looks for φ s that give the extremum to the action functional S[φ] = D L(φ,... ) dx dt. Example The Lagrangian[ for a string has the form L(z, z x, z t) = 1 2 m ( ) z 2 ( t λ z ) 2 ] x. The EL equation for the variational problem δs = 0 takes the form 2 z λ 2 z t 2 m = 0. x 2

16 Motivation In classical field theory one studies not the motion of point-particles, but rather motion of a continuum of particles (e.g. a string, a membrane, a jet of fluid etc.) Now, the trajectory of a system is described by a classical field φ = φ(x, t). The stationary action principle still holds: one starts with the Lagrangian L = L(φ,... ) and looks for φ s that give the extremum to the action functional S[φ] = D L(φ,... ) dx dt. Example The Lagrangian[ for a string has the form L(z, z x, z t) = 1 2 m ( ) z 2 ( t λ z ) 2 ] x. The EL equation for the variational problem δs = 0 takes the form 2 z λ 2 z t 2 m = 0. x 2

17 Motivation In classical field theory one studies not the motion of point-particles, but rather motion of a continuum of particles (e.g. a string, a membrane, a jet of fluid etc.) Now, the trajectory of a system is described by a classical field φ = φ(x, t). The stationary action principle still holds: one starts with the Lagrangian L = L(φ,... ) and looks for φ s that give the extremum to the action functional S[φ] = D L(φ,... ) dx dt. Example The Lagrangian[ for a string has the form L(z, z x, z t) = 1 2 m ( ) z 2 ( t λ z ) 2 ] x. The EL equation for the variational problem δs = 0 takes the form 2 z λ 2 z t 2 m = 0. x 2

18 Motivation Let F be a space of fields (for us: scalar or vector-valued functions of space-time). An observable f is a function f : F R. On the quantum level, the behavior of physical systems is no longer deterministic and we cannot use the SAP directly. Instead, we just hope to find the expectation values of observables: for f : F R, f = 1 f(φ)e i S[φ] Dφ, Z F where Z = F e i S[φ] Dφ is the normalizing factor (partition function).

19 Motivation Let F be a space of fields (for us: scalar or vector-valued functions of space-time). An observable f is a function f : F R. On the quantum level, the behavior of physical systems is no longer deterministic and we cannot use the SAP directly. Instead, we just hope to find the expectation values of observables: for f : F R, f = 1 f(φ)e i S[φ] Dφ, Z F where Z = F e i S[φ] Dφ is the normalizing factor (partition function).

20 Motivation Let F be a space of fields (for us: scalar or vector-valued functions of space-time). An observable f is a function f : F R. On the quantum level, the behavior of physical systems is no longer deterministic and we cannot use the SAP directly. Instead, we just hope to find the expectation values of observables: for f : F R, f = 1 f(φ)e i S[φ] Dφ, Z F where Z = F e i S[φ] Dφ is the normalizing factor (partition function).

21 Motivation f = 1 f(φ)e i S[φ] Dφ, Z F Intuitively, we counting contributions from all possible φ s, but the ones which are closer to the extrema of the action functional S[φ] yield greater contributions (if φ deviates a lot from a classical trajectory, then oscillations of e i S[ ] near φ will cancel each other out). Stationary phase formula Let f = f(x) be a smooth function with a unique critical point c (a, b), f (c) 0 and g = g(x) be a smooth function with vanishing derivatives at x = a, b. Then b a g(x)e i f(x) dx = 1/2 e if(c)/ I( ), where I is a smooth function g(c) on [0, ) such that I(0) = 2πe ±πi/4. f (c)

22 Motivation f = 1 f(φ)e i S[φ] Dφ, Z F Intuitively, we counting contributions from all possible φ s, but the ones which are closer to the extrema of the action functional S[φ] yield greater contributions (if φ deviates a lot from a classical trajectory, then oscillations of e i S[ ] near φ will cancel each other out). Stationary phase formula Let f = f(x) be a smooth function with a unique critical point c (a, b), f (c) 0 and g = g(x) be a smooth function with vanishing derivatives at x = a, b. Then b a g(x)e i f(x) dx = 1/2 e if(c)/ I( ), where I is a smooth function g(c) on [0, ) such that I(0) = 2πe ±πi/4. f (c)

23 Motivation The (very) bad thing about the formula f = 1 Z F f(φ)e i S[φ] Dφ: since the space of fields F is, in general, huge (infinitely-dimensional), Dφ is ill-defined. Yet, it makes sense to approximate such integrals by their finite-dimensional analogs. This leads to the problem of studying oscillating integrals of the form P (x)e i S[x] dx R n and their real counterparts R n P (x)e 1 S[x] dx

24 Motivation The (very) bad thing about the formula f = 1 Z F f(φ)e i S[φ] Dφ: since the space of fields F is, in general, huge (infinitely-dimensional), Dφ is ill-defined. Yet, it makes sense to approximate such integrals by their finite-dimensional analogs. This leads to the problem of studying oscillating integrals of the form P (x)e i S[x] dx R n and their real counterparts R n P (x)e 1 S[x] dx

25 Motivation The (very) bad thing about the formula f = 1 Z F f(φ)e i S[φ] Dφ: since the space of fields F is, in general, huge (infinitely-dimensional), Dφ is ill-defined. Yet, it makes sense to approximate such integrals by their finite-dimensional analogs. This leads to the problem of studying oscillating integrals of the form P (x)e i S[x] dx R n and their real counterparts R n P (x)e 1 S[x] dx

26 Steepest descent and stationary phase Stationary phase formula Let f = f(x) be a smooth function with a unique critical point c (a, b), f (c) 0 and g = g(x) be a smooth function with vanishing derivatives at x = a, b. Then b a g(x)e i f(x) dx = 1/2 e if(c)/ I( ), where I is a smooth function g(c) on [0, ) such that I(0) = 2πe ±πi/4. f (c) Steepest descent formula Let f = f(x) be a smooth function with a unique minimum point c (a, b), f (c) > 0 and g = g(x) be a smooth function. Then b a g(x)e 1 f(x) dx = 1/2 e f(c)/ I( ), where I is a smooth g(c) function on [0, ) such that I(0) = 2π. f (c)

27 Steepest descent and stationary phase Let B be a box region in R n and f = f(x) be a smooth function with a unique critical point c B such that... Stationary phase formula... Hf(c) 0. Then if g = g(x) is a smooth function with vanishing derivatives at the boundary of B, we have g(x)e i f(x) dx = n/2 e if(c)/ I( ), where I is a smooth function B on [0, ) such that I(0) = (2π) n/2 e ±πiσ/4 g(c) Hf(c) Steepest descent formula... Hf(c) > 0. Then if g = g(x) is a smooth function, we have g(x)e 1 f(x) dx = n/2 e f(c)/ I( ), where I is a smooth B function on [0, ) such that I(0) = (2π) n/2 g(c) Hf(c)

28 Asymptotic expansion B g(x)e 1 f(x) dx = n/2 e f(c)/ I( ) }{{}? What do we want: find a power series expansion for I( ) = A 0 + A 1 + A (Not so) bad thing: Although, I = I( ) is smooth, it is not analytic at in general: the Taylor series for I at = 0 may have the zero radius of convergence. How to fix it: We still can try to find a formal power series expansion for I. Then under some pretty mild conditions on A i s there is canonical way to sum the series up (key term: Borel summation). Thus, we can focus on finding the power series coefficients A i.

29 Asymptotic expansion B g(x)e 1 f(x) dx = n/2 e f(c)/ I( ) }{{}? What do we want: find a power series expansion for I( ) = A 0 + A 1 + A (Not so) bad thing: Although, I = I( ) is smooth, it is not analytic at in general: the Taylor series for I at = 0 may have the zero radius of convergence. How to fix it: We still can try to find a formal power series expansion for I. Then under some pretty mild conditions on A i s there is canonical way to sum the series up (key term: Borel summation). Thus, we can focus on finding the power series coefficients A i.

30 Asymptotic expansion B g(x)e 1 f(x) dx = n/2 e f(c)/ I( ) }{{}? What do we want: find a power series expansion for I( ) = A 0 + A 1 + A (Not so) bad thing: Although, I = I( ) is smooth, it is not analytic at in general: the Taylor series for I at = 0 may have the zero radius of convergence. How to fix it: We still can try to find a formal power series expansion for I. Then under some pretty mild conditions on A i s there is canonical way to sum the series up (key term: Borel summation). Thus, we can focus on finding the power series coefficients A i.

31 Asymptotic expansion B g(x)e 1 f(x) dx = n/2 e f(c)/ I( ) }{{}? What do we want: find a power series expansion for I( ) = A 0 + A 1 + A (Not so) bad thing: Although, I = I( ) is smooth, it is not analytic at in general: the Taylor series for I at = 0 may have the zero radius of convergence. How to fix it: We still can try to find a formal power series expansion for I. Then under some pretty mild conditions on A i s there is canonical way to sum the series up (key term: Borel summation). Thus, we can focus on finding the power series coefficients A i.

32 Wick s theorem Theorem [Gian-Carlo Wick] Let A be a symmetric, positively-defined d-by-d matrix and l i s be linear forms on R d (i = 1, 2,..., m). Then R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A where σ s run through all pairings on {1, 2,..., m} (that is, σ is an involution on {1, 2,..., m}). σ i

33 Wick s theorem Theorem [Gian-Carlo Wick] Let A be a symmetric, positively-defined d-by-d matrix and l i s be linear forms on R d (i = 1, 2,..., m). Then R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A where σ s run through all pairings on {1, 2,..., m} (that is, σ is an involution on {1, 2,..., m}). Example m = 0, d = 1, A = 1: R e x2 /2 dx = 2π σ i

34 Wick s theorem Theorem [Gian-Carlo Wick] Let A be a symmetric, positively-defined d-by-d matrix and l i s be linear forms on R d (i = 1, 2,..., m). Then R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A where σ s run through all pairings on {1, 2,..., m} (that is, σ is an involution on {1, 2,..., m}). Example [ ] 2 1 m = 2, d = 2, A =, l = x, l 2 = y: R xy e 2 (x2 +xy+y 2) dxdy = 2π σ i

35 Wick s theorem Theorem [Gian-Carlo Wick] Let A be a symmetric, positively-defined d-by-d matrix and l i s be linear forms on R d (i = 1, 2,..., m). Then R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A where σ s run through all pairings on {1, 2,..., m} (that is, σ is an involution on {1, 2,..., m}). Example m = 4: R d (... ) l 1, A 1 l 2 l 3, A 1 l 4 + l 1, A 1 l 3 l 2, A 1 l 4 σ i + l 1, A 1 l 4 l 2, A 1 l 4

36 Wick s theorem It is convenient to represent pairings using graphs. Namely, for a pairing σ on a set {1, 2,..., m}, consider a graph with vertices indexed by 1, 2,..., m and connect vertices i,σ(i) with an edge. Example The graph represents the pairing 1, 3, 2, 5, 4, 6.

37 Wick s theorem Example ( (2π) d/2 det A ) 1 l 1 (x)... l 4 (x)e x,ax /2 dx = R d

38 Wick s theorem: sketch of the proof The main idea: reduce everything to the one-dimensional case, where all computations can be performed in terms of the Gamma-function. Useful facts: Any real symmetric matrix A can be diagonalized. Polarization identity. Let P : V V R be a symmetric polylinear form. Then P (x 1, x 2,..., x n ) = 1 n! I {1,...,n}( 1) n I P ( x i,..., x i ). i I i I Γ ( ) 2k+1 k = (2k)! 4 k k! π

39 Wick s theorem: sketch of the proof The main idea: reduce everything to the one-dimensional case, where all computations can be performed in terms of the Gamma-function. Useful facts: Any real symmetric matrix A can be diagonalized. Polarization identity. Let P : V V R be a symmetric polylinear form. Then P (x 1, x 2,..., x n ) = 1 n! I {1,...,n}( 1) n I P ( x i,..., x i ). i I i I Γ ( ) 2k+1 k = (2k)! 4 k k! π

40 Wick s theorem: sketch of the proof The main idea: reduce everything to the one-dimensional case, where all computations can be performed in terms of the Gamma-function. Useful facts: Any real symmetric matrix A can be diagonalized. Polarization identity. Let P : V V R be a symmetric polylinear form. Then P (x 1, x 2,..., x n ) = 1 n! I {1,...,n}( 1) n I P ( x i,..., x i ). i I i I Γ ( ) 2k+1 k = (2k)! 4 k k! π

41 Wick s theorem: sketch of the proof R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A σ i We can assume that m is even. Step 1. Both sides of the desired identity are symmetric and polylinear forms with respect to l i s. Thus, it suffices to prove it for l 1 = = l m. Step 2. The desired identity is stable under a linear change of variables. Thus, we can choose a basis in such a way that A becomes diagonal. Moreover, by rescaling, we can make A = E. Then x, Ax = x x2 d Step 3. We find x 2k e x2 /2 dx = 2π (2k)! 2 k k! by substituting y = x 2 /2 and reducing it to the Gamma-function.

42 Wick s theorem: sketch of the proof R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A σ i We can assume that m is even. Step 1. Both sides of the desired identity are symmetric and polylinear forms with respect to l i s. Thus, it suffices to prove it for l 1 = = l m. Step 2. The desired identity is stable under a linear change of variables. Thus, we can choose a basis in such a way that A becomes diagonal. Moreover, by rescaling, we can make A = E. Then x, Ax = x x2 d Step 3. We find x 2k e x2 /2 dx = 2π (2k)! 2 k k! by substituting y = x 2 /2 and reducing it to the Gamma-function.

43 Wick s theorem: sketch of the proof R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A σ i We can assume that m is even. Step 1. Both sides of the desired identity are symmetric and polylinear forms with respect to l i s. Thus, it suffices to prove it for l 1 = = l m. Step 2. The desired identity is stable under a linear change of variables. Thus, we can choose a basis in such a way that A becomes diagonal. Moreover, by rescaling, we can make A = E. Then x, Ax = x x2 d Step 3. We find x 2k e x2 /2 dx = 2π (2k)! 2 k k! by substituting y = x 2 /2 and reducing it to the Gamma-function.

44 Wick s theorem: sketch of the proof R d l 1 (x)... l m (x)e x,ax /2 dx = (2π)d/2 l i, A 1 l σ(i), det A σ i We can assume that m = 2k. Step 4. What is (2k)!? It s the number of pairings on a set 2 k k! {1, 2..., 2k}.

45 Feynman s theorem Back to our main problem: find the coefficients of the asymptotic power series expansion of the integral P (x)e S(x)/ dx, B where P is a polynomial and S is a smooth function having a unique minimum critical point c B.

46 Feynman s theorem Back to our main problem: find the coefficients of the asymptotic power series expansion of the integral l 1 (x)... l m (x)e S(x)/ dx, B where S is a smooth function having a unique minimum critical point c B.

47 Feynman s theorem Back to our main problem: find the coefficients of the asymptotic power series expansion of the integral l 1 (x)... l m (x)e S(x)/ dx, B where S is a smooth function having a unique minimum critical point c B. WLOG, we can assume c = 0, S(c) = 0. Then S(x) has a (formal, at least) power series expansion S(x) = x, Ax 2 + r 3 1 r! B r(x,..., x) where B r s are symmetric polylinear forms.

48 Feynman s theorem Make a change of variables x x : l 1 (x)... l m (x)e S(x)/ dx = B m/2 l 1 (x)... l m (x)e x,ax B m/2 l 1 (x)... l m (x)e x,ax R d 2 r 3 r/2 1 Br(x,...,x) 2 r 3 r/2 1 Br(x,...,x) r! dx = r! dx + o( )

49 Feynman diagrams Definition A Feynman diagram with m external vertices is a graph (possibly, with loops and multiple edges) with m vertices of degree 1 labeled by 1, 2,..., m and finitely many unlabeled vertices of degrees 3. G 3 (m) := {isomorphism classes of Feynman diagrams with m external vertices}. Here, an isomorphism of a labeled graph is supposed to preserve the labeling.

50 Feynman diagrams Definition A Feynman diagram with m external vertices is a graph (possibly, with loops and multiple edges) with m vertices of degree 1 labeled by 1, 2,..., m and finitely many unlabeled vertices of degrees 3. G 3 (m) := {isomorphism classes of Feynman diagrams with m external vertices}. Here, an isomorphism of a labeled graph is supposed to preserve the labeling.

51 Feynman diagrams Definition A Feynman diagram with m external vertices is a graph (possibly, with loops and multiple edges) with m vertices of degree 1 labeled by 1, 2,..., m and finitely many unlabeled vertices of degrees 3. G 3 (m) := {isomorphism classes of Feynman diagrams with m external vertices}. Here, an isomorphism of a labeled graph is supposed to preserve the labeling. Example

52 Feynman s theorem Theorem l 1... l m := m/2 l 1 (x)... l m (x)e x,ax R d where = (2π)d/2 det A Γ G 3 (m) 2 r 3 r/2 1 Br(x,...,x) r! dx b(γ) Aut(Γ) F Γ(l 1,..., l m ), b(γ) = Edges of Γ internal vertices of Γ ; Aut(Γ) is the number of automorphisms of a graph Γ which leave the external vertices fixed; F Γ (l 1,..., l m ) is the Feynman amplitude of the graph Γ computed as follows.

53 Feynman s theorem Theorem l 1... l m := m/2 l 1 (x)... l m (x)e x,ax R d where = (2π)d/2 det A Γ G 3 (m) 2 r 3 r/2 1 Br(x,...,x) r! dx b(γ) Aut(Γ) F Γ(l 1,..., l m ), b(γ) = Edges of Γ internal vertices of Γ ; Aut(Γ) is the number of automorphisms of a graph Γ which leave the external vertices fixed; F Γ (l 1,..., l m ) is the Feynman amplitude of the graph Γ computed as follows.

54 Feynman s theorem Theorem l 1... l m := m/2 l 1 (x)... l m (x)e x,ax R d where = (2π)d/2 det A Γ G 3 (m) 2 r 3 r/2 1 Br(x,...,x) r! dx b(γ) Aut(Γ) F Γ(l 1,..., l m ), b(γ) = Edges of Γ internal vertices of Γ ; Aut(Γ) is the number of automorphisms of a graph Γ which leave the external vertices fixed; F Γ (l 1,..., l m ) is the Feynman amplitude of the graph Γ computed as follows.

55 Feynman s theorem Theorem l 1... l m := m/2 l 1 (x)... l m (x)e x,ax R d where = (2π)d/2 det A Γ G 3 (m) 2 r 3 r/2 1 Br(x,...,x) r! dx b(γ) Aut(Γ) F Γ(l 1,..., l m ), b(γ) = Edges of Γ internal vertices of Γ ; Aut(Γ) is the number of automorphisms of a graph Γ which leave the external vertices fixed; F Γ (l 1,..., l m ) is the Feynman amplitude of the graph Γ computed as follows.

56 Feynman amplitude of a graph For a connected graph Γ, Step 1. put the linear form l i at the i-th external vertex (i = 1,..., m); Step 2. put the polylinear form B r at each internal vertex of degree r (r = 3, 4,... ); Step 3. take contractions of these forms along the edges of Γ using the pairing, A 1. If Γ has several connected components Γ i, F Γ is defined to be the product of F Γi s.

57 Feynman amplitude of a graph For a connected graph Γ, Step 1. put the linear form l i at the i-th external vertex (i = 1,..., m); Step 2. put the polylinear form B r at each internal vertex of degree r (r = 3, 4,... ); Step 3. take contractions of these forms along the edges of Γ using the pairing, A 1. If Γ has several connected components Γ i, F Γ is defined to be the product of F Γi s.

58 Feynman amplitude of a graph For a connected graph Γ, Step 1. put the linear form l i at the i-th external vertex (i = 1,..., m); Step 2. put the polylinear form B r at each internal vertex of degree r (r = 3, 4,... ); Step 3. take contractions of these forms along the edges of Γ using the pairing, A 1. If Γ has several connected components Γ i, F Γ is defined to be the product of F Γi s.

59 Feynman amplitude of a graph For a connected graph Γ, Step 1. put the linear form l i at the i-th external vertex (i = 1,..., m); Step 2. put the polylinear form B r at each internal vertex of degree r (r = 3, 4,... ); Step 3. take contractions of these forms along the edges of Γ using the pairing, A 1. If Γ has several connected components Γ i, F Γ is defined to be the product of F Γi s.

60 Feynman theorem Example Let S(x) = x, Ax (free theory). Then l 1... l N = (2π)d/2 det A = (2π)d/2 det A Γ G 3 (N) pairings b(γ) Aut(Γ) F Γ(l 1,..., l N ) N/2 1 l i, A 1 l σ(i) That s basically the statement of Wick s theorem. i

61 Feynman amplitude of a graph Example Let d = 1, S(x) = x x3. Then x x x x = 2π Γ G 3 (4) 2 1 F Γ(x, x, x, x)

62 Feynman digrams in physics

Feynmann diagrams in a finite dimensional setting

Feynmann diagrams in a finite dimensional setting Feynmann diagrams in a finite dimensional setting Daniel Neiss Uppsala Universitet September 4, 202 Abstract This article aims to explain and justify the use of Feynmann diagrams as a computational tool

More information

Path integrals in quantum mechanics

Path integrals in quantum mechanics Path integrals in quantum mechanics Phys V3500/G8099 handout #1 References: there s a nice discussion of this material in the first chapter of L.S. Schulman, Techniques and applications of path integration.

More information

Lecture 2. Contents. 1 Fermi s Method 2. 2 Lattice Oscillators 3. 3 The Sine-Gordon Equation 8. Wednesday, August 28

Lecture 2. Contents. 1 Fermi s Method 2. 2 Lattice Oscillators 3. 3 The Sine-Gordon Equation 8. Wednesday, August 28 Lecture 2 Wednesday, August 28 Contents 1 Fermi s Method 2 2 Lattice Oscillators 3 3 The Sine-Gordon Equation 8 1 1 Fermi s Method Feynman s Quantum Electrodynamics refers on the first page of the first

More information

Euclidean path integral formalism: from quantum mechanics to quantum field theory

Euclidean path integral formalism: from quantum mechanics to quantum field theory : from quantum mechanics to quantum field theory Dr. Philippe de Forcrand Tutor: Dr. Marco Panero ETH Zürich 30th March, 2009 Introduction Real time Euclidean time Vacuum s expectation values Euclidean

More information

1 The Quantum Anharmonic Oscillator

1 The Quantum Anharmonic Oscillator 1 The Quantum Anharmonic Oscillator Perturbation theory based on Feynman diagrams can be used to calculate observables in Quantum Electrodynamics, like the anomalous magnetic moment of the electron, and

More information

HW 3 Solution Key. Classical Mechanics I HW # 3 Solution Key

HW 3 Solution Key. Classical Mechanics I HW # 3 Solution Key Classical Mechanics I HW # 3 Solution Key HW 3 Solution Key 1. (10 points) Suppose a system has a potential energy: U = A(x 2 R 2 ) 2 where A is a constant and R is me fixed distance from the origin (on

More information

Path Integrals in Quantum Field Theory C6, HT 2014

Path Integrals in Quantum Field Theory C6, HT 2014 Path Integrals in Quantum Field Theory C6, HT 01 Uli Haisch a a Rudolf Peierls Centre for Theoretical Physics University of Oxford OX1 3PN Oxford, United Kingdom Please send corrections to u.haisch1@physics.ox.ac.uk.

More information

Continuity Equations and the Energy-Momentum Tensor

Continuity Equations and the Energy-Momentum Tensor Physics 4 Lecture 8 Continuity Equations and the Energy-Momentum Tensor Lecture 8 Physics 4 Classical Mechanics II October 8th, 007 We have finished the definition of Lagrange density for a generic space-time

More information

Path integrals and the classical approximation 1 D. E. Soper 2 University of Oregon 14 November 2011

Path integrals and the classical approximation 1 D. E. Soper 2 University of Oregon 14 November 2011 Path integrals and the classical approximation D. E. Soper University of Oregon 4 November 0 I offer here some background for Sections.5 and.6 of J. J. Sakurai, Modern Quantum Mechanics. Introduction There

More information

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.

You may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator. MATHEMATICAL TRIPOS Part IB Thursday 7 June 2007 9 to 12 PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question in Section

More information

The Path Integral Formulation of Quantum Mechanics

The Path Integral Formulation of Quantum Mechanics Based on Quantum Mechanics and Path Integrals by Richard P. Feynman and Albert R. Hibbs, and Feynman s Thesis The Path Integral Formulation Vebjørn Gilberg University of Oslo July 14, 2017 Contents 1 Introduction

More information

Lecture-05 Perturbation Theory and Feynman Diagrams

Lecture-05 Perturbation Theory and Feynman Diagrams Lecture-5 Perturbation Theory and Feynman Diagrams U. Robkob, Physics-MUSC SCPY639/428 September 3, 218 From the previous lecture We end up at an expression of the 2-to-2 particle scattering S-matrix S

More information

Classical mechanics of particles and fields

Classical mechanics of particles and fields Classical mechanics of particles and fields D.V. Skryabin Department of Physics, University of Bath PACS numbers: The concise and transparent exposition of many topics covered in this unit can be found

More information

Max Margin-Classifier

Max Margin-Classifier Max Margin-Classifier Oliver Schulte - CMPT 726 Bishop PRML Ch. 7 Outline Maximum Margin Criterion Math Maximizing the Margin Non-Separable Data Kernels and Non-linear Mappings Where does the maximization

More information

Georgia Department of Education Common Core Georgia Performance Standards Framework CCGPS Advanced Algebra Unit 2

Georgia Department of Education Common Core Georgia Performance Standards Framework CCGPS Advanced Algebra Unit 2 Polynomials Patterns Task 1. To get an idea of what polynomial functions look like, we can graph the first through fifth degree polynomials with leading coefficients of 1. For each polynomial function,

More information

Physics 6010, Fall Relevant Sections in Text: Introduction

Physics 6010, Fall Relevant Sections in Text: Introduction Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational

More information

Action Principles in Mechanics, and the Transition to Quantum Mechanics

Action Principles in Mechanics, and the Transition to Quantum Mechanics Physics 5K Lecture 2 - Friday April 13, 2012 Action Principles in Mechanics, and the Transition to Quantum Mechanics Joel Primack Physics Department UCSC This lecture is about how the laws of classical

More information

GRAPH QUANTUM MECHANICS

GRAPH QUANTUM MECHANICS GRAPH QUANTUM MECHANICS PAVEL MNEV Abstract. We discuss the problem of counting paths going along the edges of a graph as a toy model for Feynman s path integral in quantum mechanics. Let Γ be a graph.

More information

LECTURE 5, FRIDAY

LECTURE 5, FRIDAY LECTURE 5, FRIDAY 20.02.04 FRANZ LEMMERMEYER Before we start with the arithmetic of elliptic curves, let us talk a little bit about multiplicities, tangents, and singular points. 1. Tangents How do we

More information

PhD in Theoretical Particle Physics Academic Year 2017/2018

PhD in Theoretical Particle Physics Academic Year 2017/2018 July 10, 017 SISSA Entrance Examination PhD in Theoretical Particle Physics Academic Year 017/018 S olve two among the four problems presented. Problem I Consider a quantum harmonic oscillator in one spatial

More information

INFINITE SEQUENCES AND SERIES

INFINITE SEQUENCES AND SERIES 11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES

More information

Course Summary Math 211

Course Summary Math 211 Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.

More information

221A Lecture Notes Steepest Descent Method

221A Lecture Notes Steepest Descent Method Gamma Function A Lecture Notes Steepest Descent Method The best way to introduce the steepest descent method is to see an example. The Stirling s formula for the behavior of the factorial n! for large

More information

arxiv: v1 [math-ph] 5 Sep 2017

arxiv: v1 [math-ph] 5 Sep 2017 Enumeration of N-rooted maps using quantum field theory K. Gopala Krishna 1, Patrick Labelle, and Vasilisa Shramchenko 3 arxiv:1709.0100v1 [math-ph] 5 Sep 017 Abstract A one-to-one correspondence is proved

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

2.3 Calculus of variations

2.3 Calculus of variations 2.3 Calculus of variations 2.3.1 Euler-Lagrange equation The action functional S[x(t)] = L(x, ẋ, t) dt (2.3.1) which maps a curve x(t) to a number, can be expanded in a Taylor series { S[x(t) + δx(t)]

More information

arxiv: v2 [math-ph] 4 Apr 2018

arxiv: v2 [math-ph] 4 Apr 2018 Enumeration of N-rooted maps using quantum field theory K. Gopala Krishna 1, Patrick Labelle, and Vasilisa Shramchenko 3 arxiv:1709.0100v [math-ph] 4 Apr 018 Abstract A one-to-one correspondence is proved

More information

Opinions on quantum mechanics. CHAPTER 6 Quantum Mechanics II. 6.1: The Schrödinger Wave Equation. Normalization and Probability

Opinions on quantum mechanics. CHAPTER 6 Quantum Mechanics II. 6.1: The Schrödinger Wave Equation. Normalization and Probability CHAPTER 6 Quantum Mechanics II 6.1 The Schrödinger Wave Equation 6. Expectation Values 6.3 Infinite Square-Well Potential 6.4 Finite Square-Well Potential 6.5 Three-Dimensional Infinite- 6.6 Simple Harmonic

More information

Polynomial Review Problems

Polynomial Review Problems Polynomial Review Problems 1. Find polynomial function formulas that could fit each of these graphs. Remember that you will need to determine the value of the leading coefficient. The point (0,-3) is on

More information

Section 3.2 Polynomial Functions and Their Graphs

Section 3.2 Polynomial Functions and Their Graphs Section 3.2 Polynomial Functions and Their Graphs EXAMPLES: P (x) = 3, Q(x) = 4x 7, R(x) = x 2 + x, S(x) = 2x 3 6x 2 10 QUESTION: Which of the following are polynomial functions? (a) f(x) = x 3 + 2x +

More information

BMT 2016 Orthogonal Polynomials 12 March Welcome to the power round! This year s topic is the theory of orthogonal polynomials.

BMT 2016 Orthogonal Polynomials 12 March Welcome to the power round! This year s topic is the theory of orthogonal polynomials. Power Round Welcome to the power round! This year s topic is the theory of orthogonal polynomials. I. You should order your papers with the answer sheet on top, and you should number papers addressing

More information

Functional differentiation

Functional differentiation Functional differentiation March 22, 2016 1 Functions vs. functionals What distinguishes a functional such as the action S [x (t] from a function f (x (t, is that f (x (t is a number for each value of

More information

Optimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30

Optimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained

More information

1.2 Functions and Their Properties Name:

1.2 Functions and Their Properties Name: 1.2 Functions and Their Properties Name: Objectives: Students will be able to represent functions numerically, algebraically, and graphically, determine the domain and range for functions, and analyze

More information

QCD on the lattice - an introduction

QCD on the lattice - an introduction QCD on the lattice - an introduction Mike Peardon School of Mathematics, Trinity College Dublin Currently on sabbatical leave at JLab HUGS 2008 - Jefferson Lab, June 3, 2008 Mike Peardon (TCD) QCD on the

More information

the EL equation for the x coordinate is easily seen to be (exercise)

the EL equation for the x coordinate is easily seen to be (exercise) Physics 6010, Fall 2016 Relevant Sections in Text: 1.3 1.6 Examples After all this formalism it is a good idea to spend some time developing a number of illustrative examples. These examples represent

More information

Local Parametrization and Puiseux Expansion

Local Parametrization and Puiseux Expansion Chapter 9 Local Parametrization and Puiseux Expansion Let us first give an example of what we want to do in this section. Example 9.0.1. Consider the plane algebraic curve C A (C) defined by the equation

More information

AP Calculus Testbank (Chapter 9) (Mr. Surowski)

AP Calculus Testbank (Chapter 9) (Mr. Surowski) AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series

More information

Curve Fitting. 1 Interpolation. 2 Composite Fitting. 1.1 Fitting f(x) 1.2 Hermite interpolation. 2.1 Parabolic and Cubic Splines

Curve Fitting. 1 Interpolation. 2 Composite Fitting. 1.1 Fitting f(x) 1.2 Hermite interpolation. 2.1 Parabolic and Cubic Splines Curve Fitting Why do we want to curve fit? In general, we fit data points to produce a smooth representation of the system whose response generated the data points We do this for a variety of reasons 1

More information

Feynman s path integral approach to quantum physics and its relativistic generalization

Feynman s path integral approach to quantum physics and its relativistic generalization Feynman s path integral approach to quantum physics and its relativistic generalization Jürgen Struckmeier j.struckmeier@gsi.de, www.gsi.de/ struck Vortrag im Rahmen des Winterseminars Aktuelle Probleme

More information

G : Statistical Mechanics

G : Statistical Mechanics G5.651: Statistical Mechanics Notes for Lecture 1 I. DERIVATION OF THE DISCRETIZED PATH INTEGRAL We begin our discussion of the Feynman path integral with the canonical ensemble. The epressions for the

More information

APPM 1350 Exam 2 Fall 2016

APPM 1350 Exam 2 Fall 2016 APPM 1350 Exam 2 Fall 2016 1. (28 pts, 7 pts each) The following four problems are not related. Be sure to simplify your answers. (a) Let f(x) tan 2 (πx). Find f (1/) (5 pts) f (x) 2π tan(πx) sec 2 (πx)

More information

Chapter 4: Interpolation and Approximation. October 28, 2005

Chapter 4: Interpolation and Approximation. October 28, 2005 Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error

More information

14 Renormalization 1: Basics

14 Renormalization 1: Basics 4 Renormalization : Basics 4. Introduction Functional integrals over fields have been mentioned briefly in Part I devoted to path integrals. In brief, time ordering properties and Gaussian properties generalize

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechanics Rajdeep Sensarma! sensarma@theory.tifr.res.in Lecture #22 Path Integrals and QM Recap of Last Class Statistical Mechanics and path integrals in imaginary time Imaginary time

More information

MATH529 Fundamentals of Optimization Constrained Optimization I

MATH529 Fundamentals of Optimization Constrained Optimization I MATH529 Fundamentals of Optimization Constrained Optimization I Marco A. Montes de Oca Mathematical Sciences, University of Delaware, USA 1 / 26 Motivating Example 2 / 26 Motivating Example min cost(b)

More information

1 Assignment 1: Nonlinear dynamics (due September

1 Assignment 1: Nonlinear dynamics (due September Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.

More information

Quantization of Scalar Field

Quantization of Scalar Field Quantization of Scalar Field Wei Wang 2017.10.12 Wei Wang(SJTU) Lectures on QFT 2017.10.12 1 / 41 Contents 1 From classical theory to quantum theory 2 Quantization of real scalar field 3 Quantization of

More information

Section Taylor and Maclaurin Series

Section Taylor and Maclaurin Series Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power

More information

MATH 1231 MATHEMATICS 1B CALCULUS. Section 5: - Power Series and Taylor Series.

MATH 1231 MATHEMATICS 1B CALCULUS. Section 5: - Power Series and Taylor Series. MATH 1231 MATHEMATICS 1B CALCULUS. Section 5: - Power Series and Taylor Series. The objective of this section is to become familiar with the theory and application of power series and Taylor series. By

More information

Quantum Field Theory Notes. Ryan D. Reece

Quantum Field Theory Notes. Ryan D. Reece Quantum Field Theory Notes Ryan D. Reece November 27, 2007 Chapter 1 Preliminaries 1.1 Overview of Special Relativity 1.1.1 Lorentz Boosts Searches in the later part 19th century for the coordinate transformation

More information

Quantum Field theories, Quivers and Word combinatorics

Quantum Field theories, Quivers and Word combinatorics Quantum Field theories, Quivers and Word combinatorics Sanjaye Ramgoolam Queen Mary, University of London 21 October 2015, QMUL-EECS Quivers as Calculators : Counting, correlators and Riemann surfaces,

More information

Combinatorial Enumeration. Jason Z. Gao Carleton University, Ottawa, Canada

Combinatorial Enumeration. Jason Z. Gao Carleton University, Ottawa, Canada Combinatorial Enumeration Jason Z. Gao Carleton University, Ottawa, Canada Counting Combinatorial Structures We are interested in counting combinatorial (discrete) structures of a given size. For example,

More information

Penalty and Barrier Methods. So we again build on our unconstrained algorithms, but in a different way.

Penalty and Barrier Methods. So we again build on our unconstrained algorithms, but in a different way. AMSC 607 / CMSC 878o Advanced Numerical Optimization Fall 2008 UNIT 3: Constrained Optimization PART 3: Penalty and Barrier Methods Dianne P. O Leary c 2008 Reference: N&S Chapter 16 Penalty and Barrier

More information

Formulas that must be memorized:

Formulas that must be memorized: Formulas that must be memorized: Position, Velocity, Acceleration Speed is increasing when v(t) and a(t) have the same signs. Speed is decreasing when v(t) and a(t) have different signs. Section I: Limits

More information

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks. Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the

More information

LSZ reduction for spin-1/2 particles

LSZ reduction for spin-1/2 particles LSZ reduction for spin-1/2 particles based on S-41 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free

More information

REVIEW REVIEW. A guess for a suitable initial state: Similarly, let s consider a final state: Summary of free theory:

REVIEW REVIEW. A guess for a suitable initial state: Similarly, let s consider a final state: Summary of free theory: LSZ reduction for spin-1/2 particles based on S-41 In order to describe scattering experiments we need to construct appropriate initial and final states and calculate scattering amplitude. Summary of free

More information

Lecture 3 Dynamics 29

Lecture 3 Dynamics 29 Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.

More information

Math221: HW# 7 solutions

Math221: HW# 7 solutions Math22: HW# 7 solutions Andy Royston November 7, 25.3.3 let x = e u. Then ln x = u, x2 = e 2u, and dx = e 2u du. Furthermore, when x =, u, and when x =, u =. Hence x 2 ln x) 3 dx = e 2u u 3 e u du) = e

More information

221A Lecture Notes Path Integral

221A Lecture Notes Path Integral A Lecture Notes Path Integral Feynman s Path Integral Formulation Feynman s formulation of quantum mechanics using the so-called path integral is arguably the most elegant. It can be stated in a single

More information

φ 3 theory on the lattice

φ 3 theory on the lattice φ 3 theory on the lattice Michael Kroyter The Open University of Israel SFT 2015 Chengdu 15-May-2015 Work in progress w. Francis Bursa Michael Kroyter (The Open University) φ 3 theory on the lattice SFT

More information

Last Update: March 1 2, 201 0

Last Update: March 1 2, 201 0 M ath 2 0 1 E S 1 W inter 2 0 1 0 Last Update: March 1 2, 201 0 S eries S olutions of Differential Equations Disclaimer: This lecture note tries to provide an alternative approach to the material in Sections

More information

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) 28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.

More information

Calculus 2502A - Advanced Calculus I Fall : Local minima and maxima

Calculus 2502A - Advanced Calculus I Fall : Local minima and maxima Calculus 50A - Advanced Calculus I Fall 014 14.7: Local minima and maxima Martin Frankland November 17, 014 In these notes, we discuss the problem of finding the local minima and maxima of a function.

More information

Separation of Variables in Linear PDE: One-Dimensional Problems

Separation of Variables in Linear PDE: One-Dimensional Problems Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,

More information

Robotics. Islam S. M. Khalil. November 15, German University in Cairo

Robotics. Islam S. M. Khalil. November 15, German University in Cairo Robotics German University in Cairo November 15, 2016 Fundamental concepts In optimal control problems the objective is to determine a function that minimizes a specified functional, i.e., the performance

More information

Announcements. Topics: Homework: - sections , 6.1 (extreme values) * Read these sections and study solved examples in your textbook!

Announcements. Topics: Homework: - sections , 6.1 (extreme values) * Read these sections and study solved examples in your textbook! Announcements Topics: - sections 5.2 5.7, 6.1 (extreme values) * Read these sections and study solved examples in your textbook! Homework: - review lecture notes thoroughly - work on practice problems

More information

Vector Derivatives and the Gradient

Vector Derivatives and the Gradient ECE 275AB Lecture 10 Fall 2008 V1.1 c K. Kreutz-Delgado, UC San Diego p. 1/1 Lecture 10 ECE 275A Vector Derivatives and the Gradient ECE 275AB Lecture 10 Fall 2008 V1.1 c K. Kreutz-Delgado, UC San Diego

More information

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7

x 3y 2z = 6 1.2) 2x 4y 3z = 8 3x + 6y + 8z = 5 x + 3y 2z + 5t = 4 1.5) 2x + 8y z + 9t = 9 3x + 5y 12z + 17t = 7 Linear Algebra and its Applications-Lab 1 1) Use Gaussian elimination to solve the following systems x 1 + x 2 2x 3 + 4x 4 = 5 1.1) 2x 1 + 2x 2 3x 3 + x 4 = 3 3x 1 + 3x 2 4x 3 2x 4 = 1 x + y + 2z = 4 1.4)

More information

Optimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints

Optimization using Calculus. Optimization of Functions of Multiple Variables subject to Equality Constraints Optimization using Calculus Optimization of Functions of Multiple Variables subject to Equality Constraints 1 Objectives Optimization of functions of multiple variables subjected to equality constraints

More information

1) The line has a slope of ) The line passes through (2, 11) and. 6) r(x) = x + 4. From memory match each equation with its graph.

1) The line has a slope of ) The line passes through (2, 11) and. 6) r(x) = x + 4. From memory match each equation with its graph. Review Test 2 Math 1314 Name Write an equation of the line satisfying the given conditions. Write the answer in standard form. 1) The line has a slope of - 2 7 and contains the point (3, 1). Use the point-slope

More information

Lecture 9: Implicit function theorem, constrained extrema and Lagrange multipliers

Lecture 9: Implicit function theorem, constrained extrema and Lagrange multipliers Lecture 9: Implicit function theorem, constrained extrema and Lagrange multipliers Rafikul Alam Department of Mathematics IIT Guwahati What does the Implicit function theorem say? Let F : R 2 R be C 1.

More information

BACKGROUND IN SYMPLECTIC GEOMETRY

BACKGROUND IN SYMPLECTIC GEOMETRY BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations

More information

Generating Functions

Generating Functions Semester 1, 2004 Generating functions Another means of organising enumeration. Two examples we have seen already. Example 1. Binomial coefficients. Let X = {1, 2,..., n} c k = # k-element subsets of X

More information

Computational Physics (6810): Session 8

Computational Physics (6810): Session 8 Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential

More information

x + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the

x + ye z2 + ze y2, y + xe z2 + ze x2, z and where T is the 1.(8pts) Find F ds where F = x + ye z + ze y, y + xe z + ze x, z and where T is the T surface in the pictures. (The two pictures are two views of the same surface.) The boundary of T is the unit circle

More information

Final reading report: BPHZ Theorem

Final reading report: BPHZ Theorem Final reading report: BPHZ Theorem Wang Yang 1 1 Department of Physics, University of California at San Diego, La Jolla, CA 92093 This reports summarizes author s reading on BPHZ theorem, which states

More information

Lagrangian for Central Potentials

Lagrangian for Central Potentials Physics 411 Lecture 2 Lagrangian for Central Potentials Lecture 2 Physics 411 Classical Mechanics II August 29th 2007 Here we will review the Lagrange formulation in preparation for the study of the central

More information

MATH 2250 Final Exam Solutions

MATH 2250 Final Exam Solutions MATH 225 Final Exam Solutions Tuesday, April 29, 28, 6: 8:PM Write your name and ID number at the top of this page. Show all your work. You may refer to one double-sided sheet of notes during the exam

More information

Stochastic Mechanics of Particles and Fields

Stochastic Mechanics of Particles and Fields Stochastic Mechanics of Particles and Fields Edward Nelson Department of Mathematics, Princeton University These slides are posted at http://math.princeton.edu/ nelson/papers/xsmpf.pdf A preliminary draft

More information

The Principle of Least Action

The Principle of Least Action The Principle of Least Action In their never-ending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made

More information

Lecture 7 Unconstrained nonlinear programming

Lecture 7 Unconstrained nonlinear programming Lecture 7 Unconstrained nonlinear programming Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University,

More information

A scaling limit from Euler to Navier-Stokes equations with random perturbation

A scaling limit from Euler to Navier-Stokes equations with random perturbation A scaling limit from Euler to Navier-Stokes equations with random perturbation Franco Flandoli, Scuola Normale Superiore of Pisa Newton Institute, October 208 Newton Institute, October 208 / Subject of

More information

Lecture 7: Sections 2.3 and 2.4 Rational and Exponential Functions. Recall that a power function has the form f(x) = x r where r is a real number.

Lecture 7: Sections 2.3 and 2.4 Rational and Exponential Functions. Recall that a power function has the form f(x) = x r where r is a real number. L7-1 Lecture 7: Sections 2.3 and 2.4 Rational and Exponential Functions Recall that a power function has the form f(x) = x r where r is a real number. f(x) = x 1/2 f(x) = x 1/3 ex. Sketch the graph of

More information

Math 1431 Final Exam Review

Math 1431 Final Exam Review Math 1431 Final Exam Review Comprehensive exam. I recommend you study all past reviews and practice exams as well. Know all rules/formulas. Make a reservation for the final exam. If you miss it, go back

More information

LEAST ACTION PRINCIPLE.

LEAST ACTION PRINCIPLE. LEAST ACTION PRINCIPLE. We continue with the notation from Pollard, section 2.1. Newton s equations for N-body motion are found in Pollard section 2.1, eq. 1.1 there. In section 2.2 they are succinctly

More information

Chapter 3 Numerical Methods

Chapter 3 Numerical Methods Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2

More information

Support Vector Machine (SVM) & Kernel CE-717: Machine Learning Sharif University of Technology. M. Soleymani Fall 2012

Support Vector Machine (SVM) & Kernel CE-717: Machine Learning Sharif University of Technology. M. Soleymani Fall 2012 Support Vector Machine (SVM) & Kernel CE-717: Machine Learning Sharif University of Technology M. Soleymani Fall 2012 Linear classifier Which classifier? x 2 x 1 2 Linear classifier Margin concept x 2

More information

We consider the problem of finding a polynomial that interpolates a given set of values:

We consider the problem of finding a polynomial that interpolates a given set of values: Chapter 5 Interpolation 5. Polynomial Interpolation We consider the problem of finding a polynomial that interpolates a given set of values: x x 0 x... x n y y 0 y... y n where the x i are all distinct.

More information

DYNAMICAL SYSTEMS

DYNAMICAL SYSTEMS 0.42 DYNAMICAL SYSTEMS Week Lecture Notes. What is a dynamical system? Probably the best way to begin this discussion is with arguably a most general and yet least helpful statement: Definition. A dynamical

More information

PAPER 51 ADVANCED QUANTUM FIELD THEORY

PAPER 51 ADVANCED QUANTUM FIELD THEORY MATHEMATICAL TRIPOS Part III Tuesday 5 June 2007 9.00 to 2.00 PAPER 5 ADVANCED QUANTUM FIELD THEORY Attempt THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY

More information

CHAPTER 6 Quantum Mechanics II

CHAPTER 6 Quantum Mechanics II CHAPTER 6 Quantum Mechanics II 6.1 The Schrödinger Wave Equation 6.2 Expectation Values 6.3 Infinite Square-Well Potential 6.4 Finite Square-Well Potential 6.5 Three-Dimensional Infinite-Potential Well

More information

Test 3 Review. y f(a) = f (a)(x a) y = f (a)(x a) + f(a) L(x) = f (a)(x a) + f(a)

Test 3 Review. y f(a) = f (a)(x a) y = f (a)(x a) + f(a) L(x) = f (a)(x a) + f(a) MATH 2250 Calculus I Eric Perkerson Test 3 Review Sections Covered: 3.11, 4.1 4.6. Topics Covered: Linearization, Extreme Values, The Mean Value Theorem, Consequences of the Mean Value Theorem, Concavity

More information

A Short Essay on Variational Calculus

A Short Essay on Variational Calculus A Short Essay on Variational Calculus Keonwook Kang, Chris Weinberger and Wei Cai Department of Mechanical Engineering, Stanford University Stanford, CA 94305-4040 May 3, 2006 Contents 1 Definition of

More information

Honors Advanced Algebra Unit 3: Polynomial Functions November 9, 2016 Task 11: Characteristics of Polynomial Functions

Honors Advanced Algebra Unit 3: Polynomial Functions November 9, 2016 Task 11: Characteristics of Polynomial Functions Honors Advanced Algebra Name Unit 3: Polynomial Functions November 9, 2016 Task 11: Characteristics of Polynomial Functions MGSE9 12.F.IF.7 Graph functions expressed symbolically and show key features

More information

Lecture for Week 2 (Secs. 1.3 and ) Functions and Limits

Lecture for Week 2 (Secs. 1.3 and ) Functions and Limits Lecture for Week 2 (Secs. 1.3 and 2.2 2.3) Functions and Limits 1 First let s review what a function is. (See Sec. 1 of Review and Preview.) The best way to think of a function is as an imaginary machine,

More information

Gaussian processes and Feynman diagrams

Gaussian processes and Feynman diagrams Gaussian processes and Feynman diagrams William G. Faris April 25, 203 Introduction These talks are about expectations of non-linear functions of Gaussian random variables. The first talk presents the

More information

CHAPTER 6 Quantum Mechanics II

CHAPTER 6 Quantum Mechanics II CHAPTER 6 Quantum Mechanics II 6.1 6.2 6.3 6.4 6.5 6.6 6.7 The Schrödinger Wave Equation Expectation Values Infinite Square-Well Potential Finite Square-Well Potential Three-Dimensional Infinite-Potential

More information

Implicit Functions, Curves and Surfaces

Implicit Functions, Curves and Surfaces Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then

More information