phy2403

Momentum of scalar field: continued.

November 9, 2018 phy2403 , , , ,

[Click here for a PDF of this post with nicer formatting]

Conservation of the field momentum.

This is a follow up to the unanswered questions I had yesterday related to the apparent time dependent terms in the previous expansion of \( P^i \) for a scalar field.

It turns out that examining the reasons that we can say that the field momentum is conserved also sheds some light on the question. \( P^i \) is not an a-priori conserved quantity, but we may use the charge conservation argument to justify this despite it not having a four-vector nature (i.e. with zero four divergence.)

The momentum \( P^i \) that we have defined is related to the conserved quantity \( T^{0\mu} \), the energy-momentum tensor, which satisfies \( 0 = \partial_\mu T^{0\mu} \) by Noether’s theorem (this was the conserved quantity associated with a spacetime translation.)

That tensor was
\begin{equation}\label{eqn:momentum:120}
T^{\mu\nu} = \partial^\mu \phi \partial^\nu \phi – g^{\mu\nu} \LL,
\end{equation}
and can be used to define the momenta
\begin{equation}\label{eqn:momentum:140}
\begin{aligned}
\int d^3 x T^{0k}
&= \int d^3 x \partial^0 \phi \partial^k \phi \\
&= \int d^3 x \pi \partial^k \phi.
\end{aligned}
\end{equation}
Charge \( Q^i = \int d^3 x j^0 \) was conserved with respect to a limiting surface argument, and we can make a similar “beer can integral” argument for \( P^i \), integrating over a large time interval \( t \in [-T, T] \) as sketched in fig. 1. That is
\begin{equation}\label{eqn:momentum:160}
\begin{aligned}
0
&=
\partial_\mu \int d^4 x T^{0\mu} \\
&=
\partial_0 \int d^4 x T^{00}
+
\partial_k \int d^4 x T^{0k} \\
&=
\partial_0 \int_{-T}^T dt \int d^3 x T^{00}
+
\partial_k \int_{-T}^T dt \int d^3 x T^{0k} \\
&=
\partial_0 \int_{-T}^T dt \int d^3 x T^{00}
+
\partial_k \int_{-T}^T dt
\inv{2} \int \frac{d^3 p }{(2 \pi)^3} p^k
\lr{
a_\Bp^\dagger a_\Bp
+ a_\Bp a_\Bp^\dagger
– a_\Bp a_{-\Bp} e^{- 2 i \omega_\Bp t}
– a_\Bp^\dagger a_{-\Bp}^\dagger e^{2 i \omega_\Bp t}
} \\
&=
\int d^3 x \evalrange{T^{00}}{-T}{T}
+
T \partial_k
\int \frac{d^3 p }{(2 \pi)^3} p^k
\lr{
a_\Bp^\dagger a_\Bp
+ a_\Bp a_\Bp^\dagger
}
-\inv{2}
\partial_k \int_{-T}^T dt
\int \frac{d^3 p }{(2 \pi)^3} p^k
\lr{
a_\Bp a_{-\Bp} e^{- 2 i \omega_\Bp t}
+ a_\Bp^\dagger a_{-\Bp}^\dagger e^{2 i \omega_\Bp t}
}.
\end{aligned}
\end{equation}

fig. 1. Cylindrical spacetime boundary.

The first integral can be said to vanish if the field energy goes to zero at the time boundaries, and the last integral reduces to
\begin{equation}\label{eqn:momentum:180}
\begin{aligned}
-\inv{2}
\partial_k \int_{-T}^T dt
\int \frac{d^3 p }{(2 \pi)^3} p^k
\lr{
a_\Bp a_{-\Bp} e^{- 2 i \omega_\Bp t}
+ a_\Bp^\dagger a_{-\Bp}^\dagger e^{2 i \omega_\Bp t}
}
&=
-\int \frac{d^3 p }{2 (2 \pi)^3} p^k
\lr{
a_\Bp a_{-\Bp} \frac{\sin( -2 \omega_\Bp T )}{-2 \omega_\Bp}
+ a_\Bp^\dagger a_{-\Bp}^\dagger \frac{\sin( 2 \omega_\Bp T )}{2 \omega_\Bp}
} \\
&=
-\int \frac{d^3 p }{2 (2 \pi)^3} p^k
\lr{
a_\Bp a_{-\Bp} + a_\Bp^\dagger a_{-\Bp}^\dagger
}
\frac{\sin( 2 \omega_\Bp T )}{2 \omega_\Bp}
.
\end{aligned}
\end{equation}
The \( \sin \) term can be interpretted as a sinc like function of \( \omega_\Bp \) which vanishes for large \( \Bp \). It’s not entirely sinc like for a massive field as \( \omega_\Bp = \sqrt{ \Bp^2 + m^2 } \), which never hits zero, as shown in fig. 2.

fig 2. sin(2 omega T)/omega

Vanishing for large \( \Bp \) doesn’t help the whole integral vanish, but we can resort to the Riemann-Lebesque lemma [1] instead and interpret this integral as one with a plain old high frequency oscillation that is presumed to vanish (i.e. the rest is well behaved enough that it can be labelled as \( L_1 \) integrable.)

We see that only the non-time dependent portion of \( \mathbf{P} \) matters from a conserved quantity point of view, and having killed off all the time dependent terms, we are left with a conservation relationship for the momenta \( \spacegrad \cdot \BP = 0 \), where \( \BP \) in normal order is just
\begin{equation}\label{eqn:momentum:200}
: \BP : = \int \frac{d^3 p}{(2 \pi)^3} \Bp a_\Bp^\dagger a_\Bp.
\end{equation}

References

[1] Wikipedia contributors. Riemann-lebesgue lemma — Wikipedia, the free encyclopedia, 2018. URL https://en.wikipedia.org/w/index.php?title=Riemann%E2%80%93Lebesgue_lemma&oldid=856778941. [Online; accessed 29-October-2018].

PHY2403: QFT I: Momentum of scalar field: time dependent terms?

November 8, 2018 phy2403 , , ,

[Click here for a PDF of this post with nicer formatting]

Way back in lecture 8, it was claimed that
\begin{equation}\label{eqn:momentum:20}
P^k = \int d^3 x \hat{\pi} \partial^k \hat{\phi} = \int \frac{d^3 p}{(2\pi)^3} p^k a_\Bp^\dagger a_\Bp.
\end{equation}
If I compute this, I get a normal ordered variation of this operator, but also get some time dependent terms. Here’s the computation (dropping hats)
\begin{equation}\label{eqn:momentum:40}
\begin{aligned}
P^k
&= \int d^3 x \hat{\pi} \partial^k \phi \\
&= \int d^3 x \partial_0 \phi \partial^k \phi \\
&= \int d^3 x \frac{d^3 p d^3 q}{(2 \pi)^6} \inv{\sqrt{2 \omega_p 2 \omega_q} }
\partial_0
\lr{
a_\Bp e^{-i p \cdot x}
+
a_\Bp^\dagger e^{i p \cdot x}
}
\partial^k
\lr{
a_\Bq e^{-i q \cdot x}
+
a_\Bq^\dagger e^{i q \cdot x}
}.
\end{aligned}
\end{equation}
The exponential derivatives are
\begin{equation}\label{eqn:momentum:60}
\begin{aligned}
\partial_0 e^{\pm i p \cdot x}
&=
\partial_0 e^{\pm i p_\mu x^\mu} \\
&=
\pm i p_0
\partial_0 e^{\pm i p \cdot x},
\end{aligned}
\end{equation}
and
\begin{equation}\label{eqn:momentum:80}
\begin{aligned}
\partial^k e^{\pm i p \cdot x}
&=
\partial^k e^{\pm i p^\mu x_\mu} \\
&=
\pm i p^k e^{\pm i p \cdot x},
\end{aligned}
\end{equation}
so
\begin{equation}\label{eqn:momentum:100}
\begin{aligned}
P^k
&=
-\int d^3 x \frac{d^3 p d^3 q}{(2 \pi)^6} \inv{\sqrt{2 \omega_p 2 \omega_q} }
p_0 q^k
\lr{
-a_\Bp e^{-i p \cdot x}
+
a_\Bp^\dagger e^{i p \cdot x}
}
\lr{
-a_\Bq e^{-i q \cdot x}
+
a_\Bq^\dagger e^{i q \cdot x}
} \\
&=
-\inv{2} \int d^3 x \frac{d^3 p d^3 q}{(2 \pi)^6} \sqrt{\frac{\omega_p}{\omega_q}} q^k
\lr{
a_\Bp a_\Bq e^{-i (p + q) \cdot x}
+ a_\Bp^\dagger a_\Bq^\dagger e^{i (p + q) \cdot x}
– a_\Bp a_\Bq^\dagger e^{i (q – p) \cdot x}
– a_\Bp^\dagger a_\Bq e^{i (p – q) \cdot x}
} \\
&=
\inv{2} \int \frac{d^3 p d^3 q}{(2 \pi)^3} \sqrt{\frac{\omega_p}{\omega_q}} q^k
\lr{
– a_\Bp a_\Bq e^{- i(\omega_\Bp + \omega_\Bq) t} \delta^3(\Bp + \Bq)
– a_\Bp^\dagger a_\Bq^\dagger e^{i(\omega_\Bp + \omega_\Bq) t} \delta^3(-\Bp – \Bq)
+ a_\Bp a_\Bq^\dagger e^{i(\omega_\Bq – \omega_\Bp) t} \delta^3(\Bp – \Bq)
+ a_\Bp^\dagger a_\Bq e^{i(\omega_\Bp – \omega_\Bq) t} \delta^3(\Bq – \Bp)
} \\
&=
\inv{2} \int \frac{d^3 p }{(2 \pi)^3} p^k
\lr{
a_\Bp^\dagger a_\Bp
+ a_\Bp a_\Bp^\dagger
– a_\Bp a_{-\Bp} e^{- 2 i \omega_\Bp t}
– a_\Bp^\dagger a_{-\Bp}^\dagger e^{2 i \omega_\Bp t}
}.
\end{aligned}
\end{equation}

What is the rationale for ignoring those time dependent terms? Does normal ordering also implicitly drop any non-paired creation/annihilation operators? If so, why?

PHY2403H Quantum Field Theory. Lecture 15b: Wick’s theorem, vacuum expectation, Feynman diagrams, phi^4 interaction, tree level diagrams, scattering, cross section, differential cross section. Taught by Prof. Erich Poppitz

November 7, 2018 phy2403 , , , , , , ,

[Here are my notes for the second half of last wednesday’s lecture]. Because the simplewick latex package is required to format these meaningfully, and that’s not available in wordpress latex, I’m not going to attempt to make a web viewable version in addition to the pdf.

PHY2403 (QFT I). Pondering the ground state bra formula.

November 5, 2018 phy2403 , ,

[Click here for a PDF of this post with nicer formatting]

In lecture 14 we found

\begin{equation}\label{eqn:qftLecture15:60}
\ket{\Omega}
=
\evalbar{
\frac{ U(t_0, -T) \ket{0} }
{
e^{-i E_0(T – t_0)} \braket{\Omega}{0}
}
}{T \rightarrow \infty(1 – i \epsilon)},
\end{equation}
and it was stated that we can also show that
\begin{equation}\label{eqn:qftLecture15:80}
\bra{\Omega}
=
\evalbar{
\frac{ \bra{0} U(T, t_0) }
{
e^{-i E_0(T – t_0)} \braket{0}{\Omega}
}
}{T \rightarrow \infty(1 – i \epsilon)}.
\end{equation}

This second statement is actually not obvious since
\begin{equation}\label{eqn:braOmega:100}
\lr{
\frac{ U(t_0, -T) \ket{0} }
{
e^{-i E_0(T – t_0)} \braket{\Omega}{0}
}
}^\dagger
=
\frac{ \bra{0} U(-T, t_0) \ket{0} }
{
e^{+i E_0(T – t_0)} \braket{0}{\Omega}
}
\end{equation}

My first thought was that I’d written down \ref{eqn:qftLecture15:80} in my notes wrong, but this is actually also consistent with [1], which our Prof is following loosely (i.e. he is explicitly filling in many of the holes in that dense little 800 page book).

The resolution of this inconsistency is that the limit point \( \infty(1 – i \epsilon) \) doesn’t work if you just conjugate, and you’d also have to conjugate that limit, so while
\ref{eqn:braOmega:100} is correct, it is only part of the story, and should really be stated as
\begin{equation}\label{eqn:braOmega:120}
\bra{\Omega}
=
\evalbar{
\frac{ \bra{0} U(-T, t_0) \ket{0} }
{
e^{+i E_0(T – t_0)} \braket{0}{\Omega}
}
}{T \rightarrow \infty(1 + i \epsilon)}.
\end{equation}
This is awkward because now our expressions for \( \bra{\Omega} \) and \( \ket{\Omega} \) approach \( T \) from different directions, and we want to evaluate both with a single limiting argument.

To resolve this, we really have to start back with the identity expansion we used in lecture 14
\begin{equation}\label{eqn:braOmega:140}
\begin{aligned}
\bra{0} e^{-i H T}
&=
\lr{
\braket{0}{\Omega}\bra{\Omega}
+ {\int\kern-1em\sum}_n \braket{0}{n} \bra{n}
}
e^{-i H T} \\
&=
\braket{0}{\Omega}\bra{\Omega}
e^{-i E_0 T}
+ {\int\kern-1em\sum}_n \braket{0}{n} \bra{n} e^{-i E_n T}.
\end{aligned}
\end{equation}
We argued (as does the text) that approaching to as \( T( 1 – i \epsilon) \) kills off the energetic states since
\begin{equation}\label{eqn:braOmega:160}
\bra{n} e^{-i E_n T}
\rightarrow
\bra{n} e^{-i E_n T} e^{-E_n T \epsilon}
\end{equation}
and the exponential damping factor is smaller for each \( E_n > E_0 \), so it can be neglected in the large \( T \) limit, leaving
\begin{equation}\label{eqn:braOmega:180}
\bra{0} e^{-i H T}
=
\lim_{T \rightarrow \infty(1 – i \epsilon)}
\braket{0}{\Omega}\bra{\Omega}.
\end{equation}
As we did for \( \ket{\Omega} \) we can shift the large time \( T \) by a small constant (this time \( -t_0 \) instead of \( t_0 \)), to give
\begin{equation}\label{eqn:braOmega:200}
\begin{aligned}
\bra{\Omega}
&=
\lim_{T \rightarrow \infty(1 – i \epsilon)}
\frac{ \bra{0} e^{-i H T} }
{
\braket{0}{\Omega} e^{-i E_0 T}
} \\
&\approx
\lim_{T \rightarrow \infty(1 – i \epsilon)}
\frac{ \bra{0} e^{-i H (T – t_0)} }
{
\braket{0}{\Omega} e^{-i E_0 (T – t_0)}
} \\
&=
\lim_{T \rightarrow \infty(1 – i \epsilon)}
\frac{ \bra{0} e^{i H_0( T – t_0)} e^{-i H (T – t_0)} }
{
\braket{0}{\Omega} e^{-i E_0 (T – t_0)}
} \\
&=
\lim_{T \rightarrow \infty(1 – i \epsilon)}
\frac{ \bra{0} U(T, t_0) }
{
\braket{0}{\Omega} e^{-i E_0 (T – t_0)}
},
\end{aligned}
\end{equation}
where the projective property \( \bra{0} e^{i H_0 \alpha} = \bra{0} \) has been used to insert a no-op (i.e. \( \bra{0} H_0 = 0 \)). This recovers the result stated in class (also: eq. (4.29) in the text.)

References

[1] Michael E Peskin and Daniel V Schroeder. An introduction to Quantum Field Theory. Westview, 1995.

Analysis of incorrect problem solution attempt (non-relativistic QFT massive scalar field with quadratic interaction)

November 1, 2018 phy2403 , ,

[Click here for a PDF of this post with nicer formatting]

An incorrect start at a field theory problem.

The following is a problem from the second PHY2403 problem set, and the start of a solution attempt I made. I’m not posting this because it shows how to do the problem, but because it was a useful problem to show that I didn’t understand a lot of the problem statement.
Understanding where I went wrong is actually pretty useful in this case.

Problem: Playing with the non-relativistic limit

Consider a real scalar relativistic field theory of mass m with \( \lambda \phi^4 \) interaction. Let there be \( N \) particles of momenta labeled by \( p_1,\cdots, p_N\), whose energies are such that they are insufficient to create any new particles. Nevertheless, the particles can scatter and exchange momenta. In what follows you will study this N-particle nonrelativistic limit in some detail.

  1. Write down the Hamiltonian of the field theory, including the interaction term, restricted to the N-particle sector of Hilbert space. (Use the creation and annihilation operator representation, i.e. write the result as sums of products of creation and annihilation operators of particles of various momenta.)
  2. Does the resulting Hamiltonian preserve particle number? Is there an associated symmetry? What is the operator that generates it?
  3. Consider now the interaction term in your reduced (to the N-particle sector of Hilbert space) Hamiltonian. How does a typical interaction term (for given configurations of momenta) act on an N-particle state? What kinds of scattering processes does it describe?
  4. What do you think is the potential, in x-space, that allows the various particles to scatter and exchange momentum? How would you describe the resulting nonrelativistic quantum system to friends who never took QFT but are well-versed in quantum mechanics?

Incorrect answer attempt.

The Lagrangian density of a massive scalar field with a \( \lambda \phi^4 \) interaction has the form
\begin{equation}\label{eqn:Nparticle:20}
\LL = \inv{2} \partial_\mu \phi \partial^\mu \phi – \inv{2} m^2 \phi^2 – \lambda \phi^4.
\end{equation}
The corresponding Hamiltonian is
\begin{equation}\label{eqn:Nparticle:40}
H = \inv{2} \int d^3x \lr{ \pi^2 + \frac{m^2}{2} (\spacegrad \phi)^2 + m^2 \phi^2 } + \lambda \int d^3 x \phi^4.
\end{equation}
In terms of creation and annihilation operators, we know the form of the non-interaction portion of the Hamiltonian, which in normal order is
\begin{equation}\label{eqn:Nparticle:60}
H_0 = \int \frac{d^3 p}{(2 \pi)^3} \omega_\Bp a_\Bp^\dagger a_\Bp,
\end{equation}
but the interaction contribution is much messier
\begin{equation}\label{eqn:Nparticle:80}
\begin{aligned}
H_{\text{int}}
&=
\lambda \int d^3 x \frac{ d^3 p d^3 k d^3 q d^3 s}{4 (2 \pi)^{12} \sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_\Bs
} }
\lr{ a_\Bp e^{-i p \cdot x} + a_\Bp e^{i p \cdot x} }
\lr{ a_\Bk e^{-i k \cdot x} + a_\Bk e^{i k \cdot x} }
\lr{ a_\Bq e^{-i q \cdot x} + a_\Bq e^{i q \cdot x} }
\lr{ a_\Bs e^{-i s \cdot x} + a_\Bs e^{i s \cdot x} } \\
&=
\lambda \int d^3 x \frac{ d^3 p d^3 k d^3 q d^3 s}{4 (2 \pi)^{12} \sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_\Bs
} }
\lr{ a_\Bp e^{-i \omega_\Bp t + i \Bp \cdot \Bx} + a_\Bp e^{i \omega_\Bp t – i \Bp \cdot \Bx} }
\lr{ a_\Bk e^{-i \omega_\Bk t + i \Bk \cdot \Bx} + a_\Bk e^{i \omega_\Bk t – i \Bk \cdot \Bx} } \\
&\quad \lr{ a_\Bq e^{-i \omega_\Bq t + i \Bq \cdot \Bx} + a_\Bq e^{i \omega_\Bq t – i \Bq \cdot \Bx} }
\lr{ a_\Bs e^{-i \omega_\Bs t + i \Bs \cdot \Bx} + a_\Bs e^{i \omega_\Bs t – i \Bs \cdot \Bx} } \\
&=
\lambda \int \frac{ d^3 p d^3 k d^3 q d^3 s}{4 (2 \pi)^{9} \sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_\Bs
} }
\Big(
a_\Bp a_\Bk a_\Bq a_\Bs e^{-i (\omega_\Bp + \omega_\Bk + \omega_\Bq + \omega_\Bs)t} \delta^3( \Bp + \Bk + \Bq + \Bs ) \\
&\qquad +
a_\Bp a_\Bk a_\Bq a_\Bs^\dagger e^{-i (\omega_\Bp + \omega_\Bk + \omega_\Bq – \omega_\Bs)t} \delta^3( \Bp + \Bk + \Bq – \Bs ) \\
&\qquad \\
&\qquad + \cdots \\
&\qquad +
a_\Bp^\dagger a_\Bk^\dagger a_\Bq^\dagger a_\Bs^\dagger e^{-i (-\omega_\Bp – \omega_\Bk – \omega_\Bq – \omega_\Bs)t} \delta^3( -\Bp – \Bk – \Bq – \Bs )
\Big) \\
&=
\lambda \int \frac{ d^3 p d^3 k d^3 q }{4 (2 \pi)^{9}
}
\Big(
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{-\Bp – \Bk – \Bq}
}}
a_\Bp a_\Bk a_\Bq a_{-\Bp -\Bk – \Bq} e^{-i (\omega_\Bp + \omega_\Bk + \omega_\Bq + \omega_{-\Bp -\Bk -\Bq})t} \\
&\qquad +
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{\Bp + \Bk + \Bq}
}}
a_\Bp a_\Bk a_\Bq a_{\Bp + \Bk + \Bq}^\dagger e^{-i (\omega_\Bp + \omega_\Bk + \omega_\Bq – \omega_{\Bp + \Bk + \Bq})t} \\
&\qquad +
\cdots \\
&\qquad +
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{-\Bp – \Bk – \Bq}
}}
a_\Bp^\dagger a_\Bk^\dagger a_\Bq^\dagger a_{-\Bp -\Bk -\Bq}^\dagger e^{-i (-\omega_\Bp – \omega_\Bk – \omega_\Bq – \omega_{-\Bp -\Bk -\Bq})t}
\Big)
\end{aligned}
\end{equation}
Assuming we can normal order these terms as in \( H_0 \), we can rewrite the interaction as
\begin{equation}\label{eqn:Nparticle:100}
\begin{aligned}
H_{\text{int}}
&=
\lambda \int \frac{ d^3 p d^3 k d^3 q }{4 (2 \pi)^{9} }
\Big(
\binom{4}{0}
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{-\Bp – \Bk – \Bq}
}}
a_\Bp a_\Bk a_\Bq a_{-\Bp -\Bk – \Bq} e^{-i (\omega_\Bp + \omega_\Bk + \omega_\Bq + \omega_{-\Bp -\Bk -\Bq})t} \\
&\qquad +
\binom{4}{1}
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{\Bp – \Bk – \Bq}
}}
a_\Bp^\dagger a_\Bk a_\Bq a_{\Bp – \Bk – \Bq} e^{-i (-\omega_\Bp + \omega_\Bk + \omega_\Bq + \omega_{\Bp – \Bk – \Bq})t} \\
&\qquad +
\binom{4}{2}
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{\Bp + \Bk – \Bq}
}}
a_\Bp^\dagger a_\Bk^\dagger a_\Bq a_{\Bp + \Bk – \Bq} e^{-i (-\omega_\Bp – \omega_\Bk + \omega_\Bq + \omega_{\Bp + \Bk – \Bq})t} \\
&\qquad +
\binom{4}{3}
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{\Bp + \Bk _ \Bq}
}}
a_\Bp^\dagger a_\Bk^\dagger a_\Bq^\dagger a_{\Bp + \Bk + \Bq} e^{-i (-\omega_\Bp – \omega_\Bk – \omega_\Bq + \omega_{\Bp + \Bk + \Bq})t} \\
&\qquad +
\binom{4}{4}
\inv{\sqrt{
\omega_\Bp \omega_\Bk \omega_\Bq \omega_{-\Bp – \Bk – \Bq}
}}
a_\Bp^\dagger a_\Bk^\dagger a_\Bq^\dagger a_{-\Bp -\Bk -\Bq}^\dagger e^{-i (-\omega_\Bp – \omega_\Bk – \omega_\Bq – \omega_{-\Bp -\Bk -\Bq})t}
\Big)
\end{aligned}
\end{equation}

If we restrict the allowed momenta to the discrete set \( \Bp \in \setlr{ \Bp_1, \Bp_2, \cdots \Bp_N} \), the total Hamiltonian including the interaction term
takes the form
\begin{equation}\label{eqn:Nparticle:120}
\begin{aligned}
\text{\(:H:\)} &=
\sum_{i = 1}^N \omega_{\Bp_i} a_{\Bp_i}^\dagger a_{\Bp_i}
+
\frac{
\lambda
}{4 }
\sum_{j,m,n = 1}^N
\Big(
\binom{4}{0}
\inv{\sqrt{
\omega_{\Bp_j} \omega_{\Bp_m} \omega_{\Bp_n} \omega_{-{\Bp_j} – {\Bp_m} – {\Bp_n}}
}}
a_{\Bp_j} a_{\Bp_m} a_{\Bp_n} a_{-\Bp -\Bk – \Bq} e^{-i (\omega_{\Bp_j} + \omega_{\Bp_m} + \omega_{\Bp_n} + \omega_{-{\Bp_j} -{\Bp_m} -{\Bp_n}})t} \\
&\qquad +
\binom{4}{1}
\inv{\sqrt{
\omega_{\Bp_j} \omega_{\Bp_m} \omega_{\Bp_n} \omega_{{\Bp_j} – {\Bp_m} – {\Bp_n}}
}}
a_{\Bp_j}^\dagger a_{\Bp_m} a_{\Bp_n} a_{{\Bp_j} – {\Bp_m} – {\Bp_n}} e^{-i (-\omega_{\Bp_j} + \omega_{\Bp_m} + \omega_{\Bp_n} + \omega_{{\Bp_j} – {\Bp_m} – {\Bp_n}})t} \\
&\qquad +
\binom{4}{2}
\inv{\sqrt{
\omega_{\Bp_j} \omega_{\Bp_m} \omega_{\Bp_n} \omega_{{\Bp_j} + {\Bp_m} – {\Bp_n}}
}}
a_{\Bp_j}^\dagger a_{\Bp_m}^\dagger a_{\Bp_n} a_{{\Bp_j} + {\Bp_m} – {\Bp_n}} e^{-i (-\omega_{\Bp_j} – \omega_{\Bp_m} + \omega_{\Bp_n} + \omega_{{\Bp_j} + {\Bp_m} – {\Bp_n}})t} \\
&\qquad +
\binom{4}{3}
\inv{\sqrt{
\omega_{\Bp_j} \omega_{\Bp_m} \omega_{\Bp_n} \omega_{{\Bp_j} + {\Bp_m} – {\Bp_n}}
}}
a_{\Bp_j}^\dagger a_{\Bp_m}^\dagger a_{\Bp_n}^\dagger a_{{\Bp_j} + {\Bp_m} + {\Bp_n}} e^{-i (-\omega_{\Bp_j} – \omega_{\Bp_m} – \omega_{\Bp_n} + \omega_{{\Bp_j} + {\Bp_m} + {\Bp_n}})t} \\
&\qquad +
\binom{4}{4}
\inv{\sqrt{
\omega_{\Bp_j} \omega_{\Bp_m} \omega_{\Bp_n} \omega_{-{\Bp_j} – {\Bp_m} – {\Bp_n}}
}}
a_{\Bp_j}^\dagger a_{\Bp_m}^\dagger a_{\Bp_n}^\dagger a_{-{\Bp_j} -{\Bp_m} -{\Bp_n}}^\dagger e^{-i (-\omega_{\Bp_j} – \omega_{\Bp_m} – \omega_{\Bp_n} – \omega_{-{\Bp_j} -{\Bp_m} -{\Bp_n}})t}
\Big)
\end{aligned}
\end{equation}

When we did the same sort of calculation for \( (\spacegrad \phi)^2 + m^2 \phi^2 \) all the time dependent terms cancelled nicely, but that isn’t obviously the case here.
However, we haven’t used the non-relativistic (low energy) constraint. That constraint can be expressed as \( \Bp^2 \ll m^2 \), in which case \( \omega_\Bp = \sqrt{ \Bp^2 + m^2 } \sim m \), the mass of each of the particles. Incorporating that into our N-particle Hamiltonian, we have
\begin{equation}\label{eqn:Nparticle:140}
\begin{aligned}
\text{\(:H:\)} &=
\sum_{i = 1}^N \omega_{\Bp_i} a_{\Bp_i}^\dagger a_{\Bp_i}
+
\frac{
\lambda
}{4 m^2 }
\sum_{j,m,n = 1}^N
\Big(
\binom{4}{0}
a_{\Bp_j} a_{\Bp_m} a_{\Bp_n} a_{-\Bp -\Bk – \Bq} e^{- 4 i m t} \\
&\qquad +
\binom{4}{1}
a_{\Bp_j}^\dagger a_{\Bp_m} a_{\Bp_n} a_{{\Bp_j} – {\Bp_m} – {\Bp_n}} e^{-3 i m t} \\
&\qquad +
\binom{4}{2}
a_{\Bp_j}^\dagger a_{\Bp_m}^\dagger a_{\Bp_n} a_{{\Bp_j} + {\Bp_m} – {\Bp_n}} \\
&\qquad +
\binom{4}{3}
a_{\Bp_j}^\dagger a_{\Bp_m}^\dagger a_{\Bp_n}^\dagger a_{{\Bp_j} + {\Bp_m} + {\Bp_n}} e^{ 3 i m t } \\
&\qquad +
\binom{4}{4}
a_{\Bp_j}^\dagger a_{\Bp_m}^\dagger a_{\Bp_n}^\dagger a_{-{\Bp_j} -{\Bp_m} -{\Bp_n}}^\dagger e^{4 i m t}
\Big).
\end{aligned}
\end{equation}

The only annoying aspect to this Hamiltonian is the \( a_{{\Bp_j} + {\Bp_m} – {\Bp_n}} \) operator in the interaction term, which is not clear to me how to interpret. That seems to imply that it is possible to create particles with linear combinations of momentum that may not be in the original set of \( N \) particle momenta. I think that this can be further fudged by invoking the non-relativistic constraint again, and decreeing that each of the uniquely indexed creation and annihilation operators are distinguishable only by index, so we can write the N-particle non-relativistic sector Hamiltonian as
\begin{equation}\label{eqn:Nparticle:170}
\text{\(:H:\)} =
\sum_{i = 1}^N \omega_{\Bp_i}
a_{i}^\dagger a_{i}
+
\frac{
3 \lambda
}{2 m^2 }
\sum_{r,s,t,u = 1}^N
a_{r}^\dagger a_{s}^\dagger a_{t} a_{u}.
\end{equation}

Commentary.

While there are a few things that are not wrong above, those correct parts are liberally mixed with a few fundamental errors.

Even before getting to the fundamential errors, there are a few minor issues too. For example, I don’t think that there is anything strictly wrong with the expansion of \ref{eqn:Nparticle:100} for example, although it appears that I confused myself by actually evaluating the delta function instead of just identifying it, something like:
\begin{equation}\label{eqn:Nparticle:190}
\begin{aligned}
H_{\text{int}}
&=
\frac{\lambda}{4} \int \frac{ d^3 p d^3 k d^3 q d^3 s}{(2 \pi)^{12} \sqrt{ \omega_\Bp \omega_\Bk \omega_\Bq \omega_\Bs } }
\Big(
\binom{4}{0}
\delta^3( \Bp + \Bk + \Bq + \Bs )
a_\Bp a_\Bk a_\Bq a_\Bs e^{-i (\omega_\Bp + \omega_\Bk + \omega_\Bq + \omega_\Bs) t} \\
&\qquad +
\binom{4}{1}
\delta^3( -\Bp + \Bk + \Bq + \Bs )
a_\Bp^\dagger a_\Bk a_\Bq a_\Bs e^{-i (-\omega_\Bp + \omega_\Bk + \omega_\Bq + \omega_\Bs)t} \\
&\qquad +
\binom{4}{2}
\delta^3( -\Bp – \Bk + \Bq + \Bs )
a_\Bp^\dagger a_\Bk^\dagger a_\Bq a_\Bs e^{-i (-\omega_\Bp – \omega_\Bk + \omega_\Bq + \omega_\Bs)t} \\
&\qquad +
\binom{4}{3}
\delta^3( -\Bp – \Bk – \Bq + \Bs )
a_\Bp^\dagger a_\Bk^\dagger a_\Bq^\dagger a_\Bs e^{-i (-\omega_\Bp – \omega_\Bk – \omega_\Bq + \omega_\Bs)t} \\
&\qquad +
\binom{4}{4}
\delta^3( -\Bp – \Bk – \Bq – \Bs )
a_\Bp^\dagger a_\Bk^\dagger a_\Bq^\dagger a_\Bs^\dagger e^{-i (-\omega_\Bp – \omega_\Bk – \omega_\Bq – \omega_\Bs)t}
\Big)
\end{aligned}
\end{equation}

So much of QFT is expressed in terms of delta functions, yet I tend to view them as something to be eliminated instead of retained. One lesson is that delta functions in QFT should be embraced, and had I done so in this case, I would have had a better chance of understanding the scattering part of the question that came later. One role of the delta function actually appears to be critical to the scattering question, as it actually encodes the constraints that lead to conservation of momentum in a “scattering process” related to the interaction term of the Hamiltonian.

I didn’t understand what was meant by a scattering process, but that was related to another more fundamental misunderstanding. That screw up was in my interpretation of what was meant by the “N particle sector of the Hilbert space.” In my reading of that phrase I tossed out “Hilbert space” as irrelevant. Part of that was somewhat reactionary, as it seems to me that Hilbert space is tossed around in so much QM in a way that makes it seem more rigorous, but we never touch on the scary mathematics that rigorously defines what a Hilbert space is. It seems to me that the phrase “Hilbert space” is often used as a pretentious way of saying “complex inner product space”. It has a more precise meaning to Mathematicians, but I don’t think that most physicists understand that meaning (and most certainly most students of QM don’t).

So, long story short, I interpreted “N particle sector of the Hilbert space” as “N particle sector” and though that meant I had to somehow discretize the system, introducing N discrete momenta into the mix. Here’s an example, using the non-interaction Hamiltonian, of exactly how I did that:
\begin{equation}\label{eqn:Nparticle:210}
\int \frac{d^3 p}{(2 \pi)^3} \omega_\Bp a_\Bp^\dagger a_\Bp \rightarrow \sum_{i = 1}^N \omega_\Bp a_\Bp^\dagger a_\Bp.
\end{equation}

This is wrong! I actually clued into part of the trouble with this, but didn’t know what the root cause was. I fudged around that my dropping the problematic side effects of having evaluated the delta function, which lead me to have terms like \( a_{\Bp + \Bq – \Bk} \) in the expansion of the interaction Hamiltonian. I couldn’t see how creation and annihilation operators that are associated with arbitrary linear combinations of the other momenta could maintain an \( N \) particle system. Those linear combinations could easily lie outside of the set of the original \( N \) particle system that I constructed in my discretization.

In office hours, Professor Poppitz had me work a small problem to illustrate the error in my ways, specifically applying the Hamiltonian to a two momentum state
\begin{equation}\label{eqn:Nparticle:230}
\begin{aligned}
H \ket{\Bp_1 \Bp_2}
&=
\lr{
\int \frac{d^3}{(2 \pi)^3}
a_\Bp^\dagger
a_\Bp
}
a_{\Bp_1}^\dagger
a_{\Bp_2}^\dagger
\ket{0} \\
&=
\int \frac{d^3}{(2 \pi)^3} \omega_\Bp
a_\Bp^\dagger
\lr{
a_{\Bp_1}^\dagger a_\Bp
+ (2 \pi) \delta^3(\Bp -\Bp_1)
}
a_{\Bp_2}^\dagger
\ket{0} \\
&=
\int \frac{d^3}{(2 \pi)^3} \omega_\Bp
a_\Bp^\dagger
a_{\Bp_1}^\dagger a_\Bp
a_{\Bp_2}^\dagger
\ket{0}
+
\omega_{\Bp_1}
a_{\Bp_1}^\dagger
a_{\Bp_2}^\dagger
\ket{0} \\
&=
\int \frac{d^3}{(2 \pi)^3} \omega_\Bp
a_\Bp^\dagger
a_{\Bp_1}^\dagger
\lr{
a_{\Bp_2}^\dagger
a_\Bp
+ (2 \pi)^3 \delta^3(\Bp – \Bp_2)
}
\ket{0}
+
\omega_{\Bp_1}\ket{\Bp_1 \Bp_2} \\
&=
\lr{ \omega_{\Bp_1} + \omega_{\Bp_2} }
\ket{\Bp_1 \Bp_2}.
\end{aligned}
\end{equation}
Observe that the Hamiltonian operates on a two momentum state, returning that state scaled by the energy associated with the sum of the momenta. Given an two particle subspace of all the possible two momentum states, perhaps with a basis like \( \setlr{\ket{\Bp_1 \Bp_2}, \ket{\Bp_3 \Bp_4}, \cdots } \), the Hamiltonian provides a mapping from that set onto itself, as it scales the states in question, but does not convert two particle states into combinations of one and three momentum states (say). This is not the case with the interaction Hamiltonian, as an operation like \( \int a^\dagger a^\dagger a^\dagger a \) doesn’t preserve a two momentum state. Example
\begin{equation}\label{eqn:Nparticle:250}
\begin{aligned}
\int
&\frac{d^3 p}{(2\pi)^3\sqrt{ 2\omega_\Bp}}
\frac{d^3 q}{(2\pi)^3\sqrt{ 2\omega_\Bq}}
\frac{d^3 r}{(2\pi)^3\sqrt{ 2\omega_\Br}}
\frac{d^3 s}{(2\pi)^3\sqrt{ 2\omega_\Bs}}
\delta^3(\Bp + \Bq + \Br -\Bs)
a_\Bp^\dagger
a_\Bq^\dagger
a_\Br^\dagger
a_\Bs
\ket{ \Bp_1 \Bp_2 } \\
&=
\int
\frac{d^3 p}{(2\pi)^3\sqrt{ 2\omega_\Bp}}
\frac{d^3 q}{(2\pi)^3\sqrt{ 2\omega_\Bq}}
\frac{d^3 r}{(2\pi)^3\sqrt{ 2\omega_\Br}}
\frac{d^3 s}{(2\pi)^3\sqrt{ 2\omega_\Bs}}
\delta^3(\Bp + \Bq + \Br -\Bs)
a_\Bp^\dagger
a_\Bq^\dagger
a_\Br^\dagger
a_\Bs
a_{\Bp_1}^\dagger
a_{\Bp_2}^\dagger
\ket{0} \\
&=
\int
\frac{d^3 p}{(2\pi)^3\sqrt{ 2\omega_\Bp}}
\frac{d^3 q}{(2\pi)^3\sqrt{ 2\omega_\Bq}}
\frac{d^3 r}{(2\pi)^3\sqrt{ 2\omega_\Br}}
\frac{d^3 s}{(2\pi)^3\sqrt{ 2\omega_\Bs}}
\delta^3(\Bp + \Bq + \Br -\Bs)
a_\Bp^\dagger
a_\Bq^\dagger
a_\Br^\dagger
\lr{
a_{\Bp_1}^\dagger
a_\Bs
+ (2 \pi)^3 \delta^3(\Bp_1 -\Bs)
}
a_{\Bp_2}^\dagger
\ket{0} \\
&=
\int
\frac{d^3 p}{(2\pi)^3\sqrt{ 2\omega_\Bp}}
\frac{d^3 q}{(2\pi)^3\sqrt{ 2\omega_\Bq}}
\frac{d^3 r}{(2\pi)^3\sqrt{ 2\omega_\Br}}
\frac{d^3 s}{(2\pi)^3\sqrt{ 2\omega_\Bs}}
\delta^3(\Bp + \Bq + \Br -\Bs)
a_\Bp^\dagger
a_\Bq^\dagger
a_\Br^\dagger
a_{\Bp_1}^\dagger
\lr{
a_{\Bp_2}^\dagger
a_\Bs
+ (2\pi)^3 \delta^3(\Bp_2 – \Bs)
}
\ket{0} \\
&\quad +
\int
\frac{d^3 p}{(2\pi)^3\sqrt{ 2\omega_\Bp}}
\frac{d^3 q}{(2\pi)^3\sqrt{ 2\omega_\Bq}}
\frac{d^3 r}{(2\pi)^3\sqrt{ 2\omega_\Br}}
\inv{\sqrt{ 2\omega_{\Bp_1}}}
\delta^3(\Bp + \Bq + \Br -\Bp_1)
a_\Bp^\dagger
a_\Bq^\dagger
a_\Br^\dagger
a_{\Bp_2}^\dagger
\ket{0} \\
&=
\inv{\sqrt{ 2\omega_{\Bp_2}}}
\lr{
\int
\frac{d^3 p}{(2\pi)^3\sqrt{ 2\omega_\Bp}}
\frac{d^3 q}{(2\pi)^3\sqrt{ 2\omega_\Bq}}
\frac{d^3 r}{(2\pi)^3\sqrt{ 2\omega_\Br}}
\delta^3(\Bp + \Bq + \Br -{\Bp_2})
}
\ket{\Bp \Bq \Br \Bp_1} \\
&\quad +
\inv{\sqrt{ 2\omega_{\Bp_1}}}
\lr{
\int
\frac{d^3 p}{(2\pi)^3\sqrt{ 2\omega_\Bp}}
\frac{d^3 q}{(2\pi)^3\sqrt{ 2\omega_\Bq}}
\frac{d^3 r}{(2\pi)^3\sqrt{ 2\omega_\Br}}
\delta^3(\Bp + \Bq + \Br -\Bp_1)
}
\ket{ \Bp \Bq \Br \Bp_2 }.
\end{aligned}
\end{equation}
Such a term from the interaction Hamiltonian maps a two momentum state to a four momentum state. Notice how the continuous representation is critical to this evaluation, as well as that the action of the non-interaction Hamiltonian on a two (or N) momentum state. Attempting any sort of discretization leaves you with an operator that cannot be evaluated.

Observe that only the \( (a^\dagger)^2 a^2 \) terms in the interaction will map \((1,2,3,\cdots,N)\)-momentum states to \((1,2,3,\cdots,N)\)-momentum states, so the language that I didn’t understand “N particle sector of the Hilbert” space was really an encoded instruction to retain only those interaction terms. I did exactly that because intuition told me to do so (and didn’t justify why I did so), but I had the wrong reasons for making that selection. I knew there was something wrong, but didn’t know what it was. What I should have done was go back to basics and root out all the aspects of the problem statement that I did not understand, and ask about those. If I had done so (in a timely fashion, and not at the last minute when I attempted this problem), then I wouldn’t have gone down a dead end path that lead to more confusion. I didn’t even get to the interesting part of this problem, which was to show the correspondence between the QFT picture and the classical QM picture. I’ll still attempt to do so, despite having lost my window to get credit for that work.