Quantum stochastic transport along chains

The spreading of a particle along a chain, and its relaxation, are central themes in statistical and quantum mechanics. One wonders what are the consequences of the interplay between coherent and stochastic transitions. This fundamental puzzle has not been addressed in the literature, though closely related themes were in the focus of the Physics literature throughout the last century, highlighting quantum versions of Brownian motion. Most recently this question has surfaced again in the context of photo-synthesis. Here we consider both an infinite tight-binding chain and a finite ring within the framework of an Ohmic master equation. With added disorder it becomes the quantum version of the Sinai-Derrida-Hatano-Nelson model, which features sliding and delocalization transitions. We highlight non-monotonic dependence of the current on the bias, and a counter-intuitive enhancement of the effective disorder due to coherent hopping.

A prototype problem in Physics is the dynamics of a particle along chain that consists of sites. If the dynamics is coherent one expects to observe ballistic motion and Bloch oscillations 1 , while for stochastic dynamics one expects to see diffusion and drift. In the presence of disorder, additional fascinating effects emerge: an Anderson localization transition in the coherent problem, and a Sinai-Derrida sliding transition in the stochastic problem. In practical applications the particle can be an exiton [2][3][4] . Past literature regarding quantum spreading in chains [5][6][7][8][9][10][11][12][13][14] , including publications that address the photo-synthesis theme [15][16][17][18][19][20][21][22][23][24] , were focused mainly on the question how noise and dissipation affect coherent transport. In a sense, our interest is in the reversed question.
In the present work we assume independent mechanisms for stochastic asymmetric (dissipative) transitions, and for coherent hamiltonian (conservative) transitions. Such setup is not common: the standard models do not allow to tune on and off the two mechanisms independently. The question arises how the two mechanisms affect each other. Can we simply "sum up" known results for stochastic transport with known results for coherent motion in noisy environment? We shall see that the answer is not trivial. The main surprises come out once we take into account the presence of disorder (see below). An optional way to phrase the question: what is the quantum version of the prototype stochastic problem that is known in the literature as random walk in random environment. As we know from the above cited works, due to disorder, the stochastic dissipative dynamics is not merely a simple minded Brownian motion. We would like to know whether coherence has any implication on the predicted disorder-related crossovers.
We consider a chain whose sites are labeled by x. The particle, or the exiton, can move from site to site (near neighbor transitions only). The transitions are determined by two major parameters: the hopping frequency (c) that controls the coherent hopping; and the fluctuations intensity (ν) that controls the environmentally-induced stochastic transitions. At finite temperature T there is also a dissipation coefficient η ν = T /(2 ) that is responsible for the asymmetry of the stochastic transitions. On top we might have bias (), on-site noisy fluctuations (γ), and different types of disorder. The model is illustrated in Fig. 1.
The dynamics is governed by a master equation for the probability matrix where the dissipators  (B) ν ∝ and γ ∝ (S)  are due to the interaction with the environment. They are responsible for the stochastic aspect of the dynamics. The Hamiltonian H c ( ) contains an on-site potential U x ( ), and a sum over hopping terms ± c x x ( /2 ) 1 . Accordingly it takes the form www.nature.com/scientificreports www.nature.com/scientificreports/ where p is the momentum operator. The unit of length is the site spacing (x is an integer), and the field is In the absence of stochastic terms, coherent transport in ordered chain leads to ballistic motion (without bias) and exhibits Bloch-oscillations (with bias). In disordered chain the spreading is suppressed due to Anderson-localization. The effect of noise and dissipation on coherent transport due to (S)  has been extensively studied. In the Caldeira-Leggett model 25,26 the interaction is with homogeneous fluctuating environment, leading to Brownian motion with Gaussian spreading. If the interaction is with non-homogeneous fluctuating environment (short spatial correlation scale) the spreading is the sum of a decaying coherent Gaussian and a scattered Stochastic Gaussian 27 . The tight binding version of this model has been studied in 28 . It has been found that the decoherence and the stochastic-like evolution are dictated by different bands of the Lindblad -spectrum that correspond, respectively, to the dephasing and to the relaxation rates in NMR studies of two-level dynamics.
In the other extreme of purely stochastic dynamics, ignoring quantum effects, the disordered model, aka random walk in random environment, has been extensively studied by Sinai, Derrida, and followers [29][30][31][32][33][34][35] . Without bias the spreading becomes sub-diffusive, while above some critical bias the drift-velocity becomes finite, aka sliding transition. Strongly related is the transition from over-damped to under-damped relaxation that has been studied for a finite-size ring geometry 36,37 . The latter involves delocalization transition that has been highlighted for non-hermitian Hamiltonians in the works of Hatano, Nelson and followers 38-45 . One should realize that the two extreme limits of coherent and stochastic spreading have to be bridged within the framework of a model that includes an  (B) term, not just an  (S) term. Furthermore, a proper modeling requires the distinction between two types of Master equations. In one extreme we have the Pauli version. Traditionally this version is justified by the secular approximation that assumes weak system-bath interaction. In the other extreme we have the Ohmic version that assumes short correlation time. The so called "singular coupling limit" can be regarded as an optional way to formalize the short correlation time assumption 46 . Clearly in the mesoscopic context it is more appropriate to adopt the Ohmic version, and regard the Pauli version of the dissipator as a formal approximation. outline. The model is presented in terms of an Ohmic master equation. The units of time are chosen such that the basic model parameters are  ν η ≡ c ( , , 1, ) and the strength of the disorder  σ . The interest is in the diffusion coefficient D, the -induced drift velocity v, the implied non equilibrium steady state (NESS) current ≡ I Lv (1/ ) for a ring of length L, and the associated Lindblad -spectrum. The latter is determined via ρ λρ = − , which provides both the relaxation-modes and the decoherence-modes. In particular we observe that the NESS current depends non-monotonically on the bias (Fig. 2a), and that surprisingly it can be enhanced by disorder (Fig. 2b). In a disordered ring, counter-intuitively, relaxation modes become over-damped if coherent transitions are switched on (Fig. 3).

Results
the model. The isolated chain is defined by the H (c) Hamiltonian Eq. (2), that describes a particle or an exiton that can hop along a one-dimensional chain whose sites are labeled by x. The field x  might be non-uniform. For the average value of the field we maintain the notation , while the random component is distributed uniformly (box distribution) within [ , ] σ σ −   . We regard each pair of neighboring sites as a two-level system (Supplementary S1). Accordingly we distinguish between two types of terms in the master equation: those that originate from temporal fluctuations of the potential (dephasing due to noisy detuning), and those that are responsible to stochastic transition between the sites (incoherent hopping). The latter are implied by the replacement is a bath operator that is characterized by fluctuation (1) Illustration of the model system. Each site of the chain is represented by a line segment positioned according to its x coordinate and potential U x ( ). Blue arrow labeled by c represents the possibility for a coherent hopping between two sites. Red arrows represent bath induced stochastic transitions between two sites. The local bath that is responsible for the latter fluctuates with intensity ν, and the induced transitions are asymmetric if η ν = T /(2 ) is non-zero (finite temperatures). Note that their ratio is exp in leading order. The green wiggle lines represent a local bath that induces fluctuations of intensity γ of the on-site potential. Without the baths it is the Anderson model for coherent transport and localization in disordered chain. In the other extreme, if only the stochastic transitions are present, it is the Sinai-Derrida model for motion in random environment. The latter exhibits a sliding transition as the bias is increased, and an associated Hatano-Nelson delocalization transition once relaxation in a closed ring is considered.
www.nature.com/scientificreports www.nature.com/scientificreports/ intensity ν, and temperature T . Hence the system bath coupling term is −W f t ( ) The baths of different bonds are uncorrelated, accordingly the bond-related dissipator takes the form . In the = c 10 case also for 8 =  . Thin-lines are a guide to the eye.
The friction terms represent the response of the bath to the rate of change of the W x . Note that for getting the conventional Fokker-Planck equation the system-bath coupling term would be x f t ( ) − , and V would become the velocity operator. Here we assume interaction with local baths that in general might have different temperatures. See Methods for some extra technical details regarding the master equation, the nature of the disorder, and the handling of the periodic boundary conditions for the ring configuration.
pauli-type dynamics. For pedagogical purpose let us consider first a uniform non-disordered ring without coherent hopping. Furthermore, let us adopt the simplified Pauli-like version of the dissipator (see Methods). Consequently the dynamics of the on-site probabilities p x xx , ρ ≡ decouples from that of the off-diagonal terms. Namely, one obtains for the probabilities a simple rate equation, where the transition rates between sites are  as expected from detailed balance considerations. It follows that the drift velocity and the diffusion coefficient are: Consequently one finds two distinct sets of modes: the stochastic-like relaxation modes that are implied by the rate equation for the probabilities, and off-diagonal decoherence modes. The latter share the same decay rate w w , where γ stands for optional extra off-diagonal decoherence due to on-site fluctuations. An evolving wavepacket (Supplementary S3) will decompose into coherent decaying component that is suppressed by factor γ − e t 0 , and an emerging stochastic component that drifts with velocity v and diffuse with coefficient D.
full ohmic treatment. The state of the particle in the standard representation is given by The master equation, Eq. (1) with Eq. (4), couples the dynamics of the on-site probabilities to that of the off-diagonal elements ρ ≠ r ( 0 ) x . The generator  can be written as a sum of several terms (Supplementary S5): Each term is a super-matrix that operates on the super vector ρ r ( ) x . The first two terms c ( )  and ( ) L E arise from the Hamiltonian Eq. (2). The  ( ) ν term arise from the first term of Eq. (4), which represent noise-induced transitions. The remaining two friction-terms (proportional to η) arise from Eq. (5), and correspond to the two terms in the Hamiltonian. A schematic representation of ρ r ( ) x and the couplings is given in Fig. 4. The coherent hopping that is generated , while ( ) L E contribute "on-site" potential. The noise operator ( ) ν ν  include the Pauli-terms that were discussed previously, and an additional term that couples the = ± r 1 elements. Together with the friction operator the spectrum for a non-disordered ring. For a non-disordered ring the super-matrix  is invariant under x-translations, and therefore we can switch to a Fourier basis where the representation is ρ r q ( ; ). Due to Bloch theorem, the matrix decompose into q-blocks in this basis. Thus in order to find the eigenvalues q s , λ and the corresponding eigenmodes we merely have to handle a one dimensional tight binding r lattice. See Methods. A representative spectrum is provided in Fig. 3a. Consider first the = q 0 eigenstates. For q 0 = the c-dependent couplings are zero. For infinite temperature ( 0 η = ) the only non-zero coupling is between r 1 = ± due to a non-Pauli term in Eq. (32). Consequently the q 0 = block contains the NESS r 0 = (which is merely the identity matrix in the standard basis), along with a pair of non-trivial decoherence modes ± , and a set of uncoupled decoherence modes r 2, 3, = ± ± … . The corresponding q s , λ eigenvalues (for 0 η = ) are: The ± modes become over-damped for small bias, while the s 1 > decoherence modes are always under-damped. Considering the q dependence of the eigenvalues λ q s , we get several bands, as illustrated in Fig. 3a. Our interest below is in the relaxation modes that are associated with q,0 λ , and determine the long time spreading. the neSS. At finite temperature ( 0 η > ) there are extra couplings that lead to a modified NESS. In leading order the NESS eigenstate is 0 the current. For non-zero field ( ≠ 0) the NESS momentum distribution is shifted. The expression for the current operator is complicated (Supplementary S2), but the net NESS current comes out a simple sum of stochastic and coherent terms: We shall further illuminate the physical significance of the second term below. In contrast with the stochastic case, the drift current might be non-monotonic in , see Fig. 2a. Furthermore, there is a convex range where the second derivative of I( )  is positive. The convexity of the current in some  range, implies a counter intuitive effect: current may become larger due to disorder. The argument goes as follows: Assume that the sample is divided into two regions, such that x  is constant in each region, but slightly smaller (larger) than  in the first (second) region. Due to the convex property it is implied that the current will be larger. Extending this argument for a general non-homogeneous (i.e. are represented by filled circles, and the off-diagonal terms by empty circles. The Lindblad generator  induces "transitions" between the elements. The blue grid lines indicate c-induced couplings. The other couplings within ≤ r 2, indicated by red and purple, are due to a local bath. For presentation purpose (to avoid a crowded set of lines) the red couplings that originate from ( ) ν ν  and EL E  η ( ) assume that the local bath is positioned at bond x′, while purple couplings that originate from η   c c ( ) assume that the bath is positioned at bond x″.
www.nature.com/scientificreports www.nature.com/scientificreports/ disordered) field, with the same average bias , we expect to observe a larger NESS current. This is indeed confirmed in Fig. 2a It is therefore enough to determine λ q,0 via second order perturbation theory with respect to the = q 0 eigenstates. To leading order in η, a lengthy calculation leads to a result that is consistent with the Einstein relation, namely  = Thus, to leading order, D is given by Eq. (8), multiplied by the expression in the square brackets in Eq. (17). We see that with coherent transitions, for zero bias, this expression takes the form The latter can be interpreted as a Drude-type result   τ = D / 2 , with relaxation time τ ν ∼ 1/ and mean free path  τ ∼ c . A similar expression has been obtained in 5,28 for a chain with noisy sites. In the other extreme of large bias Eq. (8) implies that D c  . This result, like the Drude result, can be regarded as coming from FGR transitions. But now the transitions are between Bloch site-localized states. Namely, we have hopping between neighboring sites ( ∼ 1  ), with rate of the transitions (1/τ) that is suppressed by a factor c/ 2  . The suppression factor reflects the first-order-perturbation-theory overlap of Bloch-localized wavefunctions. To summarize: we can say that D  exhibits a crossover from Drude-type transport to hopping-type transport as the field  is increased.
In fact we can proceed beyond leading order, and calculate D up to second order in η, see (Supplementary S5). Here we cite only the zero bias result: − due to the narrowing of the momentum distribution, see Methods. However the current result indicates that the leading correction is related to a different mechanism. Indeed, using a semiclassical perspective, the coupling to the bath involves a p cos( ) factor, see Methods. The zero order diffusion with rate ν arises due to stochastic term in the equation of motion for x  that involves a p sin( ) factor, Consequently, due to thermal averaging, c T (1 [1/8] , which explains, up to a factor of 2, the third term in Eq. (20). We have repeated this calculation also for a Caldeira-Leggett dissipator, and also for an (S)  dissipator. For the former the expected c T ( / ) correction to D  appears, but has a different numerical factor, while for the latter the correction comes out with an opposite sign. We can show analytically that the discrepancies are due to the modification of the correlation time 47 . Disordered ring. The so called stochastic field  T / x is responsible for the asymmetry of the incoherent transitions. Following Sinai we assume that it has a random component that is (say) box-distributed. From the works of Sinai and Derrida [30][31][32] we expect a sliding transition as  T / exceeds a critical value of order σ T ( / ) 2  . Strongly related is the delocalization transition 36,[38][39][40]43,44 for which the critical value is smaller by a numerical factor. Disregarding this factor we expect In the purely stochastic model, for   > c the relaxation is expected to be under-damped due to a delocalization transition that leads to the appearance of complex eigenvalues at the vicinity of 0 λ = . The question arises how this transition is affected by quantum coherent hopping. The naive expectation would be to witness a smaller tendency for localization in the relaxation-spectrum because we add coherent bypass that enhances the transport. But surprisingly the numerical results of Fig. 3b show that the effect goes in the opposite direction: for non-zero c, some eigenvalues become real, indicating stronger effective disorder.
Enhanced effective disorder. We turn to provide an explanation for observing enhanced effective disorder due to coherent hopping. On the basis of the non-disordered ring analysis, the relaxation modes occupy mostly the r 0, 1 = diagonals of ρ, and therefore it makes sense to exclude couplings to the higher diagonals. We verify that this does not change the qualitative picture in Fig. 3b (gray vs green symbols). The effect of the = r 1 band is to introduce virtual coherent transitions between diagonal elements (Supplementary S7). Hence we end up with an effective single-band stochastic equation with transition rates www.nature.com/scientificreports www.nature.com/scientificreports/ The disorder that is associated with x ν is hermitian, namely, it does not spoil the symmetry of the transitions, it merely implies that the we have a tight binding model with random couplings that have some dispersion w c Var ( ) . This is known as resistor network (RN) disorder. In contrast, the x η ±  term induces asymmetric transitions. This type of non-hermitian Sinai-type disorder is characterized by the dispersion Var( ) . The latter translates after non-hermitian gauge transformation to (hermitian) diagonal disorder, with ill-defined boundary conditions. The procedure to handle both types of disorder has been discussed in 36 following 39 . One defines an hermitian RN matrix by setting =  0 x  in Eq. (22). The RN matrix has a real spectrum with eigenstates that are characterized by inverse localization length that is dominated by the RN-disorder, namely, κ λ σ λ ∝ ⊥ ( ) 2 . Adding back the field  x , the eigenstates remain localized (with real eigenvalues) only in regions where localization is strong enough, that is, Supplementary S7 for an explicit expression, we deduce that the additional RN disorder is responsible for the observed numerical result.
the Wannier-Stark Ladder. We shift our attention to the full Lindblad spectrum. In the absence of coupling to the bath, the eigenstates of the Hamiltonian Eq. (2) are Bloch localized. Each eigenstate occupies a spatial region ~c/, and the corresponding eigen-energies form a ladder with spacing , that reflects the frequency of the Bloch oscillations. Weak coupling to the bath leads to damping of the Bloch oscillations. This is reflected by the Lindblad spectrum. For the = q 0 modes we have obtained Eq. (12), where we see that the eigenvalues acquire a real part, but maintain the ladder structure. But for non-zero q the c ( )  term couples the modes to the perturbation that is created at the < r 2 region by  ( ) ν , see Eqs. (31) and (32)  Regime diagram. We would like to place our results in the context of the vast quantum dissipation literature.
The prototype model of Quantum Brownian Motion (QBM), aka the Caldeira-Leggett model, involves coupling to a single bath that exerts a fluctuating homogeneous field of force. In the classical framework it leads to the standard Langevin equation and Eq. (1) becomes the standard Fokker Planck equation. In the tight-binding framework we have the identifi- , where a is the lattice constant. The standard QBM model features a single dimensionless parameter, the scaled inverse temperature β, which is the ratio between the thermal time T 1/ and damping time m/η. In the lattice problem we can define two dimensionless parameters scaled temperature (26) Accordingly β α θ = / . Note that in our model we set the units such that a 1 = , hence, disregarding 2π factor, our scaled friction parameter η is the same as α.
The standard analysis of QBM 48 reveals that quantum-implied memory effects are expressed in the regime 1  β , where a transient log(t) spreading is observed in the absence of bias, followed by diffusion. Most of the www.nature.com/scientificreports www.nature.com/scientificreports/ quantum dissipation literature, regarding the two-site spin-boson model 49 and regarding multi-site chains 50,51 , is focused in this low temperature regime, where significant deviations from the classical predictions are observed for large α of order unity. In contrast, our interest is in the α β  , 1 regime. Our η θ ( , ) regime diagram Fig. 6 is roughly divided into two regions by the line θ ∼ 1. Along this line the thermal de-Broglie wavelength of the particle is of order of the lattice constant, hence it bears formal analogy to the analysis of QBM in cosine potential 52 , where it marks the border to the regime where activation mechanism comes into action. In our tight binding model we have a single band, hence transport via thermal activation is not possible. Rather, in the 1 θ > regime, where T c  , the momentum distribution within the band is roughly flat, and the drift is dictated by Eq. (17), that is, eff where η θ η ≈ 12 eff 2 for weak field. The low temperature regime  θ 1 has not been addressed in this work, because the Ohmic master equation fails to reproduce canonical equilibration in this regime (see Methods). Still we would like to illuminate what we get (not presented) if this aspect is corrected. In the regime θ 1  the momentum distribution becomes much narrower (only low energy momenta are populated) and therefore η η ∼ eff as implied by Eq. (24). This we call Classical-like Brownian motion (CBM) regime. Once the coupling to the bath is not the simple x-coupling of the Caldeira-Leggett model, a numerical prefactor is expected (see Methods for detailed argument).

Discussion
There is a rich literature regarding Quantum Brownian motion (see 48,[53][54][55] and references within). In the condensed matter literature it is common to refer to the Caldeira-Leggett model 25,26 , where the particle is linearly coupled to a bath of harmonic oscillators that mimic an Ohmic environment. Some works study the motion of a particle in a periodic potential, possibly with bias, aka washboard potential 51,52,56 , while other refer to tight binding models 6,50,57 as in this work. The focus in those papers is mostly on non-Markovian effects: at low temperatures the fluctuations are not like "white noise", and are dominated by a high frequency cutoff c ω . Consequently the handling of long-time correlations becomes tricky. In this context the low temperature dependence of the diffusion and the mobility is modified for 1/2 α > and 1 α > . The line of study in the above models has assumed that the fluctuations that are induced by the bath are uniform in space. In some other works the dynamics of a particle that interacts with local baths has been considered. In such models the fluctuations acquire finite correlation length in space 5,[7][8][9][10][11][12][13][14]27,28 . The extreme case, as in this work, are tight binding models where the coupling is to uncorrelated baths that seat on different sites or bonds. Studies in this context assume bath that are connected at the end points 58 , or baths that act as noise source 5 . Reference 28 has analyzed the spectral properties for a chain with noisy sites, ref. 9 has considered colored noise sources strongly coupled to each site, ref. 59 has considered noisy transitions on top of coherent transitions, ref. 60 has considered transport in the presence of dephasing and disorder, ref. 11 has considered numerically transport properties of a noisy system with static disorder, and 61 has addressed some bounds in the absence of disorder. The basic question of transport in a tight-binding models has resurfaced in the context of excitation transport in photosynthetic light-harvesting complexes [15][16][17][18][19][20][21][22][23][24] .
In this article we have studied transport properties along a chain taking into account several themes that have not been combined in past studies: (a) The baths on different bonds are not correlated in space; (b) The baths are not just noise -the temperature is high but finite; (c) Without coherent hopping it is the Sinai-Derrida-Hatano-Nelson model which exhibits sliding and delocalization transitions; (d) Without baths it is a disordered chain with Anderson localization; (e) The bias might be large such that Bloch dynamics is reflected.
The "small" parameter in our analysis is the inverse temperature. The following observations have been highlighted: (1) The NESS current is the sum of stochastic and quasi-coherent terms; (2) It displays non-monotonic dependence on the bias, as shown in Fig. 2, due to crossover from Drude-type to hopping-type transport; (3) Disorder may increase the current due to convex property; (4) The interplay of stochastic and coherent transition is reflected in the Lindblad spectrum; (5) In the presence of disorder the quasi-coherent transitions enhance the localization of the relaxation modes. Thus, with regard to the Sinai-Derrida sliding transition, and the strongly related Hatano-Nelson delocalization transition, we find that adding coherent transitions "in parallel" have in some sense opposite effects: on the one hand they add bypass for the current (point (1) above), but on the other hand they enhance the tendency towards localization (point (5) above). Some of our results might be relevant to studies of optimal transport efficiency and the quantum Goldilocks effect 9 .

Methods
Master equation for disordered chain. A pedagogical presentation of the procedure for the construction of an Ohmic master equation for a two site system, and then for a chain, is presented in Supplementary S1. Each bond has a different bath, and therefore can experience different temperature and friction. Accordingly we can have disorder that originates either from the Hamiltonian (say random  x as assumed in the main text), or from having different baths (random ν x or random x η ). This extra disorder can be incorporated in a straightforward manner, and does not affect the big picture.
thermalization. The Ohmic master equation for a Brownian particle, if the coupling to the bath were −xf t ( ), is the standard Fokker-Plank equation. It leads to canonical thermal state for any friction and for any temperature. This is not the case for a discrete two level system. The agreement with the canonical result is guaranteed only to first order in η. This is reflected in Eq. (6). The same applies for a chain. Note that in this sense the Ohmic approximation is very different from the secular or Pauli approximation 62 , or specially constructed Davies Liovillian 63 , that guarantee canonical thermalization.
It is important to identify the "small parameter" that controls the deviation from canonical thermalization. The standard coupling via x induce transitions between neighboring momenta, and therefore the small parameter is Δ/T, where the level spacing Δ goes to zero in the L → ∞ limit. But for local baths, the coupling is to x x ( ) δ − α scatterers, that create transitions to all the levels within the band. Therefore the small parameter in the absence of bias is c T / . This assertion is confirmed numerically by inspection of the equilibrium momentum distribution p k ( ), see Supplementary S5. We conclude that in the regime of interest ( 1 θ > ) the Ohmic master equation can be trusted, while for lower temperatures we have to "correct" it. For the two site system the "corrected" equation is the Bloch equation, where the ratio of the rates is in agreement with detailed balance (not just in leading order in η). For a chain, we cannot justify the secular approximation, and therefore the correction procedure is ill defined.
The friction coefficient. In the Caldeira-Leggett model for Brownian particle, with interaction term −xf t ( ), the bath induced fluctuations f t ( ) are determined by a coupling constant η, and by the bath temperature T , such that at high temperatures the intensity of the fluctuations is The η parameter is defined such that the friction coefficient in Eq. (24) equals η. A straightforward generalizations 27 shows that for interaction with local baths For homogeneous distribution of local baths that have the same η, the friction coefficient becomes x-independent. In the model under consideration the coupling to the baths is W f t ( ) ∑ α α α , where α labels the sites. Disregarding commutation, it can be written as x u f t p 2 cos( ) ( ) ( ) ∑ α α α , where the u-s are site localized. It follow that the effective friction coefficient is momentum dependent, namely η η ∼ p cos( ) eff 2 . But for equilibration in the 1 θ < regime only low momenta are important hence we expect, up to numerical prefactor to observe eff η η ≈ . The failure to observe this result is due to improper thermalization, as discussed in a previous paragraph. positivity. Irrespective of η, there is another complication with the Ohmic master equation. If the temperature is low (small ν) the relaxation may lead to a sub-minimal wavepacket that violates the uncertainty principle. This reflects the observation that the Ohmic master equation is not of Lindblad form, and violates the positivity requirement. The minimal correction required to restore positivity is to couple V to an extra noise source of intensity This term is essential in the low temperature regime. We have verified numerically that the extra noise term can be neglected in the high temperature regime where our interest is focused. www.nature.com/scientificreports www.nature.com/scientificreports/ on-site dissipators. The model of 28 combines Hamiltonian term with dissipator of the  (S) type that originates from couplings via = W Q : Such dissipator leads to off-diagonal dephasing that is generated by Q Q x x ρ terms, and therefore excludes the possibility for inter-site stochastic transitions. Similar remark applies to the familiar Caldeira-Leggett model of Quantum Brownian motion 25 , where the coupling is via W x := . Namely, it is a single bath that exerts a fluctuating homogeneous force that affects equally all the sites, as in 50 . In our model the dissipation effect is local: many uncorrelated local baths. , are in agreement with FGR. For completeness, we added here a γ term that represents optional off-diagonal decoherence due to on-site noise.
Ring configuration. For numerical treatment, and for the purpose of studying relaxation dynamics, we close the chain into a ring. This means to impose periodic boundary conditions. With uniform field , one encounter a huge potential drop at the boundary. To avoid this complication we assume that the boundary bond has an infinite temperature, hence the formation of a stochastic barrier is avoided, and the circulation of the stochastic field (aka affinity) becomes  T / as desired. In the analytical treatment of a clean ring we assume that ν, and η and  in the master equation are all uniform, such that invariance under translation is regained. This cheat is valid for large ring if ρ is banded, reflecting a finite spatial correlation scale. See numerical verification in Supplementary S2. the Bloch eigenstates of a clean ring. The  super-operator in the Bloch r q ( ; ) basis, decomposes into q Blocks. Each block can be written as a sum of terms, as in Eq. (9), that operate over a one dimensional tight binding r lattice. In this representation the coherent dynamics is generated by www.nature.com/scientificreports www.nature.com/scientificreports/ the relevant dimensionless parameter is ν T / and not c T / . Figure Supplementary S4 of Supplementary S5 confirms this statement.