# Numerically exact full counting statistics of the nonequilibrium Anderson impurity model

###### Abstract

The time dependent full counting statistics of charge transport through an interacting quantum junction is evaluated from its generating function, controllably computed with the inchworm Monte Carlo method. Exact noninteracting results are reproduced; then, we continue to explore the effect of electron–electron interactions on the time-dependent charge cumulants, first-passage time distributions and -electron transfer distributions. We observe a crossover in the noise from Coulomb blockade- to Kondo-dominated physics as the temperature is decreased. In addition, we uncover long-tailed spin distributions in the Kondo regime and analyze queuing behavior caused by correlations between single electron transfer events.

## I Introduction

Mesoscopic quantum dots, molecules in junctions and other small quantum systems coupled to baths are often studied in transport experiments. It is possible to measure not only the electronic current through such systems, but also its noise and higher order moments. The complete set of moments determines the full counting statistics (FCS) of the system, which grants direct access to otherwise concealed properties: for example, the entanglement entropy, the fidelity, and thermoelectric efficiency fluctuations.

Apart from special limits, theoretical modeling of FCS in fermionic systems has so far been restricted to either noninteracting situations, or approximations whose accuracy is difficult to gauge. This difficulty is exacerbated in the experimentally relevant cases where strong electron–electron correlations and nonequilibrium effects are in play. A controlled and fully quantum description of the current and its moments is therefore sorely needed.

In this paper, we provide such a description for a model of interacting fermions in a junction, within a numerically exact inchworm quantum Monte Carlo (iQMC) calculation. By computing the generating function of lead population change (Fig. 1) we gain access to all moments of population transfer, first passage time distributions (FPTDs) and -electron transfer distributions. We analyze Fano factors and consider the shapes of the distributions, which provide remarkably detailed information about the dynamics of single tunneling events. We extract indications of queuing effects in the presence of strong interactions.

The paper proceeds as follows. Section II will motivate the full counting statistics, and Section III provide a brief overview of the main definitions and concepts used in FCS. In Section IV we describe how FCS can be obtained within iQMC, specializing to a particular model. Next, we present results for non-interacting (Sec. V) and interacting (Sec. VI) systems. The effect of interactions, temperature and voltage on the first and second charge cumulants is discussed, and we further explore how first passage times and -particle passage probabilities depend on the initial preparation. Finally, we discuss our conclusions in Section VII.

## Ii Background and Motivation

The transfer of charges across a small system, such as a molecular electronic junction or a mesoscopic quantum dot, is fundamentally a stochastic process Blanter and Büttiker (2000). Like any stochastic process, it is fully characterized by time dependent probability distributions Kampen (2007); in this case, , where is the number of charges transferred by time Braggio et al. (2006). Fluctuations in the electron current have long been known to provide experimentally accessible information not contained in the expectation value of the current Blanter and Büttiker (2000); Landauer (1998), such as dwell times Beenakker et al. (2003) and transmission probabilities Vardimon et al. (2013).

Current noise, the second moment of the current, has been of interest as a way to investigate reflection processes and traversal times in molecular junctions Selzer and Peskin (2013); Yang et al. (2014); Ochoa et al. (2015); Ridley et al. (2017) and propagation of correlated electron–hole pairs in photon-assisted transport Rychkov et al. (2005); Battista et al. (2014); Zamoum et al. (2016). A great deal of work has focused on fluctuations in periodically driven systems, including the prediction Levitov et al. (1996); Dubois et al. (2013a) and discovery Dubois et al. (2013b) of noiseless Lorentzian pulses, or “levitons”. Quantum noise has also been used to extract effective quasiparticle charges of transmitted electrons due to Cooper pair formation in superconducting junctions () Khlus (1987); Lefloch et al. (2003); Braggio et al. (2011), the quantum Hall effect () Saminadayar et al. (1997), or two-particle inelastic backscattering processes in the Kondo regime () Sela et al. (2006); Golub (2006, 2007); Zarchin et al. (2008); Yamauchi et al. (2011); Ferrier et al. (2016).

The noise-to-signal ratio, also known as the Fano factor , can be used to characterize charge transfer statistics. It takes the value when the charge transfer is Poissonian (completely uncorrelated), as generally occurs for small transmission probabilities Levitov and Reznikov (2004) and weak coupling between the molecule and bath Bomze et al. (2005). Super-Poissonian shot noise () has been observed in single Safonov et al. (2003) and double Djuric et al. (2005); Kießlich et al. (2007) quantum dot junctions and is generally associated with strong electron bunching and quasiparticle formation Blanter and Büttiker (2000) that may occur in the transient Feng et al. (2008) or steady state Zarchin et al. (2008); Ferrier et al. (2016) regimes. In addition, the Fano factor can be tuned so that it is sub-Poissonian () in quantum dot Gustavsson et al. (2006); Emary et al. (2007) and graphene nanoribbon Tworzydło et al. (2006) structures.

Instead of focusing only on the second moment of the current, it is possible to consider the full counting statistics (FCS) associated with the individual charge transfer events from an underlying moment-generating or characteristic function Kampen (2007), and recover all higher-order moments and cumulants. This idea was pioneered by Levitov and Lesovik Levitov and Lesovik (1993); Levitov et al. (1996), who evaluated the steady state FCS for noninteracting systems.

The FCS formalism for charge can be connected to entanglement entropy Klich and Levitov (2009) and fidelity Lesovik et al. (2006) in noninteracting systems, and was also generalized to a heat and work FCS giving access to efficiency fluctuations in thermoelectric junctions Esposito et al. (2015); Agarwalla et al. (2015). Parallel with the development of the field of spintronics, the FCS of spin currents has been developed in regular metallic lead junctions Di Lorenzo and Nazarov (2004, 2005); Di Lorenzo et al. (2006); Fransson and Galperin (2010) and magnetic tunnel junctions (MTJs) Tang and Wang (2014); Gaomin et al. (2017); Tang et al. (2017a), and can be applied, for example, to the detection of spin-singlet pairs in the electron counting statistics Di Lorenzo and Nazarov (2005); Schmidt et al. (2007). In the steady state, the second cumulant can be used to compute the zero-frequency current power spectrum Esposito et al. (2009); Sakano et al. (2011a), and in recent years a theory of FCS has been proposed which also gives access to frequency-dependent noise spectra Emary et al. (2007); Flindt et al. (2008); Ubbelohde et al. (2012). Theoretical Tobiska and Nazarov (2005); Saito and Utsumi (2008); Simine and Segal (2012) and experimental Pekola (2015) investigations have demonstrated that the generating function satisfies steady state fluctuation–dissipation relations equivalent to the condition for detailed balance, although recent work indicates violation of detailed balance when bound states are established with circular probability currents in the steady state Stegmann and König (2017) and modification of detailed balance when the system is undergoing transient evolution Altland et al. (2010). It has also been shown that for noninteracting electrons, the zeros of the generating function all lie on the negative real axis Abanov and Ivanov (2008, 2009); Utsumi et al. (2013), so that the position of zeros in the complex plane can shed light on interaction effects in the transport Kambly et al. (2011).

Recently, the time dependence of higher-order cumulants has become accessible to experiment Flindt et al. (2009); Fricke et al. (2010); Ubbelohde et al. (2012), motivating the study of transient cumulant generating functions and related objects which reveal more information on dynamics than the current or particle number. One such quantity, which is closely related to the generating function, is the waiting time distribution (WTD). The WTD is the probability distribution characterizing the time separating successive charge detection events Albert et al. (2011, 2012); Thomas and Flindt (2014). Experimentally, WTDs can be used to study the effect of coherence Brandes (2008) and complex internal molecular structures such as double quantum dot Wang et al. (2007); Welack et al. (2008); Schaller et al. (2009); Welack et al. (2009) or normal–superconducting junctions exhibiting Andreev reflection Rajabi et al. (2013); Droste and Governale (2016); Albert et al. (2016), as these effects can alter the traversal times for electrons crossing the system. At the level of practical quantum electronics, these dynamical effects determine the maximum operation rate for devices Elzerman et al. (2004).

Despite this flurry of theoretical activity, computing the FCS for generic interacting fermionic systems has so far proven to be challenging. In noninteracting systems, there now exist robust approaches to the study of both steady state and transient FCS, including an appealing coherent state path integral nonequilibrium Green’s function formalism (PI-NEGF) Tang and Wang (2014); Tang et al. (2014); Yu et al. (2016); Tang et al. (2017b) which can be perturbatively extended to interacting cases Souto et al. (2015); Tang et al. (2017b, c). Exact results are available for steady state FCS and noise characteristics of the AIM in the Fermi liquid regime and at the Toulouse point Komnik and Gogolin (2005); Gogolin and Komnik (2006a, b); Golub (2006, 2007); Gogolin et al. (2007). The Fermi liquid theory is, however, valid only for low voltages and temperatures Sela et al. (2006); Gogolin and Komnik (2006b); Mora et al. (2008, 2009); Sakano et al. (2011b); Oguri and Sakano (2013). Several approximate approaches to FCS in interacting systems have been successfully employed in appropriate regimes: a noncomprehensive list includes various pertubative expansions Tobiska and Nazarov (2005); Sakano et al. (2011a, b, 2012); Oguri and Sakano (2013); Dong et al. (2013); Esposito et al. (2015); Gogolin and Komnik (2006a); Schmidt et al. (2007); Schmidt and Komnik (2009); Park and Galperin (2011); Maier et al. (2011) and quantum master equation approaches Bagrets and Nazarov (2003); Brandes (2008); Flindt et al. (2008, 2010); Simine and Segal (2012); Albert et al. (2012); Kaasbjerg and Belzig (2015); Benito et al. (2016); Stegmann and König (2017); Luo et al. (2017); Kosov (2017); Wang et al. (2017). Nonperturbative numerical approaches also exist, but for models which do not capture the full complexity of interacting fermionic transport. These include hierarchical equations of motion in the spin-boson model Cerrillo et al. (2016) at high temperatures, and the density matrix renormalization group for the interacting resonant level model at its self dual point and at zero temperature Carr et al. (2011); Schmitteckert et al. (2014); Carr et al. (2015). Also, an equilibrium determinant QMC method was recently applied to cold atomic gases Humeniuk and Büchler (2017).

To date, no numerically exact method has accessed FCS in a generic nonequilibrium interacting fermionic model; nor is a method currently available which is (even in principle) suitable for arbitrary temperatures, bias voltages and interaction strengths. In fact, this is largely true even for the calculation of the second cumulant, or noise. In this paper, we provide such a numerically exact method for computing the FCS, based upon the recently developed inchworm diagrammatic quantum Monte Carlo (iQMC) method Cohen et al. (2015).

## Iii FCS: theory and general considerations

In a charge detection experiment, one studies a system in which some regions, labeled , are called “leads”. Lead is found to contain electrons on any given measurement of its total population at time . The change in charge on the lead is measured between some initial time and the measurement time . The FCS formalism aims to construct a generating function for the probability distribution characterizing , determined by the associated operator (we suppress lead index from now on). For this purpose, an auxiliary counting field is introduced, such that

(1) |

Here, is the probability that the measured number of electrons has changed by the integer after time . The distributions can be obtained from the generating function via

(2) |

and all moments and cumulants can be extracted by taking successive derivatives with respect to :

(3) | ||||

(4) |

The and contain equivalent information; in particular we note that and relate the first and second moments. Of particular interest Tang and Wang (2014); Tang et al. (2014) is the probability distribution for no charges passing across the junction, which defines the idle time probability :

(5) |

This can be used to extract the first passage time distribution :

(6) |

which describes the probability distribution of the time separating initialization of the counting experiment from the first detection of a charge transfer event Albert et al. (2012).

The counting statistics of any operator can be evaluated by way of an effective Hamiltonian transformation originally formulated by Levitov and Lesovik Levitov and Lesovik (1993). We will primarily be interested in the special case . The generating function for the moments of a change in the value associated with is given by Esposito et al. (2009)

(7) |

in terms of the modified, operator-specific propagator , which is defined by Esposito et al. (2009)

(8) |

The modified propagator itself can be written as a time-ordered integral with respect to a modified Hamiltonian:

(9) |

where

(10) |

This satisfies the relation . We express the generating function of Eq. (7) as a time ordered integral over the full Keldysh contour , composed of an “upper” branch and a “lower” branch Keldysh (1965). Using the contour representation, is given by the compact form

(11) |

where orders times later on the contour to the left.

The auxiliary field takes a different value on each branch:

(12) | ||||

When , it is immediately apparent that the generating function for all . In general, is periodic in with a periodicity of . Any value of for means that the Hamiltonian depends on the contour branch, such that time reversal symmetry is explicitly broken.

## Iv Model and Monte Carlo algorithm

### iv.1 Model and observables

We now specialize the discussion to the concrete case of the nonequilibrium Anderson impurity model (AIM) Anderson (1961). The model’s Hamiltonian can be written as follows:

(13) |

describes a small, interacting “dot” subsystem; depicts a set of large, noninteracting lead subsystems; and is a hybridization Hamiltonian coupling the dot and the leads. These three terms take the following form:

(14) | ||||

(15) | ||||

(16) |

Here, the denote dot annihilation operators labeled by a spin index , while the operators signify lead annihilation operators labeled by a spatial index , spin index and lead index . We assume two leads, denoted . The represent dot level energies, and the charging energy determines the strength of Coulomb repulsion between electrons. In this paper, we set (the particle–hole symmetric case, to which our method is not in any way limited). The hybridization parameters couple the lead states to the dot and the enumerate the energies of lead states. The system is initially prepared in the state , a decoupled equilibrium state of in which the chemical potential in the leads is given by , . is chosen to be in one of the four eigenstates of : the empty state ; the singly occupied spin up and spin down down states, , with ; and the doubly occupied dot state . These states will be denoted by the label in what follows.

The and are completely determined by the coupling density

(17) |

We choose this to be a flat band with a smooth cut-off Werner et al. (2009):

(18) |

In what follows, we set such that determines our unit of energy. We take the leads’ band cutoff to be , and their edge width to be .

At the initial time of the simulation, a coupling quench is performed and the subsystems are suddenly coupled by the introduction of into the Hamiltonian. This kind of switch-on is often referred to as the partitioned approach, and has been contrasted in the quantum transport literature to partition-free approaches (such as a voltage quench) in which the dot–lead coupling is established before a bias is added to the lead energies Stefanucci and Almbladh (2004). We note in passing that a partition-free voltage quench has also been explored within inchworm QMC Antipov et al. (2017), by augmenting the finite Keldysh branches with an imaginary time Matsubara branch, thus obtaining the full Konstantinov–Perel’ contour Konstantinov and Perel (1960).

In the present context, we are interested in the FCS for this model. We therefore set in the modified Hamiltonian of Eq. (10), where is the particle number operator in lead . This counting function describes the moments of population changes (rather than currents), but also provides access to the mean time-dependent current flowing out of the lead, , by way of a time derivative. The current noise and higher order cumulants can only be accessed at steady state. For example, the noise

(19) |

where , is related to the long time limit of the second cumulant Esposito et al. (2009). We can therefore write:

(20) | ||||

(21) |

From the first and second cumulants, one can define the time-dependent Fano factor for the population Kambly and Flindt (2013)

(22) |

In the long-time limit, this converges to the Fano factor (up to a scaling factor of ) Beenakker et al. (2003); Emary et al. (2007) for the current,

(23) |

This is the Fano factor usually observed in steady state transport experiments, where is the charge of an individual carrier and is the Poisson value of the shot noise Blanter and Büttiker (2000); Safonov et al. (2003); Beenakker et al. (2003).

### iv.2 FCS with inchworm Monte Carlo

Within the modified Hamiltonian picture of Sec. III, the problem of evaluating a generating function can be mapped onto that of evaluating the time dependence of a unit operator propagated by the modified Hamiltonian. For the AIM with as the counting field, using the Baker–Hausdorff lemma and the fact that commutes with one can show Tang and Wang (2014) that the modified Hamiltonian of Eq. (10) is equivalent to the AIM, but with lead–molecule coupling terms that acquire a branch-dependent phase:

(24) | ||||

(25) |

The generating function in Eq. (11) may then be evaluated by a standard Keldysh hybridization expansion Werner et al. (2006); Mühlbacher and Rabani (2008); Werner et al. (2009). Whereas usually in such expansions an operator corresponding to some observable is applied at the fold of the Keldysh contour, here no such operator is needed, and the full counting statistics are given directly in terms of a modified propagator:

(26) |

In this expression, is analogous to the dressed restricted propagators used in propagator noncrossing approximations Eckstein and Werner (2010); Härtle et al. (2015); Chen et al. (2017a), in bold-line QMC Gull et al. (2010, 2011a); Cohen et al. (2013, 2014a, 2014b), and in inchworm QMC Cohen et al. (2015); Antipov et al. (2017); Chen et al. (2017a, b); Dong et al. (2017), but modified to include the auxiliary field via Eq. (25):

(27) | ||||

When is expressed in the interaction picture and expanded in powers of the hybridization, it becomes a weighted sum over configurations , defined by the set of contour times at which hybridization events (including dot and lead operators) occur. Each configuration sampled at a given diagram order corresponds to the insertion of hybridization lines connecting “kink” times on the Keldysh contour where a particle of given spin is added to or removed from the dot. We illustrate a diagram of order and another of order in Fig. 2, where the full propagator is represented by a thick black line, the bare atomic propagator by a thin black line and paired hybridization events by wiggly lines.

The full propagator can then be rewritten as a sum over diagrams, which can be stochastically sampled by diagrammatic Monte Carlo algorithms Gull et al. (2011b):

(28) |

Here the first summation is over the expansion order and the second summation is over all configurations corresponding to the kink times at a given order. The are products of interacting (but purely local) atomic state propagators:

(29) | ||||

Their sign depends on the number of kinks on each branch of the contour, and . The denote a determinant given by a sum over all possible lead operator pairings . Each term in this sum is given by the product of hybridization functions at the corresponding set of time pairs, with the sign given by the permutation generating that pairing Gull et al. (2011b):

(30) |

These hybridization functions, which are dressed by the counting field in a contour-time dependent manner, must be specified on the and axes:

(31) | ||||

Here, is the Heaviside function ordering times along the Keldysh contour; are the branch indices of and , respectively; and is the lead for which FCS is being evaluated. Due to the modified coupling in Eq. (25), the phase factors in the dressed hybridization function Eq. (31) cancel when and lie on the same branch, such that only hybridization lines that cross the folding time of the Keldysh contour are modified by the counting field (see the red elements in Fig. 2). The undressed lesser and greater hybridization components,

(32) | ||||

can be expressed in terms of the level width function as defined in Eq. (17), and are parametrized only by the physical times and corresponding to the contour times and .

Exact numerical approaches to the investigation of dynamics in quantum many-body systems typically suffer from an exponential scaling of computational cost with time. In real time Monte Carlo methods, this manifests in the dynamical sign problem: an exponential growth of stochastic errors with increasing time. However, the inchworm algorithm, as applied to the real time hybridization expansion Werner et al. (2006); Mühlbacher and Rabani (2008); Werner et al. (2009); Schiró (2010); Gull et al. (2011a, b), has been shown to bypass the dynamical sign problem in a wide variety of parameters. This was not only shown for the AIM Cohen et al. (2015); Antipov et al. (2017), but also in the spin–boson model Chen et al. (2017a, b) and—within the dynamical mean field approximation—for lattice models Dong et al. (2017).

Here, the inchworm algorithm is used to efficiently sum the expansion described above, by reusing data obtained from propagators on shorter time intervals to construct propagators on longer time intervals Cohen et al. (2015). In the top panel of Fig. 3, a snapshot of an intermediate step in the inchworm procedure is illustrated for a branch-independent Hamiltonian without a counting field. Each of the and axes shown consists of an “unfolded” Keldysh contour containing times ordered according to the structure . In each inchworm step, propagators are extended along one time direction (this is called “inching”), marked by the green arrows going to the right. In dark gray, we show time points unneeded for the calculation (as we can always select one contour time to appear later on the contour than the other). The white squares correspond to equal time propagators, which can be evaluated analytically. The blue squares correspond to time arguments for which the propagator is already known from previous steps, while the light red squares correspond to arguments at which the propagator may be evaluated immediately, given the currently known propagators. The dark red squares represent values to be computed at a later steps. In light gray, we show time points which can be obtained from reflection about the axis (marked by a white dashed diagonal line), which can be seen as time reversal symmetry: the propagators on one side of this symmetry element can be obtained from those on the other via the relation

(33) |

Here, a contour time argument marked by a tilde, , denotes the same physical time as that of , but on the opposite contour.

Due to time reversal symmetry, in the original formulation of the inchworm algorithm the total number of propagators computed is , as only one quadrant of the two time plane needs to be explicitly evaluated. Extending a propagator to the right in the quadrant used is identical to extending a symmetrically placed propagator up in the mirrored quadrant, so that time reversal symmetry is also maintained in the inching direction.

In diagrammatic terms, time reversal symmetry is illustrated by in Fig. 2. In the upper part of 2, the cross-branch hybridization line corresponds to a factor of in the propagator, which depends only on physical times and not on the branch indices. When the contour branches of all times appearing in this diagram are flipped, one obtains the expansion in the lower part of Fig. 2. For example, the hybridization factor is replaced by . All other factors are similarly conjugated, and so Eq. (33) holds.

By contrast, in the presence of the counting field we evaluate the propagator for nonzero , and the two contour-reflected diagrams are modified as illustrated by the red elements in Fig. 2: the dressed cross-branch hybridization factor, denoted by the red wiggly line, is replaced by (upper diagram) and (lower diagram). Due to the unconjugated prefactor , Eq. (33) is no longer valid, and the two quadrants to either side of the line must be evaluated separately, as illustrated in the lower panel of Fig. 3. Since we no longer mirror the data or enforce the symmetry, an additional spurious numerical breaking of time reversal symmetry due to the asymmetry in inching direction can occur. To avoid this, we inch to each time point from the two possible directions and set the result to the average (illustrated by the green arrows, now going both right and up). This increases the number of simulations by an additional factor of two, so that are needed in total, but the overall quadratic scaling with the time step remains. We have found that bidirectional inching generally provides higher quality data for a given amount of computer time, even without the counting field.

## V Benchmark and validation for the non-interacting system

### v.1 Noninteracting FCS

In the absence of interactions (when ), the FCS following the coupling quench can be exactly obtained by the path integral nonequilibrium Green’s function (PI-NEGF) method, which provides an exact solution in the absence of electron–electron interactions Tang and Wang (2014). The generating function for a noninteracting AIM is given by the ratio of two Fredholm determinants:

(34) |

Here, is the matrix of two-time Green’s functions for the full molecular junction in Keldysh-rotated space Keldysh (1965); Kamenev (2011), such that the times correspond to pairs of times on the Keldysh contour with a folding point at . The Green’s function is given by

(35) |

where denotes the Green’s function of the isolated dot and is a modified Green’s function that depends upon the counting field in the left lead. When expressed in the Keldysh-rotated representation, the two hybridizations are related via

(36) |

where is the Pauli spin matrix. It is easy to show that this is equivalent to Eq. (31) when the Keldysh rotation is reversed. Further details on the implementation of Eq. (34) can be found in Refs. Tang and Wang, 2014; Tang et al., 2014; Yu et al., 2016. It was also shown in these works that the expression (34) correctly reduces to the Levitov–Lesovik formula for the generating function in the long-time limit Levitov and Lesovik (1993); Levitov et al. (1996), which can be written as

(37) | ||||

Here, denotes the Fermi function of lead and is the transmission probability for electrons to pass through the molecule.

In this work, the PI-NEGF generating function is computed for the noninteracting case without spin degeneracy, as in Ref. Tang and Wang, 2014. However, we are interested in the AIM. The generating function for the AIM can be constructed from that of the single-electron model as a product:

(38) |

The , subscripts denote the two spins. In the cases to be shown here, the system is initially decoupled from the leads with the dot in the unoccupied state , such that in the noninteracting calculations, each independent model begins with an unoccupied dot. The dot–bath coupling is then turned on at the quench time .

### v.2 QMC results for the non-interacting system

To verify the correctness of the iQMC method, we begin by performing calculations for the noninteracting case and comparing the results with exact results obtained from the PI-NEGF method. This formalism provides us with an independent verification of the new method: as the iQMC calculation is based on a hybridization (rather than an interaction) expansion, the noninteracting case is nontrivial and embodies a rigorous test of the method’s accuracy Werner et al. (2006). In Fig. 4 we present the time evolution of . The two panels on the left are generated from iQMC, while the two panels on the right are exact results. As might be expected, is satisfied to within the numerical errors. The following symmetry relations are also apparent:

(39) | ||||

We note that the time evolution of both real and imaginary parts of exhibit monotonic decay towards zero after a transient timescale on the order of .

In order to explore the correspondence between the numerical data and the exact result in more detail, we plot a series of cuts across the data in Fig 5. Each curve shown here corresponds to a different value of in the interval . The results from iQMC are shown in solid red, together with exact data in dashed black. The two sets of curves appear to be in very good agreement, with slight numerical deviations (on the order of 1%) in the iQMC results visible at long times.

The introduction of a finite voltage causes the components of the generating function to strongly oscillate as a function of time, as shown in Fig. 6. This indicates oscillations in cumulants of all orders, reflecting a universal phenomenon predicted by Flindt et al. in Ref. Flindt et al., 2009 and later observed experimentally Fricke et al. (2010). Physically, the introduction of a finite bias enhances the short-time “ringing” behavior, as has previously been observed for the lowest cumulant (the current) Härtle et al. (2015); Ridley et al. (2015). As the FCS in the presence of a voltage bias is substantially richer, it is instructive to consider a detailed comparison of specific cuts again (see Fig. 7). Remarkably, it is observable from this figure (which used less computational resources than Fig 5) that at higher voltages it is easier to converge the numerics. This is a notable property of real time iQMC Cohen and Rabani (2011); Cohen et al. (2013), which makes its regime of efficient applicability different from that of, e.g., low-energy wavefunction methods.

Knowledge of the generating functional provides access to all moments and cumulants in principle, but the numerical exercise of obtaining them in practice may not be trivial. One of the main technical difficulties stems from the fact that within iQMC, the relative error in this data is generally larger at small values. This is for two reasons: first, part of the statistical iQMC error is independent of the value of the observable. This source of error is strongly mitigated by the inchworm algorithm, but not completely removed, and results in noisy numerical derivatives. Second, there is a small systematic error due to the time discretization. This error is difficult to fully converge, and can sometimes be significant: for example, consider how the imaginary part of visibly deviates from its exact value of zero in the lower panel of Fig. 5.

Therefore, if we evaluate derivatives of in the limit of small by a finite difference formula , the value of cannot easily be taken to zero without investing unjustified amounts of computer time. Furthermore, the oscillatory form of the function hints that an increasingly small value will be needed as one propagates to longer times.

However, the dependence of the moments and cumulants on can be understood at long times by considering the Levitov–Lesovik formula for the current cumulant generating function , with given by Eq. (37). If this is inserted into the finite difference derivative with respect to a small shift , and the integrand is expanded to linear order in , one arrives at the following expressions for the logarithmic and direct derivatives:

(40) | ||||

(41) | ||||

Thus, for finite but small the first cumulant (which corresponds to Eq. (20)) exhibits no dependence on , whereas the first moment oscillates at long times as . This implies that numerical derivatives taken from will converge at a finite in the long time limit, while the error in numerical derivatives taken directly from the generating function will diverge.

To see this in practice, consider Fig. 8. Here, dot populations normalized by time () are plotted. These plots were obtained from both the logarithmic (thick black lines) and direct (dashed black lines) derivatives with a symmetric bias of . Cumulants are obtained from both derivatives for and compared with the exact result (red crosses) obtained from PI-NEGF at . This is done for high temperatures (, upper panel) and low temperatures (, lower panel). As suggested by the preceding analysis, the convergence to the limit is substantially faster when logarithmic derivatives are taken. While the direct and logarithmic derivatives are in agreement at short times, the direct derivative diverges from the exact result at long times. Without performing a full error analysis, it is difficult to determine whether the logarithmic derivative differs significantly from the exact result. At low temperatures convergence is slower, and the steady state value (equivalent to the long time limit of the left current) increases.

We repeat this analysis for the normalized second cumulants of the population in Fig. 9. The second cumulant increases linearly with time Flindt et al. (2009); Esposito et al. (2009), and is related at long times to the current noise in the left lead via Eq. (21). The cumulant extracted from the iQMC data is in excellent agreement with the exact result event for relatively large values of when the cumulant is obtained from the logarithmic derivative, but requires convergence to decreasing at smaller times when the direct derivative is used. We note that due to the symmetry of Eq. (39), it is in principle only necessary to evaluate at one value of to evaluate both the first and second cumulants. In practice we used and to eliminate some of the systematic errors due to the finite time discretization.

## Vi Results in the presence of interactions

We continue to discuss quantities of physical interest, as obtained from our calculations. In what follows, we will show noninteracting results for comparison, which could be compared to the PI-NEGF data as in the previous section. Such comparisons were performed and agreement was found, but for the sake of brevity this is no longer shown from here onwards. Furthermore, we do not perform the full (and costly) error analysis needed to rigorously define confidence intervals in the inchworm algorithm Cohen et al. (2015); Antipov et al. (2017), but roughly estimate that numerical errors are on the order of several percent.

### vi.1 FCS and Shot noise

The FCS in the presence of interactions is shown at two voltages in Fig. 1, at equilibrium and in the presence of a bias voltage. At first glance, this looks only subtly different from the noninteracting results in Figs. 4) and 6). However, these seemingly small changes in the contour plot encode entirely different physics. To see this, we will explore some of the properties that can be derived from the FCS.

First, we study the effect of interactions and temperature on the noise and Fano factor. Fig. 10 shows the evolution of the time-normalized second moment in the upper panel. In the lower panel, the population Fano factor defined in Eq. (22) is shown, at fixed voltage and for the initially empty dot state. The high (low) temperature plots are shown in red (black), and the interacting (noninteracting) results are distinguished by solid (dashed) lines. To the right extreme of the plot, the asymptotic values of the noise and Fano factor are shown. These were obtained from the coefficient of a linear fit to the first and second cumulants and , in accordance with the definition of the steady state current and noise defined in Eqs. (20) and (21). In addition, filled circles mark the Levitov–Lesovik values for the noninteracting quantum noise and Fano factor, obtained from the standard Landauer–Büttiker theory Blanter and Büttiker (2000); Ridley et al. (2017) via Eqs. (37) and (38). We note that the linear fits (which are not shown here) give reliable values to within 1% in this case, and in the noninteracting case are perfectly consistent with the Levitov–Lesovik values. The quality of the fits improves with the length of time simulated. In addition we note that asymptotic values shown to the right of Fig. 10 are independent of the initial condition.

At high temperatures (red lines), the absolute magnitude of the noise, and also the Fano factor, is enhanced with respect to the low temperature (black lines) case. This is due to the contribution to the second cumulant of thermal or Nyquist noise, which vanishes at small values of Blanter and Büttiker (2000); Ridley et al. (2017). We note that at high temperatures there is an asymptotic divergence between the solid and dashed red curves in the upper panel, as the noise appears to be suppressed by the presence of interactions. In the low temperature case, the interacting shot noise is enhanced with respect to the noninteracting value, thus we observe a crossover from noise suppression to enhancement as the temperature is reduced. At low temperature the quantum noise is mainly attributed to the presence of a finite bias Khlus (1987); Blanter and Büttiker (2000). While the thermal-to-shot noise crossover is well described by scattering theory in the noninteracting case, it is far more complex in the presence of strong interactions.

The low-temperature data corroborates well-known theoretical results by Lesovik Lesovik (1989) and others Büttiker (1990); Blanter and Büttiker (1999); Vitushinsky et al. (2008), later confirmed by experiment Iannaccone et al. (1998), that the presence of inelastic scattering processes in the junction causes excess noise at low temperatures and finite voltages. In particular, inelastic cotunneling processes have been considered, in which electrons tunnel onto and off the dot simultaneously with the creation of a virtual dot state. These leave the dot in an excited state, and have been shown to enhance the noise Sukhorukov et al. (2001); Thielmann et al. (2005); Kaasbjerg and Belzig (2015). When these virtual tunnel states involve correlated electrons of opposite spin, a Kondo singlet is formed Cronenwett et al. (1998). At high temperatures, the interaction induces Coulomb blockade on the dot, reducing the number of available transport channels. This suppresses both the current and the noise, such that the Fano factor takes a similar value to that in the noninteracting case (lower panel, red lines). The presence of a large temperature suppresses the formation of coherent states required for inelastic spin-flip and cotunneling processes to occur Furusaki and Matveev (1995); König et al. (1998); Tran et al. (2008). However, as the temperature is lowered, we enter the Kondo regime. The Kondo effect increases the inelastic spin-flip rate and facilitates spin fluctuations on the dot, enhancing the noise and resulting in an overall increase in the Fano factor compared to the noninteracting case (solid black lines in Fig. 10).

### vi.2 First passage time distribution

In Figs. 11–12 we plot the FPTD, as defined in Eq. (6), in the unbiased () and biased () cases. The inverse temperature in all plots is , and the dynamics for four different initial states of the dot are shown. Note that the initially half-occupied states are collected into the same line (shown in red), as the particle–hole symmetric parameters ensure their physical equivalence. Within our data, the two half-filled initial conditions are indeed identical to within numerical errors (not shown).

We begin with the noninteracting problem, Fig. 11. Here, a maximum occurs at of order , corresponding to the most probable first passage time , i.e. the time at which it is most likely to measure the first change in the left lead’s occupation. The first passage probability decays to zero at long times, as the likelihood that the first particle transfer has been detected at very long times becomes vanishingly small. In the unbiased case (top panel), the effect of initial condition on the first passage time distributions is minimal. In particular, the unoccupied (black) and fully-occupied (green) initial conditions are identical in the presence of particle–hole symmetry. This is not the case for the biased system (bottom panel), where electrons are driven out of the left lead. The bias voltage breaks the symmetry between the left and right leads. Therefore—since we are examining the FPTD of the left lead—the symmetry between the doubly occupied and unoccupied initial states is also broken. As the initial occupation of the dot increases in going from the to the and states, the probability density shifts to longer times. A second peak appears in the distribution, corresponding to the first tunneling event on the left lead occurring after an electron has traveled from the dot onto the right lead. This demonstrates an initial condition dependent “queuing” effect in the FPTD, which is reminiscent of classical queuing Lindley (1952).

The insets of Fig. 11 display the integral over , which by Eq. (5) is equal to the probability that when the total charge of the left lead is measured at time , it has changed from its initial value. In an unbiased system the electron flow across the terminals, as well as the probability that the electron count in the left lead has changed from its initial value, is due entirely to thermal fluctuations and the delocalization dynamics of the wavefunction, and may saturate at a value between one and zero (see upper inset). As the bias is increased, however, the active driving causes this probability to reach unity rapidly (as in the the lower inset).

Fig. 12 is identical to Fig. 11, except for the presence of a finite Coulomb repulsion energy . At very short times, when the time dependence is linear, the symmetrized interaction has essentially no effect on the FPTD. However, this changes dramatically at timescales substantially smaller than . Comparing the interacting unbiased case, shown in the upper panel, with the corresponding noninteracting result in Fig. 11, it is clear that the magnitude of the resonant peak and the value of are increased for the doubly occupied and unoccupied initial states and , whereas both quantities are decreased for the half occupied initial states . This can be understood in terms of Coulomb blockade physics: at short times, involving the transfer of the first electron, no difference between the initial conditions is observed. However, in the case, the suppression of the even charge states by the interaction causes the transfer of a second electron in the same direction to be energetically unfavorable. The opposite occurs in the and states.

Another striking difference between the interacting case and its noninteracting counterpart (top panels of Fig. 12 and Fig. 11, respectively) is the long tail appearing in the FPTD of the initially magnetized states . The effect of this tail is prevalent in the inset, where very slow relaxation to unity is observed (we show times up to to emphasize the slow relaxation). This can be attributed to the stabilization of local moments Hewson (1997) and slow spin dynamics Weiss (2012) associated with the Kondo regime, to which the system is equilibrating at these parameters. For comparison, the red dashed line in the upper inset of Fig. 12 shows dynamics for the initial states at , where the Kondo effect is suppressed by temperature; the relaxation is then faster.

In the lower panel of Fig. 12 we apply a bias voltage of to the interacting system. At such a large voltage, electrons are chiefly injected into the dot from the left lead and ejected into the right lead, simplifying the analysis of the effects of interaction on the queuing. In the case of the half occupied initial conditions , the first peak is suppressed relative to the second in the presence of interactions, because the energy penalties discussed above decrease the likelihood of electron injection from the left lead before the original electron escapes into the right lead. Similarly, the first peak for the doubly occupied initial condition becomes larger than the second, because the process where two electrons are ejected on the right before the first electron is injected on the left is energetically suppressed. Remarkably, the slow spin relaxation related to Kondo physics is substantially reduced in the inset for the initial states at this high voltage. This is consistent with claims that a small remnant of Kondo physics can survive at high voltages Wingreen and Meir (1994); Anders (2008); Cohen et al. (2014b), but requires further investigation.

### vi.3 -particle probabilities

The FPTD, related to , sheds light on the dynamical response of the system to all possible particle tunneling events. It does not distinguish between events in which different numbers of particles are transported. However, the full counting statistics contain much more information: it is possible to access for every . It is of interest to search for manifestations of correlated transport in these detailed distributions. We also consider the most probable time at which the particle number changes by , given by the maximum of the corresponding distribution. We denote this time by .

In Fig. 13, we begin by considering the probability for the particle number in the left lead to have changed by after time , for each possible initial condition and for the interacting (solid lines) and noninteracting (dashed lines) cases. We once again apply a bias voltage of , so that it is energetically favorable to move charges from the left lead onto the dot, and charges from the dot to the right lead. Regardless of the presence of interactions, the peak of the probability distributions for the different initial conditions is shifted to longer times by an increasing number of electrons initially on the dot. In the noninteracting case, this simply reflects the fact that the empty dot has more transport channels open to electron traversal events than a partially or fully filled dot. If the initial state is fully occupied, transport of the first electron is forbidden by the Fermi rule until at least one dot electron has tunneled to (with high probability) the right lead. The effect of Coulomb charging depends on the initial occupation: when the dot begins empty, interaction does not change the position of the maximum, but somewhat enhances the total probability that a single tunneling event was measured at any time, reflecting a reduction in the probability of detecting higher values of . In this case there is no initial “queue” for electron tunneling onto the dot, and the interaction makes the first electron tunneling event more energetically favorable, then suppresses the second. For all other initial states, the Coulomb repulsion shifts the distribution to longer times, as electrons on the left lead cannot enter the dot without “queuing” for electrons on the dot to first tunnel to the right lead.

In the lower panel of Fig. 13, we repeat the above analysis for the probability that at time , a change in the particle number is measured on the left lead. is shifted to longer times by the interaction for all initial conditions, because the transfer of the second electron must be preceded by that of the first, such that it always encounters some Coulomb repulsion.

The top panel of Fig. 14 shows a sequence of for , focusing on the unoccupied initial state. At short times, we see perfect agreement between the interacting and noninteracting results, as no interaction can take place before some electrons occupy the dot. Increasing the value of shifts the peaks to longer times, which makes sense physically since a change involving electrons implies all processes leading to a change of electrons have already occurred. We note that we do not show backscattering events , which have a small but finite probability even in the presence of a large voltage. In the unbiased system (not shown), backscattering events are of equal importance to the forward scattering events.

The rich, detailed information on population densities shown in Fig. 14 lends itself to a quantitative analysis of electron transfer processes. In the lower panel of the same figure, we study the variation of the distance between maximally probable times as a function of , for the interacting (black curve) and noninteracting (red curve) cases. As we propagate to longer times, this quantity can be considered a proxy for queuing effects at the single electron level at steady state: it describes the typical waiting time between subsequent tunneling events. As increases, the peak-to-peak distance stabilizes to a roughly constant value. The long-time value in the noninteracting case is of order , as might be expected for a simple rate process where the only relevant timescale is the coupling between the dot and baths. However, the peak-to-peak distance in the interacting case is significantly larger: as each electron enters the dot, it suppresses the next electron from entering by virtue of the Coulomb repulsion. Interestingly, while the point at which each distribution begins to differ significantly from zero is only weakly modified by the interaction, the width of the distributions and the weight at their tail end are almost immediately enhanced. This may indicate that fluctuations play an increasingly important role in the population transfer dynamics as the strength of the interaction is increased.

## Vii Conclusions

We presented the first numerically exact calculation of full counting statistics for a non-integrable model of interacting fermions, in this case the nonequilibrium Anderson impurity model. Using the inchworm quantum Monte Carlo method, we obtained the generating function at a variety of physical parameters ranging from the noninteracting case to the strongly correlated Kondo regime. This provides access to currents, which have been accessed before by numerically exact methods; but also to other experimentally measurable quantities which were not. This includes the current noise at steady state, the Fano factor, all higher moments and cumulants of population transfer event, and the complete time-dependent probability distributions for -electron transfer events.

The method was applied to a coupling quench, where a dot is suddenly attached to the leads at time zero and allowed to evolve to equilibrium or to a nonequilibrium steady state. After performing benchmark comparisons in the noninteracting limit to verify the accuracy of our results, we explored the effects of electron–electron interactions and nonequilibrium bias voltage on the full counting statistics and the properties derived thereof. We observed the signatures of the Coulomb blockade and Kondo effects in the noise and Fano factor, and found transient queuing effects depending on the choice of initial condition in the FPTD. By considering the individual, time-dependent probability distributions for the occurrence of tunneling events, we were further able to argue that such queuing effects persist at steady state. Rather than simply being associated with a lower effective tunneling rate, the rapid widening of the probability densities in the presence of interactions suggests that fluctuations play a significant role in the interacting dynamics.

With the availability of a reliable scheme for calculating the full counting statistics of generic interacting impurity models, a variety of important research questions can now be addressed. Some immediately relevant examples include quantum thermodynamic topics, such as the verification of quantum fluctuation–dissipation theorems and the calculation of efficiency fluctuations in quantum devices; and the evaluation of noise and Fano factors beyond the Fermi liquid regime. We observed queuing effects, but it would also be interesting to look for bunching effects in the presence of an effectively attractive interaction. As the noninteracting physics enters entirely through the hybridization function, it is also feasible to consider the counting statistics of impurities embedded in more realistic noninteracting models of materials and nanosystems, one interesting example being magnetic impurities in graphene nanoribbons da Rocha et al. (2015); Ridley and Tuovinen (2017).

Spin-dependent FCS in quantum junctions has only recently been treated in the noninteracting case Gaomin et al. (2017), and spintronic applications are anticipated. The method we have presented can easily be generalized to access spin-dependent counting statistics and multi-lead moments: one could then go beyond the analysis of Fig. 14 to consider the time-dependent probability that the left lead lost (e.g.) one spin up electron, while simultaneously the right lead gained two spin down electrons. In this context, the long-term survival of spin fluctuations at timescales where the mean magnetization has died out would embody a remarkably clear signature of strong correlation physics. Generalization to include bosonic degrees of freedom is also possible, enabling treatment of thermoelectric systems or coupling to an optical continuum. The method could further be extended to frequency dependent power spectra and related conductance spectra Rothstein et al. (2009); Gabdank et al. (2011), which can be used to study photon absorption and emission processes Engel and Loss (2004); Entin-Wohlman et al. (2007); Orth et al. (2012); Miwa et al. (2017). Additionally, noise in periodically-driven systems with interactions could now be investigated, enabling performance tuning of correlated nanoelectronic devices Strass et al. (2005).

To summarize, the full counting statistics of strongly correlated nonequilibrium quantum impurity problems can now be obtained numerically. This provides unprecedented insight into the stochastic nature of electronic transport at the resolution of single tunneling events, and allows simulation of a variety of experiments for which no systematic theory was previously available.

#### Acknowledgements

We are grateful to Michael Galperin for directing our attention to this problem. G.C. acknowledges support by the Israel Science Foundation (Grant No. 1604/16). M.R. was supported by the Raymond and Beverly Sackler Center for Computational Molecular and Materials Science, Tel Aviv University. E.G. was funded by DOE ER 46932. This research was supported by grant No. 2016087 from the United States-Israel Binational Science Foundation (BSF).

## References

- Blanter and Büttiker (2000) Y. M. Blanter and M. Büttiker, Physics reports 336, 1 (2000).
- Kampen (2007) N. V. Kampen, Stochastic processes in physics and chemistry (North Holland, 2007).
- Braggio et al. (2006) A. Braggio, J. König, and R. Fazio, Physical review letters 96, 026805 (2006).
- Landauer (1998) R. Landauer, Nature 392, 658 (1998).
- Beenakker et al. (2003) C. Beenakker, C. Schönenberger, et al., Physics Today 56, 37 (2003).
- Vardimon et al. (2013) R. Vardimon, M. Klionsky, and O. Tal, Physical Review B 88, 161404 (2013).
- Selzer and Peskin (2013) Y. Selzer and U. Peskin, The Journal of Physical Chemistry C 117, 22369 (2013).
- Yang et al. (2014) P.-Y. Yang, C.-Y. Lin, and W.-M. Zhang, Physical Review B 89, 115411 (2014).
- Ochoa et al. (2015) M. A. Ochoa, Y. Selzer, U. Peskin, and M. Galperin, The journal of physical chemistry letters 6, 470 (2015).
- Ridley et al. (2017) M. Ridley, A. MacKinnon, and L. Kantorovich, Physical Review B 95, 165440 (2017).
- Rychkov et al. (2005) V. S. Rychkov, M. L. Polianski, and M. Büttiker, Physical Review B 72, 155326 (2005).
- Battista et al. (2014) F. Battista, F. Haupt, and J. Splettstoesser, Physical Review B 90, 085418 (2014).
- Zamoum et al. (2016) R. Zamoum, M. Lavagna, and A. Crépieux, Physical Review B 93, 235449 (2016).
- Levitov et al. (1996) L. S. Levitov, H. Lee, and G. B. Lesovik, Journal of Mathematical Physics 37, 4845 (1996).
- Dubois et al. (2013a) J. Dubois, T. Jullien, C. Grenier, P. Degiovanni, P. Roulleau, and D. Glattli, Physical Review B 88, 085301 (2013a).
- Dubois et al. (2013b) J. Dubois, T. Jullien, F. Portier, P. Roche, A. Cavanna, Y. Jin, W. Wegscheider, P. Roulleau, and D. Glattli, Nature 502, 659 (2013b).
- Khlus (1987) V. Khlus, Zh. Eksp. Teor. Fiz 93, 2179 (1987).
- Lefloch et al. (2003) F. Lefloch, C. Hoffmann, M. Sanquer, and D. Quirion, Physical review letters 90, 067002 (2003).
- Braggio et al. (2011) A. Braggio, M. Governale, M. G. Pala, and J. König, Solid State Communications 151, 155 (2011).
- Saminadayar et al. (1997) L. Saminadayar, D. Glattli, Y. Jin, and B. c.-m. Etienne, Physical Review Letters 79, 2526 (1997).
- Sela et al. (2006) E. Sela, Y. Oreg, F. von Oppen, and J. Koch, Physical review letters 97, 086601 (2006).
- Golub (2006) A. Golub, Physical Review B 73, 233310 (2006).
- Golub (2007) A. Golub, Physical Review B 75, 155313 (2007).
- Zarchin et al. (2008) O. Zarchin, M. Zaffalon, M. Heiblum, D. Mahalu, and V. Umansky, Physical Review B 77, 241303 (2008).
- Yamauchi et al. (2011) Y. Yamauchi, K. Sekiguchi, K. Chida, T. Arakawa, S. Nakamura, K. Kobayashi, T. Ono, T. Fujii, and R. Sakano, Physical review letters 106, 176601 (2011).
- Ferrier et al. (2016) M. Ferrier, T. Arakawa, T. Hata, R. Fujiwara, R. Delagrange, R. Weil, R. Deblock, R. Sakano, A. Oguri, and K. Kobayashi, Nature Physics 12, 230 (2016).
- Levitov and Reznikov (2004) L. Levitov and M. Reznikov, Physical Review B 70, 115305 (2004).
- Bomze et al. (2005) Y. Bomze, G. Gershon, D. Shovkun, L. Levitov, and M. Reznikov, Physical review letters 95, 176601 (2005).
- Safonov et al. (2003) S. Safonov, A. Savchenko, D. Bagrets, O. Jouravlev, Y. V. Nazarov, E. Linfield, and D. Ritchie, Physical review letters 91, 136801 (2003).
- Djuric et al. (2005) I. Djuric, B. Dong, and H. Cui, Applied Physics Letters 87, 032105 (2005).
- Kießlich et al. (2007) G. Kießlich, E. Schöll, T. Brandes, F. Hohls, and R. Haug, Physical review letters 99, 206602 (2007).
- Feng et al. (2008) Z. Feng, J. Maciejko, J. Wang, and H. Guo, Physical Review B 77, 1 (2008).
- Gustavsson et al. (2006) S. Gustavsson, R. Leturcq, B. Simovič, R. Schleser, T. Ihn, P. Studerus, K. Ensslin, D. Driscoll, and A. Gossard, Physical review letters 96, 076605 (2006).
- Emary et al. (2007) C. Emary, D. Marcos, R. Aguado, and T. Brandes, Physical Review B 76, 161404 (2007).
- Tworzydło et al. (2006) J. Tworzydło, B. Trauzettel, M. Titov, A. Rycerz, and C. W. Beenakker, Physical Review Letters 96, 246802 (2006).
- Levitov and Lesovik (1993) L. S. Levitov and G. B. Lesovik, Pis’ma Zh. Exsp. Teor. Fiz. 58, 225 (1993).
- Klich and Levitov (2009) I. Klich and L. Levitov, Physical review letters 102, 100502 (2009).
- Lesovik et al. (2006) G. Lesovik, F. Hassler, and G. Blatter, Physical review letters 96, 106801 (2006).
- Esposito et al. (2015) M. Esposito, M. A. Ochoa, and M. Galperin, Physical Review B 91, 115417 (2015).
- Agarwalla et al. (2015) B. K. Agarwalla, J.-H. Jiang, and D. Segal, Physical Review B 92, 245418 (2015).
- Di Lorenzo and Nazarov (2004) A. Di Lorenzo and Y. V. Nazarov, Physical review letters 93, 046601 (2004).
- Di Lorenzo and Nazarov (2005) A. Di Lorenzo and Y. V. Nazarov, Physical review letters 94, 210601 (2005).
- Di Lorenzo et al. (2006) A. Di Lorenzo, G. Campagnano, and Y. V. Nazarov, Physical Review B 73, 125311 (2006).
- Fransson and Galperin (2010) J. Fransson and M. Galperin, Physical Review B 81, 075311 (2010).
- Tang and Wang (2014) G.-M. Tang and J. Wang, Physical Review B 90, 195422 (2014).
- Gaomin et al. (2017) T. Gaomin, X. Fuming, M. Shuo, and W. Jian, arXiv preprint arXiv:1712.00215 (2017).
- Tang et al. (2017a) G. Tang, X. Chen, J. Ren, and J. Wang, arXiv preprint arXiv:1705.10025 (2017a).
- Schmidt et al. (2007) T. Schmidt, A. Komnik, and A. Gogolin, Physical Review B 76, 241307 (2007).
- Esposito et al. (2009) M. Esposito, U. Harbola, and S. Mukamel, Reviews of modern physics 81, 1665 (2009).
- Sakano et al. (2011a) R. Sakano, A. Oguri, T. Kato, and S. Tarucha, Physical Review B 83, 241301 (2011a).
- Flindt et al. (2008) C. Flindt, T. Novotnỳ, A. Braggio, M. Sassetti, and A.-P. Jauho, Physical review letters 100, 150601 (2008).
- Ubbelohde et al. (2012) N. Ubbelohde, C. Fricke, C. Flindt, F. Hohls, and R. J. Haug, Nature communications 3, 612 (2012).
- Tobiska and Nazarov (2005) J. Tobiska and Y. V. Nazarov, Physical Review B 72, 235328 (2005).
- Saito and Utsumi (2008) K. Saito and Y. Utsumi, Physical Review B 78, 115429 (2008).
- Simine and Segal (2012) L. Simine and D. Segal, Physical Chemistry Chemical Physics 14, 13820 (2012).
- Pekola (2015) J. P. Pekola, Nature Physics 11, 118 (2015).
- Stegmann and König (2017) P. Stegmann and J. König, physica status solidi (b) 254 (2017).
- Altland et al. (2010) A. Altland, A. De Martino, R. Egger, and B. Narozhny, Physical Review B 82, 115323 (2010).
- Abanov and Ivanov (2008) A. Abanov and D. Ivanov, Physical review letters 100, 086602 (2008).
- Abanov and Ivanov (2009) A. G. Abanov and D. Ivanov, Physical Review B 79, 205315 (2009).
- Utsumi et al. (2013) Y. Utsumi, O. Entin-Wohlman, A. Ueda, and A. Aharony, Physical Review B 87, 115407 (2013).
- Kambly et al. (2011) D. Kambly, C. Flindt, and M. Büttiker, Physical Review B 83, 075432 (2011).
- Flindt et al. (2009) C. Flindt, C. Fricke, F. Hohls, T. Novotnỳ, K. Netočnỳ, T. Brandes, and R. J. Haug, Proceedings of the National Academy of Sciences 106, 10116 (2009).
- Fricke et al. (2010) C. Fricke, F. Hohls, C. Flindt, and R. J. Haug, Physica E: Low-dimensional Systems and Nanostructures 42, 848 (2010).
- Albert et al. (2011) M. Albert, C. Flindt, and M. Büttiker, Physical review letters 107, 086805 (2011).
- Albert et al. (2012) M. Albert, G. Haack, C. Flindt, and M. Büttiker, Physical review letters 108, 186806 (2012).
- Thomas and Flindt (2014) K. H. Thomas and C. Flindt, Physical Review B 89, 245420 (2014).
- Brandes (2008) T. Brandes, Annalen der Physik 17, 477 (2008).
- Wang et al. (2007) S.-K. Wang, H. Jiao, F. Li, X.-Q. Li, and Y. Yan, Physical Review B 76, 125416 (2007).
- Welack et al. (2008) S. Welack, M. Esposito, U. Harbola, and S. Mukamel, Physical Review B 77, 195315 (2008).
- Schaller et al. (2009) G. Schaller, G. Kießlich, and T. Brandes, Physical Review B 80, 245107 (2009).
- Welack et al. (2009) S. Welack, S. Mukamel, and Y. Yan, EPL (Europhysics Letters) 85, 57008 (2009).
- Rajabi et al. (2013) L. Rajabi, C. Pöltl, and M. Governale, Physical review letters 111, 067002 (2013).
- Droste and Governale (2016) S. Droste and M. Governale, Journal of Physics: Condensed Matter 28, 145302 (2016).
- Albert et al. (2016) M. Albert, D. Chevallier, and P. Devillard, Physica E: Low-dimensional Systems and Nanostructures 76, 209 (2016).
- Elzerman et al. (2004) J. Elzerman, R. Hanson, L. W. Van Beveren, B. Witkamp, L. Vandersypen, and L. P. Kouwenhoven, Nature 430, 431 (2004).
- Tang et al. (2014) G.-M. Tang, F. Xu, and J. Wang, Physical Review B 89, 205310 (2014).
- Yu et al. (2016) Z. Yu, G.-M. Tang, and J. Wang, Physical Review B 93, 195419 (2016).
- Tang et al. (2017b) G. Tang, Y. Xing, and J. Wang, Physical Review B 96, 075417 (2017b).
- Souto et al. (2015) R. S. Souto, R. Avriller, R. Monreal, A. Martín-Rodero, and A. L. Yeyati, Physical Review B 92, 125435 (2015).
- Tang et al. (2017c) G. Tang, Z. Yu, and J. Wang, New Journal of Physics 19, 083007 (2017c).
- Komnik and Gogolin (2005) A. Komnik and A. Gogolin, Physical review letters 94, 216601 (2005).
- Gogolin and Komnik (2006a) A. Gogolin and A. Komnik, Physical Review B 73, 195301 (2006a).
- Gogolin and Komnik (2006b) A. Gogolin and A. Komnik, Physical review letters 97, 016602 (2006b).
- Gogolin et al. (2007) A. Gogolin, R. Konik, A. W. Ludwig, and H. Saleur, Annalen der Physik 16, 678 (2007).
- Mora et al. (2008) C. Mora, X. Leyronas, and N. Regnault, Physical review letters 100, 036604 (2008).
- Mora et al. (2009) C. Mora, P. Vitushinsky, X. Leyronas, A. A. Clerk, and K. Le Hur, Physical Review B 80, 155322 (2009).
- Sakano et al. (2011b) R. Sakano, T. Fujii, and A. Oguri, Physical Review B 83, 075440 (2011b).
- Oguri and Sakano (2013) A. Oguri and R. Sakano, Journal of the Korean Physical Society 63, 423 (2013).
- Sakano et al. (2012) R. Sakano, Y. Nishikawa, A. Oguri, A. C. Hewson, and S. Tarucha, Physical review letters 108, 266401 (2012).
- Dong et al. (2013) B. Dong, G. Ding, and X. Lei, Physical Review B 88, 075414 (2013).
- Schmidt and Komnik (2009) T. Schmidt and A. Komnik, Physical Review B 80, 041307 (2009).
- Park and Galperin (2011) T.-H. Park and M. Galperin, Physical Review B 84, 205450 (2011).
- Maier et al. (2011) S. Maier, T. Schmidt, and A. Komnik, Physical Review B 83, 085401 (2011).
- Bagrets and Nazarov (2003) D. Bagrets and Y. V. Nazarov, Physical Review B 67, 085316 (2003).
- Flindt et al. (2010) C. Flindt, T. Novotnỳ, A. Braggio, and A.-P. Jauho, Physical Review B 82, 155407 (2010).
- Kaasbjerg and Belzig (2015) K. Kaasbjerg and W. Belzig, Physical Review B 91, 235413 (2015).
- Benito et al. (2016) M. Benito, M. Niklas, and S. Kohler, Physical Review B 94, 195433 (2016).
- Luo et al. (2017) J. Luo, Y. Yan, Y. Huang, L. Yu, X.-L. He, and H. Jiao, Physical Review B 95, 035154 (2017).
- Kosov (2017) D. S. Kosov, The Journal of Chemical Physics 146, 074102 (2017).
- Wang et al. (2017) C. Wang, J. Ren, and J. Cao, Physical Review A 95, 023610 (2017).
- Cerrillo et al. (2016) J. Cerrillo, M. Buser, and T. Brandes, Physical Review B 94, 214308 (2016).
- Carr et al. (2011) S. T. Carr, D. A. Bagrets, and P. Schmitteckert, Physical review letters 107, 206801 (2011).
- Schmitteckert et al. (2014) P. Schmitteckert, S. T. Carr, and H. Saleur, Physical Review B 89, 081401 (2014).
- Carr et al. (2015) S. T. Carr, P. Schmitteckert, and H. Saleur, Physica Scripta 2015, 014009 (2015).
- Humeniuk and Büchler (2017) S. Humeniuk and H. P. Büchler, Physical review letters 119, 236401 (2017).
- Cohen et al. (2015) G. Cohen, E. Gull, D. R. Reichman, and A. J. Millis, Physical review letters 115, 266802 (2015).
- Keldysh (1965) L. V. Keldysh, Soviet Physics JETP 20, 1018 (1965).
- Anderson (1961) P. W. Anderson, Physical Review 124, 41 (1961).
- Werner et al. (2009) P. Werner, T. Oka, and A. J. Millis, Physical Review B 79, 035320 (2009).
- Stefanucci and Almbladh (2004) G. Stefanucci and C.-O. Almbladh, Physical Review B 69, 195318 (2004).
- Antipov et al. (2017) A. E. Antipov, Q. Dong, J. Kleinhenz, G. Cohen, and E. Gull, Physical Review B 95, 085144 (2017).
- Konstantinov and Perel (1960) O. Konstantinov and V. Perel, Zhur. Eksptl’. i Teoret. Fiz. 39 (1960).
- Kambly and Flindt (2013) D. Kambly and C. Flindt, Journal of Computational Electronics 12, 331 (2013).
- Werner et al. (2006) P. Werner, A. Comanac, L. de’ Medici, M. Troyer, and A. J. Millis, Physical Review Letters 97, 076405 (2006).
- Mühlbacher and Rabani (2008) L. Mühlbacher and E. Rabani, Physical review letters 100, 176403 (2008).
- Eckstein and Werner (2010) M. Eckstein and P. Werner, Physical Review B 82, 115115 (2010).
- Härtle et al. (2015) R. Härtle, G. Cohen, D. Reichman, and A. Millis, Physical Review B 92, 085430 (2015).
- Chen et al. (2017a) H.-T. Chen, G. Cohen, and D. R. Reichman, The Journal of Chemical Physics 146, 054105 (2017a).
- Gull et al. (2010) E. Gull, D. R. Reichman, and A. J. Millis, Physical Review B 82, 075109 (2010).
- Gull et al. (2011a) E. Gull, D. R. Reichman, and A. J. Millis, Physical Review B 84, 085134 (2011a).
- Cohen et al. (2013) G. Cohen, E. Gull, D. R. Reichman, A. J. Millis, and E. Rabani, Physical Review B 87, 195108 (2013).
- Cohen et al. (2014a) G. Cohen, D. R. Reichman, A. J. Millis, and E. Gull, Physical Review B 89, 115139 (2014a).
- Cohen et al. (2014b) G. Cohen, E. Gull, D. R. Reichman, and A. J. Millis, Physical review letters 112, 146802 (2014b).
- Chen et al. (2017b) H.-T. Chen, G. Cohen, and D. R. Reichman, The Journal of Chemical Physics 146, 054106 (2017b).
- Dong et al. (2017) Q. Dong, I. Krivenko, J. Kleinhenz, A. E. Antipov, G. Cohen, and E. Gull, Phys. Rev. B 96, 155126 (2017).
- Gull et al. (2011b) E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. Troyer, and P. Werner, Reviews of Modern Physics 83, 349 (2011b).
- Schiró (2010) M. Schiró, Physical Review B 81, 085126 (2010).
- Kamenev (2011) A. Kamenev, Field theory of non-equilibrium systems (Cambridge University Press, 2011).
- Ridley et al. (2015) M. Ridley, A. MacKinnon, and L. Kantorovich, Physical Review B 91, 125433 (2015).
- Cohen and Rabani (2011) G. Cohen and E. Rabani, Physical Review B 84, 075150 (2011).
- Lesovik (1989) G. Lesovik, ZhETF Pisma Redaktsiiu 49, 513 (1989).
- Büttiker (1990) M. Büttiker, Physical review letters 65, 2901 (1990).
- Blanter and Büttiker (1999) Y. M. Blanter and M. Büttiker, Physical Review B 59, 10217 (1999).
- Vitushinsky et al. (2008) P. Vitushinsky, A. Clerk, and K. Le Hur, Physical review letters 100, 036603 (2008).
- Iannaccone et al. (1998) G. Iannaccone, G. Lombardi, M. Macucci, and B. Pellegrini, Physical review letters 80, 1054 (1998).
- Sukhorukov et al. (2001) E. V. Sukhorukov, G. Burkard, and D. Loss, Physical Review B 63, 125315 (2001).
- Thielmann et al. (2005) A. Thielmann, M. H. Hettler, J. König, and G. Schön, Physical review letters 95, 146806 (2005).
- Cronenwett et al. (1998) S. M. Cronenwett, T. H. Oosterkamp, and L. P. Kouwenhoven, Science 281, 540 (1998).
- Furusaki and Matveev (1995) A. Furusaki and K. Matveev, Physical Review B 52, 16676 (1995).
- König et al. (1998) J. König, H. Schoeller, and G. Schön, Physical Review B 58, 7882 (1998).
- Tran et al. (2008) T. Tran, I. Beloborodov, J. Hu, X. Lin, T. Rosenbaum, and H. Jaeger, Physical Review B 78, 075437 (2008).
- Lindley (1952) D. V. Lindley, in Mathematical Proceedings of the Cambridge Philosophical Society, Vol. 48 No. 2 (Cambridge University Press, 1952) pp. 277–289.
- Hewson (1997) A. C. Hewson, The Kondo problem to heavy fermions, Vol. 2 (Cambridge university press, 1997).
- Weiss (2012) U. Weiss, Quantum dissipative systems, Vol. 13 (World scientific, 2012).
- Wingreen and Meir (1994) N. S. Wingreen and Y. Meir, Physical review B 49, 11040 (1994).
- Anders (2008) F. B. Anders, Physical review letters 101, 066804 (2008).
- da Rocha et al. (2015) C. G. da Rocha, R. Tuovinen, R. van Leeuwen, and P. Koskinen, Nanoscale 7, 8627 (2015).
- Ridley and Tuovinen (2017) M. Ridley and R. Tuovinen, Physical Review B 96, 195429 (2017).
- Rothstein et al. (2009) E. Rothstein, O. Entin-Wohlman, and A. Aharony, Physical Review B 79, 075307 (2009).
- Gabdank et al. (2011) N. Gabdank, E. Rothstein, O. Entin-Wohlman, and A. Aharony, Physical Review B 84, 235435 (2011).
- Engel and Loss (2004) H.-A. Engel and D. Loss, Physical review letters 93, 136602 (2004).
- Entin-Wohlman et al. (2007) O. Entin-Wohlman, Y. Imry, S. Gurvitz, and A. Aharony, Physical Review B 75, 193308 (2007).
- Orth et al. (2012) C. P. Orth, D. F. Urban, and A. Komnik, Physical Review B 86, 125324 (2012).
- Miwa et al. (2017) K. Miwa, F. Chen, and M. Galperin, Scientific Reports 7, 9735 (2017).
- Strass et al. (2005) M. Strass, P. Hänggi, and S. Kohler, Physical review letters 95, 130601 (2005).