Pensieve header: An false implementation of nilpotent integration (still using exp).

Initialization

SetDirectory["C:\\drorbn\\AcademicPensieve\\Projects\\HigherRank"];
Once[<<KnotTheory`; <<Rot.m];

CCF[ℰ_] := ExpandDenominator@ExpandNumerator@Together[ℰ];
CCF[ℰ_] := Factor[ℰ];
CF[ℰ_𝔼]:=CF/@ℰ;
CF[ℰ_List]:=CF/@ℰ;
CF[sd_SeriesData] := MapAt[CF,sd,3];
CF[ℰ_] := Module[{vs=Cases[ℰ,(Subscript[(x|p), __]),∞]⋃{x,p},ps,c},
Total[CoefficientRules[Expand[ℰ], vs] /. (ps_->c_) :> CCF[c](Times@@(vs^ps))] ];

Integration

𝔼/:𝔼[A_]𝔼[B_]:=𝔼[A+B]

Unprotect[Integrate];
∫ω_.𝔼[L_]\[DifferentialD](vs_List):=Module[{n,Q,G,V,s,t,k,a,b},
n=Length@vs;
Q=-Table[((∂_(vs[[a]],vs[[b]]))L)/.Thread[vs->0],{a,n},{b,n}];
G=Inverse[Q]/2;
V=L+vs.Q.vs/2;
s=t=V;k=0;
While[0=!=t,
s+=(1/((++k)!))(t=CF@Sum[G[[a,b]](((∂_(vs[[a]],vs[[b]]))t)+((∂_(vs[[a]]))t)((∂_(vs[[b]]))t)),{a,n},{b,n}])];
PowerExpand@Factor[ω(((Det[Q](((2π))^n)))^(-1/2))]𝔼[CF@s/.Thread[vs->0]]
];
Protect[Integrate];

∫𝔼[I λ (Subscript[x, 1]^2)/2]\[DifferentialD]{(Subscript[x, 1])}

∫𝔼[-I λ (Subscript[x, 1]^2)/2]\[DifferentialD]{(Subscript[x, 1])}

∫𝔼[(I/2){(Subscript[x, 1]),(Subscript[x, 2])}.(({{a, b}, {b, c}})).{(Subscript[x, 1]),(Subscript[x, 2])}]\[DifferentialD]{(Subscript[x, 1]),(Subscript[x, 2])}

∫𝔼[-λ (Subscript[x, 1]^2)/2]\[DifferentialD]{(Subscript[x, 1])}

∫𝔼[-(Subscript[x, 1]^2)/2+ξ (Subscript[x, 1])]\[DifferentialD]{(Subscript[x, 1])}

∫𝔼[-(1/2){(Subscript[x, 1]),(Subscript[x, 2])}.(({{a, b}, {b, c}})).{(Subscript[x, 1]),(Subscript[x, 2])}+{(Subscript[ξ, 1]),(Subscript[ξ, 2])}.{(Subscript[x, 1]),(Subscript[x, 2])}]\[DifferentialD]{(Subscript[x, 1]),(Subscript[x, 2])}

I1=∫𝔼[-(1/2){(Subscript[x, 1]),(Subscript[x, 2])}.(({{a, b}, {b, c}})).{(Subscript[x, 1]),(Subscript[x, 2])}+{(Subscript[ξ, 1]),(Subscript[ξ, 2])}.{(Subscript[x, 1]),(Subscript[x, 2])}]\[DifferentialD]{(Subscript[x, 1])}

∫I1 \[DifferentialD]{(Subscript[x, 2])}

∫𝔼[-(1/2){(Subscript[y, 1]),(Subscript[y, 2])}.(({{a, b}, {b, c}})).{(Subscript[y, 1]),(Subscript[y, 2])}+{(Subscript[η, 1]),(Subscript[η, 2])}.{(Subscript[y, 1]),(Subscript[y, 2])}]\[DifferentialD]{(Subscript[y, 1]),(Subscript[y, 2])}

I1=∫𝔼[-(1/2){(Subscript[y, 1]),(Subscript[y, 2])}.(({{a, b}, {b, c}})).{(Subscript[y, 1]),(Subscript[y, 2])}+{(Subscript[η, 1]),(Subscript[η, 2])}.{(Subscript[y, 1]),(Subscript[y, 2])}]\[DifferentialD]{(Subscript[y, 1])}

∫I1 \[DifferentialD]{(Subscript[y, 2])}

∫𝔼[ξ x+η y+z(x-y)+(x^2)]\[DifferentialD]{x,z}

The (Subscript[ρ, 1]) Integrand

Adopted from pensieve://Projects/APAI/PerturbedGaussianIntegration.nb.

q[s_,i_,j_]:=(Subscript[x, i])((1-(T^s))(Subscript[p, i+1])+((T^s)-1)(Subscript[p, j+1]));
(Subscript[r, 1])[s_,i_,j_]:=s(-1+2(Subscript[p, i])(Subscript[x, i])-2(Subscript[p, j])(Subscript[x, i])+((T^s)-1)(Subscript[p, i])(Subscript[p, j])(Subscript[x, i]^2)+(1-(T^s))(Subscript[p, j]^2)(Subscript[x, i]^2)-2(Subscript[p, i])(Subscript[p, j])(Subscript[x, i])(Subscript[x, j])+2(Subscript[p, j]^2)(Subscript[x, i])(Subscript[x, j]))/2;
(Subscript[γ, 1])[φ_,k_]:=ϵ φ(1/2-(Subscript[x, k])(Subscript[p, k]));
ρ1i[s_,i_,j_]:=(T^(-s/2))𝔼[-q[s,i,j]+ϵ (Subscript[r, 1])[s,i,j]+((O[ϵ])^2)];
ρ1i[φ_,k_]:=(T^(-φ/2))𝔼[-(Subscript[x, k]) ((Subscript[p, k])-(Subscript[p, k+1]))+(Subscript[γ, 1])[φ,k]+((O[ϵ])^2)];
ρ1i[End,k_]:=𝔼[-(Subscript[x, k])(Subscript[p, k])+((O[ϵ])^2)];
ρ1i[K_]:=Module[{Cs,φ,n,c,k,ℰ},
{Cs,φ}=Rot[K]; n=Length[Cs];
ℰ=ρ1i[End,2n+1];
Do[ℰ*=ρ1i@@c,{c,Cs}];
Do[ℰ*=ρ1i[φ[[k]],k],{k,2n}];
CF@ℰ
];
ρ1vs[K_] := Union@@Table[{(Subscript[x, i]),(Subscript[p, i])},{i,2Crossings[K]+1}]

ρ1i[Knot[3,1]]

ρ1vs[Knot[3,1]]

Integration of ϵ-Series

Using Picard Iteration!

Unprotect[Integrate];
Integrate::sing="How dare you ask me to integrate a singular Gaussian!";
∫ω_.𝔼[L_SeriesData]\[DifferentialD](vs_List):=Module[{n,m,ϵ,L0,Q,Δ,G,Z,e,λ,a,b},
ϵ=L[[1]]; m=L[[5]];
n=Length@vs; L0=Normal@L/.ϵ->0;
Q=Table[(-(∂_(vs[[a]],vs[[b]]))L0)/.Thread[vs->0]/.(Subscript[(p|x), __])->0,{a,n},{b,n}];
If[(Δ=CF@Det[Q])==0,Message[Integrate::sing]; Return[]];
Z=L+vs.Q.vs/2; G=Inverse[Q];
While[
e=Normal@CF[((∂_λ)Z)-(1/2)(∑_(a=1)^n)(∑_(b=1)^n)(((∂_(vs[[a]],vs[[b]]))Z)+((∂_(vs[[a]]))Z)((∂_(vs[[b]]))Z)+((O[ϵ])^m))];
0=!=e,Z-=(∫_0^λ)e \[DifferentialD]λ
];
PowerExpand@Factor[ω(((Δ(((2π))^n)))^(-1/2))]𝔼[CF[Z/.λ->1/.Thread[vs->0]]]
];
Protect[Integrate];

∫𝔼[-(x^2)/2+δ (x^3)/6+((O[δ])^20)]\[DifferentialD]{x}

∫𝔼[-(ϕ^2)/2+ϵ (ϕ^4)/24+((O[ϵ])^16)]\[DifferentialD]{ϕ}

∫𝔼[p x+ϵ a (p^2)x+((O[ϵ])^5)]\[DifferentialD]{p,x}

∫𝔼[p x+ϵ a (p^2)x+ϵ b p (x^2)+((O[ϵ])^5)]\[DifferentialD]{p,x}

∫𝔼[p x+(p^2)x+((O[ϵ])^5)]\[DifferentialD]{p,x}

∫𝔼[(Subscript[x, 1])(Subscript[p, 1])+ϵ (Subscript[x, 1]^7)(Subscript[p, 1]^7)+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, 1]),(Subscript[p, 1])}

∫𝔼[(Subscript[x, 1])(Subscript[p, 2])+ϵ (Subscript[x, 2]^7)(Subscript[p, 1]^7)+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, 1]),(Subscript[p, 2])}

∫𝔼[(Subscript[x, 1])(Subscript[p, 2])+3(Subscript[x, 2])(Subscript[p, 1])+ϵ (Subscript[p, 2]^5)(Subscript[x, 1]^5)+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, 1]),(Subscript[x, 2]),(Subscript[p, 1]),(Subscript[p, 2])}

∫𝔼[(Subscript[x, 1])(Subscript[p, 2])+(Subscript[x, 2])(Subscript[p, 3])+(Subscript[x, 3])(Subscript[p, 1])+ϵ (Subscript[x, 1]^5)(Subscript[p, 2]^5)+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, 1]),(Subscript[x, 2]),(Subscript[x, 3]),(Subscript[p, 1]),(Subscript[p, 2]),(Subscript[p, 3])}

MatrixForm@Table[
∫𝔼[(Subscript[x, 1])(Subscript[p, 2])+(Subscript[x, 2])(Subscript[p, 3])+(Subscript[x, 3])(Subscript[p, 1])+(Subscript[ξ, i])(Subscript[x, i])+(Subscript[π, j])(Subscript[p, j])+O[ϵ]]\[DifferentialD]{(Subscript[x, 1]),(Subscript[x, 2]),(Subscript[x, 3]),(Subscript[p, 1]),(Subscript[p, 2]),(Subscript[p, 3])},
{i,3},{j,3}]

K=Knot[5,2];
∫ρ1i[K]\[DifferentialD](ρ1vs@K)

K=Knot[8,19];
∫ρ1i@K\[DifferentialD](ρ1vs@K)

Invariance Under Reidemeister 3b

lhsi=ρ1i[1,i,j]ρ1i[1,i+1,k]ρ1i[1,j+1,k+1]ρ1i[0,i]ρ1i[0,j]ρ1i[0,k]ρ1i[0,i+1]ρ1i[0,j+1]ρ1i[0,k+1];
lhs=∫lhsi\[DifferentialD]{(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[x, k+1]),(Subscript[p, i+1]),(Subscript[p, j+1]),(Subscript[p, k+1])}

rhsi=ρ1i[1,j,k]ρ1i[1,i,k+1]ρ1i[1,i+1,j+1]ρ1i[0,i]ρ1i[0,j]ρ1i[0,k]ρ1i[0,i+1]ρ1i[0,j+1]ρ1i[0,k+1];
rhs=∫lhsi\[DifferentialD]{(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[x, k+1]),(Subscript[p, i+1]),(Subscript[p, j+1]),(Subscript[p, k+1])};

lhs==rhs

Invariance Under Reidemeister 2b

lhsi=ρ1i[1,i,j]ρ1i[-1,i+1,j+1]ρ1i[0,i]ρ1i[0,j]ρ1i[0,i+1]ρ1i[0,j+1];
lhs=∫lhsi\[DifferentialD]{(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[p, i+1]),(Subscript[p, j+1])}

CF[((Subscript[p, i]) (Subscript[x, i])-(Subscript[p, j]) (Subscript[x, i])-T (Subscript[p, 2+j]) (Subscript[x, i])+(1/2) (-1+T) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2)+(1/2) (1-T) (Subscript[p, j]^2) (Subscript[x, i]^2)+(1/2) (1-T) (Subscript[p, 2+i]) (Subscript[p, 2+j]) (Subscript[x, i]^2)+(1/2) (-1+T) (Subscript[p, 2+j]^2) (Subscript[x, i]^2)+(Subscript[p, 2+j]) (Subscript[x, j])-(Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j])+(Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j])+(Subscript[p, 2+i]) (Subscript[p, 2+j]) (Subscript[x, i]) (Subscript[x, j])-(Subscript[p, 2+j]^2) (Subscript[x, i]) (Subscript[x, j]))/.(Subscript[p, k_+2])->(Subscript[p, k])]

lhsi=ρ1i[1,i,j]ρ1i[-1,i+1,j+1]ρ1i[0,i]ρ1i[0,j]ρ1i[0,i+1]ρ1i[0,j+1];
lhs=∫𝔼[(Subscript[π, i])(Subscript[p, i])+(Subscript[π, j])(Subscript[p, j])]lhsi\[DifferentialD]{(Subscript[x, i]),(Subscript[x, j]),(Subscript[p, i]),(Subscript[p, j]),(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[p, i+1]),(Subscript[p, j+1])}

lhs=∫𝔼[(Subscript[π, i])(Subscript[p, i])+(Subscript[π, j])(Subscript[p, j])+ℒ[1,i,j]+ℒ[-1,i+1,j+1]+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, i]),(Subscript[x, j]),(Subscript[p, i]),(Subscript[p, j]),(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[p, i+1]),(Subscript[p, j+1])}

rhs=∫𝔼[(Subscript[π, i])(Subscript[p, i])+(Subscript[π, j])(Subscript[p, j])+ℒ[0,i]+ℒ[0,i+1]+ℒ[0,j]+ℒ[0,j+1]+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, i]),(Subscript[x, j]),(Subscript[p, i]),(Subscript[p, j]),(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[p, i+1]),(Subscript[p, j+1])}

lhs==rhs

Invariance Under R2c

lhs=∫𝔼[(Subscript[π, i])(Subscript[p, i])+(Subscript[π, j])(Subscript[p, j])+ℒ[-1,i,j+1]+ℒ[1,i+1,j]+(Subscript[γ, 1])[-1,j+1]+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, i]),(Subscript[x, j]),(Subscript[p, i]),(Subscript[p, j]),(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[p, i+1]),(Subscript[p, j+1])}

rhs=∫𝔼[(Subscript[π, i])(Subscript[p, i])+(Subscript[π, j])(Subscript[p, j])+ℒ[0,i]+ℒ[0,i+1]+ℒ[0,j]+ℒ[-1,j+1]+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, i]),(Subscript[x, j]),(Subscript[p, i]),(Subscript[p, j]),(Subscript[x, i+1]),(Subscript[x, j+1]),(Subscript[p, i+1]),(Subscript[p, j+1])}

lhs==rhs

Invariance Under R1l

lhs=∫𝔼[(Subscript[π, i])(Subscript[p, i])+ℒ[1,i+2,i]+ℒ[1,i+1]+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, i]),(Subscript[x, i+1]),(Subscript[x, i+2]),(Subscript[p, i]),(Subscript[p, i+1]),(Subscript[p, i+2])}

rhs=∫𝔼[(Subscript[π, i])(Subscript[p, i])+ℒ[0,i]+ℒ[0,i+1]+ℒ[0,i+2]+((O[ϵ])^2)]\[DifferentialD]{(Subscript[x, i]),(Subscript[x, i+1]),(Subscript[x, i+2]),(Subscript[p, i]),(Subscript[p, i+1]),(Subscript[p, i+2])}

lhs==rhs