Wolfram Language Paclet Repository
Community-contributed installable additions to the Wolfram Language
GSberveglieri`Phi4tools`
TadSunsetSubdiagram  | |
|
| | ||||
Details and Options
Examples
(1)
Basic Examples
(1)
In[1]:=
Needs["GSberveglieri`Phi4tools`"]
Let's look at the diagram for for the theory with =3 quartic vertices, let's visualize it and print its integrand:
st
1
(2)
Γ
4
ϕ
v4
In[2]:=
[2,0,3,1],
[2,0,3,1]
 | VisualizeDiagram |
| IntegrandDiagram |
Out[2]=
,[[1]][[2]][-[1]-[2]-[3]][[3]][[1]+[2]+[3]]
1
12
The function prints the integrand in in terms of the three-dimensional components in spherical coordinates ready to be integrated.
d=3
In[3]:=
integrand2031=
@
[2,0,3,1]
 | WriteExplicit |
| IntegrandDiagram |
Out[3]=
1121+1+1+
2
[1]
ρ
2
[2]
ρ
2
[3]
ρ
2
1+++
2
(+Cos[]+Cos[])
[1]
ρ
[2]
θ
[2]
ρ
[3]
θ
[3]
ρ
2
Sin[]+CosSin[]
[2]
ρ
[2]
θ
[3]
ϕ
[3]
ρ
[3]
θ
2
[3]
ρ
2
Sin[]
[3]
θ
2
Sin
[3]
ϕ
The function returns in output the variables we have to integrate over to compute the value of this diagram. In this case we will have a six-dimensional integral to perform.
In[4]:=
 | MomVars |
Out[4]=
,,,,,
[1]
ρ
[2]
θ
[2]
ρ
[3]
θ
[3]
ρ
[3]
ϕ
Let's look at the case with the tadpole-like subdiagram substituted
In[5]:=
[2,0,3,1,"Substitutions""Analytics"],
[2,0,3,1,"Substitutions""Analytics"]
| VisualizeDiagram |
| IntegrandDiagram |
Out[5]=
,
12
In[6]:=
| WriteExplicit |
| IntegrandDiagram |
Out[6]=
-
Log
4
3
1536
3
π
This is directly the value of this Feynman diagram in .
d=3
The substitution of this tadpole-like subdiagram with sunset reduces the number of loop of the diagram by three. Let's look at an example: the diagram for for the theory with vertices =6 quartic vertices, without and with the subdiagrams substitutions.
st
1
(0)
Γ
4
ϕ
v4
In[7]:=
[0,0,6,1],
[0,0,6,1,"Substitutions""Analytics"]
| VisualizeDiagram |
| VisualizeDiagram |
Out[7]=
,
In[8]:=
| IntegrandDiagram |
| IntegrandDiagram |
Out[8]=
1
576
2
[[3]]
2
[[4]]
2
[[5]]
Out[8]=
1
576
2
The 3 substitutions reduced the number of loops from 7 to 0, as we can see printed using the function on the respective integrands
In[9]:=
| CountLoops |
 | IntegrandDiagram |
| IntegrandDiagram |
Out[9]=
{7,0}
 |
|
 |
|
""