Kondo blockade due to quantum interference in single-molecule junctions

Abstract

Molecular electronics offers unique scientific and technological possibilities, resulting from both the nanometre scale of the devices and their reproducible chemical complexity. Two fundamental yet different effects, with no classical analogue, have been demonstrated experimentally in single-molecule junctions: quantum interference due to competing electron transport pathways, and the Kondo effect due to entanglement from strong electronic interactions. Here we unify these phenomena, showing that transport through a spin-degenerate molecule can be either enhanced or blocked by Kondo correlations, depending on molecular structure, contacting geometry and applied gate voltages. An exact framework is developed, in terms of which the quantum interference properties of interacting molecular junctions can be systematically studied and understood. We prove that an exact Kondo-mediated conductance node results from destructive interference in exchange-cotunneling. Nonstandard temperature dependences and gate-tunable conductance peaks/nodes are demonstrated for prototypical molecular junctions, illustrating the intricate interplay of quantum effects beyond the single-orbital paradigm.

ARTICLE

Received 20 Jul 2016 | Accepted 6 Mar 2017 | Published 11 May 2017

Molecular electronics offers unique scientic and technological possibilities, resulting from both the nanometre scale of the devices and their reproducible chemical complexity. Two fundamental yet different effects, with no classical analogue, have been demonstrated experimentally in single-molecule junctions: quantum interference due to competing electron transport pathways, and the Kondo effect due to entanglement from strong electronic interactions. Here we unify these phenomena, showing that transport through a spin-degenerate molecule can be either enhanced or blocked by Kondo correlations, depending on molecular structure, contacting geometry and applied gate voltages. An exact framework is developed, in terms of which the quantum interference properties of interacting molecular junctions can be systematically studied and understood. We prove that an exact Kondo-mediated conductance node results from destructive interference in exchange-cotunneling. Nonstandard temperature dependences and gate-tunable conductance peaks/ nodes are demonstrated for prototypical molecular junctions, illustrating the intricate interplay of quantum effects beyond the single-orbital paradigm.

DOI: 10.1038/ncomms15210 OPEN

Kondo blockade due to quantum interference in single-molecule junctions

Andrew K. Mitchell1,2, Kim G.L. Pedersen3, Per Hedegrd4 & Jens Paaske5

1 School of Physics, University College Dublin, Dublin 4, Ireland. 2 Institute for Theoretical Physics, Utrecht University, Princetonplein 5, Utrecht 3584 CE, The Netherlands. 3 Institut fr Theorie der Statistischen Physik, RWTH Aachen University, Aachen 52074, Germany. 4 Niels Bohr Institute, University of Copenhagen, Copenhagen DK-2100, Denmark. 5 Center for Quantum Devices, Niels Bohr Institute, University of Copenhagen, Universitetsparken 5, Copenhagen DK-2100, Denmark. Correspondence and requests for materials should be addressed to A.K.M. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210

Perhaps the most important feature of nanoscale devices built from single molecules is the potential to exploit exotic quantum mechanical effects that have no

classical analogue. A prominent example is quantum interference (QI), which has already been demonstrated in a number of different molecular devices19. QI manifests as strong variations in the conductance with changes in molecular conformation, contacting or conjugation pathways or simply by tuning the back-gate voltage in a three-terminal setup. Another famous quantum phenomenon, relevant to single-molecule junctions with a spin-degenerate ground state, is the Kondo effect1016, which gives rise to a dramatic conductance enhancement below a characteristic Kondo temperature, TK. Strong electronic interactions in the molecule cause it to bind strongly to a large Kondo cloud17,18 of conduction electrons when contacted to source and drain leads. A hallmark of the Kondo effect is the proliferation of spin ips as electrons tunnel coherently through the molecule, ultimately screening its spin by formation of a many-body singlet19. In this article, we uncover the intricate interplay of these two quantum effects, nding that the combined effect of QI and Kondo physics has highly non-trivial consequences for conductance through single-molecule junctions, and can even lead to an entirely new phenomenonthe Kondo blockade.

The Kondo effect is also routinely observed in semiconductor and nanotube quantum dot devices2024, which are regarded as lead-coupled articial atoms25 and as such are often well described in terms of a single active interacting quantum orbital, tunnel-coupled to a single channel of conduction electrons comprising both source and drain leads. This Anderson impurity model (AIM) is by now rather well understood19, and a quantitative description of the Kondo peak in single quantum dots can be achieved within linear response using non-perturbative methods such as the numerical renormalization group (NRG)26,27. In particular, the conductance is a universal function of T/TK, meaning that data for different systems collapse to the same curve when rescaled in terms of their respective Kondo temperatures2023. Indeed, even for multi-orbital molecular junctions, experimental conductance lineshapes have in some cases been successfully t to the theoretical form of the AIM, suggesting that an effective single-orbital, single-channel description is valid at low temperatures15,16. However, some molecular junctions14,28,29 apparently manifest nonuniversal behaviour and unconventional gate voltage dependences of conductance and TK, hinting at new physics beyond the standard single-orbital paradigm.

The breakdown of the AIM is well-known in the context of coupled quantum dot devices3044, which can be viewed as simple articial molecular junctions due to their multi-orbital structure and the coupling to distinct source and drain channels. Already the extension to two or three orbital systems has lead to the discovery of striking phenomena such as the ferromagnetic Kondo effect3133 corresponding to a sign change of the exchange coupling, and multistage35,36 or frustrated3744 screening.

In the following, we argue that a similar kind of multi-channel, multi-orbital Kondo physics accounts for the behaviour of real molecular junctions, and can be understood as a many-body QI effect characteristic of the orbital complexity and strong electronic correlations in molecules. On entirely general grounds, we construct an effective model describing off-resonant conductance through single-molecule junctions with a spin-degenerate ground state, taking into account both interactions leading to Kondo physics, and orbital structure leading to QI. The physics of this generalized two-channel Kondo (2CK) model (including both potential scattering and exchange-cotunneling) is discussed in relation to the local density of states and observable conductance. We demonstrate how renormalized Kondo resonant conductance evolves into a novel Kondo blockade regime of suppressed

conductance due to QI (Fig. 1). As an illustration, we consider two simple relevant molecular examples, whose properties can be tuned between these limits using gate voltages to provide functionality as an efcient QI-effect transistor.

ResultsModels and mappings. The Hamiltonian describing single-molecule junctions can be decomposed as,

HHmol Hg H

leads

Hhyb: 1

Here Hmol describes the isolated molecule, and contains all information about its electronic structure and chemistry. The rst-principles characterization of molecules is itself a formidable problem when electronelectron interactions are taken into account. In practice however, the relevant molecular degrees of freedom associated with electronic transport are often effectively decoupled. This is the case for many conjugated organic molecules, where the extended p system can be treated separately in terms of an extended Hubbard model45. Reduced multi-orbital models have also been formulated using ab initio methods4649.

The leads are modelled as non-interacting conduction electrons with Hleads

Pask Ekcwaskcask where cwask creates an electron in lead a s, d (source, drain) with momentum

(or other orbital quantum number) k, and spin s m, k. The

dispersion Ek corresponds approximately to a at density of states rEr0y D E

j j

, inside a band of width 2D.

The molecule is coupled to the leads via Hhyb

Pas

tadwiascas H:c:, where cast 1a

Pk takcask is the localized orbital in lead a at the junction, and dias is a specic frontier orbital ia of the molecule, determined by the contacting geometry.

The molecule-lead hybridization is local, and specied by Gapr0 ta

j j2.

The number of electrons on the molecule, N h

Pis dwisdisi, is controlled by a gate voltage, incorporated in the model by

Hg eVg

Pis dwisdis which shifts the energy of all molecular orbitals. Deep inside the Coulomb diamond10, a substantial charging energy, EC, must be overcome to either add or remove electrons from the molecule. Provided Gs;d EC the

Hamiltonian (1) can therefore be projected onto the subspace with a xed number of electrons on the molecule. In general this requires full diagonalization of the isolated Hmol in the

Source

a

Kondo resonance

Kondo blockade

Drain

Drain Source

b

Figure 1 | Interplay between quantum interference and electronic interactions in single molecule junctions. (a) Enhanced Kondo resonant conductance; (b) Kondo blockade, where conductance precisely vanishes. Tuning between a,b by applying a back-gate voltage allows efcient manipulation of the tunnelling current.

2 NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210 ARTICLE

many-particle basis. Charging a neutral/spinless molecule by applying a back-gate voltage to add or remove an electron typically yields a net spin-12 state. For odd-integer N , we therefore

assume that the molecule hosts a spin-12 degree of freedom, S. Higher-spin molecules can arise, but are not considered here (the generalization is straightforward). To second-order in the molecule-lead coupling Hhyb, we then obtain50 an effective model of generalized 2CK form, H2CK H

leads

Hex, where

Hex

Xaa0ss012 Jaa0 S sss 0

Waa0 dss0

cy

asca0s0 : 2

Here s denotes the vector of Pauli matrices. The form of Hex is

guaranteed by spin-rotation invariance and Hermiticity. Further details of the 2CK mapping are provided in Methods. The cotunneling amplitudes form matrices in source-drain space,

J

Jss Jsd

Jsd Jdd

; 3

and are referred to as respectively exchange, and potential scattering terms. These 2CK parameters depend on the specics of molecular structure and contacting geometry in a complicated way, and must be derived from rst-principles calculations for the isolated molecule. This generalized 2CK model hosts a rich range of physics; the non-Fermi liquid critical point38 is merely a single point in its parameter space. Furthermore, any conducting molecular junction must have a Fermi liquid ground state, as demonstrated below.

Any off-resonant molecule hosting a net spin-12 is described by the above generalized 2CK model at low temperatures T EC.

The physics is robust due to the large charging energy deep in a Coulomb diamond (charge uctuations only dominate at the very edge of the Coulomb diamond29). In fact, the physics of the effective model can be regarded as exact in the renormalization group (RG) sense, despite the perturbative derivation of equation (2). Corrections to Hex obtained in higher-order perturbation theory are formally RG irrelevant, and can be safely neglected because they get smaller and asymptotically vanish on decreasing the temperature. They cannot affect the underlying physics; only the emergent energy scales can be modied (this effect is also small, since the corrections are suppressed by EC).

Experimentally relevant physical observables such as conductance can therefore be accurately extracted from the solution of the effective 2CK model (Methods). This requires sophisticated many-body techniques such as NRG, which theoretically attach the source and drain leads non-perturbatively26,27. All microscopic details of a real molecular junction are encoded in the 2CK parameters J and W, which serve as input for the NRG calculations. In particular, destructive QI produces nodes (zeros) in these parameters. Furthermore, QI nodes can be simply accessed by tuning the back-gate voltage Vg, as was

shown recently in ref. 45 for the case of conjugated organic molecules.

Importantly, two different types of QI can arise in molecular junctions due to the electronic interactions. The QI can either be of standard potential scattering type (zeros in elements of W) or of exchange type (zeros in elements of J). Potential scattering QI is analogous to that observed in non-interacting systems described by molecular orbitals. For interacting systems such as molecular junctions (which typically have large charging energies10), potential scattering QI can similarly be understood in terms of extended FeynmanDyson orbitals, which are the generalization of molecular orbitals in the many-particle basis. Information on the real-space character of these orbitals, and how QI relates to molecular structure, can be extracted from the 2CK mapping. By contrast, exchange QI has no single-particle

analogue, and cannot arise in non-interacting systems. Indeed, interactions are a basic requirement for the molecule to host a spin-12 via Coulomb blockade. The spin wavefunction is again characterized by the FeynmanDyson orbitals; depending on the molecule in question, the spin can be delocalized over the entire molecule.

In the following we uncover the effect of this QI on Kondo physics, highlighting two distinct scenarios for the resulting conductanceKondo resonance and Kondo blockade. We then go on to show that this physics is indeed realized in simple examples of molecular junctions, and can be manipulated with gate voltages.

Emergent decoupling. The generalized 2CK model can be simplied by diagonalizing the exchange term in equation (2) via the unitary transformation ca0sUa

0acas such that

UyJU

W

Wss Wsd

Wsd Wdd

Je 0

0 Jo

; Je=oJ d 4

where J 12 Jss Jdd

and d2J2 J2sd. Note that W is not

generally diagonalized by this transformation. The odd channel decouples (Jo 0) if and only if J2sdJssJdd, as is the case when

starting with a single-orbital Anderson model (see Supplementary Note 1). By contrast, real multi-orbital molecules couple to both even and odd channels (electronic propagation through the entire molecule yields J2sd JssJdd when off resonance).

However, electronic interactions play a key role here: the exchange couplings become renormalized as the temperature is reduced. A simple perturbative RG treatment hints at ow toward a two-channel strong-coupling state, since both Je and Jo initially grow. But the true low-temperature physics is much more complex, as seen in Fig. 2 from the imaginary part of the scattering T-matrix taa(o,T) pr0ImTaa(o,T) obtained by

NRG for the generalized 2CK model and plotted as a function of excitation energy o at T 0 (see Methods). The molecule spin is

ultimately always Kondo-screened by conduction electrons in the more strongly coupled even channel since Je4Jo for any nite d.

Indeed, any real molecular junction will inevitably have some degree of asymmetry in the source/drain coupling J , so that dZJ is always nite in practice. At particlehole

(ph) symmetry, the Friedel sum rule19 then guarantees that tee(0, 0) 1, characteristic of the Kondo effect. On the other

hand, Kondo correlations with the less strongly coupled odd channel are cut off on the scale of TK, and therefore too(0, 0) 0

(consistent with the optical theorem). These analytic predictions are veried by NRG results in the centre panels of Fig. 2.

In all cases the odd channel decouples on the lowest energy/ temperature scales, and the problem becomes effectively single-channel. This is an emergent phenomenon driven by interactions, not a property of the bare model. Despite the emergent decoupling of the odd channel, the Kondo effect always involves conduction electrons in both source and drain leads for any nite Jsd. From the transformation dened in equation (4), the T-matrix in the physical basis can be expressed as taao; TjUa;ej2teeo; T jUa;oj2tooo; T, such

that taa0; 0jUa;ej2 at ph symmetrysee left panels of Fig. 2.

Although the physics at T 0 is effectively single-channel, the

full temperature dependence is highly non-trivial due to the competing involvement of the odd channel (only for the oversimplied single-orbital AIM is the odd channel strictly decoupled for all T). The universal physics of the AIM is lost for d 12 Je Jo

6 J (or equivalently J2sd 6 JssJdd): conductance

lineshapes no longer exhibit scaling collapse in terms of T/TK. Indeed, Kondo screening by the even channel occurs on the scale TeK D exp 1=r0Je

, and hence depends on d. The

NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210

a

Source/drain basis

g

Even/odd basis

m

Decoupled

"Kondo cloud" Fermi liquid

1

0.5

0

1

0.5

0

1

0.5

0

Jsd /D 0.2 (1CK)

Jsd /D = 101

Jsd /D = 102

Jsd /D = 103

J /D = 104

b

c

d

e

f

h

i

j

k

l

n

o

p

q

r

t ( , 0)

1

0.5

0

1

0.5

0

1

0.5

0

1010 108 106

/D /D

Increasing even-odd frustration

Odd chain

Even chain

104 102 100 1010 108 106 104 102 100

Figure 2 | Decoupling and frustration due to the Kondo effect in molecular junctions. (af) Imaginary part of the T-matrix, characterizing the effective energy-dependent exchange, in the physical basis of source (blue) and drain (red) leads at T 0 for the effective 2CK model. taa(o, 0) is related to the

renormalized density of states in lead a at the junction. We take a representative molecule-lead coupling J 0.2D with small but nite source/drain

coupling asymmetry J

10 4D, and consider the effect of reducing the exchange-cotunneling Jsd from af. Physically, this could be achieved by

gate-tuning in the vicinity of a QI node. The frustration of Kondo screening is always eventually relieved on the lowest energy scales, below TFL minTeK; T , because dZJ is always nite in any realistic setting. (gl) Corresponding T-matrix in the even/odd (blue/red) lead basis.

(mr) Real-space competition between even/odd (blue/red) conduction electron channels, illustrated for the case where the leads are 1D quantum wires. The Kondo cloud (yellow) corresponds to the spatial region of high molecule-lead entanglement. For small dt10 3D one has a Kondo frustration cloud

embodying incipient overscreening17.

Kondo temperature itself can therefore acquire an unconventional gate voltage dependence, beyond the AIM paradigm.

For even smaller d, the Kondo effect occurs as a two-step process, with even and odd channels competing to screen the molecule spin. As in this case JeEJo, a frustration of Kondo screening sets in on the scale TeK ToK T2CKK. The incipient

frustration for T T2CKK results in only partial screening

(the molecule is overscreened, producing non-Fermi liquid signatures38,4042,51). The frustration is relieved on the much smaller scale38 T D r0d

2. The even channel eventually wins

for T T and fully Kondo-screens the molecule spin, while

the odd channel decouples. This dramatic breakdown of the single-orbital AIM paradigm is shown in Fig. 2, with the degree of even/odd frustration increasing from top to bottom. In practice, such frustration arises in a nearly symmetrical junction (small J ), tuning in the vicinity of a QI node in Jsd such that

the perturbation strength d is reduced. The rst signatures of frustration appear in conductance when T*tTeK. Only when J Jsd 0, such that d 0, does the frustration persist

down to T 0; we do not consider this unrealistic scenario in the

present work.

In real-space, the entanglement between the molecule and the leads is characterized by the Kondo cloud18a large spatial

region of extent xK vF=kBTeK penetrating both source and

drain leads (vF is the Fermi velocity). In the right panels of Fig. 2 we illustrate this for the case where the leads are 1D quantum wires; the real-space physics is then directly related to the T-matrix plotted in the left panels, as shown in ref. 17. Note that if the source/drain leads are 1D quantum wires, then the even/odd leads are also 1D quantum wires as depicted. For small d (lower panels) we have instead a Kondo frustration cloud. The frustration is only relieved at longer length scales x vF=kBT ,

beyond which the odd channel decouples.

Conductance. The current through a molecular junction is mediated by the cross terms coupling source and drain leads; the exchange and potential scattering terms Jsd and Wsd constitute two distinct conductance mechanisms. At high temperatures, the overall conductance can be understood from a simple leading-order perturbative treatment using Fermis golden rule and is simply additive45, G=G0 2pr0

2 W2sd 3J2sd

,

with

G0 2e2h 1. However, at lower temperatures, electronic

interactions lead to strong renormalization effects and rather surprising Kondo physics. Non-perturbative methods such as NRG must therefore be used to calculate the full temperature-

4 NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210 ARTICLE

a c

b

1

0.8

0.6

G (0) / G 0

T/T*

1

0.8

0.6

0.4

0.2

0 102 101 100 101 102 10 5 0 5 10

103

G (T) / G 0

G (T,V sd) / G 0

1

0.8

0.6

0.4

0.2

0

0.4

0.2

0 6 4 2

(J ss J dd ) / 2J sd T / TFL

0 2 4 6

eVsd / kB T*

Figure 3 | Kondo resonant conductance near a potential scattering quantum interference node. (a) Zero-temperature linear conductance G(0)as a function of derived 2CK parameters Jss, Jdd and Jsd at Wss Wdd Wsd 0. (b) Limiting universal conductance curves G(T/TK) and G(T/T*) in the

single-channel regime (large d, red line) and the frustrated two-channel regime (small d, black line and points), respectively. (c) Exact non-equilibrium conductance G(T, Vsd) as a function of bias voltage Vsd at various temperatures in the frustrated regime of small d, from equation (7).

dependence of conductance, as described in Methods. Note that conductance through a molecular junction cannot be obtained simply from the T-matrix (except at T 0).

The QI aspect of the problem is entirely encoded in the effective 2CK parameters, providing an enormous conceptual simplication. In particular, we identify two limiting QI scenarios relevant for conductance: Wsd 0 or Jsd 0. Exact analytic

results, supported by NRG, show that the Kondo effect survives QI in the case of Wsd 0 to give enhanced conductance at low

temperatures (Figs 1a and 3), while a Kondo-mediated QI node in the total conductance is found for Jsd 0, a Kondo blockade

(Figs 1b and 4). We demonstrate explicitly that this remarkable interplay between QI and the Kondo effect arises in two simple conjugated organic molecules on tuning gate voltages in Fig. 5.

Kondo resonance. First we focus on conductance mediated exclusively by the exchange cotunneling term Jsd, tuning to a potential scattering QI node W 0. Even though the bare Jsd is

typically small, it gets renormalized by the Kondo effect and becomes large at low temperatures. The Kondo effect therefore involves both source and drain leads (Fig. 2), leading to Kondo-enhanced conductance.

As shown in Supplementary Note 2, the fact that the odd channel decouples asymptotically implies the following exact result for the linear conductance,

GT04G0

tss0; 0tdd0; 0

T TFL ln 2 T=TeK

: T TeK;

T=T

1 : T TeK:

6

Furthermore, the abelian bosonization methods of refs 5355 can be applied to single-molecule junctions in the limit T TeK to

obtain an exact analytic expression for the full conductance crossover (Supplementary Note 3),

G T; Vsd =G0

T 2pT Re c1

1

2

T 2pT i

eVsd

2pkBT

; 7

where c1 is the trigamma function. Remarkably, this result also holds away from thermal equilibrium, at nite bias Vsd TeK.

Within linear response, equation (7) is conrmed explicitly by comparison to NRG data in Fig. 3b, while Fig. 3c shows the nonequilibrium predictions. The condition T TeK pertains to

nearly symmetric junctions, tuned near a QI node in Jsd.

Equation (7) should be regarded as a limiting scenario: conductance lineshapes for real single-molecule junctions will typically interpolate between the red and black lines of Fig. 3b.

Kondo blockade. At a QI node in the exchange-cotunneling Jsd 0, conductance through a single-molecule junction is

mediated solely by Wsd. In this case, the molecule spin is fully Kondo screened by either the source or drain lead (whichever is more strongly coupled). Only in the special but unrealistic case Jss Jdd and Jsd 0 does the frustration persist down to T 0. For

concreteness we now assume ph symmetry Wss Wdd 0, and

Jss4Jdd such that the even conduction electron channel is simply the source lead. The drain lead therefore decouples on the scale of TsK. As shown in Supplementary Note 4, one can then prove that,

GT0G0 2pr0Wsd

2 1 tss0; 0

0; 8 where tss(o, T) is the T-matrix of the source lead. The Kondo effect with the source lead, characterized by tss(0, 0) 1, therefore

exactly blocks current owing from source to drain. This is an emergent effect of interactionsat high temperatures T TK

p G0 4J2sd4J2sd Jss Jdd

2

: 5

Note that any nite interlead coupling Jsd yields unitarity conductance G G0 at T 0 in the symmetric case Jss Jdd.

The analytic result is conrmed by NRG in Fig. 3a, and further holds for all T TK; T . Equation (5) is an exact

generalization of the standard single-orbital AIM result, G0=G04JssJdd= Jss Jdd

2 4GsGd= Gs Gd

2, and reduces

to it when J2sdJssJdd.

The full temperature-dependence of conductance can also be studied with NRG. In all cases, we nd Fermi-liquid behaviour GT G0 T=TFL

2 at the lowest temperatures T TFL,

with TFL min TeK; T

(although TFL itself may have a nontrivial gate dependence). At large d J T TeK

,

the behaviour of the single-channel AIM20,52 is essentially recovered for the entire crossover (see red line, Fig. 3b). However, the universality of the AIM is lost for smaller d due

to the competing involvement of the odd screening channel. In fact, for T TeK, appreciable conductance only sets in

around T T (rather than TeK), and the entire conductance

crossover becomes a universal function of T/T*different in form from that of the AIM (see black line, Fig. 3b). The formation of the Kondo state is reected in conductance by the following limiting behaviour,

GT

NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210

0.20

Wsd /D = 0.1

Kondo

0.15

0.10

0.05

0.00

102 101 101

T / TK

G (T)/ G 0

102 103

100

Figure 4 | Kondo blockaded conductance near a quantum interference node in the exchange cotunneling. NRG results for the conductance G(T) as a function of rescaled temperature T/TK at Jsd 0 for various Wsd,

showing in all cases an overall conductance node G(0) 0. Plotted

for Jss 0.25D, Jdd 0.2D and Wss Wdd 0. Dotted lines show the

high-temperature perturbative expectation.

when tssE0, conductance is nite and the perturbative result is recovered, Gpert/G0E(2pr0Wsd)2.

The zero-temperature conductance node arising for Jsd 0 can

be understood physically, as a depletion of the local source-lead density of states at the junction, due to the Kondo effect. Conductance vanishes because locally, no source-lead states are available from which electrons can tunnel into the drain lead. From a real-space perspective17, one can think of the Kondo cloud in the source lead as being impenetrable to electronic tunnelling at low energies. The effect of this Kondo blockade is demonstrated in Fig. 4, where full NRG calculations for the conductance are shown for Jsd 0. The conductance crossover as

a function of temperature is entirely characteristic of the Kondo effect; G(T)/Gpert is a universal function of T/TK.

At low temperatures, a node in Jsd thus implies an overall conductance node, even though Wsd remains nite.

The Kondo blockade will be most cleanly observed in real single-molecule junctions that have strong molecule-lead hybridization and do not have a nearby Kondo resonance. In addition to the large Kondo temperature, the perturbative cotunneling conductance observed at high temperatures T TK

is also larger in this case, thereby increasing the contrast of the blockade on lowering the temperature.

We emphasize that the Kondo blockade is unrelated to the Fano effect34,56,57, which arises due to QI in the hybridization rather than intrinsic QI in the interacting molecule itself (Supplementary Note 4). Unlike the Kondo blockade, the Fano effect is essentially a single-channel phenomenon that does not necessitate interactions, and different (asymmetric) lineshapes result.

Gate-tunable QI in Kondo-active molecules. In real molecular junctions, the two conductance mechanisms discussed separately above (due to nite exchange Jsd and potential scattering Wsd), are typically both operative. Their mutual effect can be complicated due to renormalization from cross-terms proportional to Wsd Jsd. However, as the gate voltage Vg is tuned, both Kondo resonant and Kondo blockade regimes are often accessible due to QI nodes45 in either Jsd or Wsd. In practice, we observe that

overall conductance nodes can also be shifted away from the nodes in Jsd by marginal potential scattering Wss and Wdd (not considered above). We speculate that the conductance nodes are topological and cannot be removed by potential scattering only shifted to a different gate voltage. Precisely at the node, the low-temperature physics is universal and therefore common to all such off-resonant spin-12 molecules.

To demonstrate the gate-tunable interplay between QI and the Kondo effect in single-molecule junctions, we now consider two simple conjugated organic molecules as examples. Following ref. 45, exact diagonalization of the PariserParrPople (PPP) model58 for the sp2-hybridized p system of the molecule allows the effective 2CK model parameters to be extracted as a function of applied gate voltage (Supplementary Note 5). The 2CK model is then solved using NRG26,27, and the conductance is calculated numerically-exactly as a function of temperature. These steps are described in detail in Methods.

Figure 5 shows the conductance G(T) for junctions spanned by respectively a benzyl, (a) and an isoprene-like molecule (d), as a function of rescaled temperature T/TK at different gate voltages.

Both systems exhibit Kondo resonant and Kondo blockade physics. In panel (a), a pronounced Kondo blockade appears near Vg 0, corresponding to the midpoint of the Coulomb diamond.

Finite conductance at higher temperatures due to cotunneling Wsd is blocked at low temperatures by the Kondo effect. On increasing the gate voltage, we nd numerically that G0 eV2g,

with conductance enhancement due to renormalized Jsd (Fig. 5b). The overall conductance in this case remains rather small for all eVg analysed. We also note that the Kondo temperature varies as ln TK=D eV2g, Fig. 5c. This gate evolution of TK could be

considered as conventional from the single-orbital AIM perspective10, but the conductance itself is blockaded rather than enhanced by Kondo correlations.

However, richer physics can be accessed in junction (d). The crossovers of G(T) show perfect Kondo resonant conductance at nite eVg 2.4 eV, reaching the unitarity limit G(0) 2e2h 1.

But increasing the gate voltage slightly to Vg 2.625 eV yields

almost perfect Kondo blockade, with G(0)C0 (note the log scale). The full crossovers are entirely characteristic of the underlying correlated electron physics. Panel (e) shows the evolution of G(0) as a function of gate voltage at T 0 (and in practice for all

T TK), which exhibits nontrivial behaviour due to the

interplay between QI and the Kondo effect. The rapid switching between Kondo resonant and Kondo blockade conductance with applied gate voltage might make such systems candidates for QI-effect transistors, or other technological applications.

Finally, in panel (f), we show that the Kondo temperature also displays an unconventional gate-dependence, with TK increasing as one moves in towards eVg 0, analogous to the

effect observed experimentally in ref. 28. The Kondo temperature remains nite for all eVg, but takes its minimum value at the

Kondo resonance peak. In practice, the Kondo temperature can vary widely from system to system because it depends sensitively (exponentially) on the molecule-lead hybridization. However, Kondo temperatures up to around 30 K are commonly observed in real single-molecule junctions10.

We did not attempt an ab initio calculation of the absolute Kondo temperatures, but note that the effective bandwidth cutoff D in the effective 2CK model is essentially set by the large charging energy of the molecule. For the PPP models used for the conjugated hydrocarbons in Fig. 5, this in turn is set by the onsite Coulomb repulsion, taken to be 11 eV within the standard Ohno parametrization58. With this identication, we have TK 10 K for

the specic example shown in panel (a) at the Kondo blockade, and 0.1 K in (d). We emphasize that the Kondo blockade arises on similar temperature scales to that of the standard Kondo effect in

6 NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210 ARTICLE

a d

0.00075 eVg [eV]

eVg [eV]

1.0

1

0.1

G (T ) [2e2 h1 ]

0.50.30.2 0

0.00050

0.00025

0.00000

0.003

0.002

0.001

02 1 1 2

0 2 1 1 2

1

0.8

0.6

0.4

0.2

G (T ) [2e2 h1 ]

2.42.452.52.552.5752.62.625

0.01

0.001

0.0001

101 100 100

102 102 104

101T / TK

T / TK

T K /D

102 102

b

c

e f

102

102

G (0) [2e2 h1 ]

104

T K /D

103

G (0) [2e2 h1 ]

106

108

104

0 2

0 3 1 1 2 3

0 2

3 1 1 2 3

0

eVg [eV] eVg [eV] eVg [eV] eVg [eV]

Figure 5 | Gate-tunable Kondo resonance and Kondo blockade in simple conjugated organic molecular junctions. Single-molecule junctions basedon a benzyl (ac) and an isoprene-like molecule (df), with all carbons sp2-hybridized, were mapped to an effective 2CK model and linear conductance was calculated with NRG. (a,d) Conductance G(T) as a function of rescaled temperature T/TK for various gate voltages eVg. (a) shows Kondo blockade

G(0) 0 at eVg 0 and Kondo enhanced conductance for Vg

40.

(d) shows Kondo blockade at eVg 2.625 eV and perfect (unitarity) Kondo resonance

at eVg 2.4 eV. (b,e) G(0) as a function of gate voltage eVg at T 0; (c,f) Corresponding Kondo temperatures. Note the sensitive gate dependence of

G(0) in (e), and the corresponding unconventional non-monotonic gate dependence of the Kondo temperature in (f).

molecules, and therefore signatures should generally be observable at experimental base temperatures on gate tuning to a QI node.

The stark difference in transport properties of the two molecular junctions shown in Fig. 5 is due to differences in their QI characteristicsspecically the number and position of QI nodes in the effective 2CK parameters (Supplementary Note 5). In turn, this is related to the underlying molecular structure and contacting geometry, as explored for these alternant hydrocarbons in ref. 45. Both molecules exhibit a Kondo blockade due to a node in Jsd, but this arises at eVg 0 for the

benzyl radical in (ac), whereas there are two nodes at nite eVg for the isoprene-like molecule in (df). In general, Jsd has an odd(even) number of nodes as a function of gate in odd-membered alternant molecules if the source and drain electrodes are connected to sites of the molecule on different sublattices(the same sublattice) of the bipartite p system.

The strong Kondo resonance arising at eVg 2.4 eV for the

isoprene-like molecule is a consequence of the parity symmetry with respect to the contacting geometry, such that JssEJdd

(equation (5)). By contrast, there is no such symmetry for the benzyl molecule and Jdd happens to dominate.

Although we have exemplied the gate-tunable interplay between QI and Kondo effect with these conjugated hydrocarbon moieties, we emphasize that a Kondo blockade should be found in any off-resonant spin-12 molecule with intramolecular interactions and sufcient orbital complexity to produce a QI node in the exchange cotunneling.

DiscussionTransport through spinful Coulomb-blockaded single-molecule junctions requires a description beyond the standard single-orbital Anderson paradigm. The relevant model is instead a generalized

two-channel Kondo model, to which real molecular junctions can be exactly mapped. Experimental data for individual molecular junctions can be understood within this framework, avoiding the need for a statistical interpretation.

Quantum interference can be classied as being of either exchange or potential scattering type. Although these distinct conductance mechanisms are simply additive at high energies, where standard perturbation theory holds, the low-temperature behaviour is much richer due to electronelectron interactions which drive the Kondo effect. We show that the Kondo effect survives a quantum interference node in the potential scattering to give enhanced conductance, while a novel Kondo blockade arises in the case of an exchange cotunneling node, entirely blocking the current through the junction. This rich physics is tunable by applying a back-gate voltage, as demonstrated explicitly for two simple conjugated organic molecules, opening up the possibility of efcient Kondo-mediated quantum interference effect transistors.

The theoretical framework we present can be used to systematically study candidate molecules and help optimize the type and location of anchor groups for particular applications. Quantum chemistry techniques could be used to accomplish the Kondo model mapping for larger molecules. The effect of vibrations and dissipation (relevant at higher energies5) could also be taken into account within generalized AndersonHolstein or BoseFermi Kondo models27.

Methods

Schrieffer-Wolff transformation. We derive the effective Kondo model describing off-resonant single-molecule junctions by projecting out high-energy molecular charge uctuations from the full lead-coupled system. This is equivalent to a two-channel generalization of the standard SchriefferWolff transformation19. That is, projecting onto the subspace of Hilbert space where the number, N, of electrons on the molecule is xed. The effective Hamiltonian in this subspace has

NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210

such that in practice only two matrix elements are needed to obtain the exchange couplings Jaa and the potential scattering amplitudes Waa .

Equations (13) and (14) therefore encode all the properties of the single-molecule junction inside an N-electron Coulomb diamond. The exchange and potential scattering terms in the 2CK model are determined by the amplitudes A which, from equation (12), have contributions from both particle (p) and hole (h) processes (that is, processes involving virtual states with N 1 or N 1 electrons

on the molecule). In particular, note that all quantum interference effects are entirely encoded in the effective parameters Jaa and Waa quantum interference nodes arise if and only if molecular states are connected by particle and hole processes with equal but opposite amplitudes. As discussed in ref. 45, the appearance of such nodes can be understood in terms of the properties of the underlying FeynmanDyson orbitals.

We emphasize that the effective 2CK model is totally general, applying for any molecule with a two-fold spin-degenerate ground state, at temperatures less than the molecule charging energy so that charge uctuations on the molecule are frozen (typically the charging energy is large when deep inside a Coulomb diamond, and therefore the molecule is off-resonant). The 2CK model parameters can be obtained purely from a knowledge of the isolated molecule, and can therefore be calculated in practice using a number of established techniques (exact diagonalization, conguration interaction and so on). The low-temperature properties of the resulting 2CK model are however deeply nontrivial, requiring sophisticated many-body methods to attach the leads and account for nonperturbative renormalization effects. In the present work, we do this second step using the numerical renormalization group27.

An advantage of the effective theory is that it can also be analysed exactly on an abstract level (independently of any specic realization). This allows us to identify all the possible scenarios that could in principle arise in molecular junctions. The basic physics is arguably obfuscated rather than claried by the complexity of a full microscopic description: a brute-force method (even if that were possible) may not yield new conceptual understanding or provide general predictions beyond a case-by-case basis.

2CK parameters for the molecules presented in Fig. 5 were obtained following ref. 45; see Supplementary Figs 1 and 2.

Exact diagonalization of Hmol. In this work, we model the isolated molecule by a semi-empirical Pariser-Parr-Pople Hamiltonian59,60 for the molecular p-system:

^

Hmol

Xi;j h i

Xs"=#tijdyisdjs H:c:

X

i

the form

Heff EP H

leads

Hmol Hg HhybQ E QHQ 1QHhyb

P;

9

where P is a projection operator onto the N-electron subspace of the molecule, while Q I P projects onto the orthogonal complement. These subspaces are

connected by Hhyb, and the resolvent operator (E QHQ) 1 determines the

propagation of excited states at energy E. Note that PHleads P H

leads since P acts only on the molecular degrees of freedom, and P(Hmol Hg)P is merely a constant

and dropped in the following. For the isolated molecule, HmoljCNniENnjCNni,

where jCNni denotes the nth N-electron many-body eigenstate with energy

ENn, and EN0 is the ground state energy. jCNni spans the entire molecule and

generally has weight on all atomic/molecular basis orbitals in Hmol. So far the

treatment is exact.To second order in Hhyb, Equation (9) reduces to the effective Hamiltonian

Heff H

leads

PHhybQ EN0 Q Hmol Hg Q

1QH

hybP: 10 Virtual processes involving a given excited state jCN 1ni contribute to

Equation (10) with weight controlled by the energy denominator

hCN 1njEN0 QHmol HgQjCN 1niEN0 EN 1n eVg, which must be negative

to ensure stability of the N-electron ground state (including the electrostatic shift from the backgate described by Hg). The perturbative expansion in Hhyb is controlled by a large energy denominator, and therefore use of Equation (10) is justied deep inside the N-electron Coulomb diamond. Inserting the tunnelling Hamiltonian, Hhyb

Pask

takdwiscask H:c:, one arrives at the effective

Hamiltonian Heff H

leads

Hex, with

Hex

Xa k saksta k t akcyaksca k s

Xm ;mCNm

Aa ass;m m CNm

11

where the last line follows from the denition of local lead-electron operators cast 1a

Pk takcaks, with t2a

Pk

Xa sasta t acy asca s

Xm ;mCNm

Aa

as s;m m CNm

;

jtakj2. Here m and m0 label (degenerate)

molecular ground states with energy EN0. For odd N, the molecule often carriesa net spin-12, and so m and m0 are simply the projections Sz 12. The spin density

need not be spatially localized. We now focus on this standard case, although the generalization to arbitrary spin is straightforward when the molecular ground state for a given N is more than two-fold degenerate.

The cotunneling amplitudes can be decomposed as,

Aa ass;m mha ass;m m pa ass;m m 12

where the contributions from hole and particle propagation are given respectively by,

ha ass;m mVg

XnCNm

dy

U ni" 12

ni# 12

:

21

1

2

Xi 6 jVij ni 1 nj 1

The operator dwis creates an electron with spin s on the pz-orbital i

j i, nisdwisdis

and nini" ni#. The Coulomb interaction is given by the Ohno parametrization58 VijU=

1 ~rij

i s CN 1n

CN 1n

d

i s CNm

q 2U2=207:3 eV,

where ~rij

is the real-space distance between

two pz-orbitals i

j i and j

j i measured in ngstrm. For sp2 hybridized carbon, the

nearest neighbour overlap, tij, is tE 2.4 eV, and UE11.26 eV (ref. 61).

For suitably small molecules, equation (21) can be solved using exact diagonalization (exploiting overall conserved charge and spin) to provide the many-particle eigenstates jCNni and eigenenergies ENn. For larger molecules,

approximate methods can also be used, provided interactions are accounted for on some level. Any molecule can be addressed within our framework, provided the eigenstates and eigenenergies of the isolated molecule can be determined.

Importantly, the perturbed two-channel Kondo model derived in the previous section remains the generic Hamiltonian of interest to describe such off-resonant junctions. The diagonalization of equation (21) is required only to obtain the parameters J and W, which are then used in subsequent numerical renormalization group calculations to treat the coupling to source and drain leads. Note that the calculation of physical quantities such as conductance at lower temperatures necessitates an explicit and nonperturbative treatment of the leads, and cannot be achieved with single-particle methods or exact diagonalization alone.

However, once the generic physics of the underlying 2CK model is understood(a key goal of this paper), the transport properties and quantum interference effects of specic molecular junctions can already be rationalized and predicted from their 2CK parameters. The suitability of candidate molecules and the positions of anchor groups can therefore be efciently assessed, opening up the possibility of rational device design.

Calculation of conductance. The key experimental quantity of interest for single-molecule junction devices is the differential conductance G T; Vsd d Isd

h i=dVsd.

In this section we recap the generic framework for exact calculations of the linear response conductance G(T) G(T, Vsd-0) through a molecule, taking fully into

account renormalization effects due to electronic interactions. We then describe how the NRG27 can be used to accurately obtain G(T) for a given system described by the effective model, equation (16).

To simulate the experimental protocol, we add a time-dependent bias term to the Hamiltonian, H H2CK H0(t), with

H0t

eVsd

2 cosot

eVg EN0 E

N 1

n

i0 ; 13

pa ass;m mVg

XnCNm

d

i s CN 1n

CN 1n

dy

i s CNm

i0 : 14

The matrix elements in the numerators (referred to as Feynman-Dyson orbitals) constitute a correlated generalization of molecular orbitals, and are computed in the many-particle molecular eigenstate basis.

Since the total Hamiltonian must preserve its original spin-rotational invariance, the cotunneling amplitude must take the form

Aa ass;m m14Jaa sss smm Waa dss dmm : 15

This leads to the desired effective 2CK model (equation (2) of subsection Models and mappings):

H2CKH

leads

eVg EN0 E

N 1

n

Xa sas12Ja aS ss s Wa ads s

cy

a s cas: 16

The 2CK model parameters themselves are obtained from traces with Pauli matrices

Jaa ta t a

Xss ;mmtis sAa ass;m mtim m for ix; y; z; 17

Waa ta t a

Xss ;mmAa ass;mm 18

which, by spin-rotation invariance, further simplify to

Jaa 2ta t aAa a"#;#"2ta t a

XmAa a"";mmtzmm 19

Waa 4ta t aAa a"";"" 2Aa a"#;#"2ta t a

XmAa a"";mm 20

^

Nd

;

^

Ns

22

where ^

Na

Pk;s cwakscaks is the total number operator for lead a. We focus on

8 NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210 ARTICLE

the serial ac and dc conductance at t 0, after the system has reached an oscillating

steady state. In the zero-bias limit Vsd-0, the exact serial ac conductance at linear response follows from the Kubo formula62,

Gaco; T

e2h

2p 2 Im Ko; T

o

; 23

where K(o, T) is the Fourier transform of the retarded currentcurrent correlator,

Kt; T iyt

^

^ O0 h i

D E

T ; 24

where ^

O12

_ _ and _

Naddt

^

Na. The dc conductance is then simply,

GT lim

o!0

Gaco; T: 25

In practice, NRG is used to obtain K(o, T) numerically. The full density matrix NRG method63, established on the complete AndersSchiller basis64, provides essentially exact access to such dynamical correlation functions at any temperature T and energy scale o. We use _

Nai

to nd an expression for the current operator amenable for treatment with NRG:

i^

O JsdS ssd Wsd

Xscy sscds

^

H; ^

Na

!

H:c:; 26

where sab12

Pss cwasrss cbs .

As the ground state of any conducting molecular junction must be a Fermi liquid, the system can be viewed as a renormalized non-interacting system at T 0

(ref. 65). The zero-bias dc conductance at T 0 can therefore also be obtained from

a LandauerBttiker treatment,

GT0

2e2h 4

sd Gsd o0; T0

j j2; 27

in terms of the full retarded electronic Greens function Gabo; T $

FT

iyt

hfcast; cwbs0giT, which must be calculated non-perturbatively in the presence of

the interacting molecule. Herea1= pr0

. Note that equation (27) applies to

Fermi liquid systems only at T 0 and in the dc limit. The full temperature

dependence of G(T) must be obtained from the Kubo formula.

In practice, Gsd o; T

is obtained from the T-matrix equation19, which

describes electronic scattering in the leads due to the molecule,

Gab

o; TG0odab G0oh i2

Wab Tabo; T

;

28

where G0o is the free retarded lead electron Greens function when the

molecule is disconnected, such that Im G0o pro, and r(0) r0.

Within NRG, the T-matrix can be calculated directly51 as the retarded correlator Tabo; T $

FT

iythfaat; awb0giT. For the present problem,

the composite operators,

aa

XgWagcg" 12Jag cg"Sz cg#S

;

29

follow from equation (16) using equations of motion methods. In subsection Emergent decoupling we also present NRG results for the spectrum of the T-matrix, dened as

tabo; T pr0Im Tabo; T: 30

Note that the simple Landauer form of the MeirWingreen formula66, which relates the conductance through an interacting region to a generalized transmission function, applies only in the special case of proportionate couplings. In single-molecule junctions, the various molecular degrees of freedom couple differently to source and drain leads (which are spatially separated), and therefore this standard form of the MeirWingreen formula cannot be used, and one has to resort to using full Keldysh Greens functions (or the methods described above for linear response). The exception is when the molecule is a single orbitalthis articial limit is considered in Supplementary Note 1.

For the NRG calculations, even/odd conduction electron baths were discretized logarithmically using L 2, and Ns 15,000 states were retained at each step of the

iterative procedure. Total charge and spin projection quantum numbers were exploited to block-diagonalize the NRG Hamiltonians, and the results of Nz 2

calculations were averaged. 2CK model parameters for the molecules presented in Fig. 5 are discussed in Supplementary Note 5.

Data availability. The data supporting our ndings are available from the corresponding author on reasonable request.

References

1. Fracasso, D., Valkenier, H., Hummelen, J. C., Solomon, G. C. & Chiechi, R. C. Evidence for quantum interference in sams of arylethynylene thiolates in tunneling junctions with eutectic Ga-In (EGaIn) top-contacts. J. Am. Chem. Soc. 133, 95569563 (2011).

2. Gudon, C. M. et al. Observation of quantum interference in molecular charge transport. Nat. Nanotechnol. 7, 305309 (2012).

3. Aradhya, S. V. et al. Dissecting contact mechanics from quantum interference in singlemolecule junctions of stilbene derivatives. Nano Lett. 12, 16431647 (2012).

4. Vazquez, H. et al. Probing the conductance superposition law in single-molecule circuits with parallel paths. Nat. Nanotechnol. 7, 663667 (2012).

5. Ballmann, S. et al. Experimental evidence for quantum interference and vibrationally induced decoherence in single-molecule junctions. Phys. Rev. Lett. 109, 056801 (2012).

6. Arroyo, C. R. et al. Signatures of quantum interference effects on charge transport through a single benzene ring. Angew. Chem. Int. Ed. 52, 31523155 (2013).

7. Arroyo, C. R. et al. Quantum interference effects at room temperature in OPV-based single-molecule junctions. Nanoscale Res. Lett. 8, 234 (2013).

8. Rabache, V. et al. Direct observation of large quantum interference effect in anthraquinone solid-state junctions. J. Am. Chem. Soc. 135, 1021810221 (2013).

9. Koole, M., Thijssen, J. M., Valkenier, H., Hummelen, J. C. & van der Zant, H. S.J. Electric-eld control of interfering transport pathways in a single-molecule anthraquinone transistor. Nano Lett. 15, 55695573 (2015).10. Scott, G. D. & Natelson, D. Kondo resonances in molecular devices. ACS Nano 4, 35603579 (2010).

11. Park, J. et al. Coulomb blockade and the Kondo effect in single-atom transistors. Nature 417, 722725 (2002).

12. Liang, W., Shores, M. P., Bockrath, M., Long, J. R. & Park, H. Kondo resonance in a single-molecule transistor. Nature 417, 725729 (2002).

13. Yu, L. H. & Natelson, D. The Kondo effect in C60 single-molecule transistors. Nano Lett. 4, 7983 (2004).

14. Osorio, E. A. et al. Electronic excitations of a single molecule contacted in a three-terminal conguration. Nano Lett. 7, 33363342 (2007).

15. Parks, J. J. et al. Tuning the Kondo effect with a mechanically controllable break junction. Phys. Rev. Lett. 99, 026601 (2007).

16. Scott, G. D., Keane, Z. K., Ciszek, J. W., Tour, J. M. & Natelson, D. Universal scaling of nonequilibrium transport in the Kondo regime of single molecule devices. Phys. Rev. B 79, 165413 (2009).

17. Mitchell, A. K., Becker, M. & Bulla, R. Real-space renormalization group ow in quantum impurity systems: local moment formation and the Kondo screening cloud. Phys. Rev. B 84, 115120 (2011).

18. Srensen, E. S., Chang, M.-S., Laorencie, N. & Afeck, I. Impurity entanglement entropy and the Kondo screening cloud. J. Stat. Mech.: Theory Exp. 2007, L01001 (2007).

19. Hewson, A. C. The Kondo problem to heavy fermions (Cambridge University Press, 1997).

20. Goldhaber-Gordon, D. et al. Kondo effect in a single-electron transistor. Nature

391, 156159 (1998).

21. Cronenwett, S. M., Oosterkamp, T. H. & Kouwenhoven, L. P. A tunable Kondo effect in quantum dots. Science 281, 540544 (1998).

22. Van der Wiel, W. G. et al. The Kondo effect in the unitary limit. Science 289, 21052108 (2000).

23. Nygrd, J., Cobden, D. H. & Lindelof, P. E. Kondo physics in carbon nanotubes. Nature 408, 342346 (2000).

24. Pustilnik, M. & Glazman, L. I. Kondo effect in real quantum dots. Phys. Rev. Lett. 87, 216601 (2001).

25. Kouwenhoven, L. & Marcus, C. M. Quantum dots. Phys. World 11, 35 (1998).26. Wilson, K. G. The renormalization group: critical phenomena and the Kondo problem. Rev. Mod. Phys. 47, 773 (1975).

27. Bulla, R., Costi, T. A. & Pruschke, T. Numerical renormalization group method for quantum impurity systems. Rev. Mod. Phys. 80, 395 (2008).

28. Yu, L. H. et al. Kondo resonances and anomalous gate dependence in the electrical conductivity of singlemolecule transistors. Phys. Rev. Lett. 95, 256803 (2005).

29. Roch, N., Florens, S., Bouchiat, V., Wernsdorfer, W. & Balestro, F. Quantum phase transition in a single-molecule quantum dot. Nature 453, 633637 (2008).

30. Sasaki, S., Tamura, H., Akazaki, T. & Fujisawa, T. Fano-Kondo interplay in a side-coupled double quantum dot. Phys. Rev. Lett. 103, 266806 (2009).31. Mitchell, A. K., Jarrold, T. F. & Logan, D. E. Quantum phase transition in quantum dot trimers. Phys. Rev. B 79, 085124 (2009).

32. Mitchell, A. K., Jarrold, T. F., Galpin, M. R. & Logan, D. E. Local moment formation and Kondo screening in impurity trimers. J. Phys. Chem. B 117, 1277712786 (2013).

33. Lucignano, P., Mazzarello, R., Smogunov, A., Fabrizio, M. & Tosatti, E. Kondo conductance in an atomic nanocontact from rst principles. Nat. Mater. 8, 563567 (2009).

34.itko, R. Fano-Kondo effect in side-coupled double quantum dots at nite temperatures and the importance of two-stage Kondo screening. Phys. Rev. B 81, 115316 (2010).

35. Craig, N. J. et al. Tunable nonlocal spin control in a coupled-quantum dot system. Science 304, 565567 (2004).

36.itko, R. & Bona, J. Enhanced conductance through side-coupled double quantum dots. Phys. Rev. B 73, 035332 (2006).

NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15210

37. Jones, B. A., Varma, C. M. & Wilkins, J. W. Low-temperature properties of the two-impurity Kondo Hamiltonian. Phys. Rev. Lett. 61, 125 (1988).

38. Afeck, I. & Ludwig, A. W. W. Exact conformal-eld-theory results on the multichannel Kondo effect: single-fermion greens function, self-energy, and resistivity. Phys. Rev. B 48, 7297 (1993).

39. Bork, J. et al. A tunable two-impurity Kondo system in an atomic point contact. Nat. Phys. 7, 901906 (2011).

40. Sela, E., Mitchell, A. K. & Fritz, L. Exact crossover green function in the two-channel and two-impurity Kondo models. Phys. Rev. Lett. 106, 147202 (2011).

41. Mitchell, A. K. & Sela, E. Universal low-temperature crossover in two-channel Kondo models. Phys. Rev. B 85, 235127 (2012).

42. Mitchell, A. K., Logan, D. E. & Krishnamurthy, H. R. Two-channel Kondo physics in odd impurity chains. Phys. Rev. B 84, 035119 (2011).

43. Mitchell, A. K., Sela, E. & Logan, D. E. Two-channel Kondo physics in two-impurity Kondo models. Phys. Rev. Lett. 108, 086405 (2012).

44. Keller, A. J. et al. Universal fermi liquid crossover and quantum criticality in a mesoscopic system. Nature 526, 237240 (2015).

45. Pedersen, K. G. L. et al. Quantum interference in off-resonant transport through single molecules. Phys. Rev. B 90, 125413 (2014).

46. Karan, S. et al. Shifting the voltage drop in electron transport through a single molecule. Phys. Rev. Lett. 115, 016802 (2015).

47. Baruselli, P. P., Fabrizio, M., Smogunov, A., Requist, R. & Tosatti, E. Magnetic impurities in nanotubes: from density functional theory to Kondo many-body effects. Phys. Rev. B 88, 245426 (2013).

48. Jacob, D., Soriano, M. & Palacios, J. J. Kondo effect and spin quenching in high-spin molecules on metal substrates. Phys. Rev. B 88, 134417 (2013).

49. Korytr, R. & Lorente, N. Multi-orbital non-crossing approximation from maximally localized Wannier functions: the Kondo signature of copper phthalocyanine on Ag (100). J. Phys.: Condens. Matter 23, 355009 (2011).

50. Schrieffer, J. R. & Wolff, P. A. Relation between the Anderson and Kondo Hamiltonians. Phys. Rev. 149, 491 (1966).

51. Tth, A. I., Borda, L., von Delft, J. & Zarnd, G. Dynamical conductance in the two-channel Kondo regime of a double dot system. Phys. Rev. B 76, 155318 (2007).

52. Costi, T. A., Hewson, A. C. & Zlatic, V. Transport coefcients of the Anderson model via the numerical renormalization group. J. Phys.: Condens. Matter 6, 2519 (1994).

53. Schiller, A. & Hersheld, S. Exactly solvable nonequilibrium Kondo problem. Phys. Rev. B 51, 12896 (1995).

54. Sela, E. & Afeck, I. Nonequilibrium transport through double quantum dots: exact results near a quantum critical point. Phys. Rev. Lett. 102, 047201 (2009).

55. Mitchell, A. K., Landau, L. A., Fritz, L. & Sela, E. Universality and scaling in a charge two-channel Kondo device. Phys. Rev. Lett. 116, 157202 (2016).

56. Schiller, A. & Hersheld, S. Theory of scanning tunneling spectroscopy of a magnetic adatom on a metallic surface. Phys. Rev. B 61, 9036 (2000).57.jsghy, O., Kroha, J., Szunyogh, L. & Zawadowski, A. Theory of the fano resonance in the STM tunneling density of states due to a single Kondo impurity. Phys. Rev. Lett. 85, 2557 (2000).

58. Ohno, K. Some remarks on the Pariser-Parr-Pople method. Theor. Chim. Acta 2, 219227 (1964).

59. Pariser, R. & Parr, R. G. A semi-empirical theory of the electronic spectra and electronic structure of complex unsaturated molecules. II. J. Chem. Phys. 21, 767776 (1953).

60. Pople, J. A. Electron interaction in unsaturated hydrocarbons. Trans. Faraday Soc. 49, 13751385 (1953).

61. Soos, Z. G. & Ramasesha, S. Valence-bond theory of linear Hubbard and Pariser-Parr-Pople models. Phys. Rev. B 29, 5410 (1984).

62. Izumida, W., Sakai, O. & Shimizu, Y. Many body effects on electron tunneling through quantum dots in an Aharonov-Bohm circuit. J. Phys. Soc. Jpn 66, 717726 (1997).

63. Weichselbaum, A. & von Delft, J. Sum-rule conserving spectral functions from the numerical renormalization group. Physical Review Letters 99, 076402 (2007).

64. Anders, F. B. & Schiller, A. Real-time dynamics in quantum-impurity systems: a timedependent numerical renormalization-group approach. Phys. Rev. Lett. 95, 196801 (2005).

65. Oguri, A. & Hewson, A. C. NRG approach to the transport through a nite Hubbard chain connected to reservoirs. J. Phys. Soc. Jpn 74, 988996 (2005).

66. Meir, Y. & Wingreen, N. S. Landauer formula for the current through an interacting electron region. Phys. Rev. Lett. 68, 2512 (1992).

Acknowledgements

We thank Eran Sela and Martin Galpin for fruitful discussions. A.K.M. acknowledges funding from the D-ITP consortium, a program of the Netherlands Organisation for Scientic Research (NWO) that is funded by the Dutch Ministry of Education, Culture and Science (OCW). The Center for Quantum Devices is funded by the Danish National Research Foundation. We are grateful for use of HPC resources at the University of Cologne.

Author contributions

A.K.M. wrote the NRG code, performed NRG calculations and derived analytic results. K.G.L.P., P.H. and J.P. formulated the 2CK mapping. K.G.L.P derived effective model parameters. All authors prepared the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Mitchell, A. K. et al. Kondo blockade due to quantum interference in single-molecule junctions. Nat. Commun. 8, 15210 doi: 10.1038/ncomms15210 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

10 NATURE COMMUNICATIONS | 8:15210 | DOI: 10.1038/ncomms15210 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications