Contractions of AO basis operators in an asymmetric Fermi vacuum
The overall matrix element in an asymmetric Fermi vacuum
Single contractions
To show how a general string of operators in an asymmetric Fermi vacuum can be evaluated using contractions, we first calculate an element of the one-body density matrix

Look at the most general case, where there are

The next step requires reversing the order of summation. This rewrites any summation as

Here, we moved
We now consider the terms in square brackets. Due to the definition of the extended Thouless transformation,

We can simplify the binomail expansions due to the relation

Working out the Kronecker deltas gives us

Due to the relation between the two possible contractions of the

Multiple contractions
Let us now consider a two-body reduced co-density matrix element

We will work out both contractions individually and add them together at the end. Let us start with the first one. The first step is to apply the extended Thouless transformation

in which we can now use the contractions between the

The summations in square brackets can once again be written as binomials, which in turn can be simplified to Kronecker deltas by using the relation defined earlier. This gives us

The other contraction can be analogously derived. One has to take a phase factor in account as the two contractions cross each other. This gives

Adding the two contractions together gives us the solution for the matrix element