Quantum State of the Fermionic Carriers in a Transport Channel Connecting Particle Reservoirs

Abstract

We analyze the quantum state of fermionic carriers in a transport channel attached to a particle reservoir. The analysis is done from first principles by considering microscopic models of the reservoir and transport channel. In the case of infinite effective temperature of the reservoir we demonstrate a full agreement between the results of straightforward numerical simulations of the system dynamics and the solution of the master equation on the single-particle density matrix of the carriers in the channel. This allows us to predict the quantum state of carriers in the case where the transport channel connects two reservoirs with different chemical potentials.

1. Introduction

Electron transport in mesoscopic devices is a wide subfield in solid-state physics [1,2,3]. The studies on electron transport are aimed at control of the electron current between two or more contacts (electron reservoirs) attached to a quantum device. Recently, the problems of the same kind have been addressed for the principally different system—charge neutral atoms in laser-based devices, both experimentally [4,5,6,7], and theoretically [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22]. The advantage of the latter systems against the electron system is the perfect control over the system parameters and effective detection techniques that allow for in situ measurement of the quantum state of carriers in the device which, following Reference [4], we refer to as the transport channel connecting particle reservoirs.

On the formal level the quantum state of carriers in the transport channel is characterized by the single-particle density matrix (SPDM), the knowledge of which suffices to predict the current between reservoirs. In the present paper we analyze the SPDM of fermionic carriers from first principles with the emphasis on decoherence effect of reservoirs. Clearly, to address this problem from first principles one needs physically relevant microscopic models of the transport channel and particle reservoir. Having in mind cold atoms we model the transport channel by the tight-binding chain, which is known to adequately describe neutral atoms in deep optical lattices. (Here ‘deep’ means that the width of the ground Bloch band is smaller than the energy gap separating it from the rest of the spectrum.) As for the particle reservoir, we model it by the Two-Body Random Interaction Model (TBRIM) [23,24] that corresponds to a system of N weakly interacting spinless fermions distributed over M natural orbitals. The closed (isolated) TBRIM possesses the self-thermalization property [25,26]. This means, in particular, that one has a meaningful notion of the temperature of TBRIM, despite the fact that the system has no contact with a thermostat. It has been also shown in the recent work [27] that TBRIM retains the self-thermalization property when the system is open, which makes it an excellent model for the reservoir of fermionic particles.

The structure of the paper is as follows. After reviewing TBRIM in Section 2, we attach a finite-length tight-binding chain to this reservoir and study particles propagation across the chain in Section 3. We quantify decoherence effect of the reservoir on the carriers in the channel by the von Neumann entropy of SPDM and show that it is strictly positive. In Section 4 we compare the exact numerical results with those obtained by using the master equation on the reduced density matrix of the carriers (RDM). Finally, in the concluding Section 5 we summarise the results and give the list of open problems.

2. The Model

In this section we specify the system Hamiltonian $\widehat{H}$, which consists of the Hamiltonian of the particle reservoir ${\widehat{H}}_b$, the Hamiltonian of the transport channel ${\widehat{H}}_s$, and the coupling Hamiltonian ${\widehat{H}}_{int}$:

(1)$\widehat{H}={\widehat{H}}_b+{\widehat{H}}_s+{\widehat{H}}_{int}.$

2.1. The Particle Reservoir

We model the particle reservoir by TBRIM which describes N interacting spinless fermions distributed over M natural orbitals with the energies ${ϵ}_k$ (${ϵ}_{k+1}\ge {ϵ}_k$):

(2)$\begin{array}{c}\hfill {\widehat{H}}_b=\sum _{k=1}^M{ϵ}_k{\widehat{d}}_k^{†}{\widehat{d}}_k+{\epsilon }_b\sum _{ijkl}{V}_{ij,kl}{\widehat{d}}_i^{†}{\widehat{d}}_j^{†}{\widehat{d}}_k{\widehat{d}}_l.\end{array}$

Here operators ${\widehat{d}}_i^{†},{\widehat{d}}_i$ satisfy the usual anti-commutation relation and one-particle energies ${ϵ}_k$ and interaction constants $V_{ij,kl}$ are random (up to the obvious symmetry relations insuring hermiticity of the Hamiltonian) variables with standard deviation equal to unity. The parameter ${\epsilon }_b$ in the Hamiltonian (2) controls the strength of two-body interactions which couples every Fock state with other $K=1+N(M-N)+N(N-1)(M-N)(M-N-1)/4$ Fock states. In the paper we assume ${\epsilon }_b\ll 1$, i.e., we consider the limit of weakly interacting fermions. Yet, ${\epsilon }_b$ must larger than the critical value where TBRIM shows the transition to Quantum Chaos [28,29]. An analytical estimate for the critical interaction strength can be obtained by using the Åberg criteria [30], while numerically this transition is detected as the change of the level-spacing distribution from the Poisson distribution to the Wigner-Dyson distribution. In what follows we fix the reservoir size to $M=12$, $N=6$, and set ${\epsilon }_b=0.0085$ where the energy level statistics perfectly follows the Wigner-Dyson distribution. The chosen ${\epsilon }_b$ is approximately twice larger than the critical value. At the same time, it is small enough to speak of weakly-interacting fermions. In particular, the mean density of states, which in the case of non-interacting fermions is well approximated by the Gaussian of the width ∼$\sqrt{N}$, remains practically unaffected.

Provided the condition of Quantum Chaos is satisfied, the system (2) shows the phenomenon of self-thermalization [26]. It means that for any given eigenstate $|{\psi }_E⟩$ occupation numbers of the natural orbitals $n_k=⟨{\psi }_E|{\widehat{d}}_k^{†}{\widehat{d}}_k|{\psi }_E⟩$ obey (of course, with some fluctuations) the Fermi-Dirac distribution,

(3)$n_k=\frac{1}{e^{\beta ({ϵ}_k-\mu )}+1},$

where the inverse effective temperature $\beta$ and the chemical potential $\mu$ are uniquely determined by the eigenstate energy E and the number of particles N through the solution of the following system of two non-linear algebraic equations,

(4)$\sum _{k=1}^M\frac{1}{e^{\beta ({ϵ}_k-\mu )}+1}=N,\sum _{k=1}^M\frac{{ϵ}_k}{e^{\beta ({ϵ}_k-\mu )}+1}=E.$

Then the ground and the highest energy eigenstates of the system (2) corresponds to $\beta =\pm \infty$, while an eigenstate from the middle of the spectrum corresponds to $\beta =0$. We mention that Equations (3)–(4) also hold for the open TBRIM [27], where the number of particles changes in time.

2.2. The Transport Channel

We model the transport channel by the tight-binding chain,

(5)$\begin{array}{c}\hfill {\widehat{H}}_s=V_g\sum _{l=1}^L{\widehat{c}}_l^{†}{\widehat{c}}_l-\frac{J}{2}\left(\sum _{l=1}^{L-1}{\widehat{c}}_{l+1}^{†}{\widehat{c}}_l+h.c.\right),\end{array}$

where J is the hopping matrix element and $V_g$ has the meaning of the gate voltage. We shall characterise fermions in the channel by SPDM,

(6)$\begin{array}{c}\hfill {\rho }_{l,m}\left(t\right)=⟨\mathsf{\Psi }\left(t\right)|{\widehat{c}}_l^{†}{\widehat{c}}_m|\mathsf{\Psi }\left(t\right)⟩,\end{array}$

where $\left|\mathsf{\Psi }\right(t)⟩$ is the total wave function of the whole system defined in the Hilbert space of the dimension

(7)$\mathcal{N}=\frac{(M+L)!}{N!(M+L-N)!}.$

Since the quadratic form ${\widehat{c}}_l^{†}{\widehat{c}}_m$ in Equation (6) conserves the number of particles, the density matrix (6) can be presented as a sum of partial density matrices,

(8)$\begin{array}{c}\hfill \rho \left(t\right)=\sum _{i=1}^L{\rho }^{\left(i\right)}\left(t\right),\end{array}$

where ${\rho }^{\left(i\right)}\left(t\right)$ refer to the fixed number fermions in the channel. We note that for an isolated channel with i fermions in a pure state the matrix ${\rho }^{\left(i\right)}\left(t\right)$ has i eigenvalues equal to unity and $L-i$ eigenstates equal to zero.

2.3. The Coupling Hamiltonian

Particles from the reservoir enter the transport channel due to the coupling controlled by the Hamiltonian

(9)$\begin{array}{c}\hfill {\widehat{H}}_{int}=\epsilon \left(\sum _{k=1}^M{W}_k{\widehat{c}}_1^{†}{\widehat{d}}_k+h.c.\right),\end{array}$

where $W_k$ are random entries of the same magnitude as the interaction constants $V_{ij,kl}$ and $\epsilon$ is our control parameter. In what follows we consider the situation where initially all particles are in the reservoir, i.e.,

(10)$\begin{array}{c}\hfill |\mathsf{\Psi }(t=0)⟩=|{\psi }_E⟩\otimes |vac⟩.\end{array}$

3. System Dynamics

In this section we discuss the system dynamics governed by the Schrödinger equation with the Hamiltonian (1) for the initial condition specified in Equation (10).

3.1. Population Dynamics

Figure 1 shows the occupation numbers of the natural orbitals and the chain sites,

(11)$\begin{array}{c}{n}_k\left(t\right)=|⟨\mathsf{\Psi }\left(t\right)|{\widehat{d}}_k^{†}{\widehat{d}}_k{|\mathsf{\Psi }\left(t\right)⟩|}^2,\hfill \\ n_l\left(t\right)=|⟨\mathsf{\Psi }\left(t\right)|{\widehat{c}}_l^{†}{\widehat{c}}_l{|\mathsf{\Psi }\left(t\right)⟩|}^2,\hfill \end{array}$

as the function of time for $|{\psi }_E⟩$ from the middle of the spectrum of the system (2) where all reservoir modes are equally populated, i.e., $n_k\approx N/M$. According to Equation (3) this corresponds to the infinite temperature of the reservoir. One distinguishes two qualitatively different stages/regimes in Figure 1. During the fist stage fermionic particles propagate in the channel with the velocity determined by the hopping matrix element J in Equation (5). Reaching the boundary particles are reflected back towards the reservoir. Notice that during this stage, which we refer to as propagation stage, the number of particles in the channel monotonically increases. During the second stage, which we refer to as equilibration stage, occupation of the chain sites and natural orbitals equilibrate at $n_k=n_l=N/(M+L)$.

To get a deeper insight into the population dynamics we calculate the partial density matrices ${\rho }^{\left(i\right)}\left(t\right)$, see Equation (8). The upper panel in Figure 2 shows probabilities $P_i\left(t\right)$ to find i fermions in the channel at a given time t, which is given by the equation

(12)$\begin{array}{c}\hfill P_i\left(t\right)=\mathrm{Tr}\left[{\rho }^{\left(i\right)}\left(t\right)\right]/i.\end{array}$

Increasing the evolution time we find $P_i\left(t\right)$ to approach the value $P_i(t=\infty )={\mathcal{N}}_i/\mathcal{N}$ where ${\mathcal{N}}_i$ is dimension of the sub-space of the Hilbert space defined by the condition that there are i particles in the channel. This result proves that for infinite reservoir temperature we have complete equilibration between the system (the tight-binding chain) and the bath (TBRIM).

3.2. Decoherence Dynamics

Next we discuss decoherence effect due to the reservoir. We characterize coherence of the carriers in the transport channel by the von Neumann entropy for the normalised partial SPDMs

(13)$\begin{array}{c}\hfill {\rho }^{\left(i\right)}\left(t\right)\to \frac{{\rho }^{\left(i\right)}\left(t\right)}{P_i\left(t\right)},S_i\left(t\right)=-\mathrm{Tr}[{\rho }^{\left(i\right)}\left(t\right)log{\rho }^{\left(i\right)}\left(t\right)].\end{array}$

(The information von Neumann entropy should not be mismatch with the thermodynamic entropy given by the logarithm of the number of states.) Entropies $S_i\left(t\right)$ are depicted in the lower panel in Figure 2. It is seen that decoherence takes place immediately after the particles enter the transport channel and $S_i\left(t\right)$ quietly reach the maximally possible values ${\overline{S}}_i=-ilog(i/L)$ that corresponds to a diagonal matrix with equal matrix elements. The existence of this upper boundary is the main reason for considering the partial SPDMs, which refer to the fixed number of particles, instead of the total SPDM. In fact, the von Neumann entropy $S\left(t\right)=-\mathrm{Tr}\left[\rho \right(t)log\rho (t\left)\right]$ of the total SPDM Equation (5) depends on the mean number of particles in the chain and, thus, an increase of $S\left(t\right)$ does not necessarily indicate decoherence.

A direct consequence of the observed complete decoherence is an irreversible decay of the mean current $j\left(t\right)$,

(14)$\begin{array}{c}\hfill j\left(t\right)=\mathrm{Tr}\left[\widehat{j}\rho \left(t\right)\right],j_{l,m}=j_0\frac{{\delta }_{l,m-1}-{\delta }_{l-1,m}}{2i},\end{array}$

see red solid line in Figure 3a. We also mention that decay of the mean current is insensitive (at least, on the qualitative level) to variation of the gate voltage $V_g$, see blue dashed and dash-dotted lines in Figure 3a which correspond to $V_g=\pm 0.5$. This is in a strong contrast with the low-temperature limit, where population dynamics and the mean current crucially depend on inequality relation between the gate voltage and the Fermi energy ${ϵ}_F$ which is located at $ϵ=0$ in the considered case of half-filling $N=M/2$. Indeed, in terms of the reservoir eigenstates the low-temperature limit corresponds to $|{\psi }_E⟩$ close to the ground state, where occupation numbers $n_k$ of the natural orbitals show a pronounced step at ${ϵ}_F$. Thus, fermions cannot enter the channel if $V_g>J$, where the whole conductance band lies above the Fermi energy. The results of numerical simulation of the low-temperature limit fully confirm this conjecture, see Figure 3b. Let us also notice the enhanced residual fluctuations of the current as compared to the high-temperature limit.

4. Master Equation Approach

It is interesting to compare the results of Section 3 with solution of the master equation on the reduced density matrix $\mathcal{R}\left(t\right)={\mathrm{Tr}}_b\left[\right|\mathsf{\Psi }\left(t\right)⟩⟨\mathsf{\Psi }\left(t\right)\left|\right]$ for fermionic carriers in the transport channel. Usually, one considers the following equation:

(15)$\begin{array}{c}\hfill \frac{d\mathcal{R}}{dt}=-i[{\widehat{H}}_s,\mathcal{R}]-{\mathcal{L}}_{gain}\left(\mathcal{R}\right)-{\mathcal{L}}_{loss}\left(\mathcal{R}\right),\\ \hfill {\mathcal{L}}_{loss}\left(\mathcal{R}\right)=\frac{\gamma }{2}(1-\overline{n})({\widehat{c}}_1^{†}{\widehat{c}}_1\mathcal{R}-2{\widehat{c}}_1\mathcal{R}{\widehat{c}}_1^{†}+\mathcal{R}{\widehat{c}}_1^{†}{\widehat{c}}_1),\\ \hfill {\mathcal{L}}_{gain}\left(\mathcal{R}\right)=\frac{\gamma }{2}\overline{n}({\widehat{c}}_1{\widehat{c}}_1^{†}\mathcal{R}-2{\widehat{c}}_1^{†}\mathcal{R}{\widehat{c}}_1+\mathcal{R}{\widehat{c}}_1{\widehat{c}}_1^{†}),\end{array}$

where $\overline{n}$ is the filling factor of the reservoir and $\gamma \sim {\epsilon }^2$ is the relaxation constant. (This equation also captures the case of bosonic carries, where the prefactor $(1-\overline{n})$ in the Lindblad term ${\mathcal{L}}_{loss}$ should be replaced with $(1+\overline{n})$ and fermionic annihilation and creation operators with bosonic operators.) It should be stressed that the standard derivation of the displayed master equation assumes a number of approximations [31,32], which have to be verified [33,34]. In this sense Equation (15) implicitly refers to the high-temperature limit and is not valid in the low-temperature limit where, as it was demonstrated in the previous section, the system dynamics depends on inequality relation between the Fermi energy and the gate voltage. Also notice that Equation (15) does not involve ${ϵ}_F$ as a parameter. For this reason from now on we focus on the high-temperature limit where all required assumptions/approximations are believed to be justified.

4.1. Populations and Decoherence Dynamics

It is easy to prove that the matrix $\mathcal{R}$ in Equation (15) has the block structure where each block is associated with the fixed number of fermions in the tight-binding chain of the length L. Using these blocks we calculate the partial SPDMs,

(16)$\begin{array}{c}\hfill {\rho }_{l,m}^{\left(i\right)}\left(t\right)=\mathrm{Tr}\left[{\mathcal{R}}^{\left(i\right)}\left(t\right){\widehat{c}}_l^{†}{\widehat{c}}_m\right],\end{array}$

and then use them to calculate probabilities $P_i\left(t\right)$ to find i fermions in the transport channel and von Neumann entropies $S_i\left(t\right)$, which characterize quantum state of these fermions. The results are presented in Figure 4, which should be compared with Figure 2. We notice that, when solving Equation (15), we take into account depletion of the reservoir, i.e., the parameter $\overline{n}$ is decreased in time according to the depletion dynamics,

(17)$\begin{array}{c}\hfill \overline{n}\left(t\right)=(N-N_s\left(t\right))/M.\end{array}$

with this minor modification one finds very good agreement between the master equation approach and the exact numerical results. This agreement indicates that all assumptions/approximations used to derive Equation (15) are indeed justified. This allows us to address within the framework of the master equation the more complex problem, where the transport channel connects two high-temperature reservoirs with different filling factors.

4.2. Stationary Current between Two Reservoirs

To take into account the second reservoir the master Equation (15) on the reduced density matrix should be complimented by two additional Lindblad terms which has the same structure as the Lindblad terms in Equation (15) but involves operators ${\widehat{c}}_L^{†}$ and ${\widehat{c}}_L$ instead of the operators ${\widehat{c}}_1^{†}$ and ${\widehat{c}}_1$. From this point we redefine the parameters $\overline{n}$ and $\gamma$ as ${\overline{n}}_L$ and ${\gamma }_L$ (the left reservoir). Correspondently, the filling factor and relaxation constant of the right reservoir are denoted by ${\overline{n}}_R$ and ${\gamma }_R$ and, to be certain, we assume ${\overline{n}}_L>{\overline{n}}_R$.

The solution of the described master equation with the source and sink terms was discussed in Reference [19] for the case of bosonic carriers. Adopting the results of Reference [19] to the present case of fermionic carries we come to the following conclusions. In course of time SPDM relaxes to the three-diagonal matrix where (pure imaginary) off-diagonal elements of the matrix determine the stationary current $\overline{j}$ between the left and right reservoirs. This current is proportional to difference in the reservoir filling factors, where the proportionality coefficient A has particularly simple form in the case ${\gamma }_L={\gamma }_R\equiv \gamma$,

(18)$\begin{array}{c}\hfill A=\frac{1}{2}\frac{J\gamma }{J^2+{\gamma }^2},\end{array}$

and in the case ${\gamma }_R\ll {\gamma }_L\equiv \gamma$,

(19)$\begin{array}{c}\hfill A=\frac{{\gamma }_R\gamma }{J^2+{\gamma }^2}.\end{array}$

The latter case is of special interest for the purpose of microscopic analysis of the system dynamics. Indeed, if ${\gamma }_L\gg {\gamma }_R$ then the main source of decoherence is the left reservoir while the right reservoir barely serves as a particle sink. In the next subsection we analyze the problem where the left reservoir is modelled microscopically while the right reservoir is taken into account by using the master equation approach.

4.3. Quasi-Stationary Current

Following the discussion in the previous subsection, we consider the master equation

(20)$\begin{array}{c}\hfill \frac{d\mathcal{R}}{dt}=-i[\widehat{H},\mathcal{R}]-{\mathcal{L}}_{loss}\left(\mathcal{R}\right),\\ \hfill {\mathcal{L}}_{loss}\left(\mathcal{R}\right)=\frac{{\gamma }_R}{2}({\widehat{c}}_L^{†}{\widehat{c}}_L\mathcal{R}-2{\widehat{c}}_L\mathcal{R}{\widehat{c}}_L^{†}+\mathcal{R}{\widehat{c}}_L^{†}{\widehat{c}}_L),\end{array}$

where $\widehat{H}$ is the Hamiltonian of the left reservoir with the attached transport channel. Due to large dimension of the Hamiltonian we solve Equation (20) by using the stochastic approach [32]. Specifically, we solve the Schrödinger equation of the form [35]

(21)$\begin{array}{c}\hfill \mathrm{d}|\mathsf{\Psi }⟩=\left(-i\widehat{H}\mathrm{d}t-\frac{{\gamma }_R}{2}{\widehat{c}}_L^{†}{\widehat{c}}_L\mathrm{d}t+\sqrt{{\gamma }_R}{\widehat{c}}_L\mathrm{d}\xi \right)|\mathsf{\Psi }⟩,\end{array}$

where $\mathrm{d}\xi$ is the Wiener process with $\overline{\mathrm{d}\xi }=0$ and ${\overline{\mathrm{d}\xi }}^2=\mathrm{d}t$. Within this approach the reduced density matrix $\mathcal{R}\left(t\right)$ is found by averaging the solution of Equation (21) over different realisations of the stochastic process, i.e., $\mathcal{R}\left(t\right)=\overline{\left|\mathsf{\Psi }\right(t)⟩⟨\mathsf{\Psi }(t\left)\right|}$. The convergence of the averaging procedure is controlled against the condition $\mathrm{Tr}\left[\mathcal{R}\right(t\left)\right]=1$.

First we reproduce the result of Figure 3a. The lower solid line in Figure 5 shows the mean current in the transport channel for $\epsilon =0.1$ but slightly smaller system size $M=10$, $L=4$ (this reduces the dimension of the Hilbert space from $\mathcal{N}=\mathrm{18,564}$ to $\mathcal{N}=3003$) and ${\gamma }_R=0$. The exponential decay of the current is clearly seen. Next we set ${\gamma }_R$ to a small value ${\gamma }_R=0.04$. It is seen that $j\left(t\right)$ now decays to a finite value $\overline{j}$, i.e., we have a quasi-stationary current between the reservoirs. We stress that the observed rapid relaxation of the current to zero or finite value is exclusively due to decoherence effect of the left reservoir. In fact, matching the lower solid line to the solution of the master Equation (15) we find ${\gamma }_L=0.24$. Thus we are indeed in the regime ${\gamma }_R\ll {\gamma }_L$ where one can neglect decoherence effect of the right reservoir.

5. Conclusions

We analyzed the quantum state of fermionic carriers in the transport channel connecting two reservoirs. The analysis is done from first principles by considering a microscopic model of the reservoir (Two-Body Random Interaction Model) and the transport channel (tight-binding chain of a finite length). In the case of infinite effective temperature of the reservoirs the single-particle density matrix (SPDM) of fermions in the channel is shown to relax to a three-diagonal matrix, whose off-diagonal elements determine the stationary current between the reservoirs. We stress that relaxation of SPDM to this steady state is entirely due to decoherence effect of the reservoirs on the carriers propagating in the channel. We obtain explicit expressions for the stationary current by justifying the master equation on the reduced density matrix of the carriers, which fortunately can be solved analytically. This first-principle justification of the master equation is one of the main results of the work.

The main challenge in the context of the presented studies is the case of low reservoir temperature, where occupation numbers of its natural orbitals show a step at the Fermi energy. It is believed that in this case the quantum state of fermionic carries in the transport channel is close to the Bloch wave with $k_F=arccos[({ϵ}_F-V_g)/J]$. In the other words, the stationary SPDM has many non-zero diagonals. In principle, one can prove or disprove this conjecture numerically within the framework of the discussed microscopic model by considering a larger system size [that will reduce residual fluctuations in Figure 3b]. The other route is a generalization of the master Equation (15) onto the case of finite reservoir temperature, where it should include $k_F$ as an additional parameter.

Author Contributions

Conceptualization, A.R.K.; investigation, A.R.K. and D.N.M.; writing–original draft preparation, A.R.K.; writing–review and editing, D.N.M.

Funding

This work has been supported by Russian Science Foundation through grant N19-12-0016. We appreciate discussion with Anna A. Bychek.

Conflicts of Interest

The authors declare no conflict of interest.

Figures

Figure 1. Population dynamics of the reservoir orbitals (left to the vertical dashed line) and the lattice sites (right to the dashed line) for infinite effective temperature of the reservoir. The system size is M=12, L=6, and N=6. The value of the hopping matrix element J=0.5, the coupling constant ε=0.1.

Figure 2. Probabilities Pi to find i fermions in the transport channel (upper panel) and von Neumann entropies Si of the normalized partial single-particle density matrices (SPDMs) (lower panel) as the functions of time. The different curves correspond to i=1 (blue), 2 (green), 3 (red), 4 (light blue), and 5 (magenta). Probability to find 6 fermions is close to zero and is not shown, as well as the probability to find 0 fermions, which is given by the equation P0=1−∑i=1LPi.

Figure 3. The total current in the transport channel as the function of time in the high-temperature (a) and low-temperature (b) limits. The dash-dotted, solid, and dashed lines correspond to different values of the gate voltage Vg=−0.5,0,0.5, respectively.

Figure 4. The same quantities as in Figure 2 yet calculated by using the master equation approach. The value of the relaxation constant is adjusted to γ=0.4.

Figure 5. The mean current in the transport channel connecting two reservoirs. Dotted lines show solution of the master equation with the source (n¯L=0.5, γL=0.24) and sink (n¯R=0, γR=0,0.04,0.08,0.12 from bottom to top) terms. The solid lines are solution of the master Equation (20) for γR=0 and γR=0.04.