Transformation of the many-body Hamiltonian

We discuss the basic concept of the Schrieffer-Wolff transformation in the section Perturbative similarity transformation. We detail our proposed extended Schrieffer-Wolff transformation method in sections Symmetry specification of block-offdiagonal operators and Schrieffer-Wolff transformation as a system of linear equations for unique block-offdiagonal operators.

Perturbative similarity transformation

In principle, any Hamiltonian can be completely diagonalized by means of a particular unitary transformation $U$ . In practice, finding the particular unitary transformation $U$ often requires a complete diagonalization of the Hamiltonian to begin with. Here we recap how a perturbative similarity transformation, namely the Schrieffer-Wolff transformation, can be used to determine an approximate transformation operator $U$ , or the generator $S$ thereof, which does not diagonalize the Hamiltonian fully, but yields a block-diagonal Hamiltonian instead. These blocks consist of the orthonormal states in the Hilbert space $H$ that share the a given choice of characteristics, e.g. the local particle quantum number $n_{i}$ . We will denote the set of terms in the Hamiltonian that are already block-diagonal in the initial basis as $H_{0}$ . The remaining terms connect different blocks, i.e are block-offdiagonal, and are denoted $V$ . The complete Hamiltonian thus reads

$H = H_{0} + V .$

A unitary similarity transformation of the Hamiltonian is given by

$\tilde{H} = U^{†} H U = e^{- S} (H_{0} + V) e^{S} = e^{- S} H_{0} e^{S} + e^{- S} V e^{S},$

where $S$ is an anti-hermitian operator. One refers to it as the generator of the Schrieffer-Wolff transformation. The key problem of the Schrieffer-Wolff transformation is finding the generator $S$ such that the transformed Hamiltonian $\tilde{H}$ becomes entirely block-diagonal. In order to arrive at an equation for $S$ one makes use of the Campbell-Baker-Hausdorff formula to expand

$\tilde{H} = e^{S} H e^{- S} = m = 0 \sum \infty \frac{1}{m !} [S, H]_{m} = H + [S, H] + \frac{1}{2} [S, [S, H]] + \dots,$

where for generators $S$ satisfying $∥ S ∥ ≪ ∥ H_{0} ∥$ , with $∥ \cdot ∥$ a suitable norm, one can approximate the expression as

$\tilde{H} = H_{0} + V + [S, H_{0}] + [S, V] + O (S^{2}) .$

Considering that commutators of pairs of block-diagonal operators or pairs of block-offdiagonal operators respectively generally become block-diagonal, while the commutators of block-diagonal operators with block-offdiagonal operators become block-offdiagonal, one chooses the equation

$[S, H_{0}] = - V \Leftrightarrow [H_{0}, S] = V,$

by which one can determine the generator $S$ which removes the block-offdiagonal terms $V$ of the Hamiltonian. If a solution $S$ exists, one can use

$[S, [S, H_{0}]] = [S, - V] = - [S, V] \neq \in O (S^{2}),$

to simplify the expression for the transformed Hamiltonian

$\tilde{H} = e^{- S} (H_{0} + V) e^{S} ≃ H_{0} + \frac{1}{2} [S, V] + O (S^{2}),$

where the terms originating from $[S, V] /2$ contain the perturbative corrections arising from the consecutive application of two block-offdiagonal operators. A subsequent projection to the subspaces

$P_{n} = p \in P_{n} \sum ∣ p ⟩ ⟨ p ∣,$

yields the block-diagonal Hamiltonian

$H_{block - diagonal} = n \sum P_{n} H P_{n},$

where $n$ denotes the distinct blocks of the Hilbert space.

Symmetry specification of block-offdiagonal operators

The block-offdiagonal part of the Hamiltonian $V$ consists of a sum of block-offdiagonal terms. These in turn comprise products of individually block-offdiagonal operators $x$ . In the following we detail a method to decompose generic block-offdiagonal operators $x$ into distinct components. Each of the components exclusively connects two distinct blocks of the Hilbert space $H$ , often associated with distinct quantum numbers of a symmetry of the system. A given block-offdiagonal operator $x : H \to H$ satisfies

$[H_{0}, x] = ε z \neq = 0,$

where $ε$ denotes an arbitrary scalar and $z : H \to H$ an arbitrary operator. Let $A$ be a diagonal operator in the initial basis. It can be identical to the symmetry operator differentiating the blocks of the Hilbert space, but it is not required to be. One can use the spectrum of $A$ to expand the operator $x$ as

$x = q \sum x_{q} = q \sum β_{q} i \neq = q \prod (A - a_{i}) x,$

where the different $x_{q}$ couple the target subspace associated with the eigenvalue $a_{q}$ of $A$ to other subspaces of the Hilbert space. If the operator $A$ satisfies $[x, A] = 0,$ then both subspaces, initial and final, coupled via $x_{q}$ are specified by the eigenvalue $a_{q}$ . This is possible for the fermionic creation and annihilation operators $c_{iσ}^{†}$ and $c_{iσ}$ . The coefficients $β_{q}$ are solutions to the equation $β_{q} i \neq = q \prod (a_{q} - a_{i}) = 1 .$ The symmetry-specified block-offdiagonal operators $x_{q}$ satisfy $[x_{q}, x_{q^{'}}^{†}] \propto δ_{q q^{'}},$ where the operator $x_{q}^{†}$ denotes the hermitian conjugate of the operator $x_{q}$ .

Schrieffer-Wolff transformation as a system of linear equations for unique block-offdiagonal operators

Following the procedure outlined in the section "[Perturbative similarity transformation](#base_sw)" we separate the Hamiltonian of the system into its block-diagonal and block-offdiagonal contributions as

H = H_{0} + V .

The block-offdiagonal contribution

V

comprises each block-offdiagonal term

v

V = {v} \sum α_{v} (j \prod o^{j} i \prod x^{i}) + (j \prod o^{j} i \prod x^{i})^{†} = {v} \sum α_{v} i, j \prod o^{j} q (i) \sum x_{q (i)}^{i} + i, j \prod o^{j} q (i) \sum x_{q (i)}^{i}^{†},

where

\prod_{i} x^{i}

denotes sequences of individually block-offdiagonal operators,

\prod_{j} o^{j}

denotes sequences of individually block-diagonal operators, and

x_{q (i)}^{i}

denotes the symmetry-specified components of the operator

x^{i}

. We introduce a vector spaces

V_{0}^{h}

and

V_{0}^{a}

, for which each unique pair of hermitian, or anti-hermitian respectively, symmetry-specified operator sequences in

V

corresponds to a unique orthonormal basis vector

\overset{e}{^}_{v} = ⎩ ⎨ ⎧ (\prod_{j} o^{j} \prod_{i} x_{q (i)}^{i}) + (\prod_{j} o^{j} \prod_{i} x_{q (i)}^{i})^{†} (\prod_{j} o^{j} \prod_{i} x_{q (i)}^{i}) - (\prod_{j} o^{j} \prod_{i} x_{q (i)}^{i})^{†} \in V_{0}^{h} \in V_{0}^{a},

where

V^{h}

denotes the hermitian vector space and

V^{a}

the anti-hermitian vector space. We define the linear map

L : V_{0}^{h} \to V_{1}^{a} V_{0}^{a} \to V_{1}^{h},

where the action of the linear map on a vector

\overset{e}{^}_{v}

is given by

L \overset{e}{^}_{v} = [H_{0}, \overset{e}{^}_{v}],

and where

(L \overset{e}{^}_{v}) \in V_{1}^{a} (L \overset{e}{^}_{v}) \in V_{1}^{h} if if \overset{e}{^}_{v} \in V_{0}^{h} \overset{e}{^}_{v} \in V_{0}^{a} .

We have that

dim (V_{1}) \geq dim (V_{0}),

since

V_{1}

includes additional unique operator sequences generated by

[H_{0}, \overset{e}{^}_{v}]

. In the vector spaces

V

, we can translate the equation for the generator

S

of the Schrieffer-Wolff transformation to

L S = V \Leftrightarrow [H_{0}, S] = V,

with

S \in V_{0}^{a}

and

V \in V_{1}^{h}

. Determining the generator

S

of the Schrieffer-Wolff transformation becomes equivalent to solving the set of linear equations. In general one finds

rank (L) \leq dim (V_{0}) \leq dim (V_{1}),

and there is consequently no unique solution

S

to the set of equations. It is intuitive that the terms of

V

approximately bring about transitions between distinct eigenstates of

H_{0}

belonging to different blocks of the Hilbert space. This is reflected by

[H_{0}, \overset{e}{^}_{v}^{h}] = (Δ E_{0, v} l \prod o^{l}) \overset{e}{^}_{v}^{a},

where

Δ E_{0, v}

denotes the eigenvalue difference between the two eigenstates of

H_{0}

and

\prod_{l} o^{l}

denotes an arbitrary sequence of individually block-diagonal operators. From the previous equations we can approximate

∥ S ∥ \approx {v} \sum (\frac{α _{v}}{Δ E _{0, v}})^{2},

which highlights the necessity for a significant energy gap

Δ E_{0, v}

between the separate subspaces of the Hilbert space coupled by

V

in order for the series expansion of

\tilde{H}

O (S^{2})

to be considered a good approximation. We consider terms

v

for which

\frac{α _{v}^{2}}{∣Δ E _{0, v} ∣} \geq 1,

to be the resonant terms of

V

which should be retained in the transformed Hamiltonian

\tilde{H}

. To identify these resonant terms of

V

, we employ a singular value decomposition (SVD) of the linear map

L

and arrive at

L = U Σ W^{†} = U (Σ_{>} + Σ_{<}) W^{†} = L_{gapped} + L_{resonant},

where

Σ_{>}

comprises the singular values

σ_{i} > \frac{1}{N _{v}} {v} \sum α_{v}^{2},

where

N_{v}

denotes the number of terms

v

. We define the Moore-Penrose pseudoinverse of

L

that acts exclusively on the gapped, i.e. non-resonant, terms in

V

L_{gapped}^{+} = (W Σ_{>}^{+} U^{†}),

where

Σ_{>}^{+}

denotes the pseudoinverse of

Σ_{>}

. In the absence of a unique solution to

[H_{0}, S] = V

, the closest approximate solution is given by

S = (L_{gapped}^{+} V) \in V_{0}^{a},

and the resulting transformed Hamiltonian reads

\tilde{H} = H_{0} + V_{resonant} + [S, V_{resonant}] + \frac{1}{2} [S, V_{gapped}],

where

V_{gapped} = [H_{0}, S]

and

V_{resonant} = (V - V_{gapped})

HQS Spin Mapper Documentation

Transformation of the many-body Hamiltonian

Perturbative similarity transformation

Symmetry specification of block-offdiagonal operators

Schrieffer-Wolff transformation as a system of linear equations for unique block-offdiagonal operators