Non-Markovian full counting statistics in quantum dot molecules.

Full counting statistics of electron transport is a powerful diagnostic tool for probing the nature of quantum transport beyond what is obtainable from the average current or conductance measurement alone. In particular, the non-Markovian dynamics of quantum dot molecule plays an important role in the nonequilibrium electron tunneling processes. It is thus necessary to understand the non-Markovian full counting statistics in a quantum dot molecule. Here we study the non-Markovian full counting statistics in two typical quantum dot molecules, namely, serially coupled and side-coupled double quantum dots with high quantum coherence in a certain parameter regime. We demonstrate that the non-Markovian effect manifests itself through the quantum coherence of the quantum dot molecule system, and has a significant impact on the full counting statistics in the high quantum-coherent quantum dot molecule system, which depends on the coupling of the quantum dot molecule system with the source and drain electrodes. The results indicated that the influence of the non-Markovian effect on the full counting statistics of electron transport, which should be considered in a high quantum-coherent quantum dot molecule system, can provide a better understanding of electron transport through quantum dot molecules.

Full counting statistics1 (FCS) of electron transport through mesoscopic system has attracted considerable attention both experimentally and theoretically because it can provide a deeper insight into the nature of electron transport mechanisms, which cannot be obtained from the average current2345678910. For instance, the shot noise measurements can be used to probe the dynamical in an open double quantum dots (QDs)11, the coherent coupling between serially coupled QDs12, the evolution of the Kondo effect in a QD13, and the conduction channels of quantum conductors14. In particular, shot noise characteristics can provide information about the feature of the pseudospin Kondo effect in a laterally coupled double QDs15, the spin accumulations in a electron reservoir16, and the charge fractionalization in the ν = 2 quantum Hall edge17. In addition, the degree of entanglement of two electrons in the double QDs18, the dephasing rate in a closed QD19, the internal level structure of single molecule magnet2021 can be characterized by the super-Poissonian shot noise.

On the other hand, the quantum coherence in coupled QD system, which is characterized by the off-diagonal elements of the reduced density matrix of the QD system within the framework of the density matrix theory22, plays an important role in the electron tunneling processes and has a significant influence on electron transport2324252627282930313233. In particular, theoretical studies have demonstrated that the high-order cumulants, e.g., the shot noise, the skewness, are more sensitive to the quantum coherence than the average current in the different types of QD systems123435363738 and the quantum coherence information in a side-coupled double QD system can be extracted from the high-order current cumulants35. In fact, the non-Markovian dynamics of the QD system also plays an important role in the non-equilibrium electron tunneling processes. However, the above studies on current noise or FCS were mainly based on the different types of Markovian master equations. Although the influence of non-Markovian effect on the long-time limit of the FCS in the QD systems has received some attention333940414243444546, how the non-Markovian effect affects the FCS is still an open issue, especially the influence of the interplay between the quantum coherence and non-Markovian effect on the long-time limit of the FCS has not yet been revealed.

The aim of this report is thus to derive a non-Markovian FCS formalism based on the exact time-convolutionless (TCL) master equation and study the influences of the quantum coherence and non-Markovian effect on the FCS in QD molecule systems. It is demonstrated that the non-Markovian effect manifests itself through the quantum coherence of the considered QD molecule system, and has a significant impact on the FCS in the high quantum-coherent QD molecule system, which depends on the coupling of the considered QD molecule system with the incident and outgoing electrodes. Consequently, it is necessary to consider the influence of the non-Markovian effect on the full counting statistics of electron transport in a high quantum-coherent single-molecule system.

Results

We now study the influences of the quantum coherence and non-Markovian effect on the FCS of electronic transport through the QD molecule system. In order to facilitate discussions effectively, we consider three typical QD systems, namely, single QD without quantum coherence, serially coupled double QDs and side-coupled double QDs with high quantum coherence in a certain parameter regime (see Fig. 1). In addition, we assume the bias voltage (μ = −μ = V/2) is symmetrically entirely dropped at the QD-electrode tunnel junctions, which implies that the levels of the QDs are independent of the applied bias voltage even if the couplings are not symmetric, and choose meV as the unit of energy which corresponds to a typical experimental situation47.

Figure 1

Schematic of the two single-level QD molecules weakly coupled to two electrodes, (a) serially coupled double QDs, (b) side-coupled double QDs. Here, the two QD molecules possess high quantum coherence in the case of (Δ being the singly-occupied eigenenergy separation, k the Boltzmann constant, T the temperature of the QDs system). The hopping coupling between the two QDs, and the strength of coupling between the QDs system and the electrode α, are characterized by J and Γ, respectively.

Single quantum dot without quantum coherence

In this subsection, we consider a single QD weakly coupled to two ferromagnetic electrodes. The Hamiltonian of the considered system is described by the H = H + H + H. The QD Hamiltonian H is given bywhere creates (annihilates) an electron with spin σ and on-site energy ε (which can be tuned by a gate voltage V) in this QD system. U is the intradot Coulomb interaction between two electrons in the QD system.

The relaxation in the two ferromagnetic electrodes is assumed to be sufficiently fast, so that their electron distributions can be described by equilibrium Fermi functions. The two electrodes are thus modeled as non-interacting Fermi gases and the corresponding Hamiltonians can be expressed aswhere creates (annihilates) an electron with energy ε, spin s and momentum k in α (α = L, R) electrode, and s = + (−) denotes the majority (minority) spin states with the density of states g,. The polarization vectors p (left lead) and p (right lead) are parallel to each other, and their magnitudes are characterized by p = |p| = (g,+ − g,−)/(g,+ + g,−). The tunneling between the QD and the electrodes is described bywhere spin-up ↑ and spin-down ↓ are defined to be the majority spin and minority spin of the ferromagnet, respectively.

The QD-electrode coupling is assumed to be sufficiently weak, thus, the sequential tunneling is dominant and can be well described by the quantum master equation of reduced density matrix spanned by the eigenstates of the QD. The particle-number-resolved TCL quantum master equation for the reduced density matrix of the considered single QD is given byFor more details, see Methods section. Here, the complete basis {|0, 0〉, |↑, 0〉, |↓, 0〉, |↑, ↓〉} is chosen to describe the electronic states of this single QD system, and the single QD system parameters are chosen as , U = 5, p = 0.9 and k = 0.04.

Figure 2 shows the first four current cumulants as a function of the bias voltage for different ratios Γ/Γ describing the left-right asymmetry of the QD-electrode coupling. We found that the non-Markovian effect has no influence on the current noise behaviors of the single QD considered here, see Fig. 2. Scrutinizing Eq. (4), it is found that for the non-Markovian case the elements of the reduced density matrix are equivalent to that for the Markovian case because there are not the off-diagonal elements of the reduced density matrix. Thus, the equations of motion of the four elements of the reduced density matrix can be expressed asHere, f,+ is the Fermi function of the electrode α, and f,− = 1−f,+. The detailed procedure for calculation of the equation of motion of a reduced density matrix, see Methods section. Within the framework of the density matrix theory, the off-diagonal elements of the reduced density matrix characterize the quantum coherence of the considered QD system. Thus, the influence of the non-Markovian effect on the FCS may be associated with the quantum coherence of the considered QD system. In order to confirm this conclusion, we take serially coupled and side-coupled double QDs for illustration in the following two subsection.

Figure 2

The average current (〈I〉), shot noise (C2/C1), skewness (C3/C1) and kurtosis (C4/C1) versus bias voltage for the Morkovian and the non-Markovian case at different coupling of the single QD with two ferromagnetic electrodes, respectively.

Here, C is the zero-frequency k-order cumulant of current fluctuations. The non-Markovian effect has no influence on the first four current cumulants of the considered single QD. The single QD system parameters: , U = 5, p = 0.9 and k = 0.04.

Serially coupled double quantum dots with high quantum coherence

We now consider two serially coupled double QDs weakly connected to two metallic electrodes, see Fig. 1(a). For the sake of simplicity, the spin degree of freedom has not been considered. The double-QD is described by a spinless Hamiltonianwhere creates (annihilates) an electron with energy ε (which can be tuned by a gate voltage V) in ith QD. U is the interdot Coulomb repulsion between two electrons in the double QD system, where we consider the intradot Coulomb interaction U → ∞, so that the double-electron occupation in the same QD is prohibited. The last term of H describes the hopping coupling between the two dots with J being the hopping parameter. The two metallic electrodes are modeled as non-interacting Fermi gases and the corresponding Hamiltonians are given bywhere creates (annihilates) an electron with energy ε and momentum k in α (α = L, R) electrode. The tunneling between the double QDs and the two electrodes is described byFor the case of the weak QD-electrode coupling, the particle-number-resolved TCL quantum master equation for the reduced density matrix of the considered serially double-QD system readsHere, we can diagonalize the serially coupled double QDs Hamiltonian Hdot,2 in the basis represented by the electron occupation numbers in the QD-1 and QD-2 denoted respectively by N and N, namely, {|0, 0〉, |1, 0〉, |0, 1〉, |1, 1〉}, and obtain the corresponding four eigenstates of the considered serially coupled double QDs system48withandHere, we focus on the regime , where the hopping coupling between the two QDs strongly modifies the internal dynamics, and the off-diagonal elements of the reduced density matrix play an essential role in the electron tunneling processes23495051. In the following numerical calculations, thus, the parameters of the serially coupled double QDs system are chosen as , J = 0.001, U = 4 and k = 0.05.

When the coupling of the QD-2 with the right (drain) electrode is stronger than that of the QD-1 with the left (source) electrode, namely, Γ/Γ < 1, we plot the first four current cumulants as a function of the bias voltage for different values of the QD-2-electrode coupling Γ at Γ/Γ = 0.1 in Figs. 3(a)–3(d). We found that the non-Markovian effect has a very weak influence on the FCS. Interestingly, the high-order current cumulants the skewness and the kurtosis can still show the tiny differences, see Figs. 3(c) and 3(d). Whereas for the Γ/Γ ≥ 1 case, the non-Markovian effect has a significant impact on the FCS, see Fig. 4. Especially, for a relatively large value of the ratio Γ/Γ = 10 and the coupling of the QD-1 with the left electrode being stronger than the hoping coupling, namely, Γ/J > 1, the non-Markovian effect can induce a strong negative differential conductance (NDC) and super-Poissonian noise, see Figs. 4(e) and 4(f). In addition, in the case of Γ/Γ ≥ 1 and Γ/J > 1, the transitions of the skewness and the kurtosis from positive (negative) to negative (positive) values are observed, see the dotted line in Fig. 4(c), the dotted and dash-dot-dotted lines in Fig. 4(d), and the dash-dot-dotted line in Fig. 4(h). It is well known that the skewness and the kurtosis (both its magnitude and sign) characterize, respectively, the asymmetry of and the peakedness of the distribution around the average transferred-electron number during a time interval t, thus that provides further information for the counting statistics beyond the shot noise.

Figure 3

(a)–(d) The average current (〈I〉), shot noise (C2/C1), skewness (C3/C1) and kurtosis (C4/C1) versus bias voltage for the Morkovian and the non-Markovian case at different values of the QD-2-electrode coupling Γ with Γ/Γ = 0.1. Here, C is the zero-frequency k-order cumulant of current fluctuations. The non-Markovian effect in the Γ/Γ = 0.1 case has a weak influence on the the first four current cumulants. The serially coupled double QDs system parameters: , J = 0.001, U = 4 and k = 0.05. (e)–(h) The average current (〈I〉), shot noise (C2/C1), skewness (C3/C1) and kurtosis (C4/C1) versus bias voltage for different coupling of the serially coupled double QDs system with two metallic electrodes. Here three cases are considered, namely, (1) the Markovian and the diagonal elements of the reduced density matrix, (2) the Markovian and the off-diagonal elements of the reduced density matrix, (3) the non-Markovian and the off-diagonal elements of the reduced density matrix. The non-Markovian effect has a very weak influence on the first four current cumulants in the serially coupled double QD system with a relatively small quantum coherence. The serially coupled double QDs system parameters: , J = 1, U = 4 and k = 0.05.

Figure 4

The average current (〈I〉), shot noise (C2/C1), skewness (C3/C1) and kurtosis (C4/C1) versus bias voltage for the Morkovian and the non-Markovian case at different values of the QD-2-electrode coupling Γ.

(a)–(d) for Γ/Γ = 1, (e)–(h) for Γ/Γ = 10. Here, C is the zero-frequency k-order cumulant of current fluctuations. The non-Markovian effect in the Γ/Γ ≥ 1 case has a significant impact on the first four cumulants of transport current. The serially coupled double QDs system parameters: , J = 0.001, U = 4 and k = 0.05.

To discuss the underlying mechanisms of the current noise clearly, for the system parameters considered here, the two singly-occupied eigenstates and eigenvalues can be expressed asHere we have utilized the equations and . In this situation, the equations of motion of the six elements of the reduced density matrix are given bywhere , (Ψ is the digamma function) and . Compared with the Markovian case, it is obvious that the non-Markovian effect manifests itself through the off-diagonal elements of the reduced density matrix, namely, the quantum coherence of the considered QDs system. In Fig. 5(a), we plot the functions Φ − 0.1Φ (Γ = 0.1Γ), Φ − Φ (Γ = Γ) and 0.1Φ − Φ (Γ = 0.1Γ) as a function of bias voltage. It is clearly evident that the values of the functions Φ − 0.1Φ and Φ − Φ show significant variations with increasing bias voltage, especially in the vicinity of the bias voltages V = 2 and V = 10 because the new transport channels begin to participate in quantum transport; while 0.1Φ − Φ has a gentle variation. Consequently, the non-Markovian effects in the Γ/Γ ≥ 1 case have a remarkable impact on the FCS, see Fig. 4. Moreover, for Γ/Γ = 10 case, the non-Markovian effect has a more significant on the FCS than the Γ/Γ = 1 case, which originates from the QD-2-electrode coupling Γ is weaker than the hoping coupling J, where the electron tunneling from QD-1 can not tunnel out QD-2 very quickly and still influence the internal dynamics.

Figure 5

(a) The functions Φ − 0.1Φ (Γ = 0.1Γ), Φ − Φ (Γ = Γ) and 0.1Φ − Φ (Γ = 0.1Γ) as a function of bias voltage with U = 4 and k = 0.05. (b) The functions Φ + 0.1Φ (Γ = 0.1Γ), Φ + Φ (Γ = Γ) and 0.1Φ + Φ (Γ = 0.1Φ) as a function of bias voltage with U = 5 and k = 0.1. Here, , and Ψ is the digamma function. The variation of the value of the above mentioned function is responsible for whether the non-Markovian effect has a remarkable influence on the first four cumulants of transport current.

In order to illustrate whether the non-Markovian effect has a weak influence on the FCS in a relatively small quantum-coherent QD system, we consider the regime (J = 1), where the off-diagonal elements of the reduced density matrix have little influence on the electron tunneling processes. We find that for the J = 1 case the diagonal elements of the reduced density matrix play a major role in the electron tunneling processes, and the non-Markovian effect in this case indeed has little impact on the FCS, see Figs. 3(e)–3(h). Consequently, the influence of the non-Markovian effect on the FCS depends on the quantum coherence of the considered QD system. To prove whether this conclusion is universal or not, we take side-coupled double QDs for further illustration in the following subsection.

Side-coupled double quantum dots with high quantum coherence

We consider here a side-coupled double QDs system. In this case, the QD-1 is only weakly coupled to the two electrodes, see Fig. 1(b). The QD-electrode tunneling is thus described by

In the case of the QD-electrode weak coupling, the particle-number-resolved TCL quantum master equation for the side-coupled double QDs can be expressed asHere, the eigenstates and eigenvalues of the side-coupled double QDs system are the same as the serially coupled double QDs system. In the following numerical calculations, the parameters of the side-coupled QDs system are chosen as , J = 0.001, U = 5 and k = 0.1.

For the present side-coupled QDs system with high quantum coherence, we find that for Γ/Γ ≥ 1 case the non-Markovian effect has a more remarkable impact on the FCS than that in the serially coupled double QDs system, but the NDC does not appear, see Figs. 4 and 6. For instance, in the case of Γ/J > 1 and Γ/Γ = 1, the non-Markovian effect can further enhance the super-Poissonian shot noise, see the dotted and dash-dot-dotted lines in Fig. 6(b); and the transitions of the skewness and the kurtosis from a relatively small positive to a large negative values take place, especially for a relatively large value Γ/J the kurtosis can be further decreased to a very large negative value, see the dotted and dash-dot-dotted lines in Figs. 6(c) and 6(d). While for the Γ/J > 1 and Γ/Γ = 10 case the non-Markovian effect can enhance the shot noise to a super-Poissonian value, see the dotted and dash-dot-dotted lines in Fig. 6(f), and the transition of the kurtosis from small positive to large negative values only takes place, see the dotted and dash-dot-dotted lines in Fig. 6(h). For the system parameters considered here, namely, in the limit of , the equations of motion of the six elements of the reduced density matrix readFrom the above four equations, we find that these characteristics also originate from the quantum coherence of the side-coupled double QDs, and can also be understood in terms of the functions Φ + 0.1Φ and Φ + Φ, which have considerable variations in the vicinity of the bias voltages V = 2 and V = 12 because the new transport channels begin to enter the bias voltage window, see the solid and dashed lines in Fig. 5(b). As for the Γ/Γ < 1 case the non-Markovian effect has a slightly influence on the FCS because the function 0.1Φ + Φ has a gentle variation with increasing the bias voltage, see the dotted line in Fig. 5(b), which is the same as the serially coupled double QDs system, see Figs. 3(a)–3(d) and 7.

Figure 6

The average current (〈I〉), shot noise (C2/C1), skewness (C3/C1) and kurtosis (C4/C1) versus bias voltage for the Morkovian and the non-Markovian case at different values of the QD-1-electrode coupling Γ.

(a)–(d) for Γ/Γ = 1, (e)–(h) for Γ/Γ = 10. Here, C is the zero-frequency k-order cumulant of current fluctuations. The non-Markovian effect in the Γ/Γ ≥ 1 case has a more remarkable impact on the first four cumulants of transport current than that in the serially coupled double QDs system, but the NDC does not appear. The side-coupled double QDs system parameters: , J = 0.001, U = 5 and k = 0.1.

Figure 7

The average current (〈I〉), shot noise (C2/C1), skewness (C3/C1) and kurtosis (C4/C1) versus bias voltage for the Morkovian and the non-Markovian case at different values of the QD-1-electrode coupling Γ with Γ/Γ = 0.1.

Here, C is the zero-frequency k-order cumulant of current fluctuations. The non-Markovian effect in the Γ/Γ = 0.1 case has a slightly influence on the the first four current cumulants. The other system parameters are the same as in Fig. 6.

In addition, it should be pointed out that for Γ/Γ = 1 the non-Markovian effect has a stronger impact on the FCS than that for Γ/Γ > 1 case, which is contrary to the case of the serially coupled double QDs system. For the the side-coupled double QDs system, the quantum coherence originates from the quantum interference between the direct electron tunneling process, namely, the conduction-electron tunneling into the QD-1 and then directly tunneling out of the QD-1 onto the drain electrode, and the indirect tunneling process, namely, the conduction-electron from the source electrode first tunneling from the QD-1 to the QD-2, then tunneling back into the QD-1, and at last tunneling out of the QD-1. Thus, the fast direct tunneling process in the Γ = 10Γ case can be suppressed compared with the Γ = Γ case, which leads to the non-Markovian effect has a relatively strong impact on the FCS in the Γ/Γ = 1 case.

Discussion

We have developed a non-Markovian FCS formalism based on the exact TCL master equation, and studied the influence of the interplay between the quantum coherence and non-Markovian effect on the long-time limit of the FCS in three QD systems, namely, single QD, serially coupled double QDs and side-coupled double QDs. It is demonstrated that the non-Markovian effect manifests itself through the quantum coherence of the considered QD molecule system, and especially has a significant impact on the FCS in the high quantum-coherent QD molecule system, which depends on the coupling of the considered QD molecule system with the source and drain electrodes. For the single QD system without quantum coherence, the non-Markovian effect has no influence on the current noise properties; whereas for the serially coupled and side-coupled double QDs systems with high quantum coherence, that has a remarkable impact on the FCS when the coupling of the considered QD molecule with the incident electrode is equal to or stronger than that with the outgoing electrode. For instance, for the high quantum-coherent serially coupled double QDs system, the non-Markovian effect can induce a strong NDC and change the shot noise from the sub-Poissonian to super-Poissonian distribution in the case of and Γ > J; while for the high quantum-coherent side-coupled double QDs system, that can remarkably enhance the super-Poissonian noise or the sub-Poissonian noise for the Γ/Γ ≥ 1 case. Moreover, the non-Markovian effect can also lead to the occurrences of the skewness and kurtosis from small positive to large negative values. These results indicated that the influence of the non-Markovian effect on the long-time limit of the FCS should be considered in a highly quantum-coherent single-molecule system.

Methods

Particle-number-resolved time-convolutionless quantum master equation

We consider a general transport setup consisting of a single-level QD molecule weakly coupled to the two electrodes, see Fig. 1, which is described by the following HamiltonianHere, the first term stands for the Hamiltonians of the two electrodes, with ε being the energy dispersion, and the annihilation (creation) operators in the α electrode. The second term , which may contain vibrational or spin degrees of freedom and different types of many-body interaction, represents the QD molecule Hamiltonian, where is the creation (annihilation) operator of electrons in a quantum state denoted by μ. The third term describes the tunneling coupling between the QD molecule and the two electrodes, which is assumed to be a sum of bilinear terms that each create an electron in the QD molecule and annihilate one in the electrodes or vice versa.

The QD-electrode coupling is assumed to be sufficiently weak, so that Hhyb can be treated perturbatively. In the interaction representation, the equation of motion for the total density matrix readswithwhere and . In order to derive an exact equation of motion for the reduced density matrix ρ of the QD molecule system, it is convenient to define a super-operator according towith ρ being some fixed state of the electron electrode. Accordingly, a complementary super-operator readsFor a factorizing initial condition , , and . Using the TCL projection operator method52, one can obtain the second-order TCL master equationThe Eq. (31) is the starting point of deriving the particle-number-resolved quantum master equation. Using Eqs. (28) and (29), after some algebraic calculations we can rewrite Eq. (31) as

In order to fully describe the electron transport problem, we should record the number of electrons arriving at the drain electrode, which emitted from the source electrode and passing through the QD molecule. We follow Li and co-authors5354 and introduce the Hilbert subspace B( (n = 1, 2, …) corresponding to n electrons arriving at the drain electrode, which is spanned by the product of all many-particle states of the two isolated electrodes, and formally denoted as . Then, the entire Hilbert space of the two electrodes can be expressed as . With this classification of the electrode states, the average over states in the entire Hilbert space B in Eq. (32) should be replaced with the states in the subspace B(, and leading to a conditional TCL master equation

To proceed, two physical considerations are further implemented. (i) Instead of the conventional Born approximation for the entire density matrix , the ansatz is proposed, where being the electrode density operator associated with n electrons arriving at the drain electrode. With this ansatz for the entire density operator, tracing over the subspace B(, the Eq. (33) can be rewritten asHere we have used the orthogonality between the states in different subspaces. (ii) The extra electrons arriving at the drain electrode will flow back into the source electrode via the external closed transport circuit. Moreover, the rapid relaxation processes in the electrodes will bring the electrodes to the local thermal equilibrium states quickly, which are determined by the chemical potentials. Consequently, after the procedure done in Eq. (34), the electrode density matrices and should be replaced by . In the Schrödinger representation, the Eq. (34) can be expressed aswhere the correlation function are defined asIntroducing the following super-operatorsthen, the Eq. (35) can be rewritten as a compact formwhere . The above equation is the starting point of the non-Markovian FCS calculation.

Non-Markovian full counting statistics

In this subsection, we outline the procedure to calculate the non-Markovian FCS based on Eq. (38). The FCS can be obtained from the cumulant generating function (CGF) F (χ) which related to the probability distribution P (n, t) by5455 , where χ is the counting field. The CGF F (χ) connects with the particle-number-resolved density matrix ρ( (t) by defining . Evidently, we have e− = Tr[S (χ, t)], where the trace is over the eigenstates of the QD molecule system. Since Eq. (38) has the following form , then, S (χ, t) satisfies , where S is a column matrix, and A, C and D are three square matrices. The specific form of L can be obtained by performing a discrete Fourier transformation to the matrix element of Eq. (38). In the low frequency limit, the counting time, namely, the time of measurement is much longer than the time of tunneling through the QD molecule system. In this case, F (χ) is given by344043555657 F (χ) = −λ1 (χ) t, where λ1 (χ) is the eigenvalue of L which goes to zero for χ → 0. According to the definition of the cumulants one can express λ1 (χ) as . The low order cumulants can be calculated by the Rayleigh–Schrödinger perturbation theory in the counting parameter χ. In order to calculate the first four current cumulants we expand L to four order in χand define the two projectors40435658 and Q = Q2 = 1 − P, obeying the relations PL0 = L0P = 0 and QL0 = L0Q = L0. Here, |0〉〉 is the right eigenvector of L0, i.e., L0 |0〉〉 = 0, and is the corresponding left eigenvector. In view of L0 being singular, we also introduce the pseudoinverse according to , which is well-defined due to the inversion being performed only in the subspace spanned by Q. After a careful calculation, λ1 (χ) is given byFrom Eq. (40) we can identify the first four current cumulants:Here, it is important to emphasize that the first four cumulants C are directly related to the transport characteristics. For example, the first-order cumulant (the peak position of the distribution of transferred-electron number) gives the average current 〈I〉 = eC1/t. The zero-frequency shot noise is related to the second-order cumulant (the peak-width of the distribution) . The third-order cumulant and four-order cumulant characterize, respectively, the skewness and kurtosis of the distribution. Here, . In general, the shot noise, skewness and kurtosis are represented by the Fano factor F2 = C2/C1, F3 = C3/C1 and F4 = C4/C1, respectively.

Author Contributions

H.B.X. conceived the idea and designed the research and performed calculations. H.J.J., J.Q.L. and W.M.L. contributed to the analysis and interpretation of the results and prepared the manuscript.