# Local moment formation and Kondo screening in impurity trimers

###### Abstract

We study theoretically a triangular cluster of three magnetic impurities, hybridizing locally with conduction electrons of a metallic host. Such a cluster is the simplest to exhibit frustration – an important generic feature of many complex molecular systems in which different interactions compete. Here, low-energy doublet states of the trimer are favored by effective exchange interactions produced by strong electronic repulsion in localized impurity orbitals. Parity symmetry protects a level-crossing of such states on tuning microscopic parameters, while an avoided crossing arises in the general distorted case. On coupling to a metallic host, the behavior is shown to be immensely rich, since collective quantum many-body effects now also compete. In particular, impurity degrees of freedom are totally screened at low temperatures in a Kondo-screened Fermi liquid phase, while degenerate ground states persist in a local moment phase. Local frustration drives the quantum phase transition between the two, which may be first order or of Kosterlitz-Thouless type, depending on symmetries. Unusual mechanisms for local moment formation and Kondo screening are found, due to the orbital structure of the impurity trimer. Our results are of relevance for triple quantum dot devices. The problem is studied by a combination of analytical arguments and the numerical renormalization group.

## I Introduction and motivation

The interplay of orbital and spin degrees of freedom with electronic interactions can produce a diverse range of chemical and physical behavior. At the few-electron level of a single molecule, understanding the resulting complexity is a traditional problem in theoretical chemistry — and a challenge because strong correlations preclude an independent particle picture. At the many-electron level of clean bulk metals, by contrast, interactions are often rather unimportant, and the system adequately described by an essentially independent particle description.

Bridging between these limits is the fascinating class of ‘quantum impurity problems’, Hewson (1993) in which an interacting and, in effect, small molecular subsystem is coupled locally to the continuum of conduction electrons in a metal. The most basic example, a classic paradigm in condensed matter science,Hewson (1993); Kondo (1964); Anderson (1961) is a single magnetic impurity embedded in a metallic host. Local Coulomb repulsion favors single-occupancy of the active impurity orbital, with a local moment thus forming at high temperatures due to the free spin degree of freedom. But the Kondo effect plays a key role at low temperatures/energies: the impurity moment is screened dynamically by conduction electrons, which together form a many-body spin-singlet state (the ‘Kondo singlet’). The screening process itself can be understood in terms of a renormalization group flow, corresponding to the crossover from the local moment fixed point applicable at high temperatures, to the strong coupling fixed point describing the Kondo singlet ground state,Hewson (1993); Wilson (1975); Krishnamurthy et al. (1980); Mitchell et al. (2011a) and with physical properties exhibiting universality in terms of the crossover temperature scale, . Interestingly, there is a close connection between spin Kondo physics and dissipative tunneling relevant to electron transfer processes in chemistry, as emphasised many years ago by Jongeward and Wolynes in Ref. Jongeward and Wolynes (1983).

The Kondo effect associated with screening of a single magnetic impurity has been observed experimentally in many metals.Hewson (1993) In semi-metallic systems by contrast, where the host density of states vanishes at the Fermi level, it is well known that Kondo physics can be suppressed, and that such systems also support degenerate local moment ground states.Withoff and Fradkin (1990); Gonzalez-Buxton and Ingersent (1996); Logan and Glossop (2000); Bulla et al. (2000); Glossop and Logan (2003) Local quantum phase transitions between Kondo screened phases and local moment phases have thus attracted much interest and have been studied in detail; being sought for example in softgapped systems such as graphene,Fritz and Vojta (2013); Mitchell and Fritz (2013) d-wave superconductorsFritz and Vojta (2005) and surfaces of 3d topological insulators.Zitko (2010); Mitchell et al. (2013)

The situation is unsurprisingly richer when several impurities are present. Already in the case of two impurities, the resulting behavior can be markedly different.Jayaprakash et al. (1981); jon (); CFT (); van der Wiel et al. (2003); Chang and Chen (2009); Merino et al. (2009); Mitchell et al. (2006); Jayatilaka et al. (2011); Sela et al. (2011); Mitchell and Sela (2012); Mitchell et al. (2012); Chorley et al. (2012) If the impurities are spatially well separated, each is essentially screened independently by the Kondo effect. But an effective exchange interaction (‘RKKY’) can arise and dominate when two impurities are brought closer, with a resulting crossover to local interimpurity spin-singlet formationJayaprakash et al. (1981); jon () (indeed for strictly independent screening channels, the crossover is sharpened to a true transitionCFT (); Sela et al. (2011); Mitchell and Sela (2012); Mitchell et al. (2012)).

Multi-impurity systems, and even small molecules deposited on metallic surfaces, exhibit more complex phenomena due to competition between local interactions and the Kondo effect.Jamneala et al. (2001); Ingersent et al. (2005); Uchihashi et al. (2008); Lazarovits et al. (2005); Mitchell et al. (2009); Mitchell and Logan (2010); Wang (2007); Xiong et al. (2012); Žitko and Bonča (2008); Oguri et al. (2011); Numata et al. (2009); Seo et al. (2013) In this paper we consider a cluster of three impurities, which exhibits the frustrating effects of competing ground states. Frustration and degenerate ground states at level crossings are of course an important general feature of many complex molecular systems familiar in chemistry. Here we examine the interplay between frustration and the Kondo effect in one of its simplest realizations, to obtain a detailed understanding of the underlying physical behavior.

Real three-impurity systems, such as Cr clusters on clean Au(111) surfaces or Co clusters on Cu(100), have been studied both experimentallyJamneala et al. (2001); Uchihashi et al. (2008) and theoretically.Ingersent et al. (2005); Lazarovits et al. (2005); Kudasov and Uzdin (2002); Gotsis et al. (2006); Antal et al. (2008) But since the impurities couple to different surface sites of the host, a realistic theoretical model maps irreducibly to a three-channel problem,Ingersent et al. (2005) whose properties depend quite sensitively on details of the particular experimental realization. More general aspects of the three impurity problem can however be accessed in triple quantum dot devices,Schröer et al. (2007); Vidan et al. (2005); Rogge and Haug (2008); Gaudreau et al. (2009); Seo et al. (2013) where the metallic leads to which the dots are tunnel-coupled provide the conduction electrons. The Kondo effect has been observed in single semiconductor quantum dots coupled to metallic leads,qd: () which are often referred to as ‘artificial atoms’.Kastner (1993) By the same token, coupled dot devices behave as highly tunable ‘artificial molecules’.Jeong et al. (2001); Blick et al. (1996); Vidan et al. (2005) The primary experimental probe in all cases is of course the conductance, obtained by driving a steady current through the lead-coupled quantum dot system.

The system studied here is illustrated in Fig. 1, and consists of a triangular cluster of three quantum dots (impurities), tunnel-coupled apically to a single channel of conduction electrons (and with the metallic lead split into source and drain to allow measurement of conductance). We find both a Kondo screened phase and a free local moment phase, depending on model parameters (see Fig. 2), which themselves could be tuned in a real device by application of gate voltages.Schröer et al. (2007); Vidan et al. (2005); Rogge and Haug (2008); Gaudreau et al. (2009); Seo et al. (2013) The Kondo phase is characterized by a strong coupling Fermi liquid ground state, in which all dot degrees of freedom are screened. Zero-bias conductance through the device in this case is maximally enhanced at low temperatures due to the Kondo effect. Three distinct screening mechanisms are uncovered, depending on the relative strength of interdot tunnel-couplings, which have characteristic universal signatures in physical quantities. By contrast, the Kondo effect is totally suppressed in the local moment phase, which has a doubly-degenerate ground state and extremely low conductance.

The quantum phase transition separating the Kondo and local moment phases is studied in some detail. In the parity-symmetric case, the transition is a parity-protected level crossing.Mitchell et al. (2009) In the more general distorted case by contrast, the transition occurs at the critical end point of a line of Kondo screened states. An effective model is derived to describe this Kosterlitz-Thouless transition,Kosterlitz and Thouless (1973) which captures the characteristic vanishing of the low-energy scale, , as the transition is approached from the Kondo phase. Analytic arguments are supplemented and confirmed by exact numerics obtained using the numerical renormalization group.Wilson (1975)

## Ii Models and methods

We consider the triple quantum dot (TQD) device illustrated in Fig. 1. It consists of three equivalent and locally correlated single-level (‘Andersonian’) sites, with level energy and on-site Coulomb repulsion . Dots and are tunnel-coupled by a matrix element to form a triangular arrangement. Dot 1 is also coupled to source and drain leads (and we consider the zero-bias case where the system is in equilibrium). The Hamiltonian is decomposed as where, in standard notation,

(1a) | ||||

(1b) | ||||

(1c) | ||||

(1d) |

with the number operator for spin- electrons on dot site . describes the isolated TQD, by itself a small quasi-molecular entity that can accommodate up to 6 electrons. In the absence of electron interactions () it is in fact merely a 3-site Hückel ring, although interactions are essential, since is inversely proportional to the dot capacitance (and no quantum dot has infinite capacitance!). by contrast describes the leads, an essentially non-interacting but macroscopic metal. The hybridization term , in connecting the two subsystems via tunnel-coupling between dot 1 and the leads, ensures that the lead-coupled TQD system is an interacting, many-body problem containing macroscopic numbers of electrons. That is essential to the basic physics (despite some naive attempts to avoid it in the chemistry literature).

For equivalent leads , we define the conventional ‘local’ orbital to which dot 1 couples as,

(2) |

With this, the full TQD problem can then be mapped to an effective 1d problem by tridiagonalizing the conduction electrons,Hewson (1993); Wilson (1975) starting from the ‘zero’ orbital. One thereby obtains,

(3a) | ||||

(3b) |

with all diagonal one-electron energies if the local conduction electron density of states (DOS), , is particle-hole symmetric; with given by

(4) |

where is the free Green function for the local bath site, (correlation functions of the type are simply the Fourier transform of the usual retarded functions ).

The conduction electron and hybridization terms of the Hamiltonian enter in such quantum impurity problems through a single ‘hybridization function’,Hewson (1993) given generally by

(5) |

with , and where the real part, , follows by Hilbert transformation of the imaginary part, . As such, the free DOS and the tunnel-coupling specify completely the effect of the leads on the TQD. For simplicity, we take a wide flat DOS, which is the conventional and generic caseHewson (1993) relevant to most metallic hosts: , such that is defined inside the band, , and .

### ii.1 Numerical Renormalization Group

The Numerical Renormalization GroupWilson (1975) (NRG) is a non-perturbative technique for treating quantum impurity problems, such as the lead-coupled TQD system (for a recent review, see Ref. Bulla et al., 2008). Numerically exact thermodynamic and dynamic quantities can be calculated with NRG on essentially all relevant energy scales — from the conduction electron bandwidth (typically a few eV in bulk metals), down to the Kondo scale (which is an exponentially-smallHewson (1993) fraction of ).

The first stepWilson (1975) is to divide the conduction electron DOS into intervals whose width decreases exponentially (the points separating intervals are , with and ). The spectrum is then discretized by replacing the DOS in each interval by a single pole of the same total weight. The ‘Wilson chain’Wilson (1975); Krishnamurthy et al. (1980) is then defined by writing the conduction electron Hamiltonian as a 1d tight-binding chain of form Eq. (3a), and parameters , chosen so that the free DOS at the TQD site corresponds to the discretized spectrum. This discretized model is then diagonalized iteratively, with high-lying states discarded at each step (the truncation being justified because the strength of the hoppings decrease exponentially down the chainWilson (1975); Krishnamurthy et al. (1980)).

One thus examines the behavior of the system on progressively lower energy scales, and thermodynamics can be built up as a function of temperature.Wilson (1975); Bulla et al. (2008); Krishnamurthy et al. (1980) Dynamical quantities, such as spectral functions, can also be obtained from the full density matrix,Peters et al. (2006); Weichselbaum and von Delft (2007) calculated iteratively in the Anders-Schiller basis.Anders and Schiller (2005) For the calculations presented in this paper, we use a discretization parameter , and retain around states at each iteration.

### ii.2 Physical quantities

We consider below the contribution to thermodynamics arising from the lead-coupled TQD system, specifically the so-called excess quantitiesWilson (1975); Krishnamurthy et al. (1980) , with denoting the thermal average in the absence of the TQD itself. We focus on the entropy , uniform spin susceptibility , and ‘excess’ charge (here and refer to the spin and charge of the entire system). These quantities show characteristic signatures of the underlying fixed points, reached under renormalization on progressive reduction of the temperature/energy scale. The renormalization group flow between such fixed points shows up as crossover behavior in thermodynamics, which thus allows determination of emergent energy scales in the problem, such as the Kondo temperature, .

We also consider dynamical properties — in particular dot spectral functions (or, loosely, local densities of states), given by in terms of the local Green function for that dot, . In the zero-bias (equilibrium) limit of interest, it is naturally the local spectrum of dot 1 (itself directly tunnel-coupled to the leads) which determines the differential conductance between source and drain leads, mediated via the TQD. This is given exactlyMeir and Wingreen (1992) by

(6) |

where is the Fermi function (and is the Fermi level). Here, , where the dimensionless quantity embodies simply the relative coupling strength of dot 1 to source and drain leads; such that is maximal in the symmetric case . At zero temperature, the conductance is thus controlled by the behavior of the spectrum at the Fermi level,

(7) |

where the final identity follows from the definition of the scattering phase shift , together with the fact that all electron scattering vanishes at the Fermi level.

## Iii Results: symmetric case

We consider first the full TQD model, Eq. (1), in the parity-symmetric case and (see Fig. 1). In this limit the Hamiltonian is invariant to swapping the dot labels , meaning formally that for the permutation operator . All states can thus be classified according to parity (since ). As shown in Ref. Mitchell et al., 2009 this allows for the possibility of a (zero temperature) quantum phase transition (QPT) between two parity-distinct phases, which are of Kondo screened and local moment types. That scenario is explored in detail below, focusing primarily on the strongly correlated case . Since the most interesting physical behavior arises when the TQD is in essence triply occupied, we consider the representative case . The largest energy scale of the problem, the bandwidth , is taken as .

### iii.1 Phase diagram

Before discussing specific results for the symmetric TQD, in Fig. 2 we present the phase diagram which highlights the relevant regimes, crossover scales and phases, verified directly by full NRG calculations. The parity-symmetric TQD model supports two phases: a Kondo screened phase with strong coupling (SC) ground state, and a local moment (LM) phase with a doubly-degenerate ground state.Mitchell et al. (2009)

As shown below, however, the effective spins which remain free in the LM phase depend on the underlying parameters in the -plane. Likewise, Kondo screening proceeds by different mechanisms in the SC phase depending on relative magnitudes of the interdot tunnel-couplings. When the interdot tunnel-couplings are large compared to the interaction, (regime (a) of Fig. 2), the relevant TQD states are simple molecular orbitals (MOs); with the lowest doublet state remaining free down to when . By contrast, for each dot is essentially singly-occupied. In this case, dot 1 (connected directly to the leads) can undergo the Kondo effect on the scale of — provided the interdot tunnel-couplings are weak (). The strength of the effective coupling between the remaining dots 2 and 3 then determines the ultimate ground state of the system, with a level-crossing QPT occurring between LM and SC states as is increased through (regime (c) of Fig. 2). If by contrast the interdot couplings are stronger, , the lowest-energy TQD doublet state has weight on all three dots. This results – regime (b) of Fig. 2 – in an effective antiferromagnetic coupling to the leads if , in which case the Kondo effect is operative and the SC phase obtains; or, if , in an effective ferromagnetic coupling to the leads, such that the TQD doublet remains free down to in this LM phase.

In the following, we discuss in some detail the behavior in these three distinct regimes.

### iii.2 Molecular orbital regime

We examine first the simplest case, in which the interdot tunnel-couplings are strong, . Here, the isolated TQD states are essentially non-interacting MOs. To obtain a handle on the problem in this regime, we thus derive a low-energy effective model for the non-interacting TQD (), then incorporate the effect of interactions within lowest-order perturbation theory in .

When , a simple canonical transformation of the dot operators brings in Eq. 1 to diagonal form,

(8) |

where (). The single-particle levels (MO energies) are given by and , and for have the relative ordering . The difference in energy between the 4-electron state and the 3-electron state is thus . The lowest many-particle TQD state is thus always a 4-electron state when .

Intuition suggests however that the Coulomb repulsion will favor a 3-electron state. Treating perturbatively the interaction part of the Hamiltonian, , we thus calculate the correction to the isolated TQD states to first-order in . The resultant energy difference between the 4- and 3-electron states is then , with . Hence for as considered explicitly, , indicating that the 3-electron TQD state is indeed favored when .

Of course, the interesting behavior arises on coupling the TQD to the leads. In the relevant case of (), one can obtain an effective model valid at low temperatures/energies , by projecting the full Hamiltonian onto the 4-electron TQD state (for ) or the 3-electron state (for ). This is achieved by a Schrieffer-Wolff transformation (SWT)Schrieffer and Wolff (1966) to second order in , eliminating perturbatively excitations to higher-lying TQD states.

In the effective 4-electron sector, the spin-singlet TQD state essentially decouples from the leads. This ‘frozen impurity’ fixed point is continuously connected to the Kondo screened ground statesHewson (1993); Krishnamurthy et al. (1980) (discussed further below). More interestingly, in the 3-electron sector of the TQD, one obtains via a SWT the effective model

(9) |

with a spin- operator for the lowest (doublet) 3-electron TQD state, and the spin density of conduction electrons at the TQD. Eq. 9 is a model of Kondo form,Hewson (1993) with its underlying physics well known to depend crucially on the sign of the exchange coupling, . Importantly, we find here , wherein the ferromagnetic Kondo effect arises.Hewson (1993); Mehta et al. (2005); Koller et al. (2005) Kondo quenching of the TQD spin- is thus inoperative, with the TQD decoupling from the leads on renormalization under reduction of the temperature/energy scale, leaving asymptotically a free local moment as . The residual entropy at the corresponding LM fixed point is thus , and leading irrelevant corrections to the fixed point are non-analytic.Mehta et al. (2005); Koller et al. (2005)

Since the ground states for (LM) and (SC) cannot be continuously connected, a QPT is expected around [see arrow (a) of Fig. 2]. This physical picture is confirmed directly in Fig. 3, where NRG results for the full TQD model are presented. Specifically, we show the TQD contribution to the entropy as a function of temperature, for systems approaching the QPT between LM and SC phases. For , the 3- and 4- electron states are quasi-degenerate, and a entropy thus results. Flow to either the LM ground state with or the SC ground state with occurs on the scale of , which as such vanishes linearly with as the transition at the critical is approached from either side (and which linearity is symptomatic of a parity-protected level-crossing QPT Mitchell et al. (2009)).

The inset shows the excess chargeHewson (1993); Krishnamurthy et al. (1980) due to the TQD, , for the same systems; confirming that the transition at occurs between effective 3- and 4-electron TQD states. These results also show that precisely at the QPT (curve (g) in Fig. 3), the problem can be understood in terms of coexisting SC and LM fixed points – for example the charge , is simply an average of the triple charge characteristic of the doubly-degenerate LM fixed point and the quadruple charge of the singly-degenerate SC fixed point.

### iii.3 Spin regime: strong interimpurity coupling

More subtle behavior arises in the strongly-correlated case , where each dot is essentially singly-occupied for temperatures . However, the tunnel-couplings do generate effective exchange couplings between the dots. This can be seen by projecting (SWT) the full onto the singly-occupied manifold of TQD states, perturbatively eliminating virtual excitations to 2- and 4-electron TQD states to second-order in the tunnel-couplings , and . The resulting effective spin model is then given by

(10) |

where is a spin- operator for dot , is the conduction electron spin density at the TQD as before, and the exchange couplings are given by

(11) |

Two parity-distinct doublet ground states are then obtained for the isolated TQD, depending on whether (even parity) or (odd parity). In Ref. Mitchell et al., 2009 we derived analytically an effective low-energy model to describe this lead-coupled TQD system for but , corresponding to regime (b) in Fig. 2. In either case, or , the effective low-energy model is again a Kondo model Eq. 9, describing the residual coupling of the lowest TQD doublet state to the leads. Importantly however, for and for . The ferromagnetic Kondo effectHewson (1993) thus again arises for , with the lowest TQD doublet decoupling asymptotically at low energies. However, for the regular spin- antiferromagnetic Kondo effectHewson (1993) drives the system to a strong coupling state below the Kondo scale , in which all TQD degrees of freedom are quenched by the lead conduction electrons. The entropy is thus , corresponding to the Kondo singlet ground state.Hewson (1993); Wilson (1975); Krishnamurthy et al. (1980) On tuning through , a level-crossing QPT occurs between LM and SC phases.Mitchell et al. (2009)

One can in fact derive more generally an effective low-energy Kondo model of form Eq. 9, working directly with (i.e. without first projecting into the spin sector, as above). The effective Kondo model is obtained simply by projecting the full Hamiltonian onto the reduced TQD Hilbert space spanned by the lowest eigenstate of the numerically-diagonalized . For any the odd-parity TQD doublet couples to the leads, while the even parity doublet couples for . The QPT can of course also be realized by tuning , keeping fixed (Fig. 2), and a numerical calculation of the effective coupling so obtained, , is shown in Fig. 4. The asymptotic valuesMitchell et al. (2009) and in the LM and SC phases are recovered in the singly-occupied limit; and the discontinuous change in the sign of is seen clearly at the point , due to the level crossing of TQD states. NRG calculations on the full model (without any low-energy projections) indeed confirm that the transition occurs very close to .

### iii.4 Spin regime: weak interimpurity coupling

As suggested by arrow (c) in the phase diagram Fig. 2, a QPT between LM and SC phases arises even when the interdot couplings and are very small. In this case however, the mechanism for local moment formation and Kondo screening is rather subtle.

Since , , , the dots are still essentially singly-occupied, so the effective spin Hamiltonian Eq. 10 remains valid. In the limit where (arising when ), dot 1 decouples from dots 2 and 3, and undergoes the regular Kondo effect with the lead to which it is coupled, being screened below (with the Kondo scale as above).Hewson (1993) This must also in fact remain the case for small but non-vanishing , since the quenched Kondo singlet is already formed on the scale of . However, this fixed point contains residual unquenched degrees of freedom corresponding to the spins on dots 2 and 3. In order to analyze the stability of this fixed point, one must thus determine whether there is an effective coupling of these degrees of freedom to the remaining Fermi liquid bath states of the lead.

The effective exchange coupling acting directly between dots 2 and 3 is given simply by (see Eqs. 10, 11), such that the local singlet state is lower in energy than the local triplet by . However, there is an additional effective coupling between dots 2 and 3 due to an RKKY-type interaction mediated by the Kondo singlet formed between dot 1 and the leads. Virtual polarization of this Kondo singlet is readily shown to generate an effective ferromagnetic contribution to the coupling between the spins of dots 2 and 3; with second-order perturbation theory in the coupling within a Wilson chain formalismMitchell et al. (2011b); Mitchell et al. (2012) thereby yielding a renormalized singlet-triplet splitting, (with an undetermined positive constant). When , the 2-3 singlet thus lies lowest and decouples from the rest of the system. In consequence, at temperatures the entire system is in a singlet state, characteristic of the SC phase of regime (c) in Fig. 2.

By contrast, when the triplet formed between dots 2 and 3 lies lowest in energy. Thus in this case there are still residual TQD degrees of freedom at temperatures . To determine the stability of this state, we must again consider the effective coupling between this triplet and the rest of the system. Within the Wilson chain picture, it can be shown that the TQD triplet experiences an effective coupling to the remaining Fermi liquid bath states of the lead, mediated via the Kondo singlet formed with dot 1. The mechanism here is in fact completely analogous to that occurring in odd quantum dot chains with weak interdot coupling;Mitchell et al. (2011b) or in the asymmetric two-impurity Kondo problem.Mitchell et al. (2012) For , the effective low-energy model follows as

(12) |

where is now a spin-1 operator for the residual TQD state, is the free conduction electron Hamiltonian with the Wilson chain ‘zero orbital’ removed due to the first-stage Kondo screening involving dot 1; and is its spin density at the TQD. The key result is that the effective coupling is now antiferromagnetic, and the effective low-energy model Eq. 12 is thus a realization of the famous single-channel spin-1 Kondo model of Nozières and Blandin Nozières and Blandin (1980).

The rich physical behavior of the spin-1 Kondo modelNozières and Blandin (1980); Mehta et al. (2005); Koller et al. (2005); Mitchell et al. (2009); Logan et al. (2009); Wright et al. (2011) is thus expected on the lowest energy scales. The antiferromagnetic coupling, is renormalized upward on reduction of the temperature scale, resulting in quenching (or ‘underscreening’) of the ‘impurity’ below () Nozières and Blandin (1980); such that the ground state comprises a residual free spin- local moment, with leading RG-irrelevant corrections to the fixed point that are ferromagnetic (and non-analytic).Mehta et al. (2005); Koller et al. (2005)

In consequence, there is again a QPT on tuning , separating a Kondo screened SC phase from a LM phase (regime (c) of Fig. 2). The transition itself is expected to occur at , implying [or equivalently , see Fig. 2, arrow (c)].

To illustrate this complex behavior, Fig. 5 shows full NRG results for the TQD contribution to entropy and magnetic susceptibility for systems close to the transition. Both systems have and similar . The solid line corresponds to the LM phase, and the dashed line to the SC phase. On the temperature scale , the dots become singly occupied, and so and corresponding to three free spins-. As the temperature is reduced below , dot 1 is screened by the Kondo effect, leaving two quasi-degenerate spins on dots 2 and 3 (yielding thereby and ). For (dashed line), the residual TQD singlet state is lowest, and so and when . However, for (solid line), the crossover is first to a free residual TQD triplet state on the scale of ; then below to the ultimate stable LM fixed point describing the underscreened Kondo state, with and .

Similar behavior is observed in dynamical quantities, such as the spectral function for dot 1, the energy/frequency dependence of which is shown in Fig. 6. In the SC phase, the classic three-peak structure is observed. At high energies , the only spectral features are the Hubbard satellites, whose origin is simply dot charge fluctuations.Hewson (1993) At low energies , the narrow Kondo resonance is observed (see in particular inset (A), dashed line); with the unitarity limit reached at the Fermi level, , being characteristic of the SC fixed point. In the LM phase however (solid lines), first-stage Kondo screening of dot 1 on the scale of is followed by second stage underscreening of the residual TQD triplet state, which as above is mediated via dot 1. This results in a slow crossover on the scale of to the LM fixed point, characterized by behavior, such that only at , characteristic of the singular Fermi liquid, and reflecting the marginally irrelevant corrections to the fixed point.Mehta et al. (2005); Koller et al. (2005); Mitchell et al. (2009); Logan et al. (2009); Wright et al. (2011)

Finally, in Fig. 7 we show the zero-bias conductance through the device as a full function of temperature. Solid lines are for in the SC phase, approaching progressively closely (from (d)–(a)) the transition at ; while dashed lines are for in the LM phase. All systems have a common scale , and the data are rescaled in terms of (with chosen arbitrarily as the for which ). The phase shift in the Kondo screened SC phaseHewson (1993) implies a unitarity conductance from Eq. 7; while in the LM phaseLogan et al. (2009) gives vanishing conductance,Logan et al. (2009) . The behavior at the transition fixed point itself can again be understood here in terms of coexisting LM and SC states; yielding a conductance at the transition of . The full temperature dependence is of course rich, reflecting with decreasing renormalization group flow to the Kondo SC fixed point on the scale of , followed by flow to the transition fixed point at ; and finally flow to the stable LM or SC fixed points describing the true ground state, on the scale .

## Iv Results: distortions

We now turn to the more general case of broken parity symmetry, as occurs due to distortions of the triangular TQD structure (Fig. 1). The presence of symmetry-breaking perturbations, such as , of course precludes labeling states by a parity quantum number, since . In the symmetric case considered above, the QPT between LM and SC phases was characterized by a level crossing between parity-distinct states. We show below that the parity-broken model still supports LM and SC phases, but that the level crossing transition becomes instead a QPT of Kosterlitz-Thouless form.Kosterlitz and Thouless (1973) The schematic phase diagram, Fig. 2, still applies in the general case, but with the phase boundary now to be understood as a line of SC critical end points, with the Kondo scale vanishing as the transition is approached from the SC phase.

### iv.1 Effective models

We focus here on the strongly correlated case of primary interest, , . At temperatures , each dot becomes essentially singly-occupied, and an effective TQD spin model analogous to Eq. 10 can similarly be derived via a SWT upon perturbative elimination of virtual excitations to 2- and 4-electron TQD states. For ease of comparison with the parity-symmetric case, we retain and , defining the asymmetry in terms of the parameter . The resultant effective low-energy model then has the form

(13) |

where the effective exchange couplings , and are given in Eq. 11; and where the asymmetry enters only through the final term, with

(14) |

The low-energy manifold of the isolated TQD again comprises a pair of doublet states, denoted and (with an essentially irrelevant spin quartet occurring at higher energies). For any degree of asymmetry, these states may be expressed in terms of even- and odd-parity doublets (themselves obtained in the symmetric limit considered above). They are given by

(15a) | ||||

(15b) |

where

(16) |

The energy difference between these isolated TQD doublet states is,

(17) |

such that for all finite values of , and (or , and ). The effect of asymmetry is thus to turn the level crossing of TQD states at to an avoided crossing when . There is in otherwords a single, unique doublet TQD ground state (the state) for all model parameters in the asymmetric case. This is illustrated in the inset to Fig. 8, where we show the energies of the TQD doublet states as a function of for different values of the asymmetry (the solid line (a) is for the symmetric case ).

The implication of Eq. 17 is that below the temperature/energy scale , only the TQD doublet state is accessible. In consequence one can derive a low-energy effective model, valid for , by projecting the spin model Eq. 13 onto the ground state TQD doublet manifold spanned by . This can be done to first order in the Kondo coupling , leading to

(18) |

where is a spin- operator for the doublet TQD state ; and where the effective exchange coupling to the leads follows as

(19) |

From Eq. 16, we immediately find that is antiferromagnetic for , but is ferromagnetic, , for (and one also recovers asymptotically the results of the parity-symmetric case, , for ; and for ). The effective coupling as a function of is plotted in Fig. 8 as the solid line; points are extracted from full NRG results, and show very good agreement.

Importantly, the antiferromagnetic effective Kondo coupling vanishes continuously as [or ] for any ,

(20) |

(in contrast to the discontinuous change in the effective Kondo coupling associated with a level crossing QPT, see e.g. Fig. 4). In consequence, in the full model we expect a Kondo SC phase for (), with TQD degrees of freedom entirely quenchedHewson (1993) below a Kondo scale , which itself vanishes continuously as [or )]:

(21) |

The form of the vanishing Kondo scale as the transition is approached, Eq. 21, is that arising for a Kosterlitz-ThoulessKosterlitz and Thouless (1973) (KT) transition. The point can thus be understood as the critical end point of a line of Kondo SC fixed points. For , there is no low-energy scale; and the residual TQD doublet degree of freedom remains free down to at this LM fixed point.

### iv.2 Kosterlitz-Thouless transition

This is explored further in Figs. 9 and 10, containing NRG results for the full TQD model as the KT transition is approached from the SC phase, .

Fig. 9 shows thermodynamic quantities and as a function of temperature. On the scale , the dots become singly-occupied, and so the entropy drops to , with a corresponding tripled Curie law susceptibility, . The quartet TQD state becomes inaccessible on the scale leaving the pair of quasi-degenerate TQD doublet states (and hence ). The higher lying doublet is in turn projected out on the scale of . This is the LM fixed point, with and . The lowest TQD doublet is then screened by the Kondo effect on the scale of , which itself vanishes as the transition is approached (as evident directly from the figure). At the SC fixed point, TQD degrees of freedom are entirely quenched, giving and .

Similar behavior is observed in the spectral function of dot 1 shown in Fig. 10, together (inset) with the -dependence of the zero-bias conductance through the TQD. On the energy scale , the higher energy TQD doublet state is projected out, terminating the renormalization of its coupling to the leads. Instead, the lowest TQD doublet flows to strong coupling, with the characteristic narrow Kondo resonance evident in the spectrum on the scale of ; and reaching the unitarity limit that is characteristicHewson (1993) of the strong coupling fixed point. The Kondo resonance narrows progressively as diminishes and the QPT is approached from the SC phase; vanishing continuously ‘on the spot’ at the transition itself, such that in (and throughout) the LM phase . Corresponding behavior is naturally seen in the -dependence of the conductance, which is determined by the local spectrum as in Eq. (6). Throughout the SC phase, it is likewise enhanced at low-temperatures due to the Kondo effect, with in all cases. In the LM phase by contrast, where , the conductance vanishes – the conductance thus changing discontinuously as the transition is crossed, which appears to be a rather general signature of KT transitions in quantum dot and related systems.Hofstetter and Schoeller (2002); Garst et al. (2004); Galpin et al. (2005, 2006); Logan et al. (2009)

## V Summary and Discussion

A compact cluster of three quantum impurities, hybridizing apically with a single channel of host conduction electrons, has been shown to exhibit a rich range of physical behavior. Such a model may describe apex-coupled impurity trimers on metallic surfaces, or semiconductor TQD devices.

In the parity-symmetric case where two of the impurities/dots are in equivalent local environments, a level-crossing quantum phase transition separates a Kondo screened strong coupling (SC) phase and a free local moment (LM) phase. For strong interdot tunnel-coupling (and hence weak interactions), the transition in essence occurs between states with either 3 or 4 electrons occupying TQD molecular orbitals. For weaker couplings and strong dot electron correlations, sites of the TQD are essentially singly-occupied. Within this TQD spin-regime, the QPT between SC and LM phases again arises; although different mechanisms are uncovered, depending on the relative interdot tunneling strengths. For weak interdot couplings in particular, the Kondo effect screens a single dot, with the residual TQD degrees of freedom forming either a local singlet (SC) or an underscreened spin-1 Kondo state (LM). When parity symmetry is broken by distortions of the triangular TQD structure, we show that LM and SC phases are still supported — with the transition between them now of Kosterlitz-Thouless type, such that the low-energy Kondo scale vanishes continuously as the transition is approached from the SC phase.

Given the ubiquity of SC and LM phases, an obvious question arises: what characterizes – and distinguishes – such phases more generally? One answer is to note that the zero-bias conductance is pinned in the singly-occupied case to a unitarity value in the Kondo SC phase, but vanishes in the LM phase, . The origin of this can be traced via Eq. 7 to the scattering phase shift, , experienced by conduction electrons in either phase.

By a straightforward extension to the TQD model of the analysis given in Ref.Logan et al., 2009, a Friedel-Luttinger sum ruleLogan et al. (2009); Luttinger and Ward (1960); Langreth (1966) can be shown to relate exactly the phase shift to the excess charge () via ; where the Luttinger integral is given by

(22) |

(with and the matrices for the lead-coupled TQD Green functions and self-energies respectively). We find that the SC and LM phases are each characterized by a distinct value for the ‘topological’ quantity : from extensive NRG calculations of the Luttinger integral, and regardless of the bare underlying model parameters, we find to be characteristic of the SC phase (as indeed is well known for any Fermi liquid phaseLuttinger and Ward (1960)), and to be equally characteristic of the LM phase. The latter is also precisely as found throughout the entire LM phase of a correlated 2-level quantum dot modelLogan et al. (2009), suggesting that the Luttinger integral is indeed a universal characteristic of a LM phase – albeit that the fundamental reasons for this await an answer.

###### Acknowledgements.

We are grateful for funding from EPSRC (UK) through grant EP/1032487/1. One of us in particular (DEL) would also like to salute Peter Wolynes on the occasion of his 60 birthday, and to acknowledge with gratitude the wealth of inspiration and insight provided by him over more than 30 years. Plus a lot of fun (which is not to be underestimated in science!).## References

- Hewson (1993) A. C. Hewson, The Kondo Problem to Heavy Fermions (Cambridge University Press, Cambridge, 1993).
- Kondo (1964) J. Kondo, Prog. Theor. Phys. 32, 37 (1964).
- Anderson (1961) P. W. Anderson, Phys. Rev. 124, 41 (1961).
- Wilson (1975) K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975).
- Mitchell et al. (2011a) A. K. Mitchell, M. Becker, and R. Bulla, Phys. Rev. B 84, 115120 (2011a).
- Krishnamurthy et al. (1980) H. R. Krishnamurthy, J. W. Wilkins, and K. G. Wilson, Phys. Rev. B 21, 1003, 1044 (1980).
- Jongeward and Wolynes (1983) G. A. Jongeward and P. G. Wolynes, J. Chem. Phys. 79, 3517 (1983).
- Withoff and Fradkin (1990) D. Withoff and E. Fradkin, Phys. Rev. Lett. 64, 1835 (1990).
- Gonzalez-Buxton and Ingersent (1996) C. Gonzalez-Buxton and K. Ingersent, Phys. Rev. B 54, 15614 (1996).
- Logan and Glossop (2000) D. E. Logan and M. T. Glossop, J. Phys.: Condens. Matter 12, 985 (2000).
- Bulla et al. (2000) R. Bulla, M. T. Glossop, D. E. Logan, and T. Pruschke, J. Phys.: Condens. Matter 12, 4899 (2000).
- Glossop and Logan (2003) M. T. Glossop and D. E. Logan, J. Phys.: Condens. Matter 15, 7519 (2003).
- Fritz and Vojta (2013) L. Fritz and M. Vojta, Rep. Prog. Phys. 76, 032501 (2013).
- Mitchell and Fritz (2013) A. K. Mitchell and L. Fritz, Phys. Rev. B 88, 075104 (2013).
- Fritz and Vojta (2005) L. Fritz and M. Vojta, Phys. Rev. B 72, 212510 (2005).
- Mitchell et al. (2013) A. K. Mitchell, D. Schurucht, M. Vojta, and L. Fritz, Phys. Rev. B 87, 075430 (2013).
- Zitko (2010) R. Zitko, Phys. Rev. B 81, 241414(R) (2010).
- (18) B. A. Jones, C. M. Varma, and J. W. Wilkins, Phys. Rev. Lett. 61, 125 (1988); B. A. Jones, Physica B (Amsterdam) 171, 53 (1991).
- (19) I. Affleck and A. W. W. Ludwig, Phys. Rev. Lett. 68, 1046 (1992); I. Affleck, A. W. W. Ludwig and B. A. Jones, Phys. Rev. B 52, 9528 (1995).
- Jayatilaka et al. (2011) F. W. Jayatilaka, M. R. Galpin, and D. E. Logan, Phys. Rev. B 84, 115111 (2011).
- Sela et al. (2011) E. Sela, A. K. Mitchell, and L. Fritz, Phys. Rev. Lett. 106, 147202 (2011).
- Mitchell and Sela (2012) A. K. Mitchell and E. Sela, Phys. Rev. B 85, 235127 (2012).
- Mitchell et al. (2012) A. K. Mitchell, E. Sela, and D. E. Logan, Phys. Rev. Lett. 108, 086405 (2012).
- Jayaprakash et al. (1981) C. Jayaprakash, H. R. Krishnamurthy, and J. W. Wilkins, Phys. Rev. Lett. 47, 737 (1981).
- van der Wiel et al. (2003) W. G. van der Wiel et al., Rev. Mod. Phys. 75, 1 (2003).
- Chang and Chen (2009) A. M. Chang and J. C. Chen, Rep. Prog. Phys. 72, 096501 (2009).
- Merino et al. (2009) J. Merino, L. Borda, and P. Simon, Europhys. Lett. 85, 47002 (2009).
- Chorley et al. (2012) S. J. Chorley, M. R. Galpin, F. W. Jayatilaka, C. G. Smith, D. E. Logan, and M. R. Buitelaar, Phys. Rev. Lett. 109, 156804 (2012).
- Mitchell et al. (2006) A. K. Mitchell, M. R. Galpin, and D. E. Logan, Europhys. Lett. 76, 95 (2006).
- Jamneala et al. (2001) T. Jamneala, V. Madhavan, and M. F. Crommie, Phys. Rev. Lett. 87, 256804 (2001).
- Ingersent et al. (2005) K. Ingersent, A. W. W. Ludwig, and I. Affleck, Phys. Rev. Lett. 95, 257204 (2005).
- Mitchell et al. (2009) A. K. Mitchell, T. F. Jarrold, and D. E. Logan, Phys. Rev. B 79, 085124 (2009).
- Mitchell and Logan (2010) A. K. Mitchell and D. E. Logan, Phys. Rev. B 81, 075126 (2010).
- Uchihashi et al. (2008) T. Uchihashi, J. Zhang, J. Kröger, and R. Berndt, Phys. Rev. B 78, 033402 (2008).
- Lazarovits et al. (2005) B. Lazarovits, P. Simon, G. Zaránd, and L. Szunyogh, Phys. Rev. Lett. 95, 077202 (2005).
- Wang (2007) W. Z. Wang, Phys. Rev. B 76, 115114 (2007).
- Xiong et al. (2012) Y. C. Xiong, J. Huang, and W.-Z. Wang, J. Phys.: Condens. Matter 24, 455604 (2012).
- Žitko and Bonča (2008) R. Žitko and J. Bonča, Phys. Rev. B 77, 245112 (2008).
- Oguri et al. (2011) A. Oguri, S. Amaha, Y. Nishikawa, T. Numata, M. Shimamoto, A. C. Hewson, and S. Tarucha, Phys. Rev. B 83, 205304 (2011).
- Numata et al. (2009) T. Numata, T. Nisikawa, A. Oguri, and A. C. Hewson, Phys. Rev. B 80, 155330 (2009).
- Seo et al. (2013) M. Seo et al., Phys. Rev. Lett. 110, 046803 (2013).
- Kudasov and Uzdin (2002) Y. B. Kudasov and V. M. Uzdin, Phys. Rev. Lett. 89, 276802 (2002).
- Gotsis et al. (2006) H. J. Gotsis, N. Kioussis, and D. A. Papaconstantopoulos, Phys. Rev. B 73, 014436 (2006).
- Antal et al. (2008) A. Antal, B. Lazarovits, L. Szunyogh, B. Ujfalussy, and P. Weinberger, Phys. Rev. B 77, 174429 (2008).
- Schröer et al. (2007) D. Schröer, A. D. Greentree, L. Gaudreau, K. Eberl, L. C. L. Hollenberg, J. P. Kotthaus, and S. Ludwig, Phys. Rev. B 76, 075306 (2007).
- Vidan et al. (2005) A. Vidan, R. Westervelt, M. Stopa, M. Hanson, and A. Gossard, J. Supercond. Incorp. Novel Magn. 18, 223 (2005).
- Rogge and Haug (2008) M. C. Rogge and R. J. Haug, Phys. Rev. B 77, 193306 (2008).
- Gaudreau et al. (2009) L. Gaudreau, A. Kam, G. Granger, S. A. Studenikin, P. Zawadzki, and A. S. Sachrajda, Appl. Phys. Lett. 95, 193101 (2009).
- (49) D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. Abusch-Magder, U. Meirav, and M. A. Kastner, Nature 391, 156 (1998); S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwenhoven, Science 540, 281 (1998).
- Kastner (1993) M. A. Kastner, Physics Today 46, 24 (1993).
- Jeong et al. (2001) H. Jeong, A. M. Chang, and M. R. Melloch, Science 293, 2221 (2001).
- Blick et al. (1996) R. H. Blick, R. J. Haug, J. Weis, D. Pfannkuche, K. v. Klitzing, and K. Eberl, Phys. Rev. B 53, 7899 (1996).
- Kosterlitz and Thouless (1973) J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973).
- Bulla et al. (2008) R. Bulla, T. Costi, and T. Pruschke, Rev. Mod. Phys. 80, 395 (2008).
- Peters et al. (2006) R. Peters, T. Pruschke, and F. B. Anders, Phys. Rev. B 74, 245114 (2006).
- Weichselbaum and von Delft (2007) A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 076402 (2007).
- Anders and Schiller (2005) F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 (2005).
- Meir and Wingreen (1992) Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992).
- Schrieffer and Wolff (1966) J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (1966).
- Mehta et al. (2005) P. Mehta, N. Andrei, P. Coleman, L. Borda, and G. Zarand, Phys. Rev. B 72, 014430 (2005).
- Koller et al. (2005) W. Koller, A. C. Hewson, and D. Meyer, Phys. Rev. B 72, 045117 (2005).
- Mitchell et al. (2011b) A. K. Mitchell, D. E. Logan, and H. R. Krishnamurthy, Phys. Rev. B 84, 035119 (2011b).
- Nozières and Blandin (1980) P. Nozières and A. Blandin, J. Phys. (Paris) 41, 193 (1980).
- Logan et al. (2009) D. E. Logan, C. J. Wright, and M. R. Galpin, Phys. Rev. B 80, 125117 (2009).
- Wright et al. (2011) C. J. Wright, M. R. Galpin, and D. E. Logan, Phys. Rev. B 84, 115308 (2011).
- Hofstetter and Schoeller (2002) W. Hofstetter and H. Schoeller, Phys. Rev. Lett. 88, 016803 (2002).
- Garst et al. (2004) M. Garst, S. Kehrein, T. Pruschke, A. Rosch, and M. Vojta, Phys. Rev. B 69, 214413 (2004).
- Galpin et al. (2005) M. R. Galpin, D. E. Logan, and H. R. Krishnamurthy, Phys. Rev. Lett. 94, 186406 (2005).
- Galpin et al. (2006) M. R. Galpin, D. E. Logan, and H. R. Krishnamurthy, J. Phys.: Condens. Matter 18, 6545 (2006).
- Luttinger and Ward (1960) J. M. Luttinger and J. C. Ward, Phys. Rev. 118, 1417 (1960).
- Langreth (1966) D. C. Langreth, Phys. Rev. 150, 1516 (1966).