Mathematica Assumptions for integration [pertaining to loop corrections]

Click For Summary
SUMMARY

The discussion focuses on integrating a complex expression in Mathematica while applying assumptions for the variables b and z. The user is attempting to correctly implement assumptions using the syntax "Assuming[0 ≤ b ≤ 1 && 0 ≤ z ≤ 1, Integrate[...]]" but has encountered inconsistencies in results. It is noted that changing the assumptions from ≤ to < improves the output. Additionally, a substitution method involving known function values, such as replacing Q with Zeta[3]/3, is suggested to simplify the integration process.

PREREQUISITES
  • Familiarity with Mathematica syntax and functions
  • Understanding of integral calculus and Feynman parametrization
  • Knowledge of variable assumptions in symbolic computation
  • Experience with function substitution techniques in mathematical expressions
NEXT STEPS
  • Research "Mathematica Assumptions syntax" for proper implementation
  • Learn about "Feynman parametrization in quantum field theory" for context
  • Explore "Mathematica substitution methods" for simplifying integrals
  • Study "integrating complex expressions in Mathematica" for advanced techniques
USEFUL FOR

This discussion is beneficial for physicists, mathematicians, and researchers using Mathematica for complex integrations, particularly those working with quantum field theory and requiring precise assumptions in their calculations.

Elwin.Martin
Messages
202
Reaction score
0
Long story short I have an integral which is something like this:
(in Mathematica code)
Integrate[1/{(1 - b^2)*{((1 - y)*z + y)^2 + (1 - y)^2*(1 - z)^2} + 2*(b^2 + 1)*{(1 - y)*(1 - z)*((1 - y)*z + y)}}, {y, 0, 1} ]

Written without assumptions.

b is actually only in [0,1] and I think I can limit z to [0,1] as well. z is actually in the next integral ^^; which is from 0 to 1.

How do I add assumptions in proper syntax? I've tried two different ways and gotten VERY different results...so I'm a little concerned about my syntax and Mathematica's reference site didn't help much.

I have been using something like ", Assumptions-> 0<=b<=1, 0<=z<=1" I think...but I'm not positive about my consistency.

For those of you who don't have direct access to Mathematica or anything and don't read raw code super easily here's the integral in traditional form:
\int^1_0 \frac{dy}{{(1 - b^2)*{((1 - y)*z + y)^2 + (1 - y)^2*(1 - z)^2} + 2*(b^2 + 1)*{(1 - y)*(1 - z)*((1 - y)*z + y)}}}

If context is helpful, the integration is the result of a Feynman parametrization of a basic QCD process. Everything is scalar though, the vectors have all been dotted or otherwise taken care of and absorbed into the b term (Sort of). The rest of the integrations are a whole other mess.
 
Physics news on Phys.org
Assuming[0 ≤ b ≤ 1 && 0 ≤ z ≤ 1, Integrate[1/{(1 - b^2)*{((1 - y)*z + y)^2 + (1 - y)^2*(1 - z)^2} + 2*(b^2 + 1)*{(1 - y)*(1 - z)*((1 - y)*z + y)}}, {y, 0, 1}]]

but it looks like it is happier if you make those < rather than ≤ in all cases which I think makes sense if I look at your denominator for a few seconds.
 
Bill Simpson said:
Assuming[0 ≤ b ≤ 1 && 0 ≤ z ≤ 1, Integrate[1/{(1 - b^2)*{((1 - y)*z + y)^2 + (1 - y)^2*(1 - z)^2} + 2*(b^2 + 1)*{(1 - y)*(1 - z)*((1 - y)*z + y)}}, {y, 0, 1}]]

but it looks like it is happier if you make those < rather than ≤ in all cases which I think makes sense if I look at your denominator for a few seconds.

Thanks, I'll try that and see how it goes.

Someone else recommended I use a substitution of a known function value or something for anything that's parametized like what I have?

For example, he said that if I had some Q running from 0 to 1, I could sub {Q->Zeta[3]/3}.
Integrate.
Sub {Zeta[3]->3Q.

It came out sort of close to what I had before, but I'm not entirely sure why.
Any idea how he got this? I didn't have much of a chance to question him about it.

I'm still playing with different factorizations to make it easier for Mathematica and myself...I'll see where that takes me, too.
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 13 ·
Replies
13
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 6 ·
Replies
6
Views
4K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 9 ·
Replies
9
Views
4K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 4 ·
Replies
4
Views
2K