Dimensional Regularization and Divergences

[latentcorpse]

I'm trying question 1 in this past paper:
http://www.maths.cam.ac.uk/postgrad/mathiii/pastpapers/2009/Paper48.pdf

I'm on the bit "Show that the divergent part of this diagram can be written in the form Ap^2+B where p=p_1=-p_2 and A and B are momentum independent.

Now, as far as I can tell, this calculation is going to be exactly the same as the one on p48 of these notes:
http://www.damtp.cam.ac.uk/user/ho/Notes.pdf
(the final result being given at the top of p49) with the only difference being that when we use dimensional regularisation to analytically continue to d \in \mathbb{C}, we should set d=6-\epsilon instead of d=4-\epsilon. However, this will change the powers of some factors in the final answer, it won't change the pole structure, will it? I certainly don't see how we're going to need an Ap^2 counter term as there is no divergent piece multiplying a p^2 that I can see...

Can anyone offer some help please...

 

[fzero]

I'm trying question 1 in this past paper:
http://www.maths.cam.ac.uk/postgrad/mathiii/pastpapers/2009/Paper48.pdf

I'm on the bit "Show that the divergent part of this diagram can be written in the form Ap^2+B where p=p_1=-p_2 and A and B are momentum independent.

Now, as far as I can tell, this calculation is going to be exactly the same as the one on p48 of these notes:
http://www.damtp.cam.ac.uk/user/ho/Notes.pdf
(the final result being given at the top of p49) with the only difference being that when we use dimensional regularisation to analytically continue to d \in \mathbb{C}, we should set d=6-\epsilon instead of d=4-\epsilon. However, this will change the powers of some factors in the final answer, it won't change the pole structure, will it? I certainly don't see how we're going to need an Ap^2 counter term as there is no divergent piece multiplying a p^2 that I can see...

Can anyone offer some help please...

For d=4-\epsilon the divergence is in \Gamma(0). For d=6-\epsilon it's in \Gamma(-1).

 

[latentcorpse]

For d=4-\epsilon the divergence is in \Gamma(0). For d=6-\epsilon it's in \Gamma(-1).

So I find that \hat{\tau}_2^{(1)}(p,-p)=\frac{1}{2}g^2 \frac{1}{(4 \pi)^3} (4 \pi)^\frac{\epsilon}{2} \Gamma(-1+\frac{\epsilon}{2}) \int_0^1 d\alpha AA^{-\frac{\epsilon}{2}}
=\frac{1}{2}g^2 \frac{1}{(4 \pi)^3} (4 \pi)^\frac{\epsilon}{2} \Gamma(\frac{\epsilon}{2})\frac{1}{\frac{\epsilon}{2}-1} \int_0^1 d\alpha Ae^{-\frac{\epsilon}{2} \log{A}}
=\frac{1}{2}g^2 \frac{1}{(4 \pi)^3} (4 \pi)^\frac{\epsilon}{2} \Gamma(\frac{\epsilon}{2}+1)\frac{1}{(\frac{\epsilon}{2}-1)\frac{\epsilon}{2}} \int_0^1 d\alpha Ae^{-\frac{\epsilon}{2} \log{A}}
=\frac{1}{2}g^2 \frac{1}{(4 \pi)^3} (4 \pi)^\frac{\epsilon}{2} e^{-\gamma_E \frac{\epsilon}{2}} \frac{1}{(\frac{\epsilon}{2}-1)\frac{\epsilon}{2}} \int_0^1 d\alpha Ae^{-\frac{\epsilon}{2} \log{A}}
=\frac{1}{2}g^2 \frac{1}{(4 \pi)^3} ((4 \pi)e^{-\gamma_E})^\frac{\epsilon}{2} (-\frac{2}{\epsilon} ( 1 + \frac{\epsilon}{2} + \mathcal{O}(\epsilon^2)) \int_0^1 d\alpha A (1-\frac{\epsilon}{2} \log{A})

So this I think is my divergent piece. The p^2 term is contained in A=\alpha(1-\alpha)p^2+m^2.
What would you recommend now? Should I go through and find all pieces that are of \mathcal{O}(\frac{1}{\epsilon}) and then I should be able to read of A and B in theory, right?

 

[fzero]

Should I go through and find all pieces that are of \mathcal{O}(\frac{1}{\epsilon}) and then I should be able to read of A and B in theory, right?

Yes.

 

[latentcorpse]

Yes.

There is a bit of a complication though because, when I multiply this out I get

=-\frac{g^2}{2} \frac{1}{(4 \pi)^3} ( 4 \pi e^{-\gamma_E )^{\epsilon}{2} \int_0^1 A ( \frac{2}{\epsilon}-1+1-\frac{\epsilon}{2}) \log{A} d \alpha

Now this means I get a piece that goes as \frac{2}{\epsilon} and a piece that goes as {\epsilon}{2} (which I'm assuming we can neglect as its not divergent in the regulator limit)

So the divergent bit is =-\frac{g^2}{2} \frac{1}{(4 \pi)^3} ( 4 \pi e^{-\gamma_E )^{\epsilon}{2} \int_0^1 A \frac{2}{\epsilon} \log{A} d \alpha

But A=\alpha(1-\alpha)p^2_m^2 so I don't see how I can isolate a coefficient A of p^2 as there is a p^2 in the log as well, isn't there?

 

[fzero]

It's clear from the last line in post #3 that the \epsilon in front of the log term means that it doesn't contribute to the divergent part of the expression, so check your algebra.

 

[latentcorpse]

It's clear from the last line in post #3 that the \epsilon in front of the log term means that it doesn't contribute to the divergent part of the expression, so check your algebra.

Yeah. I made a mistake half way through and didn't correct it properly throughout.
So anyway I find A=-\frac{g^2}{2} \frac{1}{(4 \pi)^3} \frac{1}{\epsilon} and B=-g^2 \frac{m^2}{(4 \pi)^3} \frac{1}{\epsilon}

Now for the final part he asks me to calculate B using a UV cut-off. I don't understand how he wants me to do that. There is no integral in B and we need an integral to be present to use UV cut-off regularisation, don't we?

 

[fzero]

Yeah. I made a mistake half way through and didn't correct it properly throughout.
So anyway I find A=-\frac{g^2}{2} \frac{1}{(4 \pi)^3} \frac{1}{\epsilon} and B=-g^2 \frac{m^2}{(4 \pi)^3} \frac{1}{\epsilon}

Now for the final part he asks me to calculate B using a UV cut-off. I don't understand how he wants me to do that. There is no integral in B and we need an integral to be present to use UV cut-off regularisation, don't we?

B is (part of) the divergent part of an integral. Use that integral.

 

[latentcorpse]

B is (part of) the divergent part of an integral. Use that integral.

Well certainly B came from an integral over the loop momenta i.e. we had \int \frac{d^dk}{(2 \pi)^d} \dots

However, I used that integral when I introduced a gamma function during the dimensional regularisation procedure - there don't appear to be any integral signs left over i.e. the whole integral has been evaluated and has produced Ap^2+B as an answer...

 

[fzero]

Well certainly B came from an integral over the loop momenta i.e. we had \int \frac{d^dk}{(2 \pi)^d} \dots

However, I used that integral when I introduced a gamma function during the dimensional regularisation procedure - there don't appear to be any integral signs left over i.e. the whole integral has been evaluated and has produced Ap^2+B as an answer...

The problem is asking you to evaluate B using primitive cutoff regularization, not dimensional regularization.

 

[latentcorpse]

The problem is asking you to evaluate B using primitive cutoff regularization, not dimensional regularization.

Yes but B is the result of the dimensional regularisation is it not? Perhaps if it asked me to evaluate the original integral using cutoff regularisation i could try that? But B as it stands doesn't have an integral in it so it can't be evaluated any further - by any means.

So I still don't really get what I'm supposed to do because I don't know what the integral I am meant to be evaluating is?

 

[fzero]

Yes but B is the result of the dimensional regularisation is it not? Perhaps if it asked me to evaluate the original integral using cutoff regularisation i could try that? But B as it stands doesn't have an integral in it so it can't be evaluated any further - by any means.

So I still don't really get what I'm supposed to do because I don't know what the integral I am meant to be evaluating is?

The problem didn't specify how you were supposed to regulate the integral. Dimensional regularization was convenient for deriving the form of the divergence, but wasn't necessarily the only way to do it. Go back to the original integral and see if there's a limit in which you can extract B using a cutoff.

 

[latentcorpse]

The problem didn't specify how you were supposed to regulate the integral. Dimensional regularization was convenient for deriving the form of the divergence, but wasn't necessarily the only way to do it. Go back to the original integral and see if there's a limit in which you can extract B using a cutoff.

\hat{F}_2(p_1,p_2)= \frac{g^2}{2} \int_{\mathbb{R}^6} \frac{d^6k}{(2 \pi)^6} \frac{1}{(k^2+m^2)((p-k)^2+m^2)}
= \frac{g^2}{2} \frac{1}{(2 \pi)^6} \int_{S^5} d \Omega_5 \int \frac{k^5 dk}{(k^2+m^2)((p-k)^2+m^2)}
= \frac{g^2}{2^7 \pi^3} \int \frac{k^5 dk}{(k^2+m^2)((p-k)^2+m^2)} as \text{Vol}(S^5)=\pi^3
\rightarrow \frac{g^2}{2^7 \pi^3} \int_{k \leq \Lambda} \frac{k^5 dk}{(k^2+m^2)((p-k)^2+m^2)} by imposing a UV cutoff after wick rotation

Now I am wondering if there is some trick to evaluating this. the only way i can think of is introducing integrals over \alpha so that we can right \frac{1}{k^2+m^2} = \int d \alpha e^{-\alpha (k^2+m^2)} but this was from dimensional regularisation so not sure if it will be any use here!

What would you advise for evaluating this?

 

[fzero]

You don't need to compute the integral, you only need to obtain the divergent part. So you can consider the range where that comes from. Furthermore you don't have to compute the whole divergent part, only the part that is independent of p.

 

[latentcorpse]

You don't need to compute the integral, you only need to obtain the divergent part. So you can consider the range where that comes from. Furthermore you don't have to compute the whole divergent part, only the part that is independent of p.

I don't understand. I thought when we used primitive cutoff regularisation we explicitly did the integral but over a range 0 \leq k \leq \Lambda rather than 0 \leq k \leq \infty?

What should I do here then?

 

[fzero]

I don't understand. I thought when we used primitive cutoff regularisation we explicitly did the integral but over a range 0 \leq k \leq \Lambda rather than 0 \leq k \leq \infty?

What should I do here then?

Suppose that we have

I = \int \frac{d^4 k}{k^2 (k^2+m^2)}.

The divergence is in the UV, so we can introduce an intermediate scale M\gg m and write

I = \int_0^M \frac{dk k^3}{ k^2 (k^2+m^2)} + \int_M^\Lambda \frac{dk k^3}{ k^2 (k^2+m^2)} .

The divergence is completely contained in the 2nd term, which we can write as

\int_M^\Lambda \frac{dk k^3}{ k^2 (k^2+m^2)} = \int_M^\Lambda dk \frac{dk}{ k} \left( 1 - \frac{m^2}{k^2} + \cdots \right).

Now the divergence is completely contained in the leading term, since the higher order terms are finite.

 

[latentcorpse]

Suppose that we have

I = \int \frac{d^4 k}{k^2 (k^2+m^2)}.

The divergence is in the UV, so we can introduce an intermediate scale M\gg m and write

I = \int_0^M \frac{dk k^3}{ k^2 (k^2+m^2)} + \int_M^\Lambda \frac{dk k^3}{ k^2 (k^2+m^2)} .

The divergence is completely contained in the 2nd term, which we can write as

\int_M^\Lambda \frac{dk k^3}{ k^2 (k^2+m^2)} = \int_M^\Lambda dk \frac{dk}{ k} \left( 1 - \frac{m^2}{k^2} + \cdots \right).

Now the divergence is completely contained in the leading term, since the higher order terms are finite.

Why did you need to split the integral at all?

I can't see why the higher order terms should be finite?

Lastly, should I use binomial on both the terms I have in the denominator? That ((p-k)^2+m^2) term always confuses me on how to deal with it

Thanks.

 

[fzero]

Why did you need to split the integral at all?

You're right that you don't if you use the steps above. I was originally going to say that in the 2nd term we could just set m=0, but it's clearer to expand the way I ended up doing.

I can't see why the higher order terms should be finite?

Try to do the integral of the subleading term.

Lastly, should I use binomial on both the terms I have in the denominator? That ((p-k)^2+m^2) term always confuses me on how to deal with it

Thanks.

You have to be a bit careful, since (p-k)^2 depends on the angles in k space. If you had to compute the whole integral, you'd introduce Feynman parameters. However for your problem, you only have to compute the part of the divergence that is independent of p, so you can simplify things a bit before manipulating the integral.

 

[latentcorpse]

You're right that you don't if you use the steps above. I was originally going to say that in the 2nd term we could just set m=0, but it's clearer to expand the way I ended up doing.




Try to do the integral of the subleading term.

Yeah because it only becomes divergent for powers -1 or more and subleading term is -3 therefore finite.

You have to be a bit careful, since (p-k)^2 depends on the angles in k space. If you had to compute the whole integral, you'd introduce Feynman parameters. However for your problem, you only have to compute the part of the divergence that is independent of p, so you can simplify things a bit before manipulating the integral.

Well to isolate the p terms, i tried to complete the square on the denominator?

Is that right? It's not going very well:

(k^2+m^2)(p^2-2pk+k^2+m^2) = p^2k^2-2pk^3+k^4+k^2m^2 +m^2p^2-2pkm^2+k^2m^2+m^4
There's too many terms to complete the square here so I don't really know how to simplify it to get rid of the p terms?

 

[fzero]

Well to isolate the p terms, i tried to complete the square on the denominator?

Is that right? It's not going very well:

(k^2+m^2)(p^2-2pk+k^2+m^2) = p^2k^2-2pk^3+k^4+k^2m^2 +m^2p^2-2pkm^2+k^2m^2+m^4
There's too many terms to complete the square here so I don't really know how to simplify it to get rid of the p terms?

Independent of p means that you can set p=0 to evaluate the part that you want.

 

[latentcorpse]

Independent of p means that you can set p=0 to evaluate the part that you want.

OK I think I have it. So we get

\int_0^\Lambda \frac{k^5 dk}{(k^2+m^2)^2} = \int_0^\Lambda \frac{k^5 dk}{k^4(1+\frac{m^2}{k^2})^2} = \int_0^\Lambda dk k ( 1 - \frac{2m^2}{k^2} + \dots ) where all other pieces are finite

So the divergent pieces are equal to \int_0^\Lambda dk k - \frac{2m^2}{k} = \lambda - 2m^2 ( \log{\Lambda} - \log{0})

Now I'm assuming this is why you wanted me to split the integral earlier? To avoid the log(0)? Why is the log(0) appearing at all though? Isn't that a result of us having cut a corner somewhere?

In which case we get a quadratic and a logarithmic divergence of the form \frac{\Lambda^2}{2} - 2m^2 \log{\Lambda}

This doesn't quite seem to match up with the B we found earlier in post 7 though?

 

[fzero]

OK I think I have it. So we get

\int_0^\Lambda \frac{k^5 dk}{(k^2+m^2)^2} = \int_0^\Lambda \frac{k^5 dk}{k^4(1+\frac{m^2}{k^2})^2} = \int_0^\Lambda dk k ( 1 - \frac{2m^2}{k^2} + \dots ) where all other pieces are finite

So the divergent pieces are equal to \int_0^\Lambda dk k - \frac{2m^2}{k} = \lambda - 2m^2 ( \log{\Lambda} - \log{0})

Now I'm assuming this is why you wanted me to split the integral earlier? To avoid the log(0)? Why is the log(0) appearing at all though? Isn't that a result of us having cut a corner somewhere?

That divergence is fictitious. It arises from the k\rightarrow 0 part of the integral, but that's precisely where the expansion in m/k is bad. For small k, we should expand in k/m and we find no divergence at k=0.

In which case we get a quadratic and a logarithmic divergence of the form \frac{\Lambda^2}{2} - 2m^2 \log{\Lambda}

This doesn't quite seem to match up with the B we found earlier in post 7 though?

That's because you're comparing results in two different regularization schemes. In general, physical quantities should agree after renormalization, but the divergent parts are regulator dependent and won't agree.

 

[latentcorpse]

That divergence is fictitious. It arises from the k\rightarrow 0 part of the integral, but that's precisely where the expansion in m/k is bad. For small k, we should expand in k/m and we find no divergence at k=0.



That's because you're comparing results in two different regularization schemes. In general, physical quantities should agree after renormalization, but the divergent parts are regulator dependent and won't agree.

These are exactly the same poles we found in 4D i.e. we got a 1/epsilon (or a quadratic+logarithmic when we used UV cutoff). So what was the effect of higher dimensions? To change the A and B coefficients?

 

[fzero]

These are exactly the same poles we found in 4D i.e. we got a 1/epsilon (or a quadratic+logarithmic when we used UV cutoff). So what was the effect of higher dimensions? To change the A and B coefficients?

It's not the same pole, as I pointed out in post #2. In 4d there was only a log divergence, here in 6d we find a quadratic divergence.

 

[latentcorpse]

It's not the same pole, as I pointed out in post #2. In 4d there was only a log divergence, here in 6d we find a quadratic divergence.

Is the way to see that 4D only gives logarithmic because we produced a \frac{1}{\epsilon} pole by dim reg and since a \frac{1}{\epsilon} in dim reg corresponds to a \log{\frac{\Lambda}{m}} in UV cutoff we only have a log. Is that correct?

 

[fzero]

Is the way to see that 4D only gives logarithmic because we produced a \frac{1}{\epsilon} pole by dim reg and since a \frac{1}{\epsilon} in dim reg corresponds to a \log{\frac{\Lambda}{m}} in UV cutoff we only have a log. Is that correct?

Dimensional regularization of any single divergent loop integral always gives the divergence in terms of a simple pole, because the singularities of the Gamma function are simple poles. It seems like logarithmic singularities occur when the divergence is of the form \Gamma(\epsilon). Higher-order singularities will appear as \Gamma(-n+\epsilon), where n is an integer. For example, this quadratic divergence appears as \Gamma(-1), but it's still a simple pole. We can't use the power of \epsilon to determine the degree of divergence.