Finding Angular & Linear Momentum of EM Fields

[FlatLander]

Homework Statement

A point charge q is a distance a>R from the axis of an infinite solenoid (radius R, n turns per unit length, current I). Find the linear momentum and the angular momentum in the fields. (Put q on the x axis, with the solenoid along z; treat the solenoid as a nonconductor, so you don't need to worry about induced charges on its surface.)[Answer: P_em=μ₀qnIR²/2a; L_em=0]

Homework Equations

\vec{E}_q=q/4\pi\epsilon₀(1/\vec{r}²)=q/4\pi\epsilon₀(\vec{r}/r³)
\vec{B}_sol=μ₀nI\hat{z}
p_em=ε₀(\vec{E}\times\vec{B})
l_em=r\timesp_em
P_em=∫p_em d\tau
L_em=∫l_em d\tau

The Attempt at a Solution

I kind of plugged and chugged, found r²=((x-a)²+y²+z²) and \vec{r}=(x-a)\hat{x}+y\hat{y}+z\hat{z}
Plugged in for that as well. However, I eventually got to the integrations in for the P_em and realized I don't know what my limits of integration are for the volume.
I know the z is from -∞ to ∞, but I have no clue for x and y. Here is what my final line looks like so far (with me already integrating over z):
P_em=\frac{-2\mu_{0}\epsilon_{0}qnI}{4\pi\epsilon_{0}}∫\frac{(x-a)\hat{y}}{((x-a)^{2}+y^{2}} dydx

Any help would be appreciated. Thanks.

 

[FlatLander]

I figured it out. I'm too lazy to type out the solution, and since I know this could potentially help someone in the future searching for it, I'll give you the basic steps:
1. Change the integral I had in the last line to cylindrical coordinates.
[From here it's just tricky calculus]
2. Do the phi integral first and break it up into two parts
3. Here's where it gets tricky. You need to define two new variables, say A and B, where A=s^2+a^2, B=-2as.
4. Do the integrals. DO NOT do them with a calculator or other computation engine. It WILL give you wrong answer. Use an integral table or do them by hand (good luck).
5. Finding Lem is trivial from this point on, you find lem, then do the integral. Repeat steps 1-4 with the integral you get from Lem.