$\renewcommand{\bm}{\boldsymbol} \require{mediawiki-texvc}$

# Spin-Adapted DMRG Quantum Chemistry Hamiltonian¶

## Partitioning in SU(2)¶

The partitioning of Hamiltonian in left ($$L$$) and right ($$R$$) blocks is given by

$\begin{split}(\hat{H})^{[0]} =&\ \big( \hat{H}^{L} \big)^{[0]} \otimes_{[0]} \big( \hat{1}^{R} \big)^{[0]} + \big( \hat{1}^{L} \big)^{[0]} \otimes_{[0]} \big( \hat{H}^{R} \big)^{[0]} \\ &\ + \sqrt{2} \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{S}_{i}^{R} \big)^{[\frac{1}{2}]} + h.c. \right] \\ &\ + 2 \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{R} \big)^{[\frac{1}{2}]} + h.c. \right] + 2 \sum_{i\in R} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{L} \big)^{[\frac{1}{2}]} + h.c. \right] \\ &\ - \frac{1}{2} \sum_{ik\in L} \left[ \sqrt{3} \big(\hat{A}_{ik} \big)^{[1]} \otimes_{[0]} \big(\hat{P}_{ik}^{R} \big)^{[1]} + \big(\hat{A}_{ik} \big)^{[0]} \otimes_{[0]} \big(\hat{P}_{ik}^{R} \big)^{[0]} + h.c. \right] \\ &\ +\sum_{ij\in L} \left[ \big( \hat{B}_{ij} \big)^{[0]} \otimes_{[0]} \left( 2\big( \hat{Q}_{ij}^{R} \big)^{[0]} - \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[0]} \right) + \sqrt{3} \big( {\hat{B}'}_{ij} \big)^{[1]} \otimes_{[0]} \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[1]} \right]\end{split}$

where the normal and complementary operators are defined by

$\begin{split}\big( \hat{S}_{i}^{L/R} \big)^{[\frac{1}{2}]} =&\ \sum_{j\in L/R} t_{ij} \big( a_{j} \big)^{[\frac{1}{2}]} \\ \big( \hat{R}_{i}^{L/R} \big)^{[\frac{1}{2}]} =&\ \sum_{jkl\in L/R} v_{ijkl} \left[ \Big( a_{k}^\dagger \Big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \otimes_{[\frac{1}{2}]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ \big( \hat{A}_{ik} \big)^{[0/1]} =&\ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \\ \big( \hat{P}_{ik}^{R} \big)^{[0/1]} =&\ \sum_{jl\in R} v_{ijkl} \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{l} \big)^{[\frac{1}{2}]} \\ \big( \hat{B}_{ij} \big)^{[0]} =&\ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ \big( {\hat{B}'}_{ij} \big)^{[1]} =&\ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{j} \big)^{[\frac{1}{2}]}\\ \big( \hat{Q}_{ij}^{R} \big)^{[0]} =&\ \sum_{kl\in R} v_{ijkl} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \\ \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[0/1]} =&\ \sum_{kl\in R} v_{ilkj} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{l} \big)^{[\frac{1}{2}]} \\ \big( {\hat{Q}}_{ij}^{\prime \prime R} \big)^{[0]} :=&\ 2 \big( {\hat{Q}}_{ij}^{R} \big)^{[0]} - \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[0]} = \sum_{kl\in R} (2v_{ijkl} - v_{ilkj}) \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]}\end{split}$

### Derivation¶

#### CG Factors¶

From $$j_2 = 1/2$$ CG factors

$\begin{split}\bigg\langle j_1\ \left(M - \frac{1}{2} \right)\ \frac{1}{2}\ \frac{1}{2} \bigg| \left( j_1 \pm \frac{1}{2} \right)\ M \bigg\rangle =&\ \pm \sqrt{\frac{1}{2} \left( 1 \pm \frac{M}{j_1 + \frac{1}{2}} \right)} \\ \bigg\langle j_1\ \left(M + \frac{1}{2} \right)\ \frac{1}{2}\ \left( -\frac{1}{2}\right) \bigg| \left( j_1 \pm \frac{1}{2} \right)\ M \bigg\rangle =&\ \sqrt{\frac{1}{2} \left( 1 \mp \frac{M}{j_1 + \frac{1}{2}} \right)}\end{split}$

and symmetry relation

$\langle j_1\ m_1\ j_2\ m_2 |J\ M\rangle = (-1)^{j_1+j_2-J} \langle j_2\ m_2\ j_1\ m_1 |J\ M\rangle$

and

$(-1)^{j_1+\frac{1}{2}-j_1\mp\frac{1}{2}} = (-1)^{\frac{1}{2}\mp\frac{1}{2}} = \pm 1$

we have

$\begin{split}\bigg\langle \frac{1}{2}\ \frac{1}{2}\ j_1\ \left(M - \frac{1}{2} \right) \bigg| \left( j_1 \pm \frac{1}{2} \right)\ M \bigg\rangle =&\ \sqrt{\frac{1}{2} \left( 1 \pm \frac{M}{j_1 + \frac{1}{2}} \right)} \\ \bigg\langle \frac{1}{2}\ \left( -\frac{1}{2}\right)\ j_1\ \left(M + \frac{1}{2} \right) \bigg| \left( j_1 \pm \frac{1}{2} \right)\ M \bigg\rangle =&\ \pm \sqrt{\frac{1}{2} \left( 1 \mp \frac{M}{j_1 + \frac{1}{2}} \right)}\end{split}$

let $$j_1 = 1$$, we have

$\begin{split}\langle \tfrac{1}{2}\ \tfrac{1}{2}\ 1\ (M - \tfrac{1}{2}) | \tfrac{1}{2}\ M \rangle =&\ \sqrt{\tfrac{1}{2} ( 1-\frac{M}{\tfrac{3}{2}} )} \\ \langle \tfrac{1}{2}\ (-\tfrac{1}{2})\ 1\ (M + \tfrac{1}{2}) | \tfrac{1}{2}\ M \rangle =&\ -\sqrt{\tfrac{1}{2} ( 1+\frac{M}{\tfrac{3}{2}} )}\end{split}$

So the coefficients for $$[\tfrac{1}{2}] \otimes_{[\tfrac{1}{2}]} [1]$$ are

$\begin{split}[\tfrac{1}{2} + 0 = \tfrac{1}{2}] = \sqrt{\tfrac{1}{3}},\quad [-\tfrac{1}{2} + 1 = \tfrac{1}{2}] = -\sqrt{\tfrac{2}{3}} \\ [\tfrac{1}{2} + (-1) = -\tfrac{1}{2}] = \sqrt{\tfrac{2}{3}},\quad [-\tfrac{1}{2} + 0 = -\tfrac{1}{2}] = -\sqrt{\tfrac{1}{3}}\end{split}$

The coefficients for $$[1] \otimes_{[\tfrac{1}{2}]} [\tfrac{1}{2}]$$ are

$\begin{split}[0 + \tfrac{1}{2} = \tfrac{1}{2}] = -\sqrt{\tfrac{1}{3}},\quad [1 -\tfrac{1}{2} = \tfrac{1}{2}] = \sqrt{\tfrac{2}{3}} \\ [(-1) + \tfrac{1}{2} = -\tfrac{1}{2}] = -\sqrt{\tfrac{2}{3}},\quad [0 -\tfrac{1}{2} = -\tfrac{1}{2}] = \sqrt{\tfrac{1}{3}}\end{split}$

This means that the SU(2) operator exchange factor for $$[\tfrac{1}{2}] \otimes_{[\tfrac{1}{2}]} [1] \to [1] \otimes_{[\tfrac{1}{2}]} [\tfrac{1}{2}]$$ is $$-1$$. The fermion factor is $$+1$$. So the overall exchange factor for this case is $$-1$$.

#### Tensor Product Formulas¶

Singlet

$\begin{split}\big(a_p^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_q^\dagger\big)^{[1/2]} =&\ \begin{pmatrix} a_{p\alpha}^\dagger \\ a_{p\beta}^\dagger \end{pmatrix}^{[1/2]} \otimes_{[0]} \begin{pmatrix} a_{q\alpha}^\dagger \\ a_{q\beta}^\dagger \end{pmatrix}^{[1/2]} = \frac{1}{\sqrt{2}} \begin{pmatrix} a_{p\alpha}^\dagger a_{q\beta}^\dagger - a_{p\beta}^\dagger a_{q\alpha}^\dagger \end{pmatrix}^{[0]} \\ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_q\big)^{[1/2]} =&\ \begin{pmatrix} a_{p\alpha}^\dagger \\ a_{p\beta}^\dagger \end{pmatrix}^{[1/2]} \otimes_{[0]} \begin{pmatrix} -a_{q\beta} \\ a_{q\alpha} \end{pmatrix}^{[1/2]} = \frac{1}{\sqrt{2}} \begin{pmatrix} a_{p\alpha}^\dagger a_{q\alpha}+ a_{p\beta}^\dagger a_{q\beta} \end{pmatrix}^{[0]} \\ \big(a_p\big)^{[1/2]} \otimes_{[0]} \big(a_q\big)^{[1/2]} =&\ \begin{pmatrix} -a_{p\beta} \\ a_{p\alpha} \end{pmatrix}^{[1/2]} \otimes_{[0]} \begin{pmatrix} -a_{q\beta} \\ a_{q\alpha} \end{pmatrix}^{[1/2]} = \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{p\beta} a_{q\alpha} + a_{p\alpha} a_{q\beta} \end{pmatrix}^{[0]}\end{split}$

Triplet

$\begin{split}\big(a_p^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_q^\dagger\big)^{[1/2]} =&\ \begin{pmatrix} a_{p\alpha}^\dagger \\ a_{p\beta}^\dagger \end{pmatrix}^{[1/2]} \otimes_{[1]} \begin{pmatrix} a_{q\alpha}^\dagger \\ a_{q\beta}^\dagger \end{pmatrix}^{[1/2]} = \begin{pmatrix} a_{p\alpha}^\dagger a_{q\alpha}^\dagger \\ \frac{1}{\sqrt{2}} \Big( a_{p\alpha}^\dagger a_{q\beta}^\dagger + a_{p\beta}^\dagger a_{q\alpha}^\dagger \Big) \\ a_{p\beta}^\dagger a_{q\beta}^\dagger \end{pmatrix}^{[1]} \\ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_q\big)^{[1/2]} =&\ \begin{pmatrix} a_{p\alpha}^\dagger \\ a_{p\beta}^\dagger \end{pmatrix}^{[1/2]} \otimes_{[1]} \begin{pmatrix} -a_{q\beta} \\ a_{q\alpha} \end{pmatrix}^{[1/2]} = \begin{pmatrix} -a_{p\alpha}^\dagger a_{q\beta} \\ \frac{1}{\sqrt{2}} \Big( a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta} \Big) \\ a_{p\beta}^\dagger a_{q\alpha} \end{pmatrix}^{[1]} \\ \big(a_p\big)^{[1/2]} \otimes_{[1]} \big(a_q\big)^{[1/2]} =&\ \begin{pmatrix} -a_{p\beta} \\ a_{p\alpha} \end{pmatrix}^{[1/2]} \otimes_{[1]} \begin{pmatrix} -a_{q\beta} \\ a_{q\alpha} \end{pmatrix}^{[1/2]} = \begin{pmatrix} a_{p\beta} a_{q\beta} \\ -\frac{1}{\sqrt{2}} \Big( a_{p\beta} a_{q\alpha} + a_{p\alpha} a_{q\beta} \Big) \\ a_{p\alpha} a_{q\alpha} \end{pmatrix}^{[1]}\end{split}$

Doublet times singlet/triplet

$\begin{split}U^{[1/2]} = &\ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_r\big)^{[1/2]} \otimes_{[1]} \big(a_s\big)^{[1/2]} \Big] = \begin{pmatrix} a_{p\alpha}^\dagger \\ a_{p\beta}^\dagger \end{pmatrix}^{[1/2]} \otimes_{[1/2]} \begin{pmatrix} a_{r\beta} a_{s\beta} \\ -\frac{1}{\sqrt{2}} \Big( a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \Big) \\ a_{r\alpha} a_{s\alpha} \end{pmatrix}^{[1]} \\ =&\ \begin{pmatrix} -\frac{1}{\sqrt{2}}\frac{1}{\sqrt{3}} a_{p\alpha}^\dagger \Big( a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \Big) -\frac{\sqrt{2}}{\sqrt{3}} a_{p\beta}^\dagger a_{r\beta} a_{s\beta} \\ \frac{\sqrt{2}}{\sqrt{3}} a_{p\alpha}^\dagger a_{r\alpha} a_{s\alpha} +\big( -\frac{1}{\sqrt{3}}\big) \big( -\frac{1}{\sqrt{2}} \big) a_{p\beta}^\dagger \Big( a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \Big) \end{pmatrix}^{[1/2]} = \frac{1}{\sqrt{6}} \begin{pmatrix} - a_{p\alpha}^\dagger a_{r\beta} a_{s\alpha} - a_{p\alpha}^\dagger a_{r\alpha} a_{s\beta} -2 a_{p\beta}^\dagger a_{r\beta} a_{s\beta} \\ 2 a_{p\alpha}^\dagger a_{r\alpha} a_{s\alpha} +a_{p\beta}^\dagger a_{r\beta} a_{s\alpha} + a_{p\beta}^\dagger a_{r\alpha} a_{s\beta} \end{pmatrix}^{[1/2]} \\ V^{[1/2]} =&\ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_r\big)^{[1/2]} \otimes_{[0]} \big(a_s\big)^{[1/2]} \Big] = \frac{1}{\sqrt{2}} \begin{pmatrix} a_{p\alpha}^\dagger \\ a_{p\beta}^\dagger \end{pmatrix}^{[1/2]} \otimes_{[1/2]} \begin{pmatrix} -a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \end{pmatrix}^{[0]} \\ =&\ \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{p\alpha}^\dagger a_{r\beta} a_{s\alpha} + a_{p\alpha}^\dagger a_{r\alpha} a_{s\beta}\\ -a_{p\beta}^\dagger a_{r\beta} a_{s\alpha} + a_{p\beta}^\dagger a_{r\alpha} a_{s\beta}\end{pmatrix}^{[1/2]}\end{split}$

Therefore,

$\begin{split}\sqrt{3} U^{[1/2]} - V^{[1/2]} =&\ \frac{1}{\sqrt{2}} \begin{pmatrix} - a_{p\alpha}^\dagger a_{r\beta} a_{s\alpha} - a_{p\alpha}^\dagger a_{r\alpha} a_{s\beta} -2 a_{p\beta}^\dagger a_{r\beta} a_{s\beta} \\ 2 a_{p\alpha}^\dagger a_{r\alpha} a_{s\alpha} +a_{p\beta}^\dagger a_{r\beta} a_{s\alpha} + a_{p\beta}^\dagger a_{r\alpha} a_{s\beta} \end{pmatrix}^{[1/2]} - \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{p\alpha}^\dagger a_{r\beta} a_{s\alpha} + a_{p\alpha}^\dagger a_{r\alpha} a_{s\beta}\\ -a_{p\beta}^\dagger a_{r\beta} a_{s\alpha} + a_{p\beta}^\dagger a_{r\alpha} a_{s\beta}\end{pmatrix}^{[1/2]} \\ =&\ \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{p\alpha}^\dagger a_{r\beta} a_{s\alpha} - a_{p\alpha}^\dagger a_{r\alpha} a_{s\beta} -2 a_{p\beta}^\dagger a_{r\beta} a_{s\beta} +a_{p\alpha}^\dagger a_{r\beta} a_{s\alpha} - a_{p\alpha}^\dagger a_{r\alpha} a_{s\beta}\\ 2 a_{p\alpha}^\dagger a_{r\alpha} a_{s\alpha} +a_{p\beta}^\dagger a_{r\beta} a_{s\alpha} + a_{p\beta}^\dagger a_{r\alpha} a_{s\beta} +a_{p\beta}^\dagger a_{r\beta} a_{s\alpha} - a_{p\beta}^\dagger a_{r\alpha} a_{s\beta}\end{pmatrix}^{[1/2]} \\ =&\ \sqrt{2} \begin{pmatrix} - a_{p\alpha}^\dagger a_{r\alpha} a_{s\beta} - a_{p\beta}^\dagger a_{r\beta} a_{s\beta} \\ a_{p\alpha}^\dagger a_{r\alpha} a_{s\alpha} + a_{p\beta}^\dagger a_{r\beta} a_{s\alpha} \end{pmatrix}^{[1/2]}\end{split}$

Another case

$\begin{split}S^{[1/2]} = &\ \big(a_r\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_p^\dagger \big)^{[1/2]} \otimes_{[1]} \big(a_q\big)^{[1/2]} \Big] = \begin{pmatrix} -a_{r\beta} \\ a_{r\alpha} \end{pmatrix}^{[1/2]} \otimes_{[1/2]} \begin{pmatrix} -a_{p\alpha}^\dagger a_{q\beta} \\ \frac{1}{\sqrt{2}} \Big( a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta} \Big) \\ a_{p\beta}^\dagger a_{q\alpha} \end{pmatrix}^{[1]} \\ =&\ \begin{pmatrix} \frac{1}{\sqrt{2}} \frac{1}{\sqrt{3}} (-a_{r\beta}) \Big( a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta} \Big) +\frac{\sqrt{2}}{\sqrt{3}} a_{r\alpha} a_{p\alpha}^\dagger a_{q\beta} \\ -\frac{\sqrt{2}}{\sqrt{3}} a_{r\beta} a_{p\beta}^\dagger a_{q\alpha} -\frac{1}{\sqrt{2}} \frac{1}{\sqrt{3}} a_{r\alpha} \Big( a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta} \Big) \end{pmatrix}^{[1/2]} = \frac{1}{\sqrt{6}} \begin{pmatrix} -a_{r\beta} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\beta} a_{p\beta}^\dagger a_{q\beta} +2 a_{r\alpha} a_{p\alpha}^\dagger a_{q\beta}\\ -2a_{r\beta} a_{p\beta}^\dagger a_{q\alpha} -a_{r\alpha} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\alpha} a_{p\beta}^\dagger a_{q\beta} \end{pmatrix}^{[1/2]} \\ T^{[1/2]} = &\ \big(a_r\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_p^\dagger \big)^{[1/2]} \otimes_{[0]} \big(a_q\big)^{[1/2]} \Big] = \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{r\beta} \\ a_{r\alpha} \end{pmatrix}^{[1/2]} \otimes_{[1/2]} \begin{pmatrix} a_{p\alpha}^\dagger a_{q\alpha}+ a_{p\beta}^\dagger a_{q\beta} \end{pmatrix}^{[0]} \\ =&\ \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{r\beta} a_{p\alpha}^\dagger a_{q\alpha} - a_{r\beta}a_{p\beta}^\dagger a_{q\beta} \\ a_{r\alpha} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\alpha}a_{p\beta}^\dagger a_{q\beta}\end{pmatrix}^{[1/2]}\end{split}$

Therefore,

$\begin{split}\sqrt{3} S^{[1/2]} - T^{[1/2]} =&\ \frac{1}{\sqrt{6}} \begin{pmatrix} -a_{r\beta} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\beta} a_{p\beta}^\dagger a_{q\beta} +2 a_{r\alpha} a_{p\alpha}^\dagger a_{q\beta}\\ -2a_{r\beta} a_{p\beta}^\dagger a_{q\alpha} -a_{r\alpha} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\alpha} a_{p\beta}^\dagger a_{q\beta} \end{pmatrix}^{[1/2]}-\frac{1}{\sqrt{2}} \begin{pmatrix} -a_{r\beta} a_{p\alpha}^\dagger a_{q\alpha} - a_{r\beta}a_{p\beta}^\dagger a_{q\beta} \\ a_{r\alpha} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\alpha}a_{p\beta}^\dagger a_{q\beta}\end{pmatrix}^{[1/2]} \\ =&\ \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{r\beta} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\beta} a_{p\beta}^\dagger a_{q\beta} +2 a_{r\alpha} a_{p\alpha}^\dagger a_{q\beta} +a_{r\beta} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\beta}a_{p\beta}^\dagger a_{q\beta} \\ -2a_{r\beta} a_{p\beta}^\dagger a_{q\alpha} -a_{r\alpha} a_{p\alpha}^\dagger a_{q\alpha} + a_{r\alpha} a_{p\beta}^\dagger a_{q\beta} -a_{r\alpha} a_{p\alpha}^\dagger a_{q\alpha} - a_{r\alpha}a_{p\beta}^\dagger a_{q\beta} \end{pmatrix}^{[1/2]} \\ =&\ \sqrt{2} \begin{pmatrix} a_{r\beta}a_{p\beta}^\dagger a_{q\beta} +a_{r\alpha} a_{p\alpha}^\dagger a_{q\beta} \\ -a_{r\beta} a_{p\beta}^\dagger a_{q\alpha} -a_{r\alpha} a_{p\alpha}^\dagger a_{q\alpha} \end{pmatrix}^{[1/2]}\end{split}$

Triplet times triplet

$\begin{split}X^{[0]} = &\ \Big[ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_q^\dagger\big)^{[1/2]} \Big] \otimes_{[0]} \Big[ \big(a_r\big)^{[1/2]} \otimes_{[1]} \big(a_s\big)^{[1/2]} \Big] \\ =&\ \begin{pmatrix} a_{p\alpha}^\dagger a_{q\alpha}^\dagger \\ \frac{1}{\sqrt{2}} \Big( a_{p\alpha}^\dagger a_{q\beta}^\dagger + a_{p\beta}^\dagger a_{q\alpha}^\dagger \Big) \\ a_{p\beta}^\dagger a_{q\beta}^\dagger \end{pmatrix}^{[1]} \otimes_{[0]} \begin{pmatrix} a_{r\beta} a_{s\beta} \\ -\frac{1}{\sqrt{2}} \Big( a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \Big) \\ a_{r\alpha} a_{s\alpha} \end{pmatrix}^{[1]} \\ =&\ \frac{1}{\sqrt{3}} \begin{pmatrix} a_{p\alpha}^\dagger a_{q\alpha}^\dagger a_{r\alpha} s_{s\alpha} + \frac{1}{2} \Big( a_{p\alpha}^\dagger a_{q\beta}^\dagger + a_{p\beta}^\dagger a_{q\alpha}^\dagger \Big) \Big( a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \Big) + a_{p\beta}^\dagger a_{q\beta}^\dagger a_{r\beta} a_{s\beta} \end{pmatrix} \\ Y^{[0]} = &\ \Big[ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_q^\dagger\big)^{[1/2]} \Big] \otimes_{[0]} \Big[ \big(a_r\big)^{[1/2]} \otimes_{[0]} \big(a_s\big)^{[1/2]} \Big] \\ =&\ \frac{1}{\sqrt{2}} \begin{pmatrix} a_{p\alpha}^\dagger a_{q\beta}^\dagger - a_{p\beta}^\dagger a_{q\alpha}^\dagger \end{pmatrix}^{[0]} \otimes_{[0]} \frac{1}{\sqrt{2}} \begin{pmatrix} -a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \end{pmatrix}^{[0]} \\ =&\ \frac{1}{2} \Big( a_{p\alpha}^\dagger a_{q\beta}^\dagger - a_{p\beta}^\dagger a_{q\alpha}^\dagger \Big) \Big( -a_{r\beta} a_{s\alpha} + a_{r\alpha} a_{s\beta} \Big)\end{split}$

Using

$(a+b)(c+d) + (a-b)(-c+d) = (a+b)(2d) -2b(-c+d) = 2 (ad+bc)$

we have

$\begin{split}\sqrt{3} X^{[0]} + Y^{[0]} =&\ a_{p\alpha}^\dagger a_{q\alpha}^\dagger a_{r\alpha} s_{s\alpha} + a_{p\beta}^\dagger a_{q\beta}^\dagger a_{r\beta} a_{s\beta} + a_{p\alpha}^\dagger a_{q\beta}^\dagger a_{r\alpha} a_{s\beta} + a_{p\beta}^\dagger a_{q\alpha}^\dagger a_{r\beta} a_{s\alpha} \\ =&\ \sum_{\sigma\sigma'} a_{p\sigma}^\dagger a_{q\sigma'}^\dagger a_{r\sigma} s_{s\sigma'}\end{split}$

Another case

$\begin{split}Z^{[0]} = &\ \Big[ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_q\big)^{[1/2]} \Big] \otimes_{[0]} \Big[ \big(a_r^\dagger \big)^{[1/2]} \otimes_{[1]} \big(a_s\big)^{[1/2]} \Big] \\ =&\ \begin{pmatrix} -a_{p\alpha}^\dagger a_{q\beta} \\ \frac{1}{\sqrt{2}} \Big( a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta} \Big) \\ a_{p\beta}^\dagger a_{q\alpha} \end{pmatrix}^{[1]} \otimes_{[0]} \begin{pmatrix} -a_{r\alpha}^\dagger a_{s\beta} \\ \frac{1}{\sqrt{2}} \Big( a_{r\alpha}^\dagger a_{s\alpha} - a_{r\beta}^\dagger a_{s\beta} \Big) \\ a_{r\beta}^\dagger a_{s\alpha} \end{pmatrix}^{[1]} \\ =&\ \frac{1}{\sqrt{3}} \begin{pmatrix} -a_{p\alpha}^\dagger a_{q\beta} a_{r\beta}^\dagger a_{s\alpha} -\frac{1}{2} \Big( a_{p\alpha}^\dagger a_{q\alpha} - a_{p\beta}^\dagger a_{q\beta} \Big) \Big( a_{r\alpha}^\dagger a_{s\alpha} - a_{r\beta}^\dagger a_{s\beta} \Big) - a_{p\beta}^\dagger a_{q\alpha} a_{r\alpha}^\dagger a_{s\beta} \end{pmatrix} \\ W^{[0]} =&\ \Big[ \big(a_p^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_q\big)^{[1/2]} \Big] \otimes_{[0]} \Big[ \big(a_r^\dagger \big)^{[1/2]} \otimes_{[0]} \big(a_s\big)^{[1/2]} \Big] \\ =&\ \frac{1}{\sqrt{2}} \begin{pmatrix} a_{p\alpha}^\dagger a_{q\alpha}+ a_{p\beta}^\dagger a_{q\beta} \end{pmatrix}^{[0]} \otimes_{[0]} \frac{1}{\sqrt{2}} \begin{pmatrix} a_{r\alpha}^\dagger a_{s\alpha}+ a_{r\beta}^\dagger a_{s\beta} \end{pmatrix}^{[0]} \\ =&\ \frac{1}{2} \Big( a_{p\alpha}^\dagger a_{q\alpha}+ a_{p\beta}^\dagger a_{q\beta}\Big) \Big( a_{r\alpha}^\dagger a_{s\alpha}+ a_{r\beta}^\dagger a_{s\beta} \Big)\end{split}$

Using

$(a-b)(c-d) + (a+b)(c+d) = (a+b)(2c) - (2b)(c-d) = 2(ac+bd)$

we have

$\begin{split}-\sqrt{3} Z^{[0]} + W^{[0]} =&\ a_{p\alpha}^\dagger a_{q\beta} a_{r\beta}^\dagger a_{s\alpha} + a_{p\beta}^\dagger a_{q\alpha} a_{r\alpha}^\dagger a_{s\beta} + a_{p\alpha}^\dagger a_{q\alpha} a_{r\alpha}^\dagger a_{s\alpha} + a_{p\beta}^\dagger a_{q\beta} a_{r\beta}^\dagger a_{s\beta} \\ =&\ \sum_{\sigma\sigma'} a_{p\sigma}^\dagger a_{q\sigma'} a_{r\sigma'}^\dagger a_{s\sigma}\end{split}$

#### S Term¶

From second singlet formula we have

$\sqrt{2} \sum_{i\in L} \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{S}_{i}^{R} \big)^{[\frac{1}{2}]} = \sum_{i\in L} \big( t_{ij} a_{i\alpha}^\dagger a_{j\alpha} + t_{ij} a_{i\beta}^\dagger a_{j\beta} \big)$

#### R Term¶

This is the same as the S term. Note that in the expression for $$\hat{R}$$, we have a $$\otimes_{[0]}$$, this is because in the original spatial expression there is a summation over $$\sigma$$. Then there is a $$[0] \otimes_{[1/2]} [1/2]$$, which will not produce any extra coefficients.

#### AP Term¶

Using definition

$\begin{split}\big( \hat{A}_{ik} \big)^{[0/1]} =&\ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \\ \big( \hat{P}_{ik}^{R} \big)^{[0/1]} =&\ -\sum_{jl\in R} v_{ijkl} \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{l} \big)^{[\frac{1}{2}]}\end{split}$

We have

$\begin{split}&\ \sum_{ik\in L} \left[ \sqrt{3} \big(\hat{A}_{ik} \big)^{[1]} \otimes_{[0]} \big(\hat{P}_{ik}^{R} \big)^{[1]} + \big(\hat{A}_{ik} \big)^{[0]} \otimes_{[0]} \big(\hat{P}_{ik}^{R} \big)^{[0]} \right] \\ =&\ \sum_{ik\in L,jl\in R} v_{ijkl} \left[ \sqrt{3} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]}\right] \otimes_{[0]} \left[ \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] + \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]}\right] \otimes_{[0]} \left[ \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \right] \\ =&\ \sum_{ik\in L,jl\in R} v_{ijkl} \left[ \sum_{\sigma\sigma'} a_{i\sigma}^\dagger a_{k\sigma'}^\dagger a_{j\sigma} a_{l\sigma'} \right] = -\sum_{ik\in L,jl\in R,\sigma\sigma'} v_{ijkl} a_{i\sigma}^\dagger a_{k\sigma'}^\dagger a_{l\sigma'} a_{j\sigma}\end{split}$

Note that in last step, we can anticommute $$a_{l\sigma'}, a_{j\sigma}$$ because it’s assumed that in the $$\sigma$$ summation, when $$j=l$$, $$\sigma \neq \sigma'$$. Otherwise there will be two $$a$$ operators acting on the same site and the contribution is zero.

#### BQ Term¶

In spatial expression, this term is $$BQ - B'Q'$$. Now $$-\sqrt{3} Z^{[0]} + W^{[0]}$$ gives $$B'Q'$$. And $$2 W^{[0]}$$ gives $$BQ$$. Therefore,

$2 W^{[0]} - \big(-\sqrt{3} Z^{[0]} + W^{[0]}\big) = \sqrt{3} Z^{[0]} + W^{[0]}$

This looks like $$\hat{A}\hat{P}$$ term, but without $$\frac{1}{2}$$ and $$h.c.$$. But this is not correct, because the definition of $$Q, Q'$$ is not equivalent due to the index order in $$v_{ijkl}$$. So they will give different $$W^{[0]}$$. Instead we have (note that $$\big( \hat{B}_{ij} \big)^{[0]} = \big( {\hat{B}'}_{ij} \big)^{[0]}$$)

$\begin{split}&\ \sum_{ij\in L} \left[ 2\Big( \hat{B}_{ij} \Big)^{[0]} \otimes_{[0]} \Big( \hat{Q}_{ij}^{R} \Big)^{[0]} - \Big( {\hat{B}'}_{ij} \Big)^{[0]} \otimes_{[0]} \Big( {\hat{Q}'}_{ij}^{R} \Big)^{[0]} + \sqrt{3} \Big( {\hat{B}'}_{ij} \Big)^{[1]} \otimes_{[0]} \Big( {\hat{Q}'}_{ij}^{R} \Big)^{[1]} \right] \\ =&\ \sum_{ij\in L} \left[ \Big( \hat{B}_{ij} \Big)^{[0]} \otimes_{[0]} \left( \Big( 2\hat{Q}_{ij}^{R} \Big)^{[0]} - \Big( {\hat{Q}'}_{ij}^{R} \Big)^{[0]} \right) + \sqrt{3} \Big( {\hat{B}'}_{ij} \Big)^{[1]} \otimes_{[0]} \Big( {\hat{Q}'}_{ij}^{R} \Big)^{[1]} \right]\end{split}$

Note that $$B, Q$$ do not have $$[1]$$ form.

### Normal/Complementary Partitioning¶

Note that

$\sqrt{2} \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{S}_{i}^{R} \big)^{[\frac{1}{2}]} + h.c. \right] = \sqrt{2} \sum_{i\in R} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{S}_{i}^{L} \big)^{[\frac{1}{2}]} + h.c. \right]$

Therefore,

$\begin{split}&\ \sqrt{2} \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{S}_{i}^{R} \big)^{[\frac{1}{2}]} + h.c. \right] + 2 \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{R} \big)^{[\frac{1}{2}]} + h.c. \right] + 2 \sum_{i\in R} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{L} \big)^{[\frac{1}{2}]} + h.c. \right] \\ =&\ \frac{\sqrt{2}}{2} \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{S}_{i}^{R} \big)^{[\frac{1}{2}]} + h.c. \right] + \frac{\sqrt{2}}{2} \sum_{i\in R} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{S}_{i}^{L} \big)^{[\frac{1}{2}]} + h.c. \right] \\ &\ + 2 \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{R} \big)^{[\frac{1}{2}]} + h.c. \right] + 2 \sum_{i\in R} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{L} \big)^{[\frac{1}{2}]} + h.c. \right] \\ =&\ 2 \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \Big[ \big( \hat{R}_{i}^{R} \big)^{[\frac{1}{2}]} + \frac{\sqrt{2}}{4} \big( \hat{S}_{i}^{R} \big)^{[\frac{1}{2}]} \Big] + h.c. \right] + 2 \sum_{i\in R} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \Big[ \big( \hat{R}_{i}^{L} \big)^{[\frac{1}{2}]} + \frac{\sqrt{2}}{4} \big( \hat{S}_{i}^{L} \big)^{[\frac{1}{2}]} \Big] + h.c. \right]\end{split}$

So define

$\big( \hat{R}_{i}^{\prime L/R} \big)^{[\frac{1}{2}]} := \frac{\sqrt{2}}{4} \big( \hat{S}_{i}^{L} \big)^{[\frac{1}{2}]} + \big( \hat{R}_{i}^{L} \big)^{[\frac{1}{2}]} = \frac{\sqrt{2}}{4} \sum_{j\in L/R} t_{ij} \big( a_{j} \big)^{[\frac{1}{2}]} + \sum_{jkl\in L/R} v_{ijkl} \left[ \Big( a_{k}^\dagger \Big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \otimes_{[\frac{1}{2}]} \big( a_{j} \big)^{[\frac{1}{2}]}$

Here $$\frac{\sqrt{2}}{4}$$ should be understood as $$\frac{1}{2} \cdot \frac{1}{\sqrt{2}}$$. The $$\frac{1}{2}$$ is the same as spatial case, and $$\frac{1}{\sqrt{2}}$$ is because the expected $$\sqrt{2}$$ factor is not added for the $$\hat{R}$$ term.

#### Operator Exchange factors¶

Here we consider fermion and SU(2) exchange factors together. From $$j_2 = 1/2$$ CG factors

$\begin{split}\bigg\langle j_1\ \left(M - \frac{1}{2} \right)\ \frac{1}{2}\ \frac{1}{2} \bigg| \left( j_1 \pm \frac{1}{2} \right)\ M \bigg\rangle =&\ \pm \sqrt{\frac{1}{2} \left( 1 \pm \frac{M}{j_1 + \frac{1}{2}} \right)} \\ \bigg\langle j_1\ \left(M + \frac{1}{2} \right)\ \frac{1}{2}\ \left( -\frac{1}{2}\right) \bigg| \left( j_1 \pm \frac{1}{2} \right)\ M \bigg\rangle =&\ \sqrt{\frac{1}{2} \left( 1 \mp \frac{M}{j_1 + \frac{1}{2}} \right)}\end{split}$

Let $$j_1 = \frac{1}{2}$$ we have

$\begin{split}\bigg\langle \frac{1}{2}\ \left( - \frac{1}{2} \right)\ \frac{1}{2}\ \frac{1}{2} \bigg| \left( \frac{1}{2} \pm \frac{1}{2} \right)\ 0 \bigg\rangle =&\ \pm \sqrt{\frac{1}{2} } \\ \bigg\langle \frac{1}{2} \ \frac{1}{2} \ \frac{1}{2}\ \left( -\frac{1}{2}\right) \bigg| \left( \frac{1}{2} \pm \frac{1}{2} \right)\ 0 \bigg\rangle =&\ \sqrt{\frac{1}{2} }\end{split}$

The exchange factor formula is

$\begin{split}\left( \hat{X}_1^{[S_1]} \otimes_{[S]} \hat{X}_2^{[S_2]} \right)^{[S_z]} =&\ \sum_{S_{1z},S_{2z}} \hat{X}_1^{[S_1][S_{1z}]} \hat{X}_2^{[S_2][S_{2z}]} \langle SS_z| S_1S_{1z},\ S_2 S_{2z} \rangle \\ =&\ \mathrm{P}_{\mathrm{fermi}}^{\mathrm{exchange}}(N_1,N_2) \sum_{S_{1z},S_{2z}} \hat{X}_2^{[S_2][S_{2z}]} \hat{X}_1^{[S_1][S_{1z}]} \langle SS_z| S_1S_{1z},\ S_2 S_{2z} \rangle \\ =&\ \mathrm{P}_{\mathrm{fermi}}^{\mathrm{exchange}}(N_1,N_2) \frac{\langle SS_z| S_1S_{1z},\ S_2 S_{2z} \rangle} {\langle SS_z| S_2S_{2z},\ S_1 S_{1z} \rangle} \left( \hat{X}_2^{[S_2]} \otimes_{[S]} \hat{X}_1^{[S_1]} \right)^{[S_z]} \\ \hat{X}_1^{[S_1]} \otimes_{[S]} \hat{X}_2^{[S_2]} =&\ \mathrm{P}_{\mathrm{fermi}}^{\mathrm{exchange}}(N_1,N_2) \mathrm{P}_{\mathrm{SU(2)}}^{\mathrm{exchange}}(S_1, S_2, S) \hat{X}_2^{[S_2]} \otimes_{[S]} \hat{X}_1^{[S_1]}\end{split}$

For $$[1/2] \otimes_{[0]} [1/2]$$, this is

$\mathrm{P}^{\mathrm{exchange}}(\tfrac{1}{2}, \tfrac{1}{2}, 0) = (-1) \frac{\big\langle \frac{1}{2} \ \frac{1}{2} \ \frac{1}{2}\ \left( -\frac{1}{2}\right) \big| 0\ 0 \big\rangle}{\big\langle \frac{1}{2} \ \left( -\frac{1}{2}\right) \ \frac{1}{2}\ \frac{1}{2} \big| 0\ 0 \big\rangle} = (-1) \frac{\sqrt{\frac{1}{2}}}{-\sqrt{\frac{1}{2}}} = 1$

For $$[1/2] \otimes_{[1]} [1/2]$$, this is

$\mathrm{P}^{\mathrm{exchange}}(\tfrac{1}{2}, \tfrac{1}{2}, 1) = (-1) \frac{\big\langle \frac{1}{2} \ \frac{1}{2} \ \frac{1}{2}\ \left( -\frac{1}{2}\right) \big| 1\ 0 \big\rangle}{\big\langle \frac{1}{2} \ \left( -\frac{1}{2}\right) \ \frac{1}{2}\ \frac{1}{2} \big| 1\ 0 \big\rangle} = (-1) \frac{\sqrt{\frac{1}{2}}}{\sqrt{\frac{1}{2}}} = -1$

From CG factors

$\langle 1\ m_1 \ 1 \ (-m_1) | 0 \ 0 \rangle = \frac{(-1)^{1-m_1}}{\sqrt{3}}$

we have

$\mathrm{P}^{\mathrm{exchange}}(1, 1, 0) = (+1) \frac{\big\langle 1 \ 1 \ \ 1\ -1 \big| 0\ 0 \big\rangle}{\big\langle 1 \ -1 \ 1\ 1 \big| 0\ 0 \big\rangle} = (+1) \frac{\frac{(-1)^{0}}{\sqrt{3}}}{\frac{(-1)^{2}}{\sqrt{3}}} = 1$

we have

$\begin{split}(\hat{H})^{[0], NC} =&\ \big( \hat{H}^{L} \big)^{[0]} \otimes_{[0]} \big( \hat{1}^{R} \big)^{[0]} + \big( \hat{1}^{L} \big)^{[0]} \otimes_{[0]} \big( \hat{H}^{R} \big)^{[0]} \\ &\ + 2 \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{\prime R} \big)^{[\frac{1}{2}]} + \big( a_{i}\big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{\prime R\dagger} \big)^{[\frac{1}{2}]} \right] + 2 \sum_{i\in R} \left[ \big( \hat{R}_{i}^{\prime L\dagger} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{i} \big)^{[\frac{1}{2}]} + \big( \hat{R}_{i}^{\prime L} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{i}^\dagger \big)^{[\frac{1}{2}]}\right] \\ &\ - \frac{1}{2} \sum_{ik\in L} \left[ \big(\hat{A}_{ik} \big)^{[0]} \otimes_{[0]} \big(\hat{P}_{ik}^{R} \big)^{[0]} + \sqrt{3} \big(\hat{A}_{ik} \big)^{[1]} \otimes_{[0]} \big(\hat{P}_{ik}^{R} \big)^{[1]} + \big(\hat{A}_{ik}^\dagger \big)^{[0]} \otimes_{[0]} \big(\hat{P}_{ik}^{R\dagger} \big)^{[0]} + \sqrt{3} \big(\hat{A}_{ik}^\dagger \big)^{[1]} \otimes_{[0]} \big(\hat{P}_{ik}^{R\dagger} \big)^{[1]} \right] \\ &\ +\sum_{ij\in L} \left[ \big( \hat{B}_{ij} \big)^{[0]} \otimes_{[0]} \big( {\hat{Q}}_{ij}^{\prime\prime R} \big)^{[0]} + \sqrt{3} \big( {\hat{B}'}_{ij} \big)^{[1]} \otimes_{[0]} \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[1]} \right]\end{split}$

With this normal/complementary partitioning, the operators required in left block are

$\big\{ \big( \hat{H}^L \big)^{[0]}, \big( \hat{1}^{L} \big)^{[0]}, \big( a_{i}^\dagger \big)^{[\frac{1}{2}]}, \big( a_{i} \big)^{[\frac{1}{2}]}, \big( \hat{R}_{k}^{\prime L\dagger} \big)^{[\frac{1}{2}]}, \big( \hat{R}_{k}^{\prime L} \big)^{[\frac{1}{2}]}, \big(\hat{A}_{ij} \big)^{[0]}, \big(\hat{A}_{ij} \big)^{[1]}, \big(\hat{A}_{ij}^\dagger \big)^{[0]}, \big(\hat{A}_{ij}^\dagger \big)^{[1]}, \big( \hat{B}_{ij} \big)^{[0]}, \big( {\hat{B}'}_{ij} \big)^{[1]} \big\}\quad (i,j\in L, k\in R)$

The operators required in right block are

$\big\{ \big( \hat{1}^{R} \big)^{[0]}, \big( \hat{H}^{R} \big)^{[0]}, \big( \hat{R}_{i}^{\prime R} \big)^{[\frac{1}{2}]}, \big( \hat{R}_{i}^{\prime R\dagger} \big)^{[\frac{1}{2}]}, \big( a_{k} \big)^{[\frac{1}{2}]}, \big( a_{k}^\dagger \big)^{[\frac{1}{2}]}, \big(\hat{P}_{ij}^{R} \big)^{[0]}, \big(\hat{P}_{ij}^{R} \big)^{[1]}, \big(\hat{P}_{ij}^{R\dagger} \big)^{[0]}, \big(\hat{P}_{ij}^{R\dagger} \big)^{[1]}, \big( {\hat{Q}}_{ij}^{\prime\prime R} \big)^{[0]}, \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[1]} \big\}\quad (i,j\in L, k\in R)$

Assuming that there are $$K$$ sites in total, and $$K_L/K_R$$ sites in left/right block (optimally, $$K_L \le K_R$$), the total number of operators (and also the number of terms in Hamiltonian with partition) in left or right block is

$N_{NC} = 1 + 1 + 2K_L + 2K_R + 4K_L^2 + 2K_L^2 = 6K_L^2 + 2K + 2$

### Complementary/Normal Partitioning¶

Note that due the CG factors, exchange any $$\otimes_{[0]}$$ product will not produce extra sign.

$\begin{split}(\hat{H})^{[0], CN} =&\ \big( \hat{H}^{L} \big)^{[0]} \otimes_{[0]} \big( \hat{1}^{R} \big)^{[0]} + \big( \hat{1}^{L} \big)^{[0]} \otimes_{[0]} \big( \hat{H}^{R} \big)^{[0]} \\ &\ + 2 \sum_{i\in L} \left[ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{\prime R} \big)^{[\frac{1}{2}]} + \big( a_{i}\big)^{[\frac{1}{2}]} \otimes_{[0]} \big( \hat{R}_{i}^{\prime R\dagger} \big)^{[\frac{1}{2}]} \right] + 2 \sum_{i\in R} \left[ \big( \hat{R}_{i}^{\prime L\dagger} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{i} \big)^{[\frac{1}{2}]} + \big( \hat{R}_{i}^{\prime L} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{i}^\dagger \big)^{[\frac{1}{2}]}\right] \\ &\ - \frac{1}{2} \sum_{jl\in R} \left[ \big(\hat{P}_{jl}^{L} \big)^{[0]} \otimes_{[0]} \big(\hat{A}_{jl} \big)^{[0]} + \sqrt{3} \big(\hat{P}_{jl}^{L} \big)^{[1]} \otimes_{[0]} \big(\hat{A}_{jl} \big)^{[1]} + \big(\hat{P}_{jl}^{L\dagger} \big)^{[0]} \otimes_{[0]} \big(\hat{A}_{jl}^\dagger \big)^{[0]} + \sqrt{3} \big(\hat{P}_{jl}^{L\dagger} \big)^{[1]} \otimes_{[0]} \big(\hat{A}_{jl}^\dagger \big)^{[1]} \right] \\ &\ +\sum_{kl\in R} \left[ \big( {\hat{Q}}_{kl}^{\prime\prime L} \big)^{[0]} \otimes_{[0]} \big( \hat{B}_{kl} \big)^{[0]} + \sqrt{3} \big( {\hat{Q}}_{kl}^{\prime L} \big)^{[1]} \otimes_{[0]} \big( {\hat{B}'}_{kl} \big)^{[1]} \right]\end{split}$

Now the operators required in left block are

$\big\{ \big( \hat{H}^L \big)^{[0]}, \big( \hat{1}^{L} \big)^{[0]}, \big( a_{i}^\dagger \big)^{[\frac{1}{2}]}, \big( a_{i} \big)^{[\frac{1}{2}]}, \big( \hat{R}_{k}^{\prime L\dagger} \big)^{[\frac{1}{2}]}, \big( \hat{R}_{k}^{\prime L} \big)^{[\frac{1}{2}]}, \big(\hat{P}_{kl}^{L} \big)^{[0]}, \big(\hat{P}_{kl}^{L} \big)^{[1]}, \big(\hat{P}_{kl}^{L\dagger} \big)^{[0]}, \big(\hat{P}_{kl}^{L\dagger} \big)^{[1]}, \big( {\hat{Q}}_{kl}^{\prime\prime L} \big)^{[0]}, \big( {\hat{Q}}_{kl}^{\prime L} \big)^{[1]} \big\}\quad (k,l\in R, i\in L)$

The operators required in right block are

$\big\{ \big( \hat{1}^{R} \big)^{[0]}, \big( \hat{H}^{R} \big)^{[0]}, \big( \hat{R}_{i}^{\prime R} \big)^{[\frac{1}{2}]}, \big( \hat{R}_{i}^{\prime R\dagger} \big)^{[\frac{1}{2}]}, \big( a_{k} \big)^{[\frac{1}{2}]}, \big( a_{k}^\dagger \big)^{[\frac{1}{2}]}, \big(\hat{A}_{kl} \big)^{[0]}, \big(\hat{A}_{kl} \big)^{[1]}, \big(\hat{A}_{kl}^\dagger \big)^{[0]}, \big(\hat{A}_{kl}^\dagger \big)^{[1]}, \big( \hat{B}_{kl} \big)^{[0]}, \big( {\hat{B}'}_{kl} \big)^{[1]} \big\}\quad (k,l\in R, i\in L)$

The total number of operators (and also the number of terms in Hamiltonian with partition) in left or right block is

$N_{CN} = 1 + 1 + 2K_L + 2K_R + 4K_R^2 + 2K_R^2 = 6K_R^2 + 2K + 2$

## Blocking¶

The enlarged left/right block is denoted as $$L*/R*$$. Make sure that all $$L$$ operators are to the left of $$*$$ operators. (The exchange factor for this is -1 for doublet $$\otimes$$ triplet and +1 doublet $$\otimes$$ singlet.)

First we have

$\begin{split}\big( \hat{R}_{i}^{L/R} \big)^{[1/2]} =&\ \sum_{jkl\in L/R} v_{ijkl} \left[ \big( a_{k}^\dagger \big)^{[1/2]} \otimes_{[0]} \big( a_{l} \big)^{[1/2]} \right] \otimes_{[1/2]} \big( a_{j} \big)^{[1/2]} \\ =&\ \frac{1}{\sqrt{2}} \sum_{jkl\in L/R} v_{ijkl} \begin{pmatrix} a_{k\alpha}^\dagger a_{l\alpha}+ a_{k\beta}^\dagger a_{l\beta} \end{pmatrix}^{[0]} \otimes_{[1/2]} \big( a_{j} \big)^{[1/2]} \\ =&\ \frac{1}{\sqrt{2}} \sum_{jkl\in L/R} v_{ijkl} \begin{pmatrix} -a_{k\alpha}^\dagger a_{l\alpha}a_{j\beta} - a_{k\beta}^\dagger a_{l\beta}a_{j\beta} \\ a_{k\alpha}^\dagger a_{l\alpha}a_{j\alpha}+ a_{k\beta}^\dagger a_{l\beta}a_{j\alpha} \end{pmatrix}^{[1/2]}\end{split}$

From the formula $$\sqrt{3} U^{[1/2]} - V^{[1/2]}$$ we have

$\big( \hat{R}_{i}^{L/R} \big)^{[1/2]} = \frac{\sqrt{3}}{2} \sum_{jkl\in L/R} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \Big] - \frac{1}{2} \sum_{jkl\in L/R} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \Big]$

From the formula $$\sqrt{3} S^{[1/2]} - T^{[1/2]}$$ we have (for $$k\neq l$$)

$\big( \hat{R}_{i}^{L/R} \big)^{[1/2]} = \frac{\sqrt{3}}{2} \sum_{jkl\in L/R} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \Big] - \frac{1}{2} \sum_{jkl\in L/R} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \Big[ \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \Big]$

We have

$\begin{split}\big( \hat{R}_{i}^{\prime L*} \big)^{[1/2]} =&\ \big( \hat{R}_{i}^{\prime L} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \\ &\ + \sum_{j \in L} \left[ \sum_{kl\in *} v_{ijkl} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \otimes_{[\frac{1}{2}]} \big( a_{j} \big)^{[\frac{1}{2}]} + \sum_{j \in *} \left[ \sum_{kl\in L} v_{ijkl} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \otimes_{[\frac{1}{2}]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ &\ - \frac{1}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \\ &\ - \frac{1}{2} \sum_{k \in *} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{k \in *} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right]\\ &\ - \frac{1}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right]\\ &\ - \frac{1}{2} \sum_{l\in *} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in L} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{l\in *} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in L} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \\ =&\ \big( \hat{R}_{i}^{\prime L} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \\ &\ + \sum_{j \in L} \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[\frac{1}{2}]} \left[ \sum_{kl\in *} v_{ijkl} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] + \sum_{j \in *} \left[ \sum_{kl\in L} v_{ijkl} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \otimes_{[\frac{1}{2}]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ &\ - \frac{1}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \\ &\ - \frac{1}{2} \sum_{k \in *} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{k \in *} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} \\ &\ - \frac{1}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right]\\ &\ - \frac{1}{2} \sum_{l\in *} \left[ \sum_{jk\in L} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_l\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{l\in *} \left[ \sum_{jk\in L} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_l\big)^{[1/2]}\end{split}$

After reordering of terms

$\begin{split}\big( \hat{R}_{i}^{\prime L*} \big)^{[1/2]} =&\ \big( \hat{R}_{i}^{\prime L} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \\ &\ - \frac{1}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \\ &\ + \sum_{j \in L} \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[\frac{1}{2}]} \left[ \sum_{kl\in *} v_{ijkl} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \\ &\ - \frac{1}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right]\\ &\ - \frac{1}{2} \sum_{k \in *} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{k \in *} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} \\ &\ + \sum_{j \in *} \left[ \sum_{kl\in L} v_{ijkl} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \right] \otimes_{[\frac{1}{2}]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ &\ - \frac{1}{2} \sum_{l\in *} \left[ \sum_{jk\in L} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_l\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{l\in *} \left[ \sum_{jk\in L} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_l\big)^{[1/2]} \\ =&\ \big( \hat{R}_{i}^{\prime L} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \\ &\ - \frac{1}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jl\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \\ &\ + \frac{1}{2} \sum_{j\in L} \big(a_j\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{kl\in *} (2 v_{ijkl} - v_{ilkj}) \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_l\big)^{[1/2]} \right] +\frac{\sqrt{3}}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \left[ \sum_{jk\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right]\\ &\ - \frac{1}{2} \sum_{k \in *} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[0]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{k \in *} \left[ \sum_{jl\in L} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} \\ &\ + \frac{1}{2} \sum_{j\in *} \left[ \sum_{kl\in L} (2v_{ijkl} - v_{ilkj}) \big(a_k^\dagger\big)^{[1/2]} \otimes_{[0]} \big(a_l\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_j\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{l\in *} \left[ \sum_{jk\in L} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1]} \big(a_j\big)^{[1/2]} \right] \otimes_{[1/2]} \big(a_l\big)^{[1/2]}\end{split}$

By definition (The overall exchange factor for $$[1/2] \otimes_{[0]} [1/2]$$ is 1, and for $$[1/2] \otimes_{[1]} [1/2]$$ is -1)

$\begin{split}\big( \hat{A}_{ik} \big)^{[0/1]} =&\ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \\ \big( \hat{A}_{ik}^\dagger \big)^{[0]} =&\ \big( a_{i} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{k} \big)^{[\frac{1}{2}]} = \big( a_{k} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{i} \big)^{[\frac{1}{2}]} \\ \big( \hat{A}_{ik}^\dagger \big)^{[1]} =&\ -\big( a_{i} \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{k} \big)^{[\frac{1}{2}]} = \big( a_{k} \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{i} \big)^{[\frac{1}{2}]} \\ \big( \hat{P}_{ik}^{R} \big)^{[0/1]} =&\ \sum_{jl\in R} v_{ijkl} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ \big( \hat{B}_{ij} \big)^{[0]} =&\ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ \big( {\hat{B}'}_{ij} \big)^{[1]} =&\ \big( a_{i}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{j} \big)^{[\frac{1}{2}]}\\ \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[1]} =&\ \sum_{kl\in R} v_{ilkj} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{l} \big)^{[\frac{1}{2}]} \\ \big( {\hat{Q}}_{ij}^{\prime \prime R} \big)^{[0]} =&\ \sum_{kl\in R} (2v_{ijkl} - v_{ilkj}) \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]}\end{split}$

we have

$\begin{split}\big( \hat{R}_{i}^{\prime L*,NC} \big)^{[1/2]} =&\ \big( \hat{R}_{i}^{\prime L} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \\ &\ - \frac{1}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{P}_{ik}^{*} \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{k \in L} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{P}_{ik}^{*} \big)^{[1]} \\ &\ + \frac{1}{2} \sum_{j\in L} \big(a_j\big)^{[1/2]} \otimes_{[1/2]} \big( {\hat{Q}}_{ij}^{\prime \prime *} \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{l\in L} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \big( {\hat{Q}}_{il}^{\prime *} \big)^{[1]}\\ &\ - \frac{1}{2} \sum_{k \in *,jl\in L} v_{ijkl} \big( \hat{A}_{jl}^\dagger \big)^{[0]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{k \in *,jl\in L} v_{ijkl} \big( \hat{A}_{jl}^\dagger \big)^{[1]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} \\ &\ + \frac{1}{2} \sum_{j\in *,kl\in L} (2v_{ijkl} - v_{ilkj}) \big( \hat{B}_{kl} \big)^{[0]} \otimes_{[1/2]} \big(a_j\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{l\in *,jk\in L} v_{ijkl} \big( {\hat{B}'}_{kj} \big)^{[1]} \otimes_{[1/2]} \big(a_l\big)^{[1/2]} \\ \big( \hat{R}_{i}^{\prime L*,CN} \big)^{[1/2]} =&\ \big( \hat{R}_{i}^{\prime L} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \\ &\ - \frac{1}{2} \sum_{k \in L,jl\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{A}_{jl}^\dagger \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{k \in L,jl\in *} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{A}_{jl}^\dagger \big)^{[1]} \\ &\ + \frac{1}{2} \sum_{j\in L,kl\in *} (2 v_{ijkl} - v_{ilkj}) \big(a_j\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{B}_{kl} \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{l\in L,jk\in *} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \big( {\hat{B}'}_{kj} \big)^{[1]} \\ &\ - \frac{1}{2} \sum_{k \in *} \big( \hat{P}_{ik}^{L} \big)^{[0]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{k \in *} \big( \hat{P}_{ik}^{L} \big)^{[1]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} \\ &\ + \frac{1}{2} \sum_{j\in *} \big( {\hat{Q}}_{ij}^{\prime \prime L} \big)^{[0]} \otimes_{[1/2]} \big(a_j\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{l\in *} \big( {\hat{Q}}_{il}^{ \prime L} \big)^{[1]} \otimes_{[1/2]} \big(a_l\big)^{[1/2]}\end{split}$

To generate symmetrized $$P$$, we need to change the $$A$$ line to the following

$- \frac{1}{4} \sum_{k \in *,jl\in L} (v_{ijkl} + v_{ilkj}) \big( \hat{A}_{jl}^\dagger \big)^{[0]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{4} \sum_{k \in *,jl\in L} (v_{ijkl} - v_{ilkj}) \big( \hat{A}_{jl}^\dagger \big)^{[1]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]}$

Similarly,

$\begin{split}\big( \hat{R}_{i}^{\prime R*,NC} \big)^{[1/2]} =&\ \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^R \big)^{[0]} + \big( \hat{1}^* \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime R} \big)^{[1/2]} \\ &\ - \frac{1}{2} \sum_{k \in *} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{P}_{ik}^{R} \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{k \in *} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{P}_{ik}^{R} \big)^{[1]} \\ &\ + \frac{1}{2} \sum_{j\in *} \big(a_j\big)^{[1/2]} \otimes_{[1/2]} \big( {\hat{Q}}_{ij}^{\prime \prime R} \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{l\in *} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \big( {\hat{Q}}_{il}^{\prime R} \big)^{[1]}\\ &\ - \frac{1}{2} \sum_{k \in R,jl\in *} v_{ijkl} \big( \hat{A}_{jl}^\dagger \big)^{[0]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{k \in R,jl\in *} v_{ijkl} \big( \hat{A}_{jl}^\dagger \big)^{[1]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} \\ &\ + \frac{1}{2} \sum_{j\in R,kl\in *} (2v_{ijkl} - v_{ilkj}) \big( \hat{B}_{kl} \big)^{[0]} \otimes_{[1/2]} \big(a_j\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{l\in R,jk\in *} v_{ijkl} \big( {\hat{B}'}_{kj} \big)^{[1]} \otimes_{[1/2]} \big(a_l\big)^{[1/2]} \\ \big( \hat{R}_{i}^{\prime R*,CN} \big)^{[1/2]} =&\ \big( \hat{R}_{i}^{\prime *} \big)^{[1/2]} \otimes_{[1/2]} \big( \hat{1}^R \big)^{[0]} + \big( \hat{1}^* \big)^{[0]} \otimes_{[1/2]} \big( \hat{R}_{i}^{\prime R} \big)^{[1/2]} \\ &\ - \frac{1}{2} \sum_{k \in *,jl\in R} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{A}_{jl}^\dagger \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{k \in *,jl\in R} v_{ijkl} \big(a_k^\dagger\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{A}_{jl}^\dagger \big)^{[1]} \\ &\ + \frac{1}{2} \sum_{j\in *,kl\in R} (2 v_{ijkl} - v_{ilkj}) \big(a_j\big)^{[1/2]} \otimes_{[1/2]} \big( \hat{B}_{kl} \big)^{[0]} +\frac{\sqrt{3}}{2} \sum_{l\in *,jk\in R} v_{ijkl} \big(a_l\big)^{[1/2]} \otimes_{[1/2]} \big( {\hat{B}'}_{kj} \big)^{[1]} \\ &\ - \frac{1}{2} \sum_{k \in R} \big( \hat{P}_{ik}^{*} \big)^{[0]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{k \in R} \big( \hat{P}_{ik}^{*} \big)^{[1]} \otimes_{[1/2]} \big(a_k^\dagger\big)^{[1/2]} \\ &\ + \frac{1}{2} \sum_{j\in R} \big( {\hat{Q}}_{ij}^{\prime \prime *} \big)^{[0]} \otimes_{[1/2]} \big(a_j\big)^{[1/2]} -\frac{\sqrt{3}}{2} \sum_{l\in R} \big( {\hat{Q}}_{il}^{ \prime *} \big)^{[1]} \otimes_{[1/2]} \big(a_l\big)^{[1/2]}\end{split}$

Number of terms

$\begin{split}N_{R',NC} =&\ (2 + 4K_L + 4K_L^2) K_R + (2 + 4 + 4K_R) K_L = 4K_L^2 K_R + 8K_L K_R + 2K + 4 K_L \\ N_{R',CN} =&\ (2 + 4K_L + 4) K_R + (2 + 4K_R^2 + 4K_R) K_L = 4K_R^2 K_L + 8K_R K_L + 2K + 4 K_R\end{split}$

Blocking of other complementary operators is straightforward

$\begin{split}\big( \hat{P}_{ik}^{L*,CN} \big)^{[0/1]} =&\ \big( \hat{P}_{ik}^{L} \big)^{[0/1]} \otimes_{[0/1]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[0/1]} \big( \hat{P}_{ik}^{*} \big)^{[0/1]} + \sum_{j \in L, l \in *} v_{ijkl} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{j} \big)^{[\frac{1}{2}]} + \sum_{j \in *, l \in L} v_{ijkl} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ =&\ \big( \hat{P}_{ik}^{L} \big)^{[0/1]} \otimes_{[0/1]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[0/1]} \big( \hat{P}_{ik}^{*} \big)^{[0/1]} \pm \sum_{j \in L, l \in *} v_{ijkl} \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{l} \big)^{[\frac{1}{2}]} + \sum_{j \in *, l \in L} v_{ijkl} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{j} \big)^{[\frac{1}{2}]} \\ \big( \hat{P}_{ik}^{R*,NC} \big)^{[0/1]} =&\ \big( \hat{P}_{ik}^{*} \big)^{[0/1]} \otimes_{[0/1]} \big( \hat{1}^R \big)^{[0]} + \big( \hat{1}^* \big)^{[0]} \otimes_{[0/1]} \big( \hat{P}_{ik}^{R} \big)^{[0/1]} \pm \sum_{j \in *, l \in R} v_{ijkl} \big( a_{j} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{l} \big)^{[\frac{1}{2}]} + \sum_{j \in R, l \in *} v_{ijkl} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{j} \big)^{[\frac{1}{2}]}\end{split}$

and

$\begin{split}\big( {\hat{Q}}_{ij}^{\prime \prime L*,CN} \big)^{[0]} =&\ \big( {\hat{Q}}_{ij}^{\prime \prime L} \big)^{[0]} \otimes_{[0]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[0]} \big( {\hat{Q}}_{ij}^{\prime \prime *} \big)^{[0]} + \sum_{k\in L, l \in *} (2v_{ijkl} - v_{ilkj}) \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} + \sum_{k\in *, l \in L} (2v_{ijkl} - v_{ilkj}) \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} \\ =&\ \big( {\hat{Q}}_{ij}^{\prime \prime L} \big)^{[0]} \otimes_{[0]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[0]} \big( {\hat{Q}}_{ij}^{\prime \prime *} \big)^{[0]} + \sum_{k\in L, l \in *} (2v_{ijkl} - v_{ilkj}) \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} + \sum_{k\in *, l \in L} (2v_{ijkl} - v_{ilkj}) \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \\ \big( {\hat{Q}}_{ij}^{\prime \prime R*,NC} \big)^{[0]} =&\ \big( {\hat{Q}}_{ij}^{\prime \prime *} \big)^{[0]} \otimes_{[0]} \big( \hat{1}^R \big)^{[0]} + \big( \hat{1}^* \big)^{[0]} \otimes_{[0]} \big( {\hat{Q}}_{ij}^{\prime \prime R} \big)^{[0]} + \sum_{k\in *, l \in R} (2v_{ijkl} - v_{ilkj}) \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} + \sum_{k\in R, l \in *} (2v_{ijkl} - v_{ilkj}) \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]}\end{split}$

and

$\begin{split}\big( {\hat{Q}}_{ij}^{\prime L*,CN} \big)^{[1]} =&\ \big( {\hat{Q}}_{ij}^{\prime L} \big)^{[1]} \otimes_{[1]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1]} \big( {\hat{Q}}_{ij}^{\prime *} \big)^{[1]} + \sum_{k\in L, l \in *} v_{ilkj} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{l} \big)^{[\frac{1}{2}]} + \sum_{k\in *, l \in L} v_{ilkj} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{l} \big)^{[\frac{1}{2}]} \\ =&\ \big( {\hat{Q}}_{ij}^{\prime L} \big)^{[1]} \otimes_{[1]} \big( \hat{1}^* \big)^{[0]} + \big( \hat{1}^L \big)^{[0]} \otimes_{[1]} \big( {\hat{Q}}_{ij}^{\prime *} \big)^{[1]} + \sum_{k\in L, l \in *} v_{ilkj} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{l} \big)^{[\frac{1}{2}]} - \sum_{k\in *, l \in L} v_{ilkj} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \\ \big( {\hat{Q}}_{ij}^{\prime R*,CN} \big)^{[1]} =&\ \big( {\hat{Q}}_{ij}^{\prime *} \big)^{[1]} \otimes_{[1]} \big( \hat{1}^R \big)^{[0]} + \big( \hat{1}^* \big)^{[0]} \otimes_{[1]} \big( {\hat{Q}}_{ij}^{\prime R} \big)^{[1]} + \sum_{k\in *, l \in R} v_{ilkj} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{l} \big)^{[\frac{1}{2}]} - \sum_{k\in R, l \in *} v_{ilkj} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]}\end{split}$

## Middle-Site Transformation¶

$\begin{split}\big( \hat{P}_{ik}^{L,NC\to CN} \big)^{[0/1]} =&\ \sum_{jl\in L} v_{ijkl} \big( a_{l} \big)^{[\frac{1}{2}]} \otimes_{[0/1]} \big( a_{j} \big)^{[\frac{1}{2}]} = \sum_{jl\in L} v_{ijkl} \big( \hat{A}_{jl}^\dagger \big)^{[0/1]} \\ \big( {\hat{Q}}_{ij}^{\prime \prime L,NC\to CN} \big)^{[0]} =&\ \sum_{kl\in R} (2v_{ijkl} - v_{ilkj}) \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[0]} \big( a_{l} \big)^{[\frac{1}{2}]} = \sum_{kl\in R} (2v_{ijkl} - v_{ilkj}) \big( \hat{B}_{kl} \big)^{[0]} \\ \big( {\hat{Q}}_{ij}^{\prime L,NC\to CN} \big)^{[1]} =&\ \sum_{kl\in R} v_{ilkj} \big( a_{k}^\dagger \big)^{[\frac{1}{2}]} \otimes_{[1]} \big( a_{l} \big)^{[\frac{1}{2}]} = \sum_{kl\in R} v_{ilkj} \big( {\hat{B}'}_{kl} \big)^{[1]}\end{split}$