Physics:Schrieffer–Wolff transformation

Short description: In physics, an operator version of second-order perturbation theory

In quantum mechanics, the Schrieffer–Wolff transformation is a unitary transformation used to determine an effective (often low-energy) Hamiltonian by decoupling weakly interacting subspaces. Using a pertubative approach, the transformation can be constructed such that the interaction between the two subspaces vanishes up to the desired order in the perturbation. The transformation also perturbatively diagonalizes the system Hamiltonian to first order in the interaction. In this, the Schrieffer–Wolff transformation is an operator version of second-order perturbation theory. The Schrieffer–Wolff transformation is often used to project out the high energy excitations of a given quantum many-body Hamiltonian in order to obtain an effective low energy model. The Schrieffer–Wolff transformation thus provides a controlled perturbative way to study the strong coupling regime of quantum-many body Hamiltonians.

Although commonly attributed to the paper in which the Kondo model was obtained from the Anderson impurity model by J.R. Schrieffer and P.A. Wolff., Joaquin Mazdak Luttinger and Walter Kohn used this method in an earlier work about non-periodic k·p perturbation theory. Using the Schrieffer–Wolff transformation, the high energy charge excitations present in Anderson impurity model are projected out and a low energy effective Hamiltonian is obtained which has only virtual charge fluctuations. For the Anderson impurity model case, the Schrieffer–Wolff transformation showed that the Kondo model lies in the strong coupling regime of the Anderson impurity model.

Derivation

Consider a quantum system evolving under the time-independent Hamiltonian operator [math]\displaystyle{ H }[/math] of the form:[math]\displaystyle{ H = H_0 + V }[/math]where [math]\displaystyle{ H_0 }[/math] is a Hamiltonian with known eigenstates [math]\displaystyle{ |m\rangle }[/math] and corresponding eigenvalues [math]\displaystyle{ E_m }[/math], and where [math]\displaystyle{ V }[/math] is a small perturbation. Moreover, it is assumed without loss of generality that [math]\displaystyle{ V }[/math] is purely off-diagonal in the eigenbasis of [math]\displaystyle{ H_0 }[/math], i.e., [math]\displaystyle{ \langle m | V |m\rangle = 0 }[/math] for all [math]\displaystyle{ m }[/math]. Indeed, this situation can always be arranged by absorbing the diagonal elements of [math]\displaystyle{ V }[/math] into [math]\displaystyle{ H_0 }[/math], thus modifying its eigenvalues to [math]\displaystyle{ E'_m = E_m + \langle m|V| m\rangle }[/math].

The Schrieffer–Wolff transformation is a unitary transformation which expresses the Hamiltonian in a basis (the "dressed" basis) where it is diagonal to first order in the perturbation [math]\displaystyle{ V }[/math]. This unitary transformation is conventionally written as:[math]\displaystyle{ H' = e^{S} H e^{-S} }[/math]When [math]\displaystyle{ V }[/math] is small, the generator [math]\displaystyle{ S }[/math] of the transformation will likewise be small. The transformation can then be expanded in [math]\displaystyle{ S }[/math] using the Baker-Campbell-Haussdorf formula[math]\displaystyle{ H' = H + [S,H] + \frac{1}{2}[S,[S,H]]+ \dots }[/math]Here, [math]\displaystyle{ [A,B] }[/math] is the commutator between operators [math]\displaystyle{ A }[/math] and [math]\displaystyle{ B }[/math]. In terms of [math]\displaystyle{ H_0 }[/math] and [math]\displaystyle{ V }[/math], the transformation becomes[math]\displaystyle{ H' = H_0 + V + [S,H_0] + [S,V] + \frac{1}{2}[S,[S,H_0]] + \frac{1}{2}[S,[S,V]] + \dots }[/math]The Hamiltonian can be made diagonal to first order in [math]\displaystyle{ V }[/math] by choosing the generator [math]\displaystyle{ S }[/math] such that[math]\displaystyle{ V + [S, H_0] = 0 }[/math]This equation always has a definite solution under the assumption that [math]\displaystyle{ V }[/math] is off-diagonal in the eigenbasis of [math]\displaystyle{ H_0 }[/math]. Substituting this choice in the previous transformation yields:[math]\displaystyle{ H' = H_0 + \frac{1}{2}[S,V] + O(V^3) }[/math]This expression is the standard form of the Schrieffer–Wolff transformation. Note that all the operators on the right-hand side are now expressed in a new basis "dressed" by the interaction [math]\displaystyle{ V }[/math] to first order.

In the general case, the difficult step of the transformation is to find an explicit expression for the generator [math]\displaystyle{ S }[/math]. Once this is done, it is straightforward to compute the Schrieffer-Wolff Hamiltonian by computing the commutator [math]\displaystyle{ [S,V] }[/math]. The Hamiltonian can then be projected on any subspace of interest to obtain an effective projected Hamiltonian for that subspace. In order for the transformation to be accurate, the eliminated subspaces must be energetically well separated from the subspace of interest, meaning that the strength of the interaction [math]\displaystyle{ V }[/math] must be much smaller than the energy difference between the subspaces. This is the same regime of validity as in standard second-order perturbation theory.

Particular case

This section will illustrate how to practically compute the Schrieffer-Wolff (SW) transformation in the particular case of an unperturbed Hamiltonian that is block-diagonal.

But first, to properly compute anything, it is important to understand what is actually happening during the whole procedure. The SW transformation [math]\displaystyle{ W=e^S }[/math] being unitary, it does not change the amount of information or the complexity of the Hamiltonian. The resulting shuffle of the matrix elements creates, however, a hierarchy in the information (e.g. eigenvalues), that can be used afterward for a projection in the relevant sector. In addition, when the off-diagonal elements coupling the blocks are much smaller than the typical unperturbed energy scales, a perturbative expansion is allowed to simplify the problem.

Consider now, for concreteness, the full Hamiltonian [math]\displaystyle{ H=H_0+V }[/math] with an unperturbed part [math]\displaystyle{ H_0 }[/math] made of independent blocks [math]\displaystyle{ H_0^i }[/math]. In physics, and in the original motivation for the SW transformation, it is desired that each block corresponds to a distinct energy scale. In particular, all degenerate energy levels should belong to the same block. This well-split Hamiltonian is our starting point [math]\displaystyle{ H_0 }[/math]. A perturbative coupling [math]\displaystyle{ V }[/math] takes now on a specific meaning: the typical matrix element coupling different sectors must be much smaller than the eigenvalue differences between those sectors. The SW transformation will modify each block [math]\displaystyle{ H_0^i }[/math] into an effective Hamiltonian [math]\displaystyle{ H_{eff}^i }[/math] incorporating ("integrating out") the effects of the other blocks via the perturbation [math]\displaystyle{ V }[/math]. In the end, it is sufficient to look at the sector of interest (called a projection) and to work with the chosen effective Hamiltonian to compute, for instance, eigenvalues and eigenvectors. In physics, this would generate effective low- (or high-)energy Hamiltonians.

As mentioned in the previous section, the difficult step is the computation of the generator [math]\displaystyle{ S }[/math] of the SW transformation. To obtain results comparable to second-order perturbation theory, it is enough to solve the equation [math]\displaystyle{ [H_0,S]=V }[/math] (see Derivation). A simple trick in two steps is available when [math]\displaystyle{ H_0 }[/math] is block-diagonal.

The first step consists of finding the unitary transformation [math]\displaystyle{ U }[/math] diagonalizing [math]\displaystyle{ H_0 }[/math]. Since each block [math]\displaystyle{ H_0^i }[/math] can be diagonalized with a unitary transformation [math]\displaystyle{ U^i }[/math] (this is the matrix of right-eigenvectors of [math]\displaystyle{ H_0^i }[/math]), it is enough to build [math]\displaystyle{ U=\operatorname{diag}(U^i) }[/math], composed of the smaller rotations [math]\displaystyle{ U^i }[/math] on its diagonal, to transform [math]\displaystyle{ H_0 }[/math] into a purely diagonal matrix [math]\displaystyle{ D_0=\operatorname{diag}(d_i) }[/math].

The application of [math]\displaystyle{ U }[/math] to the whole matrix [math]\displaystyle{ H }[/math] yields then [math]\displaystyle{ H_D=U^{-1}HU= D_0 + V' }[/math] with a transformed perturbation [math]\displaystyle{ V'=U^{-1}VU }[/math], which remains off-diagonal. In this new form, the second step to compute [math]\displaystyle{ S }[/math] becomes very simple, since we obtain an explicit expression, in components: [math]\displaystyle{ S_{ij}=\frac{V'_{ij}}{d_i-d_j} }[/math] where [math]\displaystyle{ d_i }[/math] denotes the [math]\displaystyle{ i }[/math]th element on the diagonal of [math]\displaystyle{ D_0 }[/math]. The reason for this comes from the observation that, for any matrix [math]\displaystyle{ A=(A_{ij}) }[/math], and diagonal matrix [math]\displaystyle{ D=\operatorname{diag}(d_i) }[/math], we have the relation [math]\displaystyle{ [D,A]_{ij}=(d_i-d_j)A_{ij} }[/math]. Since the generator for [math]\displaystyle{ H_D }[/math] is defined by [math]\displaystyle{ [D_0,S]=V' }[/math], the above formula follows immediately. As expected, the associated operator [math]\displaystyle{ W=e^{S} }[/math] is unitary (it satisfies [math]\displaystyle{ W^\dagger=W^{-1} }[/math]) because the denominator of [math]\displaystyle{ S }[/math] changes sign when transposed, and [math]\displaystyle{ V }[/math] is Hermitian.

Using the last formula in the derivation, the second-order Schrieffer-Wolff-transformed Hamiltonian [math]\displaystyle{ H'=e^{S}H_D e^{-S} }[/math] has now an explicit form as a function of its elementary terms [math]\displaystyle{ D_0 }[/math] and [math]\displaystyle{ V' }[/math]: [math]\displaystyle{ H'_{ij} = d_i \delta_{ij} + \frac{1}{2} \sum_k V'_{ik} V'_{kj} \left( \frac{1}{d_i-d_k} + \frac{1}{d_j-d_k} \right) + O(V'^{3}) }[/math]

The "dressed" states have an energy [math]\displaystyle{ E'_n = d_n + \sum_k \frac{V'_{nk} V'_{kn}}{d_n-d_k} }[/math] following the recipe for first-order (non-degenerate) perturbation theory. This is applicable since the SW transformation is based on the approximation [math]\displaystyle{ |V_{ij}|\ll |d_i-d_j| }[/math] Note that the unitary rotation [math]\displaystyle{ U }[/math] does not affect the eigenvalues, meaning that [math]\displaystyle{ E'_n }[/math] is also a meaningful approximation for the original Hamiltonian [math]\displaystyle{ H }[/math].

The "dressed" states themselves can be derived, in first-order perturbation theory too, as [math]\displaystyle{ \psi'_{n} = \psi_{R,n} + \frac{1}{2}\sum_{k\neq n}\sum_j \psi_{R,k} V'_{kj} V'_{jn} \frac{2 d_j - d_k - d_n}{(d_j-d_n)(d_j-d_k)(d_n-d_k)} }[/math] Notice the index [math]\displaystyle{ R }[/math] to the unperturbed eigenstate [math]\displaystyle{ \psi_R }[/math] of [math]\displaystyle{ D_0 }[/math] to recall the current rotated basis of [math]\displaystyle{ H_D }[/math]. To express the eigenstates in the natural basis [math]\displaystyle{ \psi }[/math] of [math]\displaystyle{ H }[/math] itself, it is necessary to perform the unitary transformation [math]\displaystyle{ \psi_R \to U^{-1}\psi }[/math].

Further reading

• Phillips, Philip (2012). Advanced Solid State Physics (Second ed.). New York: Cambridge University Press. pp. 109–114. ISBN 978-1-107-49346-9. od