Vtyjkyuk

I am trying to do a six dimensional integration involving functions of two vectors $vecp$ and $vecp'$ in spherical coordinates. So $vecp$ is reprsented by $p,theta,phi$, and $vecp'$ is represented by $p',theta',phi'$. The $theta$'s range from 0 to $pi$, and the $phi$'s range from 0 to $2pi$. The angle between the two vectors is given by $cosgamma = costhetacostheta' + sinthetasintheta'cos(phi-phi')$. I cannot take any of these to vectors to be the $z$-axis, since the $z$-axis is fixed by a third vector.

Suppose the function that I have to integrate is
$$
fraccos(phi-phi')[p^2 + p'^2 + 2pp'cosgamma + M^2]^2
$$
where $M$ is a real constant. Then the integrals over $phi$ and $phi'$ would be of the form
$$
int_0^2pi dphi' int_0^2pi dphi fraccos(phi-phi')[a + bcos(phi - phi')]^2
$$
where I have used $a$ and $b$ to represent the terms that are independent of $phi$ and $phi'$.
If I now do the following substitution
$$
sin(phi - phi') = u
$$
then the integration becomes
$$
int_0^2pi dphi' int_-sinphi'^-sinphi' dufrac1[a+bsqrt1-u^2]^2
$$
which is just zero.

There are a bunch of other functions which involve $cos(phi-phi')$ that I am trying to integrate. I was expecting most of these to be nonzero. But from what I see it appears to me that if it involves $cos(phi-phi')$ then for the given range of $phi$ the integral will always be zero.

Am I doing something wrong here?

The basic problem is that your substitution is not one to one in the entire range of the integral.

You can correct the situation by splitting up the integral using intervals on which it is one to one

note that $ cos theta = -sqrt1-sin^2theta$ when $ theta in (frac pi 2, frac 3pi 2)$

In General
$$ I= int_0^2picostheta f(cos theta)dtheta
\= int_0^frac pi 2costheta f(cos theta)dtheta +int_frac pi 2^frac 3pi 2costheta f(cos theta)dtheta + int_frac 3pi 2^2picostheta f(cos theta)dtheta $$
now we can make the sine substitution in all three integrals, being careful to use $costheta = -sqrt1-u^2 $ in the second integral
$$
I=int _0^1f(sqrt1-u^2)du+int _1^-1f(-sqrt1-u^2)du
+int _-1^0f(sqrt1-u^2)du
\=int_-1^1Bigg[ f(sqrt1-u^2)-f(-sqrt1-u^2) Bigg]du
$$