SnegSimplify

◼ SnegSimplify[expr] attempts to rewrite an expression in terms of higher level sneg functions.

◼ SnegSimplify handles expressions featuring number, Hubbard, spin and hopping operators.

◼ See also: sumSimplifyKD, sumSimplify, SimplifyKD.

Further examples

Load sneg and declare c and d to be fermionic operators.

<< sneg`sneg`;
snegfermionoperators[c, d];

number[c[]]
SnegSimplify[%]

hubbard[c[]]
SnegSimplify[%]

hop[c[], d[]]
SnegSimplify[%]

spinx[c[]]
SnegSimplify[%]
spiny[c[]]
SnegSimplify[%]
spinz[c[]]
SnegSimplify[%]

spinx[c[]]·spiny[d[]]
SnegSimplify[%]

spinspin[c[], d[]]
SnegSimplify[%]
spinspinxy[c[], d[]]
SnegSimplify[%] // Expand
spinspinz[c[], d[]]
SnegSimplify[%]

Required factorizations are performed automatically

2 number[c[]]
SnegSimplify[%]
-2/3 number[c[]] // Expand
SnegSimplify[%]

Three-site Hubbard Hamiltonian

n = 3;
H = hamiltonian[Hubbard, c[], n, U, t, ϵ]
SnegSimplify[H]