Pensieve header: Normally Ordered Exponentials at 1-Co.

"NO” for “Normal Order”.

DP[P_, x_->(Subscript[D, α_]),y_->(Subscript[D, β_])][f_] := Total[CoefficientRules[P, {x,y}] /. ({m_,n_}->c_) :> c D[f,{α,m},{β,n}]]

ϵ /: (ϵ^(n_))/;n≥1:=0;
Λ[b_, c_,u_,w_,α_,β_,δ_,ν_]:=2 c w α δ ν+2 c u β δ ν+2 c u w (δ^2) ν+u (w^2) α (δ^2) (ν^3)-(1/2) b (u^2) (w^2) (δ^4) (ν^3)+(1/2) (w^2) (α^2) δ (2+b δ) (ν^3)-(u^2) w β (δ^2) (1+2 b δ) (ν^3)-(1/2) (u^2) (β^2) δ (2+3 b δ) (ν^3)+2 c (δ+α β ν)-2 b u w (δ^2) (ν^2) (δ+α β ν)+w α (ν^2) (2 δ+α β ν)-u β (1+2 b δ) (ν^2) (2 δ+α β ν)-(1/2) b ν (2 (δ^2)+4 α β δ ν+(α^2) (β^2) (ν^2));

Simp[P_.𝔼[Q_]] := Together[P]𝔼[Together[Q]];
(*P1_.𝔼[Q1_] ≡ P2_.𝔼[Q2_] := Simplify[P1==P2∧Q1==Q2];*)
P1_.𝔼[Q1_] ≡ P2_.𝔼[Q2_] := Together[P1-P2]==0∧Together[Q1-Q2]==0;

𝔼 /: 𝔼[Q1_]𝔼[Q2_] := Simp@𝔼[Q1+Q2]

NO[(Subscript[u, i_]),(Subscript[c, j_]),k_][P_.𝔼[Q_]] := Simp@Module[{q(*,α,β,Θ*)},
q=(E^(-α))β (Subscript[u, k])+α (Subscript[c, k])+Θ;
(E^(-q))DP[P,(Subscript[c, j])->(Subscript[D, α]),(Subscript[u, i])->(Subscript[D, β])][(E^q)]𝔼[q] /. {α->Together[(∂_(Subscript[c, j]))Q],β->(∂_(Subscript[u, i]))Q,Θ->(Q/.(Subscript[c, j])|(Subscript[u, i])->0)}
];

NO[(Subscript[w, i_]),(Subscript[c, j_]),k_][P_.𝔼[Q_]] := Simp@Module[{q(*,α,β,Θ*)},
q=(E^α)β (Subscript[w, k])+α (Subscript[c, k])+Θ;
(E^(-q))DP[P,(Subscript[c, j])->(Subscript[D, α]),(Subscript[w, i])->(Subscript[D, β])][(E^q)]𝔼[q] /. {α->Together[(∂_(Subscript[c, j]))Q],β->(∂_(Subscript[w, i]))Q,Θ->(Q/.(Subscript[c, j])|(Subscript[w, i])->0)}
];

NO[(Subscript[w, i_]),(Subscript[u, j_]),k_][P_.𝔼[Q_]] := Simp@Module[
{α0=(∂_(Subscript[w, i]))Q/.(Subscript[u, j])->0, β0=(∂_(Subscript[u, j]))Q/.(Subscript[w, i])->0, δ0=(∂_((Subscript[w, i]),(Subscript[u, j])))Q,Θ0=Q/.(Subscript[w, i])|(Subscript[u, j])->0, q(*,α,β,δ, Θ, ν*)},
q=-(Subscript[b, k])ν α β+ν β (Subscript[u, k])+ν δ (Subscript[u, k])(Subscript[w, k])+ν α (Subscript[w, k])+Θ;
(E^(-q))DP[P,(Subscript[w, i])->(Subscript[D, α]),(Subscript[u, j])->(Subscript[D, β])][ν(1+ϵ ν Λ[(Subscript[b, k]),(Subscript[c, k]), (Subscript[u, k]), (Subscript[w, k]), α,β,δ,ν]) (E^q)]𝔼[q]/. {α->α0,β->β0,δ->δ0, Θ->Θ0, ν->(((1+(Subscript[b, k]) δ0))^(-1))}
];

m[i_,j_,kk_][Z_] :=Simp@Module[{x,y},
Z // ReplaceAll[(Subscript[b, i|j])->(Subscript[b, kk])]// NO[(Subscript[w, i]),(Subscript[c, j]),x] // NO[(Subscript[w, x]),(Subscript[u, j]),y] // ReplaceAll[{(Subscript[c, x|y])->(Subscript[c, x]),(Subscript[w, j])->(Subscript[w, y])}] // NO[(Subscript[u, i]),(Subscript[c, x]),x] // ReplaceAll[{(Subscript[c, i|x])->(Subscript[c, kk]),(Subscript[u, x|y])->(Subscript[u, kk]),(Subscript[w, y])->(Subscript[w, kk]),(Subscript[b, x|y])->(Subscript[b, kk])}]
]

Rp[i_,j_]:=𝔼[(Subscript[b, i])(Subscript[c, j])+(((E^(Subscript[b, i]))-1)/(Subscript[b, i]))(Subscript[u, i])(Subscript[w, j])]; Rm[i_,j_]:=𝔼[-(Subscript[b, i])(Subscript[c, j])+(((E^(-(Subscript[b, i])))-1)/(Subscript[b, i]))(Subscript[u, i])(Subscript[w, j])];

Q0=𝔼[ (Subscript[u, 1]) (Subscript[w, 1])+ (Subscript[u, 2]) (Subscript[w, 3])];
t1=Q0 // m[1,2,4]

NO[(Subscript[u, 1]),(Subscript[c, 1]),1][𝔼[(Subscript[b, 2])(Subscript[c, 1])+(Subscript[w, 2])(Subscript[u, 1])+(Subscript[u, 3])(Subscript[w, 4])]]

Q0=𝔼[(Subscript[u, 1]) (Subscript[w, 1])+ (Subscript[u, 3]) (Subscript[w, 2])];

SeedRandom[0];
Q0=𝔼[Sum[RandomInteger[{-2,2}](Subscript[c, i]),{i,3}]+Sum[RandomInteger[{-2,2}](Subscript[u, i])(Subscript[w, j]),{i,3},{j,3}]]

Q0=𝔼[Sum[(Subscript[β, i])(Subscript[c, i]),{i,3}]+Sum[(Subscript[α, i,j])(Subscript[u, i])(Subscript[w, j]),{i,3},{j,3}]]

Q0=𝔼[-2 (Subscript[b, 1]) (Subscript[c, 1])-(Subscript[b, 2]) (Subscript[c, 1])+2 (Subscript[b, 3]) (Subscript[c, 1])+(Subscript[b, 2]) (Subscript[c, 2])-(Subscript[b, 3]) (Subscript[c, 2])+2 (Subscript[b, 1]) (Subscript[c, 3])+(Subscript[b, 2]) (Subscript[c, 3])-(Subscript[b, 3]) (Subscript[c, 3])-(Subscript[u, 3]) (Subscript[w, 1])+2 (Subscript[u, 2]) (Subscript[w, 2])+(Subscript[u, 3]) (Subscript[w, 2])+(Subscript[u, 1]) (Subscript[w, 3])-(Subscript[u, 2]) (Subscript[w, 3])]

t1=Q0 // m[1,2,4]//m[4,3,5]
t2=Q0 // m[2,3,4]//m[1,4,5]
t3=(t1≡t2)

Rp[1,2]Rm[4,3]

Rp[1,2]Rm[4,3] // m[1,4,1] // m[2,3,2]

Rm[1,2]Rp[4,3] // m[4,1,4]//m[2,3,2]

t1=Rp[1,2]Rp[3,4]Rp[5,6] // m[3,5,a]//m[1,6,b]//m[2,4,c]

t2=Rp[1,2]Rp[3,4]Rp[5,6]//m[1,3,a]//m[2,5,b]//m[4,6,c]

t1≡t2

t3=Rm[12,1]Rm[2,7]Rm[8,3]Rm[4,11]Rp[16,5]Rp[6,13]Rp[14,9]Rp[10,15]

t3//m[1,2,1]

t3//m[1,2,1]//m[1,3,1]

Do[t3=t3//m[1,kk,1], {kk, 2,16}]; t3

Rp[1,2]Rp[3,4]Rp[5,6] // m[1,4,4]//m[4,5,5]//m[5,2,2]//m[2,3,3]//m[3,6,6]//m[6,1,1]