Includes only up to I.4. Reorganization pending.

#### I.2

##### I.2.1

We need to use as the new Hamiltonian, and we can start by replacing it in :

.

In general and don’t commute, but as is an infinitesimal, we can assume that they **do** commute:

(using (21))

Applying the same operations as the text we get

and, when going to the continuum limit,

.

##### I.2.2

Let’s start by checking (18):

(17)

If we apply *n* times to the left side we get

.

And, if we do the same thing to the right side, we obtain

.

Dividing this result by (17), we get:

,

that matches with (18).

Now we can start with (24). Let’s copy equation (22):

.

If we derive the left side with respect to and we get:

,

and if we do the same with the right side we are left with

If we not set to zero, we can find the following equality:

.

This can turn quite messy for higher derivatives, but we can use the series expansion of the exponential to save a lot of work: if we do derivatives, only the terms with variables can survive after we set to zero. Taking this into account, let’s try to apply this to a fourth derivative:

.

Zee asks to do a sixth derivative:

We can concentrate in the first term, as the second term is the same that was previously done:

Combining this with the previous result:

Ordering the terms:

#### I.3

##### I.3.1

The propagator is

.

As the integral is Lorentz invariant, we can choose , :

.

If we express the 3-momentum volume element as

,

we get

.

Integrating out and reordering factors:

.

Integrating over :

.

Completing the integral, following Michael Kachelrieß and using :

From Mathematica:

From Wikipedia: is exponentially decaying.

##### I.3.2

The propagator must be given by

.

Setting to zero, we get:

,

and, using ,

.

As in the case in the book, we have two poles: one at and another at . But in this case we don’t have the exponential divergence problem, so we can choose any contour. Let’s try both to check that they give the same result.

Counterclockwise, upper half plane, simple pole at :

.

Clockwise, lower half plane, simple pole at :

.

Yes, they match. Adding this result to the bigger integral:

.

Using Wolfram Alpha we get

,

where is a modified Bessel function that decays exponentially.

##### I.3.3

We need to show that the advanced propagator,

,

only propagates to the future. Let’s start as Zee does, integrating over by the method of contours and, first of all , getting the position of the poles:

Using and in an open sense and defining , we get:

.

If , we have and, as is small, we can linearize the function around the old pole:

.

Replacing into the original expression and taking into account that is a generical term for “small positive value”:

,

so the two poles are in the upper half of the complex plane.

When closing the contour we need to be sure that the integral over the contour converges, so the contour must be closed over the side that makes the exponential small (because an exponential divergence would overwhelm the decay). As the integral is

,

positive values of give us the requirement of values of with positive imaginary component to ensure an exponent with a negative real component. So in this case we would catch the poles inside the contour and get a nonzero result. In the other case, m the contour would have to be closed in the lower half plane, missing the poles and getting a zero result.

The retarded propagator is symmetrical, only switching the position of the poles.

#### I.4

##### W(J(k))

Let’s start with a derivation inside the chapter:

Moving all the differentials together, removing the multiple integral signs and reordering the exponentials:

.

Doing the last integration (observe that a factor is given by the integration of the exponential):

.

Combining the exponentials and integrating:

##### Energy

Now let’s do the energy integral without consulting the appendix 😀

.

Let’s define and do the substitution . By spherical symmetry, we can use as the “north pole”:

(adds an odd function)

Closing the contour (it’s not so easy to prove that the complex section doesn’t contribute but it’s true) through the upper half plane (a positive imaginary part of gives a negative exponent) we pick up the pole at :

.

Replacing this value for the integral:

##### I.4.1

We can start with the 2-dimensional form of (6):

Expressing the angular part of the integral as a Bessel Function:

Putting this result in the previous expression:

Integrating:

,

where is a modified Bessel function of the second kind and . If the energy would seem constant, but agreement with electrostatic theory would seem to require a logarithmic potential and the integral certainly is divergent for . Let’s use a reference at to avoid this problems:

Using Mathematica:

Let’s analyze now the *D*-dimensional form of (6):

We can express the denominator as an integral over an auxiliary variable:

Setting aside the dependent terms:

Substituting :

Doing the *n*-dimensional gaussian integral:

Using Wolfram Alpha:

Putting the result of the integral in the expression for the energy:

Checking for values of D with known results:

Using Wolfram Alpha:

Hi, I want to know how to use Mathematica to calculate the last integration in 1.3.1. Can you teach me how? Many thanks!

You can just evaluate the integral, though you will get a cleaner result if you restrict the domain of and :