## Landauer-Büttiker and Thouless conductance

L. Bruneau , V. Jakši´ c 1 2 , Y. Last 3 , C.-A. Pillet 4

1 Département de Mathématiques and UMR 8088 CNRS and Université de Cergy-Pontoise 95000 Cergy-Pontoise, France

2 Department of Mathematics and Statistics McGill University 805 Sherbrooke Street West Montreal, QC, H3A 2K6, Canada

3 Institute of Mathematics The Hebrew University 91904 Jerusalem, Israel

4 Université de Toulon, CNRS, CPT, UMR 7332, 83957 La Garde, France Aix-Marseille Université, CNRS, CPT, UMR 7332, 13288 Marseille, France

Dedicated to the memory of Markus Büttiker

Abstract. In the independent electron approximation, the average (energy/charge/entropy) current flowing through a finite sample S connected to two electronic reservoirs can be computed by scattering theoretic arguments which lead to the famous Landauer-Büttiker formula. Another well known formula has been proposed by Thouless on the basis of a scaling argument. The Thouless formula relates the conductance of the sample to the width of the spectral bands of the infinite crystal obtained by periodic juxtaposition of S . In this spirit, we define Landauer-Büttiker crystalline currents by extending the Landauer-Büttiker formula to a setup where the sample S is replaced by a periodic structure whose unit cell is S . We argue that these crystalline currents are closely related to the Thouless currents. For example, the crystalline heat current is bounded above by the Thouless heat current, and this bound saturates iff the coupling between the reservoirs and the sample is reflectionless. Our analysis leads to a rigorous derivation of the Thouless formula from the first principles of quantum statistical mechanics.

## 1 Introduction

Recent rigorous formulations of the Landauer-Büttiker formula [AJPP, BSP, CJM, N] in the framework of nonequilibrium quantum statistical mechanics have opened the way to the mathematical study of a variety of related transport phenomena in quantum mechanics [GJW, JLPa, JLPi]. This paper is a continuation of this line of research. Our main goal here is to provide a mathematically rigorous proof of the celebrated Thouless conductance formula.

This paper is organized as follows. In Section 1.1 we review the Electronic Black Box Model and the corresponding Landauer-Büttiker formula. These topics have been discussed from both technical and pedagogical point of view in [AJPP, BSP] and the reader may consult these works for additional information. The notions of Thouless energy and Thouless conductance are reviewed in Section 1.2. Our main results are stated in Section 2. The proofs are given in Section 3.

Acknowledgment. The research of V.J. was partly supported by NSERC. The research of Y.L. was partly supported by The Israel Science Foundation (Grant No. 1105/10) and by Grant No. 2010348 from the United States-Israel Binational Science Foundation (BSF), Jerusalem, Israel. A part of this work has been done during a visit of L.B. to McGill University supported by NSERC. Another part was done during the visit of V.J. to The Hebrew University supported by ISF.

## 1.1 The Electronic Black Box Model and the Landauer-Büttiker formula

The Electronic Black Box Model (abbreviated EBBM) describes a finite sample S connected to two infinitely extended electronic reservoirs R l/r (where l/r stands for left/right) in the independent electron approximation. The coupling between the sample and the reservoirs allows for a flow of energy/charge/entropy through the joint system R l + S + R r .

We shall restrict our attention to one-dimensional samples in the tight binding approximation. Thus, the sample S is a free Fermi gas with one-particle Hilbert space h S = glyph[lscript] 2 ( Z L ) , where Z L = [1 , L ] ∩ Z is a finite lattice. Its Hamiltonian h S is a Jacobi matrix with parameters { J x } 1 ≤ x&lt;L , { λ x } 1 ≤ ≤ x L ,

$$( h _ { \mathcal { S } } u ) ( x ) = J _ { x } u ( x + 1 ) + J _ { x - 1 } u ( x - 1 ) + \lambda _ { x } u ( x ), \quad x \in Z _ { L }, \quad \ \ ( 1. 1 )$$

and Dirichlet boundary condition u (0) = u L ( +1) = 0 . The reservoir R l/r is a free Fermi gas with one-particle Hilbert space h l/r and one-particle Hamiltonian h l/r . The one-particle Hilbert space and Hamiltonian of the composite system R l + S + R r are

$$\mathfrak { h } = \mathfrak { h } _ { l } \oplus \mathfrak { h } _ { \mathcal { S } } \oplus \mathfrak { h } _ { r },$$

$$h _ { 0 } = h _ { l } \oplus h _ { \mathcal { S } } \oplus h _ { r }.$$

The coupling of the sample with the reservoir R l/r is realized by the hopping Hamiltonian

$$v _ { l / r } = | \chi _ { l / r } \rangle \langle \psi _ { l / r } | + | \psi _ { l / r } \rangle \langle \chi _ { l / r } |,$$

where χ l/r ∈ h l/r is a unit vector while ψ l = δ 1 and ψ r = δ L are Kronecker delta functions in h S . The one-particle Hamiltonian of the coupled system is

$$h = h _ { 0 } + \kappa v = h _ { 0 } + \kappa ( v _ { l } + v _ { r } ),$$

glyph[negationslash]

where κ = 0 is the coupling strength. As observed in [BJP], for the purposes of discussing transport properties of the coupled system R l + S + R r one may assume, without loss of generality, that χ l/r is a cyclic vector for h l/r . Hence, passing to the spectral representation we may assume that h l/r acts as multiplication by E on

$$\mathfrak { h } _ { l / r } = L ^ { 2 } ( \mathbb { R }, \mathrm d \nu _ { l / r } ( E ) ),$$

where ν l/r is the spectral measure of h l/r associated to χ l/r . Moreover, in this representation one has χ l/r ( E ) = 1 and ν l/r ( R ) = 1 .

Let Γ ( ) -h be the fermionic Fock space over h and denote by a f ( ) /a ∗ ( f ) the annihilation/creation operator on Γ ( ) -h associated to f ∈ h . The algebra of observables of the EBBM is the C ∗ -algebra O of all bounded operators on Γ ( ) -h generated by the identity I and { a ∗ ( f ) a g ( ) | f, g ∈ h } .

Let H = dΓ( h ) be the second quantized Hamiltonian. The map

$$\tau ^ { t } ( a ^ { * } ( f ) a ( g ) ) = e ^ { i t H } a ^ { * } ( f ) a ( g ) e ^ { - i t H } = a ^ { * } ( e ^ { i t h } f ) a ( e ^ { i t h } g ),$$

extends to a strongly continuous group of ∗ -automorphisms of O and the pair ( O , τ t ) is a C ∗ -dynamical system. The states of the EBBM are normalized positive linear functionals on O .

Since we are dealing with independent electrons, quasi-free states on O will be of particular relevance. Let ρ be a self-adjoint operator of h such that 0 ≤ ρ ≤ I . The gauge-invariant quasi-free state of density ρ is the unique state ω ρ on O satisfying

$$\omega _ { \rho } ( a ^ { * } ( f _ { 1 } ) \cdots a ^ { * } ( f _ { n } ) a ( g _ { m } ) \cdots a ( g _ { 1 } ) ) = \delta _ { n, m } \det \{ \langle g _ { i }, \rho f _ { j } \rangle \},$$

for any integers n, m and all f , . . . , f 1 n , g 1 , . . . , g m ∈ h .

The initial state ω 0 of the EBBM is the gauge-invariant quasi-free state of density

$$\rho _ { 0 } = \rho _ { l } \oplus \rho _ { \mathcal { S } } \oplus \rho _ { r },$$

where ρ l/r is the operator of multiplication by the Fermi-Dirac density

$$\rho _ { l / r } ( E ) = \frac { 1 } { 1 + \mathrm e ^ { \beta _ { l / r } ( E - \mu _ { l / r } ) } }.$$

In other words, the reservoir R l/r is initially in thermal equilibrium at inverse temperature β l/r &gt; 0 and chemical potential µ l/r ∈ R . None of our results depends on the particular choice of the initial density ρ S of the sample.

The EBBM is the quantum dynamical system ( O , τ t , ω 0 ) .

The basic questions regarding this system concern its behavior in the large time limit t → ∞ . To deal with this limit we observe that the initial state ω 0 = ω ρ 0 evolves in time as

$$\omega _ { t } = \omega _ { 0 } \circ \tau ^ { t } = \omega _ { \rho _ { t } },$$

where the density at time t is given by

$$\rho _ { t } = e ^ { - i t h } \rho _ { 0 } e ^ { i t h }.$$

One expects that ω t → ω + as t → ∞ , where the non-equilibrium steady state ω + carries energy/charge/entropy current induced by the initial temperature/chemical potential differential. These currents can be expressed by Landauer-Büttiker formulas involving transmission properties of individual electrons evolving under the dynamics generated by the one-particle Hamiltonian h . These heuristics is mathematically formalized as follows.

The observables describing energy and charge current out of R l/r are

$$\Phi _ { l / r } = \mathrm d \Gamma ( - \mathrm i [ h, h _ { l / r } ] ) = \kappa \left ( a ^ { * } ( \mathrm i h _ { l / r } \chi _ { l / r } ) a ( \psi _ { l / r } ) + a ^ { * } ( \psi _ { l / r } ) a ( \mathrm i h _ { l / r } \chi _ { l / r } ) \right ),$$

$$\mathcal { I } _ { l / r } = \mathrm d \Gamma ( - \mathrm i [ h, 1 _ { l / r } ] ) = \kappa \left ( a ^ { * } ( \mathrm i \chi _ { l / r } ) a ( \psi _ { l / r } ) + a ^ { * } ( \psi _ { l / r } ) a ( \mathrm i \chi _ { l / r } ) \right ).$$

The entropy current associated to the heat flux dissipated into the reservoirs is

$$\mathcal { J } = - \beta _ { l } ( \Phi _ { l } - \mu _ { l } \mathcal { I } _ { l } ) - \beta _ { r } ( \Phi _ { r } - \mu _ { r } \mathcal { I } _ { r } ).$$

The transmittance of the sample is the function 1

$$\mathbb { R } \ni E \mapsto \mathcal { T } ( E ) = 4 \kappa ^ { 4 } | \langle \psi _ { l }, ( h - E - i 0 ) ^ { - 1 } \psi _ { r } \rangle | ^ { 2 } \, \text{Im} \, F _ { l } ( E ) \, \text{Im} \, F _ { r } ( E ), \quad \quad ( 1. 2 )$$

↦

where

$$F _ { l / r } ( E ) = \langle \chi _ { l / r }, ( h _ { l / r } - E - \mathrm i 0 ) ^ { - 1 } \chi _ { l / r } \rangle. & & ( 1. 3 ) \\ \omega _ { \cdot \cdot \cdot } / E \searrow \infty \ e n d \, ^ { 2 }$$

Note that Im F l/r ( E ) ≥ 0 and 2

$$\{ E \in \mathbb { R } \, | \, \mathcal { T } ( E ) > 0 \} = \Sigma _ { l } \cap \Sigma _ { r },$$

where

$$\Sigma _ { l / r } = \left \{ E \in \mathbb { R } \, | \, \text{Im} \, F _ { l / r } ( E ) > 0 \right \} \\ \text{not set to the absolute} \, \text{contain} \, \text{contain} \, \text{$\epsilon_{k}$} \,.$$

is the essential support of the absolutely continuous spectrum of h l/r .

For E ∈ R we set

$$\zeta _ { l / r } ( E ) = \beta _ { l / r } ( E - \mu _ { l / r } ), \quad \Delta _ { l / r } ( E ) = \rho _ { l / r } ( E ) - \rho _ { r / l } ( E ),$$

and

$$\varsigma ( E ) = ( \zeta _ { r } ( E ) - \zeta _ { l } ( E ) ) \Delta _ { l } ( E ) = ( \zeta _ { l } ( E ) - \zeta _ { r } ( E ) ) \Delta _ { r } ( E ). \\ \intertext { t i c r } \frac { \lambda } { R } = \frac { R } { R } \text{ and }.. \text{ } -.. \text{ } \text{ then } \text{$\sim E\uplus$ combination in both that it is attribute, no one}.$$

Our basic assumption in the following is:

Note that if β l = β r and µ l = µ r then ς ( E ) vanishes identically but that it is strictly positive for all E ∈ R otherwise.

Assumption A. The one-particle Hamiltonian h has no singular continuous spectrum.

The starting point of this paper is the following result [AJPP].

Theorem 1.1 Suppose that Assumption A holds. Then for all A ∈ O the limit

$$\omega _ { + } ( A ) = \lim _ { t \rightarrow \infty } \frac { 1 } { t } \int _ { 0 } ^ { t } \omega _ { s } ( A ) d s,$$

exists and defines a state ω + on O . Moreover, the following formulas for the average steady currents hold:

$$\langle \Phi _ { l / r } \rangle _ { + } = \omega _ { + } ( \Phi _ { l / r } ) = \frac { 1 } { 2 \pi } \int _ { \mathbb { R } } \mathcal { T } ( E ) E \Delta _ { l / r } ( E ) \mathrm d E,$$

$$\langle \mathcal { I } _ { l / r } \rangle _ { + } = \omega _ { + } ( \mathcal { I } _ { l / r } ) = \frac { 1 } { 2 \pi } \int _ { \mathbb { R } } \mathcal { T } ( E ) \Delta _ { l / r } ( E ) \mathrm d E,$$

$$\langle \mathcal { J } \rangle _ { + } = \omega _ { + } ( \mathcal { J } ) = \frac { 1 } { 2 \pi } \int _ { \mathbb { R } } \mathcal { T } ( E ) \varsigma ( E ) \text{d} E.$$

1 This function is defined for Lebesgue a.e. E ∈ R .

2 This identity is understood modulo sets of Lebesgue measure zero.

We finish this section with a number of remarks regarding this result.

Remark 1. Theorem 1.1 deals with the simplest non-trivial setting in the study of electronic transport in the independent electron approximation. Various generalizations of this setting and of Theorem 1.1 can be found in [AJPP, N, BSP].

Remark 2. We have chosen units in such a way that the electronic charge e , the reduced Planck constant glyph[planckover2pi1] , and the Boltzmann constant k B are unity. The energy, charge, and entropy currents in the above formulas are expressed in units of 1 / glyph[planckover2pi1] , e/ glyph[planckover2pi1] and k B / glyph[planckover2pi1] . Note also that these formulas do not include the spin degeneracy which should be accounted for by a factor 2 .

Remark 3. If Assumption A is replaced by the stronger assumption that h has no singular spectrum, then

$$\omega _ { + } ( A ) = \lim _ { t \rightarrow \infty } \omega _ { t } ( A ),$$

holds for all A ∈ O .

$$\rho _ { \beta, \mu } = \frac { 1 } { 1 + \mathrm e ^ { \beta ( h - \mu ) } }.$$

Remark 4. If β l = β r = β and µ l = µ r = µ , then ω + is the thermal equilibrium state of the coupled system R l + S + R r for the given intensive thermodynamic parameters, i.e., the quasi-free state of density glyph[negationslash]

In this case, all currents vanish in average. In the following we shall exclude this possibility and assume that β l = β r or/and µ l = µ r . The state ω + is then a non-equilibrium steady state of the EBBM and (1.6) are the Landauer-Büttiker formulas for the steady state currents.

glyph[negationslash]

Remark 5. Apart from the choice of the intensive thermodynamic parameters β l/r and µ l/r , the transmission coefficient T ( E ) completely determines the steady state currents. It follows from (1.4) that the steady state currents are non-vanishing iff the Lebesgue measure | Σ l ∩ Σ r | of the set Σ l ∩ Σ r is strictly positive, i.e., iff there exists an open scattering channel between the left and the right reservoir. Note the obvious conservation laws

$$\omega _ { + } ( \Phi _ { l } ) + \omega _ { + } ( \Phi _ { r } ) = 0, \quad \omega _ { + } ( \mathcal { I } _ { l } ) + \omega _ { + } ( \mathcal { I } _ { r } ) = 0.$$

Entropy balance in the steady state ω + implies that

$$\langle \mathcal { J } \rangle _ { + } = - \beta _ { l } ( \langle \Phi _ { l } \rangle _ { + } - \mu _ { l } \langle \mathcal { I } _ { l } \rangle _ { + } ) - \beta _ { r } ( \langle \Phi _ { r } \rangle _ { + } - \mu _ { r } \langle \mathcal { I } _ { r } \rangle _ { + } ), \\ \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{cases}$$

coincides with the rate of entropy production in S [AJPP]. It follows from the above observation that 〈J 〉 + ≥ 0 , and that the inequality is strict iff | Σ l ∩ Σ r | &gt; 0 .

Remark 6. The proof of Theorem 1.1 is based on the scattering theory of the pair ( h, h 0 ) and elucidates the physical meaning of T ( E ) . It follows from the trace class scattering theory that the wave operators 3

$$w _ { \pm } = s - \lim _ { t \to \pm \infty } e ^ { i t h } e ^ { - i t h _ { 0 } } 1 _ { \text{ac} } ( h _ { 0 } ),$$

3 1 ac ( h 0 ) is the spectral projection on the absolutely continuous part of the spectrum of h 0 .

exist and are complete. The scattering matrix s = w w ∗ + -is a unitary operator on

$$\mathfrak { h } _ { \text{ac} } ( h _ { 0 } ) = \text{Ran} \, 1 _ { \text{ac} } ( h _ { 0 } ) = \text{Ran} \, 1 _ { \text{ac} } ( h _ { l } ) \oplus \text{Ran} \, 1 _ { \text{ac} } ( h _ { r } ),$$

and acts as the operator of multiplication by a unitary 2 × 2 matrix

$$s ( E ) = \begin{bmatrix} s _ { l l } ( E ) & s _ { l r } ( E ) \\ s _ { r l } ( E ) & s _ { r r } ( E ) \end{bmatrix}.$$

One then has

$$\mathcal { T } ( E ) = | s _ { l r } ( E ) | ^ { 2 } = | s _ { r l } ( E ) | ^ { 2 }, \\ \mathbb { u } \colon = \mathbb { m } \cdots \dots \dots \dots \dots \dots \dots \dots$$

i.e., the transmittance T ( E ) is the transmission probability between the left and right reservoir at energy E .

Remark 7. An additional insight into T ( E ) can be obtained by taking into account the spatial structure of the reservoirs. Suppose that the reservoirs are non-trivial in the sense that the measures ν l/r have non-vanishing absolutely continuous component. The standard orthogonal polynomial construction (see Theorem I.2.4 in [Si]) provides a unitary operator U : h → glyph[lscript] 2 ( Z ) such that the following holds:

$$( i ) \ U \mathfrak { h } _ { l } = \ell ^ { 2 } ( \left ] - \infty, 0 \right ] \cap \mathbb { Z } ), U \mathfrak { h } _ { \mathcal { S } } = \mathfrak { h } _ { \mathcal { S } } \text{ and } U \mathfrak { h } _ { r } = \ell ^ { 2 } ( [ L + 1, \infty [ \cap \mathbb { Z } ).$$

- (ii) There is a Jacobi matrix on glyph[lscript] 2 (] - ∞ , 0] ∩ Z ) , with parameters { J x } x&lt; 0 , { λ x } x ≤ 0 and Dirichlet boundary condition such that for u ∈ glyph[lscript] 2 (] -∞ , 0] ∩ Z )

$$( U h _ { l } U ^ { * } u ) ( x ) = J _ { x } u ( x + 1 ) + J _ { x - 1 } u ( x - 1 ) + \lambda _ { x } u ( x ), \quad ( x \leq 0, u ( 1 ) = 0 ).$$

(iii) Uh U S ∗ = h S .

- (iv) There is a Jacobi matrix on glyph[lscript] 2 ([ L + 1 , ∞∩ [ Z ) , with parameters { J x } x&gt;L , { λ x } x&gt;L and Dirichlet boundary condition such that for u ∈ glyph[lscript] 2 ([ L +1 , ∞∩ [ Z )

$$( U h _ { r } U ^ { * } u ) ( x ) = J _ { x } u ( x + 1 ) + J _ { x - 1 } u ( x - 1 ) + \lambda _ { x } u ( x ), \quad ( x > L, u ( L ) = 0 ).$$

$$( \nu ) \ U \chi _ { l } = \delta _ { 0 } \text{ and } U \chi _ { r } = \delta _ { L + 1 }.$$

$$( v i ) \text{ if } J _ { 0 } = J _ { L } = \kappa, \text{ then for } u \in \ell ^ { 2 } ( \mathbb { Z } )$$

$$( U h U ^ { * } u ) ( x ) = J _ { x } u ( x + 1 ) + J _ { x - 1 } u ( x - 1 ) + \lambda _ { x } u ( x ).$$

It follows that our EBBM is unitarily equivalent to a Jacobi matrix EBBM on glyph[lscript] 2 ( Z ) .

For z ∈ C + let u l/r ( · , z ) be the unique solution of the equation

$$J _ { x } u _ { l / r } ( x + 1, z ) + J _ { x - 1 } u _ { l / r } ( x - 1, z ) + \lambda _ { x } u _ { l / r } ( x, z ) = z u _ { l / r } ( x, z ),$$

that is square summable at ±∞ and normalized by u l/r (0 , z ) = 1 . For all x ∈ Z and Lebesgue a.e. E ∈ R the limit

$$\lim _ { \epsilon \downarrow 0 } u _ { l / r } ( x, E + \text{i} \epsilon ) = u _ { l / r } ( x, E ),$$

exists, is finite, and solves (1.7) with z = E . For Lebesgue a.e. E ∈ Σ r the solution u r ( E ) is not a multiple of a real solution and so u r ( E ) is also a solution of (1.7) linearly independent of u r ( E ) . Hence, for Lebesgue a.e. E ∈ Σ r one has

$$u _ { l } ( E ) = \alpha ( E ) \overline { u _ { r } ( E ) } + \beta ( E ) u _ { r } ( E ).$$

The spectral reflection probability of [GNP, GS] is

$$R _ { r } ( E ) = \left | \frac { \beta ( E ) } { \alpha ( E ) } \right | ^ { 2 }.$$

glyph[negationslash]

One extends R r to R by setting R r ( E ) = 1 for E ∈ Σ r . R E l ( ) is defined analogously. A simple computation (see Section 5 in [JLPa]) gives

$$R _ { r } ( E ) = R _ { l } ( E ) = | s _ { l l } ( E ) | ^ { 2 } = | s _ { r r } ( E ) | ^ { 2 }.$$

Hence,

$$\mathcal { T } ( E ) = 1 - R _ { r } ( E ) = 1 - R _ { l } ( E ), \\ \colon \dots \dots \colon \colon \dots \colon \dots \colon \dots \colon \colon \dots \colon \colon \dots \colon \colon \colon \dots \colon \colon \colon \colon \dots \colon \colon \colon \colon \dots \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{matrix}$$

is the spectral transmission probability at the energy E .

## 1.2 The Thouless energy and the Thouless conductance formula

We start with a review of the notions of Thouless energy and conductance as discussed in the physics literature. We follow [La] and our main goal is to extract mathematically well-defined quantities that correspond to these heuristic notions.

## 1.2.1 Heuristics

Let δt be the typical time spent in S by an electron on its journey from one reservoir to the other. The Mandelstam-Tamm time-energy uncertainty relation δE δt glyph[greaterorsimilar] 1 ([MT], see also [FP, Bu]) imposes a lower bound on the energy spread δE of the electron wave function. This sets the Thouless energy scale E Th ∼ δE [Th]. Assuming a diffusive behavior, one has L 2 ∼ Dδt where L is the sample size and D the diffusion constant, and hence

$$E _ { \text{Th} } \gtrsim \frac { D } { L ^ { 2 } }.$$

Einstein's relation

$$\sigma = D \varrho \lesssim L ^ { 2 } E _ { \text{Th} } \varrho,$$

further links D to the conductivity σ of the sample and its density of states glyph[rho1] . Finally, glyph[rho1] relates to the typical energy level spacing ∆ E of the sample as

$$\varrho L \Delta E \sim 1.$$

It follows that

$$\sigma \lesssim L \frac { E _ { \text{Th} } } { \Delta E },$$

and hence that the conductance g = L -1 σ of the one-dimensional sample satisfies

$$g \lesssim g _ { \text{Th} } = \frac { E _ { \text{Th} } } { \Delta E }.$$

The quantity at the r.h.s. of this formula is known as the Thouless conductance. Although the above derivation is extremely heuristic in nature, the Thouless conductance g Th and the closely related Thouless energy E Th are widely accepted by both theoretical and experimental physicists and play an important role in the scaling theory of localization [AALR].

The above argument and the resulting inequality on the l.h.s. of (1.8) suggest that we should consider g Th as an upper bound on the conductance of the sample S . Indeed, for such a microscopic system the conductance is not an intrinsic property of the sample, but depends also on the reservoirs and the nature of the coupling, which determines the reservoir's ability to feed the available energy levels of the sample. Thus, one expects that saturation of the Thouless bound (1.8) occurs for optimal feeding of the sample by the reservoirs, a property of the joint system R l + S + R r which we shall try to elucidate in the remaining part of this paper (see Remarks 1 and 4 after Theorem 2.3).

## 1.2.2 The crystal model

There is a simple way to make the sample transparent to incoming electrons and hence to ensure its optimal feeding: it suffices to implement the reservoirs in such a way that the joint system R l + S + R r is periodic. We shall call crystalline EBBM the model obtained by repeating the sample so as to obtain a periodic crystal with unit cell S (a construction closely related to the scaling argument of [ET]).

Consider the periodic Jacobi matrix h crystal on glyph[lscript] 2 ( Z ) obtained by extending the Jacobi parameters { λ x } 1 ≤ ≤ x L and { J x } 1 ≤ x&lt;L of the sample Hamiltonian h S to the entire lattice Z by setting J L = κ S and

$$J _ { x + n L } = J _ { x }, \quad \lambda _ { x + n L } = \lambda _ { x },$$

glyph[negationslash]

for any n ∈ Z and x ∈ Z L . The internal coupling constant κ S is a priori an arbitrary parameter, except for the obvious constraint κ S = 0 . In practice it will be determined by the physics of the problem. In a model where the sample parameters J x are independent copies of a random

variable, κ S will be another instance of this variable. If h S is a discrete Schrödinger operator with J x = J for all x ∈ Z L the choice of κ S = J appears natural.

In the crystaline EBBM model, the single particle Hilbert spaces of the reservoirs are h l = glyph[lscript] 2 (] -∞ , 0] ∩ Z ) and h r = glyph[lscript] 2 ([ L +1 , ∞∩ [ Z ) , the corresponding single particle Hamiltonians are Jacobi matrices with parameters {{ J x } x&lt; 0 , { λ x } x&lt; 0 } , {{ J x } x&gt;L , { λ x } x&gt;L } and Dirichlet boundary condition, χ l = δ 0 , χ r = δ L +1 , and the coupling constant is set to κ = κ S . The one particle Hamiltonian of the coupled system is h crystal .

We emphasize that the crystallization of the sample is specified by the pair ( S , κ S ) and not by S alone.

The Bloch-Floquet decomposition of h crystal reads

$$\ell ^ { 2 } ( \mathbb { Z } ) = \int _ { \mathfrak { B } _ { L } } ^ { \oplus } \mathfrak { h } _ { \mathcal { S } } \, \mathrm d k, \quad h _ { \text{crystal} } = \int _ { \mathfrak { B } _ { L } } ^ { \oplus } h ( k ) \, \mathrm d k,$$

where B L = [ -π/L, π/L ] is the first Brillouin zone of the crystal and h k ( ) is the self-adjoint operator on h S obtained from (1.1) by replacing Dirichlet by Bloch boundary conditions

$$u ( 0 ) = \mathrm e ^ { - \mathrm i k L } u ( L ), \quad u ( L + 1 ) = \mathrm e ^ { \mathrm i k L } u ( 1 ).$$

The spectrum of h k ( ) consists of L eigenvalues ε 1 ( k ) ≤ · · · ≤ ε L ( k ) which are even functions of k , real analytic and strictly monotone on ]0 , π/L [ . Moreover

$$\varepsilon _ { L } ( 0 ) > \varepsilon _ { L } ( \pi / L ) \geq \varepsilon _ { L - 1 } ( \pi / L ) > \varepsilon _ { L - 1 } ( 0 ) \geq \varepsilon _ { L - 2 } ( 0 ) > \cdots$$

Thus, the spectrum of h crystal is

$$\text{sp} ( h _ { \text{crystal} } ) = \bigcup _ { k \in \mathcal { B } _ { L } } \text{sp} ( h ( k ) ) = \bigcup _ { j = 1 } ^ { L } B _ { j },$$

where B j is a closed interval with boundary points ε j (0) and ε j ( π/L ) (see Theorem 5.3.4 in [Si] and Figure 1).

For simplicity, suppose that the reservoirs are at zero temperature, i.e., β l = β r = ∞ , and assume that µ l &lt; µ r . In this case the function ∆ r is the characteristic function of the interval [ µ , µ l r ] and the Landauer-Büttiker formula for the steady charge current out of the right reservoir becomes

$$\omega _ { + } ( \mathcal { J } _ { r } ) = \frac { 1 } { 2 \pi } \int _ { \mu _ { l } } ^ { \mu _ { r } } \mathcal { T } ( E ) \mathrm d E.$$

Since the transmittance of a unit cell in a perfect crystal is the characteristic function of its spectrum, we get

$$\omega _ { + } ( \mathcal { J } _ { r } ) = \frac { 1 } { 2 \pi } | \text{sp} ( h _ { \text{crystal} } ) \cap [ \mu _ { l }, \mu _ { r } ] |,$$

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0185v1-figure-a01fae93-ac6a-440c-b12c-b1fbf237b933.jpg)

-

Figure 1: Eigenvalues of h k ( ) and the band spectrum of h crystal .

from which we infer that the mean conductance of the sample on the energy window I = [ µ , µ l r ] is given by

$$g ( I ) = \frac { 1 } { 2 \pi } \frac { | \text{sp} ( h _ { \text{crystal} } ) \cap I | } { | I | }. \\ \cdot \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

We shall now argue that this expression can be interpreted as the Thouless conductance associated to the energy window I . In order for the level spacing ∆ E to be well defined, the interval I should contain several bands B j of h crystal . The energy uncertainty within a single band B j ⊂ I is of the order of the band width | B j | = | ε j ( π/L ) -ε j (0) | which coincides with the variation of the eigenvalue ε j ( k ) as the Bloch boundary condition changes from periodic to anti-periodic (see Figure 2). A rough but convenient estimate of the energy uncertainty within the window I is given by the arithmetic mean

$$\delta E \sim \frac { \sum _ { B _ { j } \subset I } | B _ { j } | } { \sum _ { B _ { j } \subset I } 1 } \sim \frac { | \text{sp} ( h _ { \text{crystal} } ) \cap I | } { \sum _ { B _ { j } \subset I } 1 }.$$

On the other hand, the mean level spacing within I is given by

$$\Delta E \sim \frac { | I | } { \sum _ { B _ { j } \subset I } 1 },$$

and the Thouless conductance becomes

$$g _ { \text{Th} } \sim \frac { \delta E } { \Delta E } \sim \frac { | \text{sp} ( h _ { \text{crystal} } ) \cap I | } { | I | }.$$

Figure 2: Energy uncertainty δE j and level spacing ∆ E j for the j -th band in the window I .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0185v1-figure-11fbb480-4078-4ed3-8f0b-dfd99eddf1ef.jpg)

According to the previous arguments, we shall define the Thouless conductance of the pair ( S , κ S ) for the energy window I ⊂ R by where h crystal is the periodization of h S with J L = κ S .

$$g _ { \text{Th} } ( I ) = \frac { 1 } { 2 \pi } \frac { | \text{sp} ( h _ { \text{crystal} } ) \cap I | } { | I | }, \\ \intertext { \text{periodization of $h_{\text{c} }$ with $J_{r} = \kappa_{\text{c} }$}.$$

Note that if κ S = 0 , then sp( h crystal ) = sp( h S ) and g Th ( I ) = 0 for all intervals I . More generally, the width of the L Bloch bands satisfies | B j | = O ( κ S ) as κ S → 0 , and

$$g _ { \text{Th} } ( I ) = \mathcal { O } ( \kappa _ { \mathcal { S } } )$$

in this limit.

## 2 The crystalline limit

To further elaborate on the connection between Thouless conductance and the Landauer-Büttiker formula, we shall now consider the approximation of h crystal by finite repetitions of the sample S connected to arbitrary reservoirs.

Let h crystal be as in the previous section. Given a positive integer N , let h ( N ) S be the restriction of h crystal to the finite lattice Z NL = [1 , NL ] ∩ Z with Dirichlet boundary condition. Hence h ( N ) S is a Jacobi matrix acting on h ( N ) S = glyph[lscript] 2 ( Z NL ) whose Jacobi parameters satisfy

$$J _ { x + n L } = J _ { x }, \quad \lambda _ { x + n L } = \lambda _ { x }, \quad ( x \in Z _ { L }, n = 0, 1, \dots, N - 1 ), \\ \cdot \quad \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \, \tau \}$$

where { J x } 1 ≤ x&lt;L and { λ x } 1 ≤ ≤ x L are the Jacobi parameters of the original sample Hamiltonian h S and J L = κ S . The pairs ( h ( N ) S , h ( N ) S ) define a sequence of sample systems which are coupled to the reservoirs R l/r as in Section 1.1. The reservoirs' single particle Hilbert spaces and Hamiltonains ( h l/r , h l/r ) , the vectors χ l/r , and the coupling strength κ do not depend on N , and one takes ψ l = δ 1 , ψ r = δ NL . We assume that the one-particle Hamiltonian h ( N ) of the coupled systems satisfies Assumption A for all N 4 and we denote by ω ( N ) + , Φ ( N ) l/r , I ( N ) l/r , J ( N ) the respective NESS and flux observables. We are interested in the large N limit of the charge, energy and entropy steady currents (see Figure 3).

4 Besides the crystaline ones, see Section 1.2.2, a concrete example of reservoirs where this is the case is h l/r = glyph[lscript] 2 ( Z + ) , h l/r = -k ∆ , k &gt; 0 . For other examples and general results regarding this point we refer the reader to [GJW].

Figure 3: The EBBM described by the Hamiltonian h ( N ) for N = 7 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0185v1-figure-db93dec7-4245-43ae-8f59-cb64848a4c84.jpg)

Let h ( ) l crystal and h ( r ) crystal be the restrictions of h crystal to glyph[lscript] 2 (( -∞ , 0] ∩ Z ) and glyph[lscript] 2 ([1 , ∞ ∩ ) Z ) with Dirichlet boundary conditions. Denote by

$$m _ { l } ( E ) & = \langle \delta _ { 0 }, ( h ^ { ( l ) } _ { \text{crystal} } - E - \text{i0} ) ^ { - 1 } \delta _ { 0 } \rangle, \\ m _ { r } ( E ) & = \langle \delta _ { 1 }, ( h ^ { ( r ) } _ { \text{crystal} } - E - \text{i0} ) ^ { - 1 } \delta _ { 1 } \rangle, \\.. \. \. \ \hat { \ }. \ \hat { \ }. \ \hat { \ }. \ \hat { \ }. \ \hat { \ }. \ \hat { \ }. \ \hat { \ }. \ \hat { \ }. \.$$

the respective Weyl m -functions. One easily shows that Im m l/r ( E ) &gt; 0 for Lebesgue a.e. E ∈ sp( h crystal ) . We set T ∞ ( E ) = 0 for E ∈ R \ (sp( h crystal ) ∩ Σ l ∩ Σ ) r and

$$\mathcal { T } _ { \infty } ( E ) = \left [ 1 + \frac { 1 } { 4 } \left ( \frac { | \kappa _ { S } ^ { 2 } m _ { r } ( E ) - \kappa ^ { 2 } F _ { r } ( E ) | ^ { 2 } } { \ln { ( \kappa _ { S } ^ { 2 } m _ { r } ( E ) ) } \ln { ( \kappa ^ { 2 } F _ { r } ( E ) ) } } + \frac { | \kappa _ { S } ^ { 2 } m _ { l } ( E ) - \kappa ^ { 2 } F _ { l } ( E ) | ^ { 2 } } { \ln { ( \kappa _ { S } ^ { 2 } m _ { l } ( E ) ) } \ln { ( \kappa ^ { 2 } F _ { l } ( E ) ) } } \right ) \right ] ^ { - 1 } \ ( 2. 2 )$$

for E ∈ sp( h crystal ) ∩ Σ l ∩ Σ r . Obviously, 0 ≤ T ∞ ( E ) ≤ 1 for Lebesgue a.e. E ∈ R .

Let

$$\mathcal { T } _ { N } ( E ) = 4 \kappa ^ { 4 } | \langle \psi _ { l }, ( h ^ { ( N ) } - E - \mathrm i 0 ) ^ { - 1 } \psi _ { r } \rangle | ^ { 2 } \, \text{Im} \, F _ { l } ( E ) \, \text{Im} \, F _ { r } ( E ), \\ \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \end{$$

be the transmittance of the N -fold repeated pair ( S , κ S ) . Our main technical result is:

Theorem 2.1 For any f ∈ L 1 ( R ) one has

$$\lim _ { N \to \infty } \int \mathcal { T } _ { N } ( E ) f ( E ) \, \mathrm d E = \int \mathcal { T } _ { \infty } ( E ) f ( E ) \, \mathrm d E.$$

The proof of this theorem is given in Section 3. As an immediate consequence, one has

## Theorem 2.2

$$\langle$$

$$\text{Theorem 2.2} \\ \langle \Phi _ { l / r } \rangle _ { \infty } & = \lim _ { N \to \infty } \omega _ { + } ^ { ( N ) } ( \Phi _ { l / r } ^ { ( N ) } ) = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } \mathcal { T } _ { \infty } ( E ) E \Delta _ { l / r } ( E ) \, \text{d} E, \\ \langle \mathcal { I } _ { l / r } \rangle _ { \infty } & = \lim _ { N \to \infty } \omega _ { + } ^ { ( N ) } ( \mathcal { I } _ { l / r } ^ { ( N ) } ) = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } \mathcal { T } _ { \infty } ( E ) \Delta _ { l / r } ( E ) \, \text{d} E, \\ \langle \mathcal { J } \rangle _ { \infty } & = \lim _ { N \to \infty } \omega _ { + } ^ { ( N ) } ( \mathcal { J } ^ { ( N ) } ) = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } \mathcal { T } _ { \infty } ( E ) \varsigma ( E ) \, \text{d} E.$$

Obviously, the conservation laws

$$\langle \Phi _ { l } \rangle _ { \infty } + \langle \Phi _ { r } \rangle _ { \infty } & = 0, \\ \langle \mathcal { I } _ { l } \rangle _ { \infty } + \langle \mathcal { I } _ { r } \rangle _ { \infty } & = 0, \\ \cdot \.$$

hold, as well as the entropy balance relation

$$\langle \mathcal { J } \rangle _ { \infty } = - \beta _ { l } ( \langle \Phi _ { l } \rangle _ { \infty } - \mu _ { l } \langle \mathcal { I } _ { l } \rangle _ { \infty } ) - \beta _ { r } ( \langle \Phi _ { r } \rangle _ { \infty } - \mu _ { r } \langle \mathcal { I } _ { r } \rangle _ { \infty } ).$$

## 2.1 Thouless conductance revisited

Given a finite sample described by h S , h S , κ S and its periodization h crystal , the Thouless currents associated to chemical potentials µ l/r and inverse temperatures β l/r are defined by setting T ∞ ( E ) = 1 in the formulas (2.4), i.e., by assuming that the transport between the reservoirs is reflectionless. The Thouless current formulas are:

$$\langle \Phi _ { l / r } \rangle _ { \text{Th} } = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } E \Delta _ { l / r } ( E ) \text{d} E,$$

$$\langle \mathcal { I } _ { l / r } \rangle _ { \text{Th} } = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } \Delta _ { l / r } ( E ) \text{d} E,$$

$$\langle \Phi _ { l / r } \rangle _ { \text{Th} } & = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } E \Delta _ { l / r } ( E ) \text{d} E, \\ \langle \mathcal { I } _ { l / r } \rangle _ { \text{Th} } & = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } \Delta _ { l / r } ( E ) \text{d} E, \\ \langle \mathcal { J } \rangle _ { \text{Th} } & = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) } \varsigma ( E ) \text{d} E. \\ \text{intervation now}$$

One has again the conservation laws

$$\langle \Phi _ { l } \rangle _ { \text{Th} } + \langle \Phi _ { r } \rangle _ { \text{Th} } & = 0, \\ \langle \mathcal { I } _ { l } \rangle _ { \text{Th} } + \langle \mathcal { I } _ { r } \rangle _ { \text{Th} } & = 0,$$

and

$$\langle \mathcal { J } \rangle _ { \text{Th} } = - \beta _ { l } ( \langle \Phi _ { l } \rangle _ { \text{Th} } - \mu _ { l } \langle \mathcal { I } _ { l } \rangle _ { \text{Th} } ) - \beta _ { r } ( \langle \Phi _ { r } \rangle _ { \text{Th} } - \mu _ { r } \langle \mathcal { I } _ { r } \rangle _ { \text{Th} } ). \\ \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots.$$

Recall the definition of crystalline EBBM given in Section 1.2.2. The following theorem is a direct consequence of the definition of Thouless currents (2.5) and of Theorem 2.2.

## Theorem 2.3 (1)

$$\langle \mathcal { J } \rangle _ { \text{Th} } = \sup \, \langle \mathcal { J } \rangle _ { \infty }, \\ \dots \dots \dots \dots$$

where the supremum is taken over all realizations of the reservoirs. Moreover, this supremum is achieved if the EBBM is crystalline.

$$( 2 ) \ I f \mu _ { l } = \mu _ { r } = \mu, \, \mu \leq \inf \text{sp} ( h _ { \text{crystal} } ), \, a n d \, \beta _ { l } > \beta _ { r }, \, \text{then}$$

$$\langle \Phi _ { r } \rangle _ { T h } = \sup \, \langle \Phi _ { r } \rangle _ { \infty },$$

and the supremum is achieved if the EBBM is crystalline. If β l &lt; β r , then this result holds for 〈 Φ l 〉 . If µ ≥ sup sp( h crystal ) , the same results hold with exchange of l and r .

$$( 3 ) \ I f \, \beta _ { l } = \beta _ { r } \, a n d \, \mu _ { l } < \mu _ { r }, \, t h e n$$

$$\langle \mathcal { I } _ { r } \rangle _ { \text{Th} } = \sup \, \langle \mathcal { I } _ { r } \rangle _ { \infty }, \\.. \,. \,. \,. \,. \,. \,......$$

and the supremum is achieved if the EBBM is crystalline. If µ &gt; µ l r , then this result holds for 〈I 〉 l .

Remark 1. All suprema in the previous theorem are achieved iff the transport between the reservoirs is reflectionless, that is, iff T ∞ ( E ) = 1 for Lebesgue a.e. E ∈ sp( h crystal ) . The crystalline EBBM provide such reservoirs. To elucidate this point further, note that T ∞ ( E ) = 1 for Lebesgue a.e. E ∈ sp( h crystal ) iff σ h ( crystal ) ⊂ Σ l ∩ Σ r and

$$\kappa _ { \mathcal { S } } ^ { 2 } m _ { r / l } ( E ) = \kappa ^ { 2 } F _ { r / l } ( E ),$$

for Lebesgue a.e. E ∈ sp( h crystal ) . Since sp( h crystal ) has positive Lebesgue measure, the theory of boundary values of analytic functions (see [Ja] or any book on harmonic analysis) yields that for any z = E +i glyph[epsilon1] with glyph[epsilon1] &gt; 0 ,

$$\kappa _ { \mathcal { S } } ^ { 2 } \langle \delta _ { 0 }, ( h _ { \text{crystal} } ^ { ( l ) } - z ) ^ { - 1 } \delta _ { 0 } \rangle = \kappa ^ { 2 } \langle \chi _ { l }, ( h _ { l } - z ) ^ { - 1 } \chi _ { l } \rangle,$$

$$\kappa _ { \mathcal { S } } ^ { 2 } \langle \delta _ { 1 }, ( h _ { \text{crystal} } ^ { ( r ) } - z ) ^ { - 1 } \delta _ { 1 } \rangle = \kappa ^ { 2 } \langle \chi _ { r }, ( h _ { r } - z ) ^ { - 1 } \chi _ { r } \rangle.$$

These relations imply

$$\kappa _ { \mathcal { S } } ^ { 2 } \nu _ { \text{crystal} } ^ { ( l / r ) } = \kappa ^ { 2 } \nu _ { l / r },$$

where ν ( l/r ) crystal is the spectral measure of h ( l/r ) crystal associated to δ /δ 0 1 . Since ν ( l/r ) crystal and ν l/r are probability measures, we conclude that κ 2 S = κ 2 and ν ( l/r ) crystal = ν l/r . Thus, the transport between the reservoirs is reflectionless iff κ 2 = κ 2 S and h r/l is unitarily equivalent to h ( r/l ) crystal . In other words, all suprema in the previous theorem are achieved iff the EBBM is unitarily equivalent to a crystalline EBMM up to (for the transport purposes) irrelevant choice of the sign of κ .

Remark 2. Part (2) holds whenever E ∆ l/r ( E ) has a definite sign on sp( h crystal ) and, similarly, Part (3) holds whenever ∆ l/r ( E ) has a definite sign on sp( h crystal ) . In the lack of a definite sign one cannot expect a variational characterization of Thouless currents in terms of the crystalline Landauer-Büttiker currents. Note that

$$\Delta _ { l / r } ( E ) = \frac { \sinh ( ( \zeta _ { r / l } ( E ) - \zeta _ { l / r } ( E ) ) / 2 ) } { 2 \cosh ( \zeta _ { l } ( E ) / 2 ) \cosh ( \zeta _ { r } ( E ) / 2 ) },$$

so that the sign of ∆ l/r ( E ) is the same as the sign of

$$\zeta _ { r / l } ( E ) - \zeta _ { l / r } ( E ) = ( \beta _ { r / l } - \beta _ { l / r } ) E - ( \beta _ { r / l } \mu _ { l / r } - \beta _ { l / r } \mu _ { l / r } ). \\ \dots \quad \cdots \quad \cdots \quad \cdots \quad \cdots.$$

In the non-trivial case β l = β r , ∆ l/r ( E ) changes the sign at precisely one point glyph[negationslash]

$$E _ { c } = \frac { \beta _ { r / l } \mu _ { l / r } - \beta _ { l / r } \mu _ { l / r } } { \beta _ { r / l } - \beta _ { l / r } }. \\ \twoheadrightarrow { \sim } \left ( h \right ) \quad \cdots \, \colon = \, \colon \, \text{sub-lab} \, \text{to-show} \,.$$

If additional information about sp( h crystal ) is available, the above fact can be used to obtain further relations between the crystalline Landauer-Büttiker currents (2.4) and Thouless currents (2.5).

Remark 3. If β l = β r = ∞ , µ &lt; µ l r , and I = [ µ , µ l r ] , then

$$\langle \Phi _ { r } \rangle _ { \text{Th} } = \frac { 1 } { 2 \pi } \int _ { \text{sp} ( h _ { \text{crystal} } ) \cap I } E \text{d} E,$$

$$\langle \mathcal { I } _ { r } \rangle _ { \text{Th} } = \frac { 1 } { 2 \pi } | \text{sp} ( h _ { \text{crystal} } ) \cap I |.$$

Thus, the Thouless formula (1.9) indeed describes the maximal conductance at zero temperature for the given potential interval I :

$$g _ { \text{Th} } ( I ) = \frac { \langle \mathcal { I } _ { r } \rangle _ { \text{Th} } } { \mu _ { r } - \mu _ { l } } = \sup \frac { \langle \mathcal { I } _ { r } \rangle _ { \infty } } { \mu _ { r } - \mu _ { l } }. \\ \dots \dots \dots \dots \dots \dots \dots \dots$$

Remark 4. Since the crystalline Landauer-Büttiker formulas (2.4) are derived from the first principles of quantum statistical mechanics, Theorem 2.3 can be considered a rigorous quantum statistical derivation of the Thouless energy formula. This derivation also identifies the heuristic notion of "optimal feeding" of electrons with reflectionless transport between the reservoirs.

## 3 Proof of Theorem 2.1

## 3.1 Sample transmittance and Green matrix

We first connect the transmittance (1.2) to the sample's Green function. Recall that F l/r are defined in (1.3) and denote by F E ( ) the 2 × 2 diagonal matrix with entries F l ( E ) and F r ( E ) . We also introduce the 2 × 2 Green matrices G ( N ) S ( z ) and G ( N ) ( z ) with entries

$$G _ { S, a b } ^ { ( N ) } ( z ) = \langle \psi _ { a }, ( h _ { \mathcal { S } } ^ { ( N ) } - z ) ^ { - 1 } \psi _ { b } \rangle,$$

$$G _ { a b } ^ { ( N ) } ( z ) = \langle \psi _ { a }, ( h ^ { ( N ) } - z ) ^ { - 1 } \psi _ { b } \rangle,$$

where a, b ∈ { l, r } (recall that ψ l = δ 1 and ψ r = δ NL ). As usual, we write G ( N ) ab ( E ) = G ( N ) ab ( E +i0) .

The full Green matrix G ( N ) and the sample Green matrix G ( N ) S are related by (see Lemma 2.1 in [BJP])

$$G _ { \mathcal { S } } ^ { ( N ) } ( E ) = ( I - \kappa ^ { 2 } G _ { \mathcal { S } } ^ { ( N ) } ( E ) F ( E ) ) G ^ { ( N ) } ( E ),$$

from which we deduce

$$G _ { l r } ^ { ( N ) } ( E ) = \frac { G _ { \mathcal { S }, l r } ^ { ( N ) } ( E ) } { \det ( I - \kappa ^ { 2 } G _ { \mathcal { S } } ^ { ( N ) } ( E ) F ( E ) ) }. \quad \quad \ \ ( 3. 1 ) \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

Combined with (1.2) and (1.6) this allows us the expression of the transmittance of the sample and hence the steady currents in terms of the Green matrix G ( N ) S .

## 3.2 Green and transfer matrix

Our next step is to relate the sample's Green matrix to the transfer matrix of the periodic Jacobi matrix h crystal .

Following [Si] the transfer matrix at energy E is defined by

$$T _ { n } ( E ) = A _ { n } ( E ) \cdots A _ { 1 } ( E ),$$

where

$$A _ { x } ( E ) = J _ { x } ^ { - 1 } \left [ \begin{array} { c c } E - \lambda _ { x } & - 1 \\ J _ { x } ^ { 2 } & 0 \end{array} \right ].$$

Note that det A x ( E ) = 1 for any x and hence det T n ( E ) = 1 as well. A function u satisfies the finite difference equation h crystal u = Eu if and only if for any x one has

$$A _ { x } \left [ \begin{array} { c } u ( x ) \\ J _ { x - 1 } u ( x - 1 ) \end{array} \right ] = \left [ \begin{array} { c } u ( x + 1 ) \\ J _ { x } u ( x ) \end{array} \right ].$$

Lemma 3.1 For any x, y, u, v ∈ C one has

$$G _ { \mathcal { S } } ^ { ( N ) } ( E ) \left [ \begin{array} { c } x \\ y \end{array} \right ] = \left [ \begin{array} { c } u \\ v \end{array} \right ] \Longleftrightarrow T _ { N L } ( E ) \left [ \begin{array} { c } u \\ - x \end{array} \right ] = \left [ \begin{array} { c } - \kappa _ { \mathcal { S } } ^ { - 1 } y \\ \kappa _ { \mathcal { S } } v \end{array} \right ].$$

↦

In other words, the matrix P : ( x, y, u, v ) → ( u, -x, -κ -1 S y, κ S v ) maps the graph of G ( N ) S ( E ) to that of T NL ( E ) .

Remark. This lemma is a slight generalization of Lemma 2.2 in [BJP]. We include the proof for the reader convenience.

Proof. Fix N and E ∈ R \ sp( h ( N ) S ) . For f ∈ glyph[lscript] 2 ( Z NL ) , the function

$$u ( x ) = \langle \delta _ { x }, ( h _ { \mathcal { S } } ^ { ( N ) } - E ) ^ { - 1 } f \rangle$$

satisfies the finite difference equation

$$( h _ { \text{crystal} } - E ) u = f,$$

with boundary conditions u (0) = u NL ( +1) = 0 . Using the transfer matrix

$$T ( x, y ) = A _ { x } A _ { x - 1 } \cdots A _ { y + 1 },$$

the solution of the initial value problem for (3.2) can be written as

$$\left [ \begin{array} { c } u ( x + 1 ) \\ J _ { x } u ( x ) \end{array} \right ] = T ( x, 0 ) \left [ \begin{array} { c } u ( 1 ) \\ J _ { 0 } u ( 0 ) \end{array} \right ] - \sum _ { y = 1 } ^ { x } T ( x, y - 1 ) \left [ \begin{array} { c } 0 \\ f ( y ) \end{array} \right ].$$

Setting x = NL and taking the boundary conditions and J NL = J 0 = κ S into account yields

$$\left [ \begin{array} { c } 0 \\ \kappa _ { \mathcal { S } } u ( N L ) \end{array} \right ] = T _ { N L } ( E ) \left [ \begin{array} { c } u ( 1 ) \\ 0 \end{array} \right ] - \sum _ { y = 1 } ^ { N L } T ( N L, y - 1 ) \left [ \begin{array} { c } 0 \\ f ( y ) \end{array} \right ],$$

which is an equation for the unknown u (1) and u NL ( ) . Setting f = δ 1 and f = δ NL , we obtain the following equations for the entries of the matrix G ( N ) S ( E ) :

$$T _ { N L } ( E ) & \left [ \begin{array} { c } G _ { \mathcal { S }, l l } ^ { ( N ) } ( E ) \\ - 1 \end{array} \right ] = \left [ \begin{array} { c } 0 \\ \kappa _ { \mathcal { S } } G _ { \mathcal { S }, r l } ^ { ( N ) } ( E ) \end{array} \right ], \\ _ { - \quad \Phi \Phi } \Phi _ { - \quad \Phi \Phi } \Phi _ { - \quad \Phi \Phi }$$

$$T _ { N L } ( E ) \left [ \begin{array} { c } G _ { \mathcal { S }, l r } ^ { ( N ) } ( E ) \\ 0 \end{array} \right ] = \left [ \begin{array} { c } - \kappa _ { \mathcal { S } } ^ { - 1 } \\ \kappa _ { \mathcal { S } } G _ { \mathcal { S }, r r } ^ { ( N ) } ( E ) \end{array} \right ].$$

Thus, the two linearly independent vectors

$$\left \{ [ G _ { \mathcal { S }, l l } ^ { ( N ) } ( E ), - 1, 0, \kappa _ { \mathcal { S } } G _ { \mathcal { S }, r l } ^ { ( N ) } ( E ) ] ^ { T }, [ G _ { \mathcal { S }, l r } ^ { ( N ) } ( E ), 0, - \kappa _ { \mathcal { S } } ^ { - 1 }, \kappa _ { \mathcal { S } } G _ { \mathcal { S }, r r } ^ { ( N ) } ( E ) ] ^ { T } \right \},$$

span the graph of T NL ( E ) . One easily checks that they are the images by the matrix P of the two vectors [1 , 0 , G ( N ) S ,ll ( E , G ) ( N ) S ,rl ( E )] T and [0 , 1 , G ( N ) S ,lr ( E , G ) ( N ) S ,rr ( E )] T which span the graph of G ( N ) S ( E ) . ✷

## 3.3 Transfer matrix eigenvalues and eigenfunctions

Since the Jacobi matrix h crystal is periodic, one has T NL ( E ) = T L ( E ) N for any N . The eigenvalues and eigenvectors of the one-period transfer matrix T L ( E ) will thus play an important role. Since det T L ( E ) = 1 , we may write its eigenvalues as α E ( ) and α E ( ) -1 . It is a standard result (see, e.g., [RS4, Si]) that these eigenvalues are either complex conjugated (and hence of modulus 1) or real. The first case occurs iff | tr T L ( E ) | ≤ 2 which is in turn equivalent to E ∈ sp( h crystal ) . We denote by

$$\Psi _ { \pm } ( E ) = \left [ \begin{array} { c } \phi _ { \pm } ( E ) \\ \kappa _ { \mathcal { S } } \psi _ { \pm } ( E ) \end{array} \right ]$$

an eigenvector of T L ( E ) associated to its eigenvalue α E ( ) ± 1 with the following conventions:

(N1) When E ∈ sp( h crystal ) we chose | α E ( ) | -1 &lt; 1 &lt; α E | ( ) | and real eigenvectors.

glyph[negationslash]

(N2) When E ∈ sp( h crystal ) we chose Ψ ( -E ) = Ψ ( + E ) . Since the two eigenvectors are linearly independent, they can be normalized by φ ± ( E ) = 1 , which further implies Im ψ ± ( E ) = 0 . We then select α E ( ) such that Im ψ + ( E &gt; ) 0 .

glyph[negationslash]

Using Lemma 3.1 and T NL ( E ) = T L ( E ) N , we get

$$G _ { S } ^ { ( N ) } ( E ) & = \frac { - \kappa _ { S } ^ { - 1 } } { D _ { N } ( E ) } \left [ \begin{array} { c c } \phi _ { + } ( E ) \phi _ { - } ( E ) ( \alpha ( E ) ^ { N } - \alpha ( E ) ^ { - N } ) & \phi _ { + } ( E ) \psi _ { - } ( E ) - \phi _ { - } ( E ) \psi _ { + } ( E ) \\ \phi _ { + } ( E ) \psi _ { - } ( E ) - \phi _ { - } ( E ) \psi _ { + } ( E ) & \psi _ { + } ( E ) \psi _ { - } ( E ) ( \alpha ( E ) ^ { N } - \alpha ( E ) ^ { - N } ) \end{array} \right ] \\ \dots \dots$$

where

$$D _ { N } ( E ) = \alpha ( E ) ^ { N } \phi _ { + } ( E ) \psi _ { - } ( E ) - \alpha ( E ) ^ { - N } \phi _ { - } ( E ) \psi _ { + } ( E ).$$

An elementary calculation yields

$$\det ( I - \kappa ^ { 2 } G _ { S } ^ { ( N ) } ( E ) F ( E ) ) = \frac { \alpha ( E ) ^ { N } \widetilde { \phi } _ { + } ( E ) \widetilde { \psi } _ { - } ( E ) - \alpha ( E ) ^ { - N } \widetilde { \phi } _ { - } ( E ) \widetilde { \psi } _ { + } ( E ) } { \alpha ( E ) ^ { N } \phi _ { + } ( E ) \psi _ { - } ( E ) - \alpha ( E ) ^ { - N } \phi _ { - } ( E ) \psi _ { + } ( E ) }, \\ \intertext { a n o }$$

where

$$\widetilde { \psi } _ { \pm } ( E ) = \psi _ { \pm } ( E ) + \eta ^ { 2 } \kappa _ { \mathcal { S } } \phi _ { \pm } ( E ) F _ { l } ( E ), \quad \widetilde { \phi } _ { \pm } ( E ) = \phi _ { \pm } ( E ) + \eta ^ { 2 } \kappa _ { \mathcal { S } } \psi _ { \pm } ( E ) F _ { r } ( E ), \quad ( 3. 3 )$$

and we have set

$$\eta = \frac { \kappa } { \kappa _ { \mathcal { S } } }.$$

Inserting the last relations into (3.1) leads to the following expression for the off-diagonal element of the full Green matrix

$$G _ { l r } ^ { ( N ) } ( E ) = - \kappa _ { \mathcal { S } } ^ { - 1 } \frac { \phi _ { + } ( E ) \psi _ { - } ( E ) - \phi _ { - } ( E ) \psi _ { + } ( E ) } { \alpha ( E ) ^ { N } \widetilde { \phi } _ { + } ( E ) \widetilde { \psi } _ { - } ( E ) - \alpha ( E ) ^ { - N } \widetilde { \phi } _ { - } ( E ) \widetilde { \psi } _ { + } ( E ) }.$$

## 3.4 The large N limit

We proceed to evaluate the large N limit of T N ( E ) in distributional sense.

We write the left hand side of (2.3) as T N, 1 ( f ) + T N, 2 ( f ) with

$$\mathcal { T } _ { N, 1 } ( f ) = \int _ { \text{sp} ( h _ { \text{crystal} } ) \cap \Sigma _ { l } \cap \Sigma _ { r } } \mathcal { T } _ { N } ( E ) f ( E ) \, \mathrm d E,$$

$$\mathcal { T } _ { N, 2 } ( f ) = \int _ { \mathbb { R } \langle \text{sp} ( h _ { \text{crystal} } ) \cap \Sigma _ { l } \cap \Sigma _ { r } ) } \mathcal { T } _ { N } ( E ) f ( E ) \, \mathrm d E.$$

To deal with T N, 2 ( f ) we prove

Lemma 3.2 For almost all E ∈ R \ (sp( h crystal ) ∩ Σ l ∩ Σ ) r one has

$$\lim _ { N \to \infty } \mathcal { T } _ { N } ( E ) = 0.$$

Proof. Combining (1.2) and (3.4), we get

$$\mathcal { T } _ { N } ( E ) = 4 \kappa ^ { 2 } \eta ^ { 2 } \left | \frac { \phi _ { + } ( E ) \psi _ { - } ( E ) - \phi _ { - } ( E ) \psi _ { + } ( E ) } { \alpha ( E ) ^ { N } \widetilde { \phi } _ { + } ( E ) \widetilde { \psi } _ { - } ( E ) - \alpha ( E ) ^ { - N } \widetilde { \phi } _ { - } ( E ) \widetilde { \psi } _ { + } ( E ) } \right | ^ { 2 } \, \text{Im} \, F _ { l } ( E ) \, \text{Im} \, F _ { r } ( E ). \\ \intertext { a n d } \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \, \dots \end{array}$$

It follows from (1.4) that T N ( E ) = 0 for almost all E ∈ R \ (Σ l ∩ Σ ) r and all N . Thus, it suffices to consider E ∈ (Σ l ∩ Σ ) r \ sp( h crystal ) .

glyph[negationslash]

For such E , the eigenvectors Ψ ( ± E ) are real and | α E ( ) | &gt; 1 by Condition (N1). Moreover, Im F l/r ( E ) &gt; 0 by (1.5). We claim that this implies ˜ φ + ( E ψ ) ˜ -( E ) = 0 from which the result clearly follows. To prove this claim, we argue by contradiction. Assume that ˜ φ + ( E ) = 0 (the case ˜ ψ -( E ) = 0 is similar). It follows from (3.3) that

$$\text{Im} \, \widetilde { \phi } _ { + } ( E ) = \eta ^ { 2 } \kappa _ { \mathcal { S } } \psi _ { + } ( E ) \text{Im} \, F _ { r } ( E ) = 0,$$

and hence ψ + ( E ) = 0 . Using (3.3) again we get φ + ( E ) = ˜ φ + ( E ) = 0 . We conclude that Ψ ( + E ) = 0 , a contradiction. ✷

Since 0 ≤ T N ( E ) ≤ 1 and f ∈ L 1 ( R ) , Lemma 3.2 and the dominated convergence theorem yield that T N, 2 ( f ) → 0 as N →∞ . It thus remains to analyze T N, 1 ( f ) . To this end, we relate the eigenvectors of T L to the Weyl m -functions (2.1).

Lemma 3.3 For any E ∈ sp( h crystal ) , we have

$$\psi _ { + } ( E ) = - \frac { 1 } { \kappa _ { \mathcal { S } } m _ { r } ( E ) }, \quad \psi _ { - } ( E ) = - \kappa _ { \mathcal { S } } m _ { l } ( E ).$$

Proof. Let us write the one-period transfer matrix as

$$T _ { L } ( E ) = \left [ \begin{array} { c c } a ( E ) & b ( E ) \\ c ( E ) & d ( E ) \end{array} \right ].$$

Since det T L ( E ) = 1 , we can write the discriminant of the quadratic equation

$$c ( E ) z ^ { 2 } + ( a ( E ) - d ( E ) ) z - b ( E ) = 0,$$

as (tr T L ( E )) 2 -4 , which is negative for E ∈ sp( h crystal ) . Thus, (3.6) has two complex conjugate solutions. The solution with positive imaginary part is m E r ( ) and the other one is 1 /κ m 2 S l ( E ) . We refer the reader to Section 5.2 in [Si] for a proof of these facts.

Using the representation (3.5) and taking Condition (N2) into account, one easily shows that the eigenvalues and eigenvectors of T L ( E ) are determined by α E ( ) = e i θ E ( ) and

$$\kappa _ { \mathcal { S } } \psi _ { \pm } ( E ) = \frac { 1 } { b ( E ) } \left ( \mathrm e ^ { \pm i \theta ( E ) } - a ( E ) \right ) = \frac { 1 } { b ( E ) } \left ( \frac { d ( E ) - a ( E ) } { 2 } \pm \mathrm i \sin \theta ( E ) \right ),$$

where θ E ( ) ∈ -] π, π [ is defined by

$$\cos \theta ( E ) = \frac { 1 } { 2 } \text{tr} T _ { L } ( E ), \quad \text{sign} ( \theta ( E ) ) = \text{sign} ( b ( E ) ).$$

Note that the product of the two solutions of (3.6) is

$$- \, \frac { b ( E ) } { c ( E ) } = \frac { m _ { r } ( E ) } { \kappa _ { \mathcal { S } } ^ { 2 } m _ { l } ( E ) } = | m _ { r } ( E ) | ^ { 2 } > 0,$$

so that b E ( ) and c E ( ) have opposite signs. It follows that the solutions of (3.6) are

$$m _ { r } ( E ) = \frac { 1 } { c ( E ) } \left ( \frac { d ( E ) - a ( E ) } { 2 } - \text{i} \sin \theta ( E ) \right ),$$

$$\frac { 1 } { \kappa _ { \mathcal { S } } ^ { 2 } m _ { l } ( E ) } = \frac { 1 } { c ( E ) } \left ( \frac { d ( E ) - a ( E ) } { 2 } + \text{i} \sin \theta ( E ) \right ).$$

Comparing these relations with (3.7) and using (3.9) yield the result.

✷

To formulate our next result, let the real functions r E ( ) ≥ 0 and ϑ E ( ) be defined by the following polar decomposition

$$\frac { m _ { l } ( E ) - \eta ^ { 2 } F _ { l } ( E ) } { \overline { m } _ { l } ( E ) - \eta ^ { 2 } F _ { l } ( E ) } \, \frac { m _ { r } ( E ) - \eta ^ { 2 } F _ { r } ( E ) } { \overline { m } _ { r } ( E ) - \eta ^ { 2 } F _ { r } ( E ) } \, \frac { \overline { m } _ { r } ( E ) } { m _ { r } ( E ) } = r ( E ) \mathrm e ^ { i \vartheta ( E ) }.$$

Lemma 3.4 Let I ⊂ sp( h crystal ) ∩ Σ l ∩ Σ r . If there exists δ &lt; 1 such that r E ( ) ≤ δ for almost all E ∈ I , then

$$\lim _ { N \to \infty } \int _ { I } \mathcal { T } _ { N } ( E ) f ( E ) \, \mathrm d E = \int _ { I } \mathcal { T } _ { \infty } ( E ) f ( E ) \, \mathrm d E,$$

for any f ∈ L 1 ( R ) .

Proof. Combining Lemma 3.3 with (1.2) and (3.4) we can write the transmittance of the N -fold repeated sample as

$$\mathcal { T } _ { N } ( E ) & = \frac { 1 6 \eta ^ { 4 } \text{Im} \, m _ { r } ( E ) \text{Im} \, m _ { l } ( E ) } { | \overline { m } _ { r } ( E ) - \eta ^ { 2 } F _ { r } ( E ) | ^ { 2 } | \overline { m } _ { l } ( E ) - \eta ^ { 2 } F _ { l } ( E ) | ^ { 2 } } \frac { \text{Im} \, F _ { l } ( E ) \text{Im} \, F _ { r } ( E ) } { | 1 - r ( E ) \text{e} ^ { i ( 2 N \theta ( E ) + \vartheta ( E ) ) | ^ { 2 } } }, \quad ( 3. 1 1 ) \\ \cdot \quad \lambda \equiv \cdot \cdot \cdot \cdot \cdot \cdot \cdot \lambda \equiv \cdot \cdot \cdot \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot \cdot \cdot \cdot \lambda \cdot$$

where θ E ( ) is defined by (3.8). Expanding the right hand side of (3.11) in powers of r E ( ) , one obtains

$$\mathcal { T } _ { N } ( E ) = \mathcal { T } _ { \infty } ( E ) \sum _ { k \in \mathbb { Z } } r ( E ) ^ { | k | } \text{e} ^ { \text{i} k ( 2 N \theta ( E ) + \vartheta ( E ) ) },$$

where T ∞ ( E ) is given by (2.2). Since r E ( ) ≤ δ &lt; 1 on I , this expansion is uniformly convergent for E ∈ I , and we have

$$\lim _ { N \to \infty } \int _ { I } \mathcal { T } _ { N } ( E ) f ( E ) \, \mathrm d E = \lim _ { N \to \infty } \sum _ { k \in \mathbb { Z } } \int _ { I } \mathcal { T } _ { \infty } ( E ) r ( E ) ^ { | k | } \mathrm e ^ { i k ( 2 N \theta ( E ) + \vartheta ( E ) ) } f ( E ) \, \mathrm d E.$$

Since the function θ E ( ) is strictly monotone in each band of h crystal (see, e.g., Sections 5.3-5.4 in [Si]), the Riemann-Lebesgue lemma yields the result. ✷

For k &gt; 0 , set

$$I _ { k } = \{ E \in \text{sp} ( h _ { \text{crystal} } ) \cap \Sigma _ { l } \cap \Sigma _ { r } \, | \, \text{Im} \, F _ { l / r } ( E ) \geq 1 / k \text{ and } \, | \text{tr} \, T _ { L } ( E ) | < 2 - 1 / k \},$$

and I ′ k = (sp( h crystal ) ∩ Σ l ∩ Σ ) r \ I k . Obviously,

$$\lim _ { k \to \infty } | I _ { k } ^ { \prime } | = 0.$$

It follows from the proof of Lemma 3.3 that there is glyph[epsilon1] k &gt; 0 such that Im m l/r ( E ) ≥ glyph[epsilon1] k for almost every E ∈ I k . One easily concludes that there exists δ k &lt; 1 such that

$$& \left | \frac { m _ { l / r } ( E ) - \eta ^ { 2 } F _ { l / r } ( E ) } { \overline { m } _ { l / r } ( E ) - \eta ^ { 2 } F _ { l / r } ( E ) } \right | \leq \delta _ { k }, \\ \ast \quad \left | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right | \right |$$

holds for almost every E ∈ I k . Thus, for such E ,

$$r ( E ) = \left | \frac { m _ { l } ( E ) - \eta ^ { 2 } F _ { l } ( E ) } { \overline { m } _ { l } ( E ) - \eta ^ { 2 } F _ { l } ( E ) } \right | \left | \frac { m _ { r } ( E ) - \eta ^ { 2 } F _ { r } ( E ) } { \overline { m } _ { r } ( E ) - \eta ^ { 2 } F _ { r } ( E ) } \right | \leq \delta _ { k } ^ { 2 } < 1.$$

## Writing

$$\lambda \left | \mathcal { T } _ { N, 1 } ( f ) - \int \mathcal { T } _ { \infty } ( E ) f ( E ) \, \mathrm d E \right | & \leq \left | \int _ { I _ { k } } ( \mathcal { T } _ { N } ( E ) - \mathcal { T } _ { \infty } ( E ) ) f ( E ) \, \mathrm d E \right | \\ & \quad + \int _ { I _ { k } ^ { \prime } } | ( \mathcal { T } _ { N } ( E ) - \mathcal { T } _ { \infty } ( E ) ) f ( E ) | \, \mathrm d E, \\ \mu \dots \dots \dots \nu \dots \nu \dots \nu \dots \nu$$

and applying Lemma 3.4 to I k , we get

$$\lim _ { N \to \infty } \sup \left | \mathcal { T } _ { N, 1 } ( f ) - \int \mathcal { T } _ { \infty } ( E ) f ( E ) \, \mathrm d E \right | \leq 2 \int _ { I ^ { \prime } _ { k } } | f ( E ) | \, \mathrm d E.$$

This estimate and (3.12) yield

$$\lim _ { N \to \infty } \mathcal { T } _ { N, 1 } ( f ) = \int \mathcal { T } _ { \infty } ( E ) f ( E ) \, \mathrm d E$$

and Theorem 2.1 follows.

## References

[AALR] Abrahams, E., Anderson, P.W., Licciardello, D.C., and Ramakrishnan, T.V.: Scaling theory of localization: Absence of quantum diffusion in two dimensions. Phys. Rev. Lett. 42 , 673-676 (1979).

- [AJPP] Aschbacher, W., Jakši´ c, V., Pautrat, Y., and Pillet, C.-A.: Transport properties of quasi-free fermions. J. Math. Phys. 48 , 032101 (2007).

[BSP] Ben Sâad, R., and Pillet, C-A.: A geometric approach to the Landauer-Büttiker formula. J. Math. Phys. 55 , 075202 (2014).

- [Bu] Busch, P.: The time-energy uncertainty relation. In Time in Quantum Mechanics, second edition. J.G. Muga, R. Sala Mayato and Í.L. Egusquiza editors. Springer, Berlin, 2008.
- [BJP] Bruneau, L., Jakši´ c, V ., and Pillet, C.A.: Landauer-Büttiker formula and Schrödinger conjecture. Commun. Math. Phys., 319 , 501-513 (2013).

[CJM] Cornean, H.D., Jensen, A., and Moldoveanu, V.: A rigorous proof of the LandauerBüttiker formula. J. Math. Phys. 46 , 042106 (2005).

- [ET] Edwards, J.T., and Thouless, D.J.: Numerical studies of localization in disordered systems. J. Phys. C: Solid State Phys. 5 , 807-820 (1972).

| [FP]   | Fröhlich, J., and Pfeifer, P.: Generalized time-energy uncertainty relations and bounds on lifetimes of resonances. Rev. Mod. Phys. 67 , 759-779 (1995).                                  |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [GJW]  | Grech, P., Jakši´, c V., and Westrich, M.: The spectral structure of the electronic black box Hamiltonian. Lett. Math. Phys. 103 , 1135-1147 (2013).                                      |
| [GNP]  | Gesztesy, F., Nowell, R., and Potz, W.: One-dimensional scattering theory for quan- tum systems with nontrivial spatial asymptotics. Diff. Integral Eqs. 10 , 521-546 (1997).             |
| [GS]   | Gesztesy, F., and Simon, B.: Inverse spectral analysis with partial information on the potential I. The case of an a.c. component in the spectrum. Helv. Phys. Acta 70 , 66-71 (1997).    |
| [Ja]   | Jakši´, c V.: Topics in spectral theory. Open Quantum Systems I. The Hamiltonian Approach. Lecture Notes in Mathematics 1880 , 235-312, Springer, 2006.                                   |
| [JLPa] | Jakši´, c V., Landon, B., and Panati, A.: A note on reflectionless Jacobi matrices. To appear in Commun. Math. Phys.                                                                      |
| [JLPi] | Jakši´, c V., Landon, B., and Pillet, C.-A.: Entropic fluctuations in XY chains and reflectionless Jacobi matrices. Ann. Henri Poincaré, 14 , 1775-1800 (2013).                           |
| [La]   | Last, Y.: Conductance and spectral properties. Ph.D. Thesis, Technion (1994).                                                                                                             |
| [N]    | Nenciu, G.: Independent electrons model for open quantum systems: Landauer- Büttiker formula and strict positivity of the entropy production. J. Math. Phys. 48 , 033302 (2007).          |
| [MT]   | Mandelstam, L., and Tamm, I.: The uncertainty relation between energy and time in non-relativistic quantum mechanics. J. Phys. (Moscow) 9 , 249-254 (1945).                               |
| [RS4]  | Reed, M., and Simon, B.: Methods of Modern Mathematical Physics IV. Analysis of Operators. Academic Press, New York, 1978.                                                                |
| [Si]   | Simon, B.: Szegö's Theorem and Its Descendants. Spectral theory for L 2 Perturba- tions of Orthogonal Polynomials. M.B. Porter Lectures. Princeton University Press, Princeton, NJ, 2011. |
| [Th]   | Thouless, D.J.: Maximum metallic resistance in thin wires. Phys. Rev. Lett. 39 , 1167-1169 (1977).                                                                                        |
Laurent Bruneau Vojkan Jaksic Yoram Last Claude-Alain Pillet
Landauer-Büttiker and Thouless conductance
In the independent electron approximation, the average
(energy/charge/entropy) current flowing through a finite sample S connected to
two electronic reservoirs can be computed by scattering theoretic arguments
which lead to the famous Landauer-B\"uttiker formula. Another well known
formula has been proposed by Thouless on the basis of a scaling argument. The
Thouless formula relates the conductance of the sample to the width of the
spectral bands of the infinite crystal obtained by periodic juxtaposition of S.
In this spirit, we define Landauer-B\"uttiker crystalline currents by extending
the Landauer-B\"uttiker formula to a setup where the sample S is replaced by a
periodic structure whose unit cell is S. We argue that these crystalline
currents are closely related to the Thouless currents. For example, the
crystalline heat current is bounded above by the Thouless heat current, and
this bound saturates iff the coupling between the reservoirs and the sample is
reflectionless. Our analysis leads to a rigorous derivation of the Thouless
formula from the first principles of quantum statistical mechanics.