Collect[-(1/2) ϵ (Subscript[ca, 1,2])+ϵ (Subscript[p, i]) (Subscript[x, i]) (Subscript[ca, 1,2])-ϵ (Subscript[p, j]) (Subscript[x, i]) (Subscript[ca, 1,2])+(1/2) ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10])-(1/2) T ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10])-(1/2) ϵ (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10])+(1/2) T ϵ (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10])+ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10])-ϵ (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10]),ϵ,CF]

Collect[-(1/2) ϵ (Subscript[ca, 1,2])-ϵ (Subscript[p, k]) (Subscript[x, k]) (Subscript[ca, 1,10]),ϵ,CF]

{Collect[-(1/2) ϵ (Subscript[ca, 1,2])+ϵ (Subscript[p, i]) (Subscript[x, i]) (Subscript[ca, 1,2])-ϵ (Subscript[p, j]) (Subscript[x, i]) (Subscript[ca, 1,2])-(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[x, i]) (Subscript[ca, 1,2]^2)+(1/2) (ϵ^2) (Subscript[p, j]) (Subscript[x, i]) (Subscript[ca, 1,2]^2)+(1/2) ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10])-(1/2) T ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10])-(1/2) ϵ (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10])+(1/2) T ϵ (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10])+ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10])-ϵ (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10])-(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])+T (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])+(1/2) (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])-T (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])-2 (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])+2 (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])+(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10]^2)-(1/2) T (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10]^2)-(1/2) (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10]^2)+(1/2) T (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10]^2)-(1/3) (ϵ^2) (Subscript[p, i]^2) (Subscript[p, j]) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)+(1/3) T (ϵ^2) (Subscript[p, i]^2) (Subscript[p, j]) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)+(5/6) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)-(2/3) T (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)-(1/6) (T^2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)-(1/2) (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)+(1/3) T (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)+(1/6) (T^2) (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)+(ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10]^2)-(ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10]^2)-(1/2) (ϵ^2) (Subscript[p, i]^2) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)+(ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)+(1/2) T (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)-(1/2) (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)-(1/2) T (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)-(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]^2) (Subscript[ca, 1,10]^2)+(1/2) (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]) (Subscript[x, j]^2) (Subscript[ca, 1,10]^2)+(ϵ^2) (Subscript[ca, 2,1])-2 (ϵ^2) (Subscript[p, i]) (Subscript[x, i]) (Subscript[ca, 2,1])+2 (ϵ^2) (Subscript[p, j]) (Subscript[x, i]) (Subscript[ca, 2,1])-(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[cb, 2,10])+(1/2) T (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[cb, 2,10])+(1/2) (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[cb, 2,10])-(1/2) T (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[cb, 2,10])-(ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j]) (Subscript[cb, 2,10])+(ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]) (Subscript[cb, 2,10]),ϵ,CF],
Collect[(1/2) ϵ (Subscript[ca, 1,2])-ϵ (Subscript[p, i]) (Subscript[x, i]) (Subscript[ca, 1,2])+ϵ (Subscript[p, j]) (Subscript[x, i]) (Subscript[ca, 1,2])-(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[x, i]) (Subscript[ca, 1,2]^2)+(1/2) (ϵ^2) (Subscript[p, j]) (Subscript[x, i]) (Subscript[ca, 1,2]^2)-(1/2) ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10])+((ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,10]))/(2 T))+(1/2) ϵ (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10])-((ϵ (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,10]))/(2 T))-ϵ (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10])+ϵ (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]) (Subscript[ca, 1,10])+(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])-(1/2) (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])-(1/3) (ϵ^2) (Subscript[p, i]^2) (Subscript[p, j]) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)+(((ϵ^2) (Subscript[p, i]^2) (Subscript[p, j]) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2))/(3 T))+(5/6) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)-(((ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2))/(6 (T^2)))-((2 (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2))/(3 T))-(1/2) (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2)+(((ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2))/(6 (T^2)))+(((ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^3) (Subscript[ca, 1,10]^2))/(3 T))-(1/2) (ϵ^2) (Subscript[p, i]^2) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)+(ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)+(((ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2))/(2 T))-(1/2) (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2)-(((ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]^2) (Subscript[x, j]) (Subscript[ca, 1,10]^2))/(2 T))-(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]^2) (Subscript[ca, 1,10]^2)+(1/2) (ϵ^2) (Subscript[p, j]^3) (Subscript[x, i]) (Subscript[x, j]^2) (Subscript[ca, 1,10]^2)-(ϵ^2) (Subscript[ca, 2,1])+2 (ϵ^2) (Subscript[p, i]) (Subscript[x, i]) (Subscript[ca, 2,1])-2 (ϵ^2) (Subscript[p, j]) (Subscript[x, i]) (Subscript[ca, 2,1])+(1/2) (ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[cb, 2,10])-(((ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]^2) (Subscript[cb, 2,10]))/(2 T))-(1/2) (ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[cb, 2,10])+(((ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]^2) (Subscript[cb, 2,10]))/(2 T))+(ϵ^2) (Subscript[p, i]) (Subscript[p, j]) (Subscript[x, i]) (Subscript[x, j]) (Subscript[cb, 2,10])-(ϵ^2) (Subscript[p, j]^2) (Subscript[x, i]) (Subscript[x, j]) (Subscript[cb, 2,10]),ϵ,CF]
} /. {(Subscript[ca, 2,1])->0,(Subscript[cb, 2,10])->3/2}

Column[{
Collect[(1/2) ϵ (Subscript[ca, 1,2])+ϵ (Subscript[p, k]) (Subscript[x, k]) (Subscript[ca, 1,10])-(ϵ^2) (Subscript[p, k]) (Subscript[x, k]) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])-(ϵ^2) (Subscript[ca, 2,1])-(ϵ^2) (Subscript[p, k]) (Subscript[x, k]) (Subscript[cb, 2,10]),ϵ,CF],
Collect[-(1/2) ϵ (Subscript[ca, 1,2])-ϵ (Subscript[p, k]) (Subscript[x, k]) (Subscript[ca, 1,10])+(ϵ^2) (Subscript[p, k]) (Subscript[x, k]) (Subscript[ca, 1,2]) (Subscript[ca, 1,10])-(ϵ^2) (Subscript[p, k]) (Subscript[x, k]) (Subscript[ca, 1,10]^2)+(ϵ^2) (Subscript[ca, 2,1])+(ϵ^2) (Subscript[p, k]) (Subscript[x, k]) (Subscript[cb, 2,10]),ϵ,CF]
}/. {(Subscript[ca, 2,1])->0,(Subscript[cb, 2,10])->3/2}]