Coherent population trapping by dark state formation in a carbon nanotube quantum dot

Illumination of atoms by resonant lasers can pump electrons into a coherent superposition of hyperfine levels which can no longer absorb the light. Such superposition is known as a dark state, because fluorescent light emission is then suppressed. Here we report an all-electric analogue of this destructive interference effect in a carbon nanotube quantum dot. The dark states are a coherent superposition of valley (angular momentum) states which are decoupled from either the drain or the source leads. Their emergence is visible in asymmetric current−voltage characteristics, with missing current steps and current suppression which depend on the polarity of the applied source-drain bias. Our results demonstrate coherent-population trapping by all-electric means in an artificial atom. Transport in quantum systems is complex and can be suppressed by coherent superposition of the involved states. Here, the authors find all-electronic suppression of transport in a carbon nanotube originating from coherent population trapping and give criteria for the presence of such a dark state.

C oherent population trapping (CPT) occurs in three-level, Lambda-type atomic systems coupled to two quasiresonant electromagnetic modes 1,2 . The light beams coherently excite the atom from two low-lying states to a common excited state. By proper detuning of the lasers, the system has a finite probability to decay into a coherent superposition of the low-lying states which is decoupled from the light, a so-called dark state (DS). In the stationary limit, the DS is occupied with probability one and the CPT is perfect.
Quantum dots offer the possibility to engineer artificial atoms and molecules by proper circuit design, and hence to probe CPT in effective Lambda-systems. Proposals [3][4][5] have considered microwave irradiated double quantum dot analogs of the seminal experiment 6 . Since localization of the electrons in the DS also implies a vanishing current through the double quantum dot, this allows the electrical detection of CPT by recording variations of the current as the microwaves parameters are tuned.
All-electrical realizations of CPT are also desirable. Several proposals have predicted trapping in a coherent superposition of spatially localized states under appropriate bias conditions, 7-11 whose signatures have been observed in vertical double dot experiments 12,13 . A second possibility is, similar to the original optical experiment, CPT in a linear superposition of energy eigenstates, predicted for highly symmetric triple quantum dot setups 14,15 or for molecular junctions with intrinsic orbital degeneracies [16][17][18] .
In this work we report the experimental realization of the second all-electrical realization of CPT in a single carbon nanotube (CNT) quantum dot. The effect requires the presence of quasi orbitally degenerate states which can form coherent superpositions. Under given conditions, tunneling events into and out of the dot successively trap the system in a DS, i.e., a coherent superposition of the degenerate levels which is decoupled from one of the leads. This in turn yields a characteristic current suppression as the bias voltage or the gate voltage are tuned. The situation is illustrated in the quantum dot setup of Fig. 1a for the case of a positive electrochemical potential drop between left and right leads. The coherent superposition of two degenerate states results in a coupled state (CS) and a DS which is decoupled from the right lead. This allows electrons to enter the DS from the left while preventing them to leave it to any of the two leads. CPT occurs and current is suppressed. For opposite bias no suppression takes place. There is a close analogy between the experiment generating dark states in atoms and the one discussed here: In both cases the coupling to an external drive (the laser fields for the atom, the DC bias for the quantum dot) enables inelastic transitions among three energy eigenstates of the isolated system. The spontaneous formation of a coherent superposition of system eigenstates (the dark state) suppresses these transitions, despite the presence of the driving. Here we demonstrate that such a situation has been realized in a CNTbased quantum dot. We notice that while in the optical setup it is the coherent drive that defines the linear combination yielding the DS 1 , in the transport setup the effect is more subtle. As seen in Eq. (2) below, the linear combination is set by the phase of the tunneling matrices which couple the quantum dot to the leads.

Results
Dark states in CNTs. Similar to graphene, CNTs posses an orbital (valley) degree of freedom, arising from the two inequivalent Dirac points K, K′ in the Brillouin zone. In CNTs of the zig-zag class, such orbital degree of freedom corresponds to the longitudinal orbital momentum h' z , accounting for clockwise ' z ¼ À' ð Þ or anticlockwise ' z ¼ ' ð Þ rotations along the tube waist 19 (see Fig. 1b). The finite length of the CNT quantum dot results in discrete longitudinal momentum states, with quantum number m, and hence in a shell structure like for the atomic bound states. In the absence of spin−orbit coupling [20][21][22] and valley mixing 21,[23][24][25] , each shell consists of four degenerate bound states with spin, σ = ↑, ↓, and angular momentum, ' z ¼ ± ', degrees of freedom. Thus each shell can accommodate up to N = 4 electrons. As shown in the following, for a CNT quantum dot near the N = 0 ↔ N = 1 resonance, a coherent superposition of angular momentum states can form which is decoupled from one of the two leads, and hence is a DS for the quantum dot for appropriate polarity of the applied bias voltage (see Fig. 1c). Due to the particle-hole symmetry of the many-body spectrum, CPT is expected near the N = 3 ↔ N = 4 resonance if the voltage polarity is reversed.
Experimental signatures of CPT. Measurements are performed on a suspended CNT grown on top of prepatterned leads. Such CNTs are usually called ultraclean owing to their low level of impurities 26 . In Fig. 2a we show the experimentally measured differential conductance G of our ultraclean CNT quantum dot as a function of the applied bias voltage V b and of a back-gate voltage V g . Coulomb diamonds are clearly visible, with a characteristic fourfold periodicity, a signature of the successive filling of CNT shells with four electrons each. Noticeably, three almost identical diamonds are followed by a larger one. The width of a Coulomb diamond is a measure of the energy required to fill the CNT with an extra electron, which accounts for charging effects and single-particle properties. For the small diamonds only a charging energy has to be paid, which indicates almost degenerate states within a shell (negligibly small spin−orbit coupling and valley mixing), as well as a small exchange energy for the middle diamond. For the larger diamonds the shell is full (N = 0,4,8,…), such that filling the CNT with an extra electron requires to pay a a in out Fig. 1 All-electronic dark states in a carbon nanotube quantum dot. a Quantum dot setup to probe coherent population trapping in a dark state (DS). For the chosen bias voltage polarity, an electron can enter a DS from the left lead but it cannot leave it by tunneling to the right lead. Precession allows population transfer from the dark state to the coupled state (CS). b A CNT features degenerate angular momentum states. Linear combinations of angular momentum states for single-electron occupancy (N = 1) of a given nanotube shell can give rise to CS and DS. The latter is represented by the eight-shaped orbital in panel (c). c Effective Lambda-system for coherent population trapping in a CNT quantum dot. Because an electron in the DS cannot leave the carbon nanotube, the state N = 0, corresponding to an empty shell configuration, can no longer be reached. CNT carbon nanotube charging energy and the mean inter-shell spacing. By closer inspection of the current voltage characteristics in shell 1 and shell 2, we observe signatures of current suppression in the form of faint Coulomb diamond edges and negative differential conductance (NDC) for the 0 ↔ 1 transition, as indicated by the red and blue arrows, respectively. The same pattern occurs also at the 3 ↔ 4 transition for opposite bias polarity.
The most frequent reason for NDC is weakly conducting (slow) excited states, resulting from strong asymmetries in the tunneling coupling of the system to the leads. This mechanism is improbable here, as it implies very different amplitudes in the coupling between time reversal partners (the different valley states) and the same lead. The strong correlation of the NDC with the conductance suppression at the edges of both the N + 1 and N + 3 Coulomb diamonds, and the absence of the latter in the N + 2 diamonds, suggests that CPT with formation of electronic dark state provides the correct interpretation of our experiments. To support such a statement, we have compared experimental gate traces and bias traces with theoretical calculations for the stationary current of a CNT. Figure 2b shows experimental gate traces for V b = ±3.045 mV (green/purple line in Fig. 2a); the numerical calculations are depicted in Fig. 2c. The parameters used can be found in Table 1. The same outcome is seen in the theoretical traces. For further inspection, we have focused on the 0 ↔ 1 and 3 ↔ 4 transitions of shell 1. The respective experimental stability diagrams are shown on an enlarged gate voltage scale in Fig. 3a, b, while the corresponding numerical results are depicted in Fig. 3c, d. We observe similar behavior as in the experiment within a range of 3 mV around zero bias. At larger voltages, besides additional diagonal lines due to excited states, the experimental curves display a horizontal line not understood at present. Experimental and theoretical gate traces at fixed V b are shown together in Fig. 3e, f.
At the 0 ↔ 1 resonance both the experimental and theoretical gate traces show a rectangular shaped current at positive bias, typical of quantum dot behavior in the sequential tunneling regime; at negative bias, however, the current first increases and then gradually decreases as the gate increases, indicating trapping of a single electron. At the 3 ↔ 4 resonance similar current shapes are observed for opposite bias voltage polarity and upon gate voltage mirroring, a signature of trapping of a single hole. The dependence of the current on the bias voltage is analyzed in more detail in Fig. 4. We show the current for the 0 ↔ 1 transition at V g = 11.771 V and V g = 11.781 V in Fig. 4a, b, respectively; these positions are marked in Fig. 3c by vertical lines and the corresponding symbols. The NDC and the faint (almost missing) resonant line are highlighted by a blue and red arrow, respectively.
Again, the agreement between theory and experiment is remarkable. As discussed below, all the characteristic features observed in Figs. 2-4 can be explained in terms of CPT in a DS, combined with a precessional motion which transfers population between the dark and the coupled state, as sketched in Fig. 1a.
Orbital degeneracy and CNT spectrum. The general ingredients to describe charge transport across a quantum dot in the   2 Experimental signatures of coherent population trapping. a Experimental differential conductance as a function of back-gate voltage V g and bias voltage V b . Twelve consecutive Coulomb diamonds can be assigned to three carbon nanotube shells which get progressively occupied with N = 1 to N = 4 electrons. Current suppression (negative differential conductance) and faint Coulomb diamond borders are observed, indicated by blue and red arrows, respectively. b Current vs. gate voltage for the two values V b = ±3.045 mV of the bias voltages corresponding to the green/purple lines in panel (a). Solid lines provide here a guide to the eye. Current suppression associated to coherent population trapping is indicated by red arrows. c Numerically evaluated stationary current qualitatively reproducing the experiment. The parameters used in the simulation are shown in Table 1   Table 1 Parameters Numerical parameters to fit all three shells of the experiment. Only Γ L and Δϕ vary with the shell; all other parameters are the same for all shells sequential tunneling regime are a tunneling Hamiltonian coupling the dot to the electrodes, and the electronic spectrum of the isolated system in the energy range set by the electrochemical potentials of the lead electrodes. The presence of orbital degeneracies (or quasi-degeneracies) is decisive for the occurrence of CPT. As discussed above, the single particle energy spectrum of a CNT of finite length is fully characterized by a shell quantum number m and the pair ðσ; ' z Þ, accounting for the spin and orbital degrees of freedom. Curvature-induced spin−orbit coupling and valley mixing remove the intra-shell degeneracy. The amplitude of the spin−orbit coupling is largest near the Dirac point and of the order of a fraction of meV. However, it strongly decreases for states away from the bottom of the CNT conduction band 23 , which is the case for the gate voltage range in which CPT is seen in our experiment. Similarly, valley mixing due to disorder is strongly suppressed in ultraclean CNTs, and is forbidden by symmetry in CNTs of the zig-zag class 25 , which suggests that we have measured such kind of tube in our experiment. We neglect both perturbations in the following. The Hamiltonian we used to describe the CNT quantum dot is thus given in Eq. (5) and it considers three shells only m = 0, 1, 2. The only relevant parameters of the Hamiltonian are the inter-shell spacing ε 0 , the interaction U accounting for charging effects and the exchange interaction J. A many-body state with N particles is thus characterized by its total energy E, the total angular momentum L z , and the total spin quantum numbers S and S z , i.e., it has the form |N, E; S, S z , L z 〉. In our three-shell model, we have fixed the energy E 0 and the particle number of the configuration with the shell m = 0 completely full, and the upper two shells m = 1, 2 completely empty. The N = 0 groundstate |0, E 0 ; 0, 0, 0〉 ≡ |0〉 is depicted in Fig. 5. The N = 1 groundstate is fourfold degenerate. A basis is the quadruplet of states 1; E 1 ; 1 2 ; σ; ' z È É obtained by adding one electron with quantum numbers ðσ; ' z Þ on shell m = 1. These states are also graphically shown in Fig. 5, where we used the short-cut notation jσ; ' z i. Further examples of many-body states with N = 1, 2 electrons are in Supplementary Note 1. Here we exemplarily focus on CPT at the 0 ↔ 1 resonance, which involves the N = 0 and the N = 1 groundstates.
Tunneling matrix. The tunneling Hamiltonian allows for transitions between dot states with different particle number. Its form, given in Eq. (6), is rather standard. The associated complex tunneling amplitude t αkm' z accounts for the overlap between an electron wave function in lead α, characterized by the momentum k, and a CNT wave function for shell m and angular momentum ' z in the contact region. We assume that, due to the nanotube curvature, tunneling is local and occurs only through those CNTs atoms closest to the leads. In this case, as discussed in the Methods, the rate matrix Γ m α , defined element-wise by is in general nondiagonal in the angular momentum basis. Figure 6 shows the situation in which a finite amount of CNT atoms is locally contacted to the leads. In a realistic setup the leads are not flat and only few CNT atoms have a relevant tunneling coupling. For a single atom contact, or in the more general surface Γ-point approximation discussed in the Supplementary Note 3, it takes the simple form Γ m where the phase ϕ m α describes a global property of contact α for shell m. Both approximations can be gradually relaxed giving rise to an attenuation of the destructive interference. Exp Theory Theory 11.79 12.  Dynamics and coherent population trapping. The dynamical quantity of interest here is the stationary current I. It follows from the stationary reduced density operatorρ 1 ¼ lim t!1ρ ðtÞ and the current operatorÎ α at lead α according to I :¼ I L ¼ Tr CNT fÎ Lρ 1 g ¼ ÀI R . Notice that this requires to take the trace over the full spectrum of the isolated CNT. For weak tunneling coupling, the dynamical equations forρ are easily obtained from the Liouville−von Neumann equation for the total density operator by treating the tunneling Hamiltonian as a perturbation 16,17 . For general values of the gate and bias voltages such equations have to be solved numerically. Analytical solutions are possible when the system is tuned near a transition involving only N and N + 1 particles groundstates, which is the case of interest here.
Let us consider the 0 ↔ 1 transition. In the one-dimensional N = 0 subspace the density matrix is a number, ρ 0 . In the fourdimensional N = 1 subspace, it is block-diagonal in spin (since spin is conserved during tunneling) but not in angular momentum. The contributions from different spin configurations can be summed up in the dynamical equations yielding a set of coupled equations for ρ 0 and a 2 × 2 matrix ρ 1 15 . Away from the exact resonance (i.e. from the border of the Coulomb diamond), one finds for positive electrochemical potential drop eV b ) k B T, with Γ m α ¼ Γ α and the coherence matrices R α . Equation (1) describes relaxation governed by the terms Γ α R α as well as precession through the Lamb Shift contribution α ω α R α . The latter originates from virtual processes from the system to the leads 14,18,27 and its effect will be discussed later. Due to the nondiagonal form of the R α in the angular momentum basis, also the stationary density matrix ρ 1 1 is not diagonal there. Its diagonalization yields the stationary eigenstates, and the associated eigenvalues define the occupation probabilities of the eigenstates. To proceed, we assume for the coherence matrices the simple form R α' z ' ′ z ¼ e iϕ α ð' z À' ′ z Þ and ϕ R ≠ ϕ L . In this case, eigenstates of ρ 1 1 are a decoupled state (DS) at the right lead and its associated orthogonal coupled state (CS), with eigenvalues 1 and 0, respectively. Hence, Eq. (1) predicts CPT with probability one in the DS, and consequently a vanishing stationary current (see Fig. 1a). Accounting also for the spin degree of freedom, the dark/coupled states at the right lead are the linear combinations jDS; σ; Ri where R indicates the right lead. Equation (2) also implies that | DS, σ; R〉 satisfies the equation ϕ Rd À'σ Þ is the electron annihilator involved in the local tunneling process at the right lead. Notice the explicit dependence on the tunneling phase acquired upon tunneling at the lead R. If ϕ R = ϕ L , the DS in Eq. (2) is also decoupled at lead L; transport solely occurs through the CS and CPT cannot occur. Thus CPT requires ϕ R ≠ ϕ L .
The DS wave function in Eq. (2) is explicitly shown in Fig. 6 on the example of a (12, 0) CNT. We have assumed the angular coordinate of the contact atoms at the right and left lead to be rotated by a small angle θ = π/24. The DS has a node at the contact positions at the right lead but not at the left lead. The corresponding CS is shown in Supplementary Fig. 1 and has finite weight at both contacts.
Let us now turn to the impact of the Lamb shift term in Eq. (1). It introduces a precession of the Bloch vector in the CS/DS basis of Eq. (2) with population transfer between dark and coupled states. The frequencies ω L , ω R are given in Eq. (11). For the situation indicated in Fig. 1a, ω L ≠ 0 allows the electrons in a DS to precess into the coupled state and from there to escape, yielding a small stationary current I ¼ 4eΓ L ρ 1 0 . We find from Eq. (1) the expression for Δϕ = (ϕ L − ϕ R )/2 ≠ 0. Because the precession frequencies strongly depend on gate and bias voltage, also the effectiveness of CPT does, as clearly seen from the gradual variation of the current in Fig. 4. Here, the traces in panels (a) and (b) correspond to two distinct values of V g . Clearly, at negative bias voltages (corresponding to positive eV b ), CPT is more pronounced in panel (b) than in panel (a), indicating a smaller ω L . As the bias polarity is changed, also the role of the precession frequencies is exchanged. The fact that at positive bias the current shows standard Coulomb steps, giving rise to rectification, is because of large ω R (compared to ω L ) for both of the chosen V g values. This rectification is a consequence of the asymmetry in the coupling to the leads. Equal couplings would restore the left−right symmetry, the rectification would vanish but the interference blocking would remain. The dependence of the frequencies ω L/R on V b for the same parameters used in the simulation of Fig. 4a is explicitly shown in Supplementary Fig. 2. In Supplementary Fig. 3, we additionally show the dependence of the current when additional inelastic relaxation processes are added in the numerical simulations.
We observe that, due to the particle-hole symmetry of the spectrum with respect to half-filling (shell with N = 2), the equations for the reduced density operator near the 3 ↔ 4 transition immediately follow from Eq. (1) upon exchanging R with L and V b with −V b . This implies that the CPT features at the 3 ↔ 4 transition can be obtained from the ones at the 0 ↔ 1 resonance by mirroring of both V b and V g . This is clearly observed in Fig. 3.
Similar to the N = 0 groundstate, also the N = 2 groundstate is a singlet for positive exchange interaction J (see Eq. (7) of the Supplementary Note 2). However, CPT is not seen at the 1 ↔ 2 and 2 ↔ 3 resonances. Whether or not interference occurs depends on the strength of the exchange coupling, since J/2 is the separation to the first excited doublet of states which do not form dark states. For J ' hΓ ' eV b the excited doublet is soon into the transport window and no CPT is seen. The N = 2 states and a special realization of a DS for the case J = 0 are further discussed in Supplementary Note 1.

Discussion
The results presented so far show a remarkable quantitative agreement between the experimental data and the theoretical predictions, strongly supporting the claim that the observed current suppression features are due to CPT. A natural question is how robust CPT is, and under which conditions can it be observed in other CNT-based quantum dots. According to our model, the effect is quite generic, as the main requirements are: first, the presence of a symmetry S of the system yielding degenerate energy states (for weak symmetry breaking the level splitting should be smaller than the tunneling broadening Γ = Γ L + Γ R ); second, the tunneling matrices being not diagonal in the basis associated to the symmetry S and with modulus of the offdiagonal elements of the coherence matrices R α close to one; third, the strong Coulomb interaction enforcing single electron tunneling. The requirements above are simultaneously met for ndoped (electron conduction) CNTs of the zig-zag class, whose bound states have angular momentum (valley) degeneracy. For this nanotube class, symmetry breaking perturbations are spin −orbit coupling and valley mixing. While the former is an intrinsic property of the CNT and decreases away from the band gap, valley mixing is due to disorder or to perturbations which break the rotational symmetry. Suspended ultraclean CNTs 26 show very weak disorder and symmetry breaking mostly occurs due to the presence of contact leads in the nonsuspended portion of the tube. In a realistic CNT-device the curvature of the nanotube and some roughness of the contacts causes tunneling to occur locally through few single carbon atoms. This ensures on the one hand that the tunneling matrix is not diagonal in the angular momentum basis, and on the other that the tunneling is a small perturbation and hence that valley mixing is small. Furthermore, weak tunneling makes it easier to reach the sequential tunneling regime, which is typically observed for CNTs in the electron conduction regime 19 . We notice that the second requirement among the ones listed above rules out the possibility that the excited states blocking recently reported in a CNT-based quantum dot with broken fourfold degeneracy is due to the CPT discussed here 28 . While the three above-mentioned properties provide a sufficient criterion for the presence of dark states, they are not intended as a necessary one. There are indeed various setups where transport can be suppressed due to the formation of a coherent superposition of quantum mechanical states. However, it is quite different whether the coherent superposition occurs among localized position eigenstates 7,12,13 , or degenerate energy eigenstates as in our setup. For the system proposed in ref. 7 , the DS is an energy eigenstate of the isolated system whose formation does not require phase coherence during the tunneling process. We also comment on the interference phenomena observed in the linear conductance of molecular junctions in the strong coupling regime 29,30 . Since energy conservation is required only within an accuracy set by ΔE % hΓ, in the strong coupling regime, where Γ is the largest scale in the problem, the CPT based on linear superpositions of quasi-degenerate states of the isolated system loses its significance. Rather, interference based on a coherent superposition of different trajectories along the molecule becomes possible 31 . Similarly, the zero-bias interference discussed for offresonant transport through single molecules 32 require energynonconserving cotunneling transitions, and are thus distinct from the interference blockade phenomena presented in this work. We finally comment on the canyon of conductance observed in an InSb nanowire quantum dot in magnetic field 33 . In this experiment, the conductance suppression is visible throughout the entire gate voltage−magnetic field region and it encompasses both the cotunneling as well as the sequential tunneling regime. As a consequence of the strong coupling of the degenerate levels to the leads, Γ > k B T, higher order virtual processes (charge fluctuations) strongly influence the position of the conductance canyon, especially in the sequential tunneling regime. A detailed study of current suppression accounting for level broadening in the framework of the so-called 2vN approximation is found in ref. 34 .
In conclusion, we have demonstrated experimentally and theoretically that CPT with dark state formation can be realized in a CNT quantum dot by all-electrical means. In the case considered here, the orbital degeneracies result from the interplay of the tubular nanotube geometry and the underlying graphene honeycomb lattice. However, the phenomenon is rather generic and is expected to occur in other highly symmetric quantum dot systems, e.g. in organic molecules in scanning tunneling microscope setups 35 , or in symmetric triangular quantum dots 15 , in the weak tunneling regime. In other setups degeneracies can be tuned by an external magnetic field 12,13 . Trapping in a dark state can be achieved by simple tuning of the bias or gate voltage.

Methods
Sample details. On a highly p-doped Si substrate with a 500 nm thermally grown SiO 2 layer, ring-like electrode structures were defined by electron beam lithography and evaporation of 20 nm Re and 40 nm Co. A catalyst for the CVD process was deposited in the center of the electrode ring structure to increase the chance of a CNT connecting the contacts. The CNT growth was performed as the last sample fabrication step, a procedure that is aimed at yielding ultraclean devices 26 . Neither magnetoresistance nor magneto-optical Kerr effect (MOKE) measurements showed magnetic behavior of the Re/Co contacts after CVD, ruling out spin valve effects with a similar I-V characteristics 27 . The current as well as the differential conductance were measured using an Ithaco DL1211 I-V-converter.
Model. The single-electron transistor setup is modeled by the total Hamiltonian H ¼Ĥ CNT þĤ leads þĤ tun , describing the CNT quantum dot (QD) weakly coupled to leads by a tunneling HamiltonianĤ tun .
The CNT-QD consists of a set of longitudinal modes, called shells, with shell index m. The electrons have both spin σ ∈ {↑, ↓} and valley ' z 2 f'; À'g degrees of freedom and, since spin−orbit coupling and valley mixing are negligible in the experiment, the single-particle spectrum of a single shell is assumed to be fourfold degenerate. We include a charging term U and the exchange interaction J. The latter strongly depends on the CNT chirality and radius and was assumed to be 10 μeV. The CNT-QD Hamiltonian reads 19,36 where in the numerical calculations only three shells (m ∈ {0, 1, 2}) are considered (see Fig. 5). The gate voltage V g applied with a level arm α g ensures particle-hole symmetry of Eq. (5) with respect to shell m = 1 for eα g V g = ε 0 + 6U + J/4. Further, the occupation operatorn m' z ¼ P σd y m' z σdm'z σ , and the spin operator b S m' z ¼ 1 2 P σσ′d y m' z σ s σσ′dm' z σ′ are defined in terms of annihilation (creation) operatorsd ðyÞ m' z σ of an electron in shell m with angular momentum ' z and spin projection σ.
The CNT Hamiltonian can be diagonalized analytically by using the basis corresponding to the eigenstates of the total particle numberN ¼ P m' zn m' z , total spinŜ 2 ¼ P m' zŜ 2 m' z , total spin projectionŜ z ¼ 1 2 P m' z σ σd y m' z σdm'z σ , and total angular momentum operatorL z ¼ P m' z ' znm' z . Accordingly, many-body states are uniquely defined by the vector set {|N, E; S, S z , L z 〉}. In our three-shell model, the N = 0 groundstate corresponds to the shell m = 0 being completely full, and the N = 1 groundstate is given by the quadruplet 1; ε 0 ; 1 2 ; σ; ' z È É , as shown in Fig. 5. Examples of many-body states with occupation N = 2 can be found in Supplementary Note 1.
The electrons in the leads are considered as fermionic reservoirs of noninteracting electrons at temperature T and chemical potentials μ L = μ 0 + eηV b , μ R = μ 0 + e(η − 1)V b for the left (L) and right (R) lead, respectively, with V b the applied bias voltage. The parameter 0 < η < 1 accounts for an asymmetric bias drop at the two leads. For the investigated setup a good fit to the data was obtained for almost symmetric potential drop (η = 0.55). For simplicity, the underlying Hamiltonian is chosen to be the one of a noninteracting electron gaŝ H leads ¼ P αkσ ε kĉ y αkσĉαkσ , whereĉ ðyÞ αkσ annihilates (creates) an electron in lead α ∈ {L, R} with momentum k and spin projection σ.
The leads are weakly coupled to the CNT via the tunneling Hamiltonian Thus,Ĥ tun ¼ P αĤtun;α removes (adds) an electron from the left/right lead (α = L/R) with momentum k, energy ε k and spin σ and it adds (removes) an electron in the dot with the same spin σ in the state ðm; ' z ; σÞ.
The tunneling Hamiltonian Eq. (6) rules the dynamics of the coupled CNTleads system. For weak coupling, it can be treated as a perturbation. Expanding the Liouville von Neumann equation for the total density operator up to second order in H tun , and taking the trace over the reservoirs, the equations for the reduced density operator are obtained 16,17 . They are ruled by a tunneling kernel, which, according to Eq. (6), has the form and is in general nondiagonal in the angular momentum basis. Rather Γ m α ¼ Γ m α R α , where the coherence matrices R α are hermitian. Explicitly, R α' z ' z ¼ 1 while the off-diagonal elements satisfy 0 jR α' z ;À' z j 1. As discussed in Supplementary Note 3, the R α become exactly diagonal in the limit in which all the atoms of the CNT are equally coupled to the leads. In the opposite case of a single atom contact along the CNT circumference, or in the more general so-called surface Γ-point approximation, see Supplementary Note 3, one finds the optimal coherence, and the coherence matrices at lead α = L,R assume the form R α' z ' ′ z ¼ e iϕ α ð' z À' ′ z Þ used in the main manuscript. Then, in the dark and coupled states basis of Eq. (2), R R becomes diagonal but not R L . Explicitly, The good agreement found between theory and experiment suggests that in our setup the coherence matrices have off-diagonal elements with modulus close to one in the angular momentum basis, i.e. tunneling occurs only at few atomic positions where the CNT is closest to the leads.
We observe that in some nanotube-based devices the quantum dot region is defined electrostatically. This yields a strong variation of the tunneling coupling and charging energy, and in turn of the current, with the gate voltage 19 . In contrast, in our experiment we did not observe strong variations of the conductance as the gate voltage was swept from the p-doped region across the band gap down to the deep n-doped regime considered here. This supports our assumption that the bottleneck for transport is given by the contact region between the metallic leads and the CNT.
Tunneling dynamics at the 0 ↔ 1 resonance. For the parameters used in this work, the single particle groundstate is fourfold degenerate, while the two particles groundstate is only at a distance J=2 % hΓ from the doublet of first excited states, as discussed in Supplementary Note 1. Thus, the secular approximation is nonvalid when transitions involving the low-lying two-particle states are considered. To account for the nonsecular contributions, the equations for the reduced density matrix have thus been derived along the lines discussed in ref. 17 . In addition, a relaxation term was added to account for inelastic processes due to e.g. phonons. Such equations have been then solved numerically within the three-shell approximation for the CNT spectrum discussed above, and by including all excitations up to a cut-off energy 1.5ε 0 .
Here we focus on the dynamics near the 0 ↔ 1 resonance, where only the 0 and 1 particle groundstates play a role and an analytical treatment is possible. Furthermore, the secular approximation holds true yielding for the stationary reduced operator the equation where L is the Liouville superoperator that contains the system internal dynamics in a commutator structure together with the Lamb shift contributionĤ LS , a tunneling part L tun and a relaxation part L rel . For positive potential drop eV b , the equation of motion for the ρ 0 and ρ 1 subblocks, where the contribution over spin configurations has been summed, reads where ρ 1,th = 1/2 is the thermal density matrix for the one-electron subblock. The dynamics is thus fully governed by the coherence matrices R α . They enter in the tunneling terms Γ α R α as well as in the angular Lamb shift Hamiltonian matrix The latter introduces a precession with frequencies 17 , where we chose an offset in the gate voltage such that the resonance is exactly at where ψ is the digamma function. These precession frequencies clearly depend on the gate voltage and, via the chemical potentials, also on the bias voltage. This dependence is further analyzed in Supplementary Note 2. Notice also the dependence of ω α on the charging energy U and exchange interaction J, that results from virtual processes involving also the two electron subspaces. The above equations are general. To proceed, we assume full coherence, such that R α' z ' ′ z ¼ e iϕ α ' z À' ′ z ð Þ . Then, when expressed in matrix form, Eq. (10) yields the set of coupled equations for populations and coherences À Γ R 2 e 2iϕ R þ iω À Á ðρ 'À' Þ Ã À Γ rel 2 ρ '' À ρ À'À' À Á ; _ ρ À'À' ¼ ÀΓ R ρ À'À' þ 2Γ L ρ 0 À Γ R 2 e À2iϕ R À iω Ã À Á ρ 'À' À Γ R 2 e 2iϕ R þ iω À Á ðρ 'À' Þ Ã À Γ rel 2 ρ À'À' À ρ '' withω ¼ ω L e 2iϕ L þ ω R e 2iϕ R . These equations can be solved in the stationary limit with _ ρ 0 ¼ 0 and _ ρ ';' ′ ¼ 0, together with the condition Trρ ¼ ρ '' þ ρ À'À' þ ρ 0 ¼ 1. If ϕ L = ϕ R and Γ rel = 0 the dark state is completely decoupled from the dynamics and therefore the stationary solution is not uniquely defined but rather depends on the initial state. Any finite relaxation rate or phase difference solves this problem. Since the general expression for the current is quite lengthy, we focus on two limiting cases. For vanishing relaxation rate but Δϕ = (ϕ L − ϕ R )/2 ≠ 0, we obtain the current IðΓ rel ¼ 0Þ ¼ 4eΓ L ρ 0 1 given in Eq. (4) of the main manuscript. As mentioned before, this expression is only valid for a non-zero phase difference. For negative bias the current is the same upon exchanging L ↔ R and the overall sign. This expression shows that the current is completely suppressed for Δϕ = π/2 despite finite Lamb shift. The second limit that can be analyzed is a finite relaxation rate and Δϕ = 0. Interestingly, the dependence of the current on the relaxation rate drops completely and we obtain Notice that this same expression for the current holds true for transport through fourfold degenerate levels in the absence of interference. The behavior of NATURE COMMUNICATIONS | https://doi.org/10.1038/s41467-018-08112-x ARTICLE NATURE COMMUNICATIONS | (2019) 10:381 | https://doi.org/10.1038/s41467-018-08112-x | www.nature.com/naturecommunications the current as a function of the phase difference Δϕ is further analyzed in Supplementary Note 2.
Parameters. All parameters used for numerical calculation can be found in Table 1. Only a few of them can be considered as free fitting parameters. In fact, the size of the Coulomb diamonds, the slopes of the resonant lines corresponding to transitions between ground states, and the height and broadening of the Coulomb oscillations allow us to fix most of them (i.e. U, ε 0 , T, η and the average tunneling rate (Γ L + Γ R )/2). This analysis leaves four fitting parameters (Γ L /Γ R , Δϕ, J and Γ rel ). However, J and Γ rel are not sensitive parameters, as long as both are not much larger than the tunneling rates. The crucial fitting parameters are instead the difference of the tunneling phases Δϕ and the ratio Γ L /Γ R of the tunneling rates. They control respectively the strength of the interference and the rectification in the I -V characteristics. Differently from the first set of parameters, the values of these last two were not extracted from experimental data, but rather tuned freely until convergence of the numerically calculated transport characteristics to the experimental ones.