Papers
Topics
Authors
Recent
Search
2000 character limit reached

Three-Impurity Anderson Model

Updated 5 July 2026
  • The Three-Impurity Anderson Model comprises two variants: a three-orbital impurity with multipolar Kondo exchange and a trimer cluster with effective inter-impurity hopping and limited screening channels.
  • It employs a Schrieffer–Wolff transformation to derive a generalized Kondo Hamiltonian that integrates spin, orbital, and quadrupolar interactions, resulting in a two-stage screening process.
  • The trimer variant exhibits frustration-induced non-Fermi-liquid behavior and Kosterlitz–Thouless criticality, emphasizing the impact of effective hybridization rank on impurity screening.

Searching arXiv for the cited papers to ground the article in the current record. The Three-Impurity Anderson Model denotes, in the literature considered here, two closely related but distinct Anderson-impurity settings. In one sense, it refers to a three-orbital Anderson impurity model with a single correlated impurity carrying three orbital degrees of freedom and coupled to a conduction bath; in another, it denotes a three-site multi-impurity cluster or trimer, in which three correlated orbitals or moments are embedded in a host. Both realizations are governed by the same general ingredients—local Coulomb interactions, hybridization to itinerant electrons, and low-energy screening—but they organize the impurity sector differently and therefore generate different infrared phenomena. Representative results include multipolar Kondo exchange and two-stage screening in the three-orbital Kanamori problem, and channel-limited screening, frustration-induced non-Fermi-liquid behavior, and Kosterlitz–Thouless criticality in trimer geometries (Horvat et al., 2016, Eickhoff et al., 2020).

1. Terminological scope and model variants

A persistent ambiguity in usage arises because “three-impurity” may refer either to three correlated sites or to one impurity with three orbitals. The three-orbital formulation studied in "Low-energy physics of three-orbital impurity model with Kanamori interaction" is a single-impurity Anderson model whose impurity interaction has orbital SO(3)SO(3) and spin SU(2)SU(2) symmetry and is intended to describe systems with partially occupied t2gt_{2g} shells (Horvat et al., 2016). By contrast, "Strongly correlated multi-impurity models: The crossover from a single-impurity problem to lattice models" treats a genuine multi-impurity setting in which several correlated orbitals form a cluster coupled to a host, with the trimer as the minimal nontrivial case (Eickhoff et al., 2020).

The distinction matters because the low-energy degrees of freedom are not the same. In the three-orbital single-impurity problem, the impurity carries spin, orbital dipole, and orbital quadrupole moments, and the effective Kondo problem is correspondingly multipolar. In the trimer, the decisive objects are the effective inter-impurity hoppings generated by the real part of the hybridization self-energy and the number of independent screening channels determined by the rank of its imaginary part. A common misconception is therefore to treat the two problems as interchangeable. The literature summarized here instead supports a narrower conclusion: they belong to the same broad Anderson-impurity family, but their microscopic content and infrared mechanisms are different.

2. Three-orbital impurity with Kanamori interaction

The three-orbital Anderson impurity model is defined by

H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},

with a flat conduction bath, impurity-bath hybridization, and a local interaction of Kanamori form. In the notation of the cited work,

Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},

Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},

and

Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.

Here NdN_d is the impurity occupancy, S\mathbf{S} the total impurity spin, L\mathbf{L} the total impurity orbital angular momentum, SU(2)SU(2)0 the Hubbard repulsion, and SU(2)SU(2)1 Hund’s coupling. The conduction band is taken as flat with density of states SU(2)SU(2)2, and the hybridization width is SU(2)SU(2)3 (Horvat et al., 2016).

The orbital sector has SU(2)SU(2)4 symmetry rather than bare SU(2)SU(2)5 symmetry. For SU(2)SU(2)6 orbitals, the natural impurity multipoles are the spin SU(2)SU(2)7, the orbital dipole SU(2)SU(2)8, and the orbital quadrupoles SU(2)SU(2)9. The orbital generators t2gt_{2g}0 span the eight generators of t2gt_{2g}1, even though the bare Hamiltonian has only orbital t2gt_{2g}2 symmetry. The same work introduces an interpolation parameter t2gt_{2g}3 connecting the Dworin–Narath model t2gt_{2g}4 and the Kanamori model t2gt_{2g}5,

t2gt_{2g}6

which is useful for isolating the effect of reduced orbital symmetry.

The physically emphasized regime is the two-electron sector t2gt_{2g}7, relevant to Hund’s metals. In that sector the atomic ground state has

t2gt_{2g}8

i.e. the Hund’s-rule high-spin state. This regime is singled out because the local moment is relatively large and its screening is strongly suppressed, producing low coherence scales and strong-correlation effects.

3. Schrieffer–Wolff reduction and multipolar Kondo structure

The low-energy physics of the three-orbital Kanamori problem is derived by a Schrieffer–Wolff transformation that integrates out charge fluctuations and projects onto the impurity ground-state multiplet. The general form given in the paper is

t2gt_{2g}9

where H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},0 projects onto the impurity ground multiplet with H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},1 electrons and H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},2 onto virtual charge-excited multiplets (Horvat et al., 2016).

For the Kanamori case this yields a generalized Kondo Hamiltonian

H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},3

The bath operators entering this expression are the local bath charge density H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},4, the bath spin H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},5, the bath orbital dipole H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},6, and the bath orbital quadrupole H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},7. The low-energy problem is therefore not a purely spin-exchange Kondo model; it contains spin, orbital, quadrupolar, and mixed spin-orbital exchange.

A central quantitative result is the hierarchy of bare couplings in the H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},8, H=Hbath+Hhyb+Himp,H = H_{\mathrm{bath}} + H_{\mathrm{hyb}} + H_{\mathrm{imp}},9, Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},0 ground state. The spin-spin coupling Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},1 is much smaller than the orbital couplings for most parameters and can even become ferromagnetic near the Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},2 valence-fluctuation regime, whereas Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},3, Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},4, Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},5, and Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},6 are comparatively large. The hierarchy is summarized as

Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},7

Because the Kanamori interaction has only orbital Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},8 symmetry, the orbital-dipole and orbital-quadrupole couplings are split: Hbath=k,m,σϵkckmσckmσ,H_{\mathrm{bath}} = \sum_{k,m,\sigma} \epsilon_k\, c^\dagger_{km\sigma} c_{km\sigma},9 with Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},0 and vanishing when the orbital Hund term is absent. This splitting distinguishes the bare Kanamori problem from the Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},1-symmetric Dworin–Narath limit.

4. Renormalization, dynamical Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},2, and two-stage screening

The perturbative renormalization-group analysis is formulated in terms of one-loop scaling equations

Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},3

Among the explicit equations given are

Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},4

Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},5

Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},6

together with corresponding equations for Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},7, Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},8, and Hhyb=Vm,σcmσdmσ+h.c.,H_{\mathrm{hyb}} = V \sum_{m,\sigma} c^\dagger_{m\sigma} d_{m\sigma} + \mathrm{h.c.},9 (Horvat et al., 2016). The qualitative implication emphasized in the paper is that couplings related to orbital degrees of freedom renormalize faster than those associated primarily with spin degrees of freedom.

This flow structure produces a two-stage screening process. Orbital moments are screened first at a higher scale, while spin moments are screened later at a much lower scale. The orbital Kondo temperature is larger, whereas the spin Kondo temperature is often an order of magnitude smaller; for Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.0 above the zero-Hund Kondo scale, the separation becomes pronounced. The slow quenching of the local spin moment is thus traced both to the smallness or even ferromagnetic sign of Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.1 and to the faster renormalization of the orbital sector.

A major conceptual result is that the distinction between orbital Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.2 and orbital Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.3 symmetry becomes irrelevant in the infrared. Even when the bare couplings are split, the renormalization-group flow drives

Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.4

so that the splittings of coupling constants flow to zero. The paper concludes that the Kanamori and Dworin–Narath problems therefore describe the same low-energy physics and flow to the same low-energy Fermi-liquid fixed point.

The perturbative picture is corroborated by numerical renormalization group (NRG) calculations. At low temperatures, the finite-size spectra for Kanamori and Dworin–Narath become identical, confirming the same Fermi-liquid fixed point. The spin susceptibility stays larger and saturates later than the orbital susceptibility, directly exhibiting slower spin screening. The occupancy dependence is also significant: near Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.5, spin screening becomes especially suppressed because the local atomic state is Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.6, whereas near Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.7, Hund’s coupling becomes unimportant.

5. Low-energy mapping of the genuine three-impurity trimer

For the genuine three-impurity Anderson model, the starting point is a multi-impurity Hamiltonian

Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.8

with a host conduction band, a correlated impurity cluster, and hybridization between them (Eickhoff et al., 2020). The crucial object is the matrix hybridization function

Himp=12(U3J)Nd(Nd1)2JS2J2L2+ϵ0Nd.H_{\mathrm{imp}} = \frac{1}{2}(U-3J)N_d(N_d-1) - 2J\,\mathbf{S}^2 - \frac{J}{2}\mathbf{L}^2 + \epsilon_0 N_d.9

which encodes the influence of the host conduction band on the impurity cluster. In the low-energy, wide-band limit, the central decomposition is

NdN_d0

The real part is absorbed into the cluster Hamiltonian as an effective hopping matrix,

NdN_d1

and thereby generates delocalization and antiferromagnetic RKKY-like exchange within the cluster. The paper’s formulation is explicit that this mechanism arises directly from the conduction-band structure, without introducing ad hoc Heisenberg couplings. The imaginary part

NdN_d2

defines the hybridization matrix governing charge fluctuations and screening. After diagonalization, its nonzero eigenvalues are interpreted as the couplings of transformed cluster orbitals to effective conduction channels.

This leads to a classification of multi-impurity models into two kinds. In the first kind,

NdN_d3

so there are enough independent screening channels to fully quench all impurity moments by a multi-stage Kondo process. In the second kind,

NdN_d4

so Kondo screening alone is insufficient and the cluster must form its singlet or residual low-energy state through inter-impurity correlations. The paper identifies this distinction as a rigorous replacement for a phenomenological exhaustion criterion and uses it to organize the low-energy physics of finite impurity arrays.

For the trimer with NdN_d5, this framework is already rich enough to accommodate magnetic frustration, competition between Kondo screening and intersite exchange, underscreened intermediate-spin regimes, and frustration-induced non-Fermi-liquid fixed points. The decisive control parameters are no longer merely the bare microscopic couplings, but the interplay between the effective hopping generated by NdN_d6 and the number of effective screening channels encoded in NdN_d7.

6. Frustration, non-Fermi-liquid regimes, and KT criticality in trimers

The trimer acquires especially interesting structure in two settings emphasized by the cited work. In the idealized equilateral-triangle or NdN_d8-symmetric case, the hybridization matrix is constrained to two independent parameters, NdN_d9 and S\mathbf{S}0, and the rank remains S\mathbf{S}1 in the parameter regime studied, so this is formally a first-kind trimer (Eickhoff et al., 2020). The low-energy behavior then depends sensitively on the magnitude of the effective hopping. Small S\mathbf{S}2 gives a fully screened Fermi liquid, intermediate S\mathbf{S}3 gives a frustrated Kondo NFL regime, and large S\mathbf{S}4 destroys the NFL fixed point and restores a Fermi liquid through a crossover.

Within this S\mathbf{S}5-symmetric setting, the paper reproduces two known NFL phases: the Frustrated Kondo regime, which is stable against moderate particle-hole asymmetry, and the Isospin Kondo regime, which is unstable to particle-hole asymmetry. The mapped Anderson formulation makes clear that these regimes are fragile because the same effective hopping that generates frustration also breaks the symmetry needed to protect the isospin Kondo fixed point.

The more channel-limited situation arises for a finite trimer cluster in 1d. There the paper states that

S\mathbf{S}6

so the three-impurity problem is generically of the second kind. In this regime the effective bath cannot screen all three impurity spins independently. At or near half filling, the trimer may enter an S\mathbf{S}7 low-temperature phase, but as the filling is shifted toward the band edges, the cluster develops progressively larger ferromagnetic alignment and unscreened moments. The reported sequence for S\mathbf{S}8 is a low-temperature regime with singlet behavior, then a transition to a regime where the cluster behaves like an effective spin-S\mathbf{S}9, and ultimately, near the band edge, an impurity entropy

L\mathbf{L}0

corresponding to a residual triplet-like degeneracy in the effective low-energy description.

The physical origin of trimer frustration is traced to competition among the effective hoppings L\mathbf{L}1 generated by L\mathbf{L}2. At a special filling these hoppings become comparable, producing maximal magnetic frustration. The isolated cluster then shows a crossing of two different doublets; in the full problem this appears as a cusp in the low-energy crossover scale. Near that point, the NRG data are consistent with Kosterlitz–Thouless-type quantum phase transitions. The crossover scale is extracted from the impurity entropy through

L\mathbf{L}3

and near the critical point obeys

L\mathbf{L}4

The interpretation given is that an effective collective spin-L\mathbf{L}5 becomes decoupled from the continuum because ferromagnetic cluster correlations suppress its final-stage screening.

7. Conceptual synthesis and significance

Taken together, these results show that the phrase Three-Impurity Anderson Model names a small but conceptually diverse class of impurity problems rather than a single universal Hamiltonian. In the three-orbital single-impurity realization, the low-energy description is a multipolar Kondo model in which orbital and quadrupolar exchanges dominate over spin exchange, producing two-stage screening and a strongly suppressed spin scale. In the trimer realization, the low-energy description is a cluster-plus-effective-baths model in which the real part of the host self-energy generates inter-impurity exchange while the rank of the imaginary part limits the available screening channels (Horvat et al., 2016, Eickhoff et al., 2020).

The two strands of the literature also identify different symmetry mechanisms. In the Kanamori problem, the central result is dynamical restoration of orbital L\mathbf{L}6 at low energies despite a bare L\mathbf{L}7 interaction. In the trimer, by contrast, symmetry can stabilize or destabilize non-Fermi-liquid regimes, but the decisive distinction is between first-kind and second-kind models, i.e. whether the number of effective screening channels is sufficient to quench the cluster moments.

A plausible implication is that precise usage should distinguish between three-orbital Hund impurity physics and three-site cluster impurity physics. The former is organized by multipolar exchange and separated orbital and spin Kondo scales; the latter by frustration, effective hopping, and channel-counting constraints. Within that narrower terminology, the literature summarized here presents the trimer as the minimal setting where frustration, limited screening, and collective-moment formation can generate NFL or KT behavior, and the three-orbital Kanamori model as the canonical impurity realization of Hund’s-metal screening with dynamically emergent L\mathbf{L}8 infrared structure.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Three-Impurity Anderson Model.