Resistance saturation in semi-conducting polyacetylene molecular wires

Realizing the promises of molecular electronic devices requires an understanding of transport on the nanoscale. Here, we consider a Su-Schrieffer-Heeger model for semi-conducting trans-polyacetylene molecular wires in which we endow charge carriers with a finite lifetime. The aim of this exercise is two-fold: (i) the simplicity of the model allows an insightful numerical and analytical comparison of the Landauer and Kubo linear-response formalism; (ii) we distill the prototypical characteristics of charge transport through gapped mesoscopic systems and compare these to bulk semiconductors. We find that both techniques yield a residual differential conductance at low temperatures for contacted polyacetylene chains of arbitrary length—in line with the resistivity saturation in some correlated narrow-gap semiconductors. Quantitative agreement, however, is limited to not too long molecules. Indeed, while the Landauer transmission is suppressed exponentially with the system size, the Kubo response only decays hyperbolically. Our findings inform the choice of transport methodologies for the ab initio modelling of molecular devices.


Introduction
The last decades have witnessed tremendous experimental progress in the field of nanoelectronics, pushing investigation of transport properties toward the ultimate scale of single-molecule architectures [1][2][3][4].At the nanoscale, semiclassical approaches break down as the nature of electron transport is inherently quantum mechanical [5,6], showing phenomena ranging from conductance quantization [7] and quantum interference [2,8], to the anomalous quantum Hall effect [9].At the length scale of individual molecules in, e.g., a mechanical break-junction setup [2][3][4], electron transport is essentially ballistic.However, in complex quantum devices, both the coupling to the environment and to internal degrees of freedom of the molecular bridge entail sources of incoherence for the charge carriers.A fundamental theoretical understanding of the mechanisms of quantum transport and the role of incoherence effects across multiple length scales is challenging [4,10].At the same time, it is of pivotal importance to the interpretation of experiments [3] and to eventually harvest the potential of quantum technologies for next-generation nanoelectronic devices.
Among the most prominent methodologies to describe transport phenomena on the nanoscale are the Landauer-Büttiker formalism and Kubo linear response [11][12][13][14].In fact, it has been established that the Landauer formula [12] itself can be derived from linear response [15][16][17].It was motivated that the omission of higher order terms (vertex corrections) in the Landauer transmission is admissible even for interacting electrons when, at zero temperature, inelastic processes subside [18].In general, however, the current has to be recast in terms of an effective transmission kernel, where the Green's function is dressed by a many-body selfenergy and the molecule-lead coupling is renormalized by vertex corrections [19][20][21][22].While in practice evaluating such vertex corrections is challenging and thus not often attempted, numerical results seem to suggest [19] that offresonant, near-equilibrium transport can still be qualitatively described in terms of quasi-particle scattering.
Using the Kubo formula, it was recently shown [23,24] that finite lifetimes of charge carriers can lead to non-trivial transport phenomena in bulk correlated narrow-gap semiconductors [25].From this perspective, it is interesting to ask how quasi-particle incoherence manifests in transport through semiconducting nanoscopic systems.The similarities of the Landauer and the Kubo formalism for coherent transport further raise the question how differently the charge carriers' finite lifetimes are reflected in the two approaches.
To answer these questions, we investigate in this paper the transport through a simple semiconducting molecule, transpolyacetylene.For the charge carriers, we impose an ad hoc finite lifetime that serves to mimic scattering events taking place within the molecule, e.g., from electron-electron scattering, disorder or electron-phonon coupling.We find quasi-particle incoherence to lead to significant, qualitative changes in the molecule's conductive properties.In particular, we find the resistance to saturate below a surprisingly large characteristic temperature.Further, we reveal marked differences between the Kubo and the Landauer results.In fact, we motivate-numerically and analytically-that the Kubo formalism is more general, in the sense that it contains the Landauer conductance as its small-scattering rate limit.Our results advocate that future investigations of semiconducting organic systems-and devices build from them-need to be mindful of the influence of incoherence onto transport properties-even at room temperature.

System and model
Polyacetylene (PA) is an organic conjugated polymer with chemical formula (CH) n .Of particular interest is the structure of trans-polyacetylene which dimerizes due to electron-phonon coupling through a Peierls mechanism.Thus, trans-PA consists of alternating single ( C−C ) and double ( C=C ) carbon bonds.Experimental investigations suggest bond-lengths d C−C = 1.44 Å and d C=C = 1.36Å [26,27], with a relative angle ≈ 122 • , yielding a unit cell length of d = 2a = a 1 + a 2 ≈ 2.45Å, where a is the lattice spacing of the undimerized chain, whereas a 1 and a 2 are the projection of the single-and double-bonds along the chain, respectively.Following the seminal work of Su, Schrieffer, and Heeger [28], a single dimerized PA chain with N unit cells can be described within an effective tight-binding model for -electrons with nearest-neighbor hopping where c ( †) i is the annihilation (creation) operator for an electron on a C-p z atomic orbital (AO) i of spin , while t 0 = 2.5eV is the hopping integral for the undimerized chain, = 4.1eV∕ Å is the electron-phonon coupling con- stant, and u i = (−1) i u 0 is the displacement describing the dimerization, with u 0 = 0.04 Å.With these parameters, the dimerized hoppings are given by t 1 = t C−C = 2.172eV and t 2 = t C=C = 2.828eV , which yields a Peierls gap Δ ∞ = 2(t 2 − t 1 ) = 8 u 0 = 1.312eV [28].
The unit cell and the bond structure of PA, including the lattice spacing and the hopping, are shown in Fig. 1a, b.The Peierls gap between the highest occupied (HOMO) and the lowest unoccupied (LUMO) molecular orbitals is found to decay as 1/N towards the bulk Peierls gap Δ ∞ , see Fig. 1c.The representative distribution of molecular orbital (MO) eigenenergies are shown in Fig. 1d for a N = 32 PA chain.The spectral function of this PA chain, shown Fig. 1e, closely resembles the one of the infinitely-long PA chain (periodic bulk).For the components projected onto individual C's AO in Fig. 1f, those towards the centre of the chain indeed resemble the spectral function of an infinitely long (dimerized) one-dimensional chain, delimited by its characteristic van Hove singularities.The spectrum close to the edges of the chain, however, assumes a shape that is semi-circular outside the hybridization gap.Analytical arguments (for the undimerized chain) supporting the numerical spectra are discussed in Appendix 1.All projected spectral functions are even, i.e., ℑG ii ( ) = ℑG ii (− ) , which is a manifestation of the particle- hole symmetry of the model.Moreover, since the lattice is bipartite, the SSH Hamiltonian is invariant under a sublattice (chiral) symmetry which maps AOs i → 2N + 1 − i.
The SSH model also raised interest in the community for its topological properties [29][30][31].Here, we are taking the intra-cell and inter-cell hopping to be t 2 and t 1 , respectively, with t 2 > t 1 .With this choice, the system is topologically trivial.However, PA can dimerize in another pattern, in which single-and double-bonds are interchanged, leading to a nonzero winding number and a topologically nontrivial state.Then, the SSH model supports zero-energy edge states and can host domain walls, giving rise to solitons [29,30].As a consequence, for varying termination, the PA chain displays dramatically different transport properties [32,33], a complexity beyond the scope of this work. (1)

Electron transport theory
In a typical transport (two-terminal) setup, the system is divided into three regions, in which a central (or scattering) region ( C ) bridges left ( L ) and right ( R ) electrodes.The electrodes are connected at infinity to particle reservoirs which account for the dissipation necessary to establish a steady state.Each reservoir injects thermalized carriers which are distributed according to the Fermi-Dirac statistics at the equilibrium chemical potential L and R , corresponding to a bias voltage V b = ( L − R )∕e , whereas the transport in the central region is phase-coherent.Then, the Hamiltonian has the following form where H C describes the scattering region, H L and H R describe the L and R lead, respectively, and the off-diagonal elements account for the coupling between these regions.Typically, it is required that the two leads are coupled only through the scattering region, i.e., H LR = 0.1 (2) , Since the leads are infinitely extended, a Hamiltonian formulation of the scattering problem is not numerically feasible.Hence the leads are instead accounted for by embedding them into the scattering region.The retarded Green's function of the scattering region then reads where the embedding self-energy of lead is defined as in terms of the corresponding surface Green's function There are different possible routes to take into account the electrodes [34,35], and also the AO beyond the ones of the -system [36].For the sake of simplicity, in the following, we assume that the molecule-lead coupling is restricted to the AO at the edges of the PA chain ( and r, respectively), and we employ a wide-band approximation [35] for the leads.Then, the embedding self-energy of each lead reduces to an imaginary constant, i.e., ( L ) = − Γ L ∕2 and ( R ) rr = − Γ R ∕2.
In the following, we recap the Landauer and Kubo theory for electron transport within the Green's function formalism.Making the formulae explicit for the SSH model will facilitate to highlight similarities and differences between the approaches.

Landauer
Within the Landauer formalism, the electron transmission function through the scattering region is given by where encloses the spectral information of lead .For a deeper understanding of the channel structure of the electron transport, it is useful to make the trace in the Landauer formula (7) explicit, yielding Considering the approximation on the structure of the leads from above, (Γ L ) ii � ∝ i i � and (Γ R ) jj � ∝ jr j � r , this reduces to Hence, the Landauer transmission is controlled by the (nonlocal) Green's function that links the and r AOs, corresponding to the outer C atoms which are connected to the L and R lead, respectively.
Although a derivation from linear response is possible [15][16][17]19], the Landauer transmission, Eq. ( 10), is de facto agnostic to the applied perturbation as it includes no information on external couplings.Indeed, Landauer describes the amplitude (absolute square) for the process of adding an electron on one end of the molecule (AO r) and removing it from the other (AO ), In this tunneling process (which is independent of spin σ), inequivalent bonds or effects of dimerization only manifest implicitly, via their effect onto the hopping amplitudes in the Hamiltonian.
Finally, we recall that, in the presence of many-body correlations within the scattering region, the Landauer transmission, Eq. ( 7) is approximate.Besides considering the many-body Green's function, one also needs to include a vertex correction [19,22,41].Following Ferretti et al. [41], the corrected expression reads where and C = −2ℑ C is defined in terms of the many-body self- energy of the scattering region.Considering the matrix structure of the leads' and many-body self-energies, it is easy to verify that, including the vertex corrections, Eq. ( 10) becomes For all interpretational purposes, the vertex corrections correspond to an asymmetric renormalization of the moleculelead coupling, introducing an effective many-body electrode in which the electrons undergo some scattering processes, before being re-injected into the system [22].

Kubo
Within the Kubo formalism, the formula for the optical transmission along the x direction (when neglecting vertex corrections [71][72][73][74][75]) is given by [76] (11) where L = Nd − a 1 is the length of a PA chain containing N unit cells,2 the Green's function including the embedding self-energy of the leads is given by Eq. ( 3) and the matrix elements of the velocity are given by For lattice models, such as the SSH Hamiltonian Eq. ( 1), these velocities are typically derived from the Peierls substitution approach [78][79][80], in which the vector potential couples to the atomic positions (projected onto the direction of transport).
Making the trace explicit as above, and expressing the velocity in terms of the hopping, yields For dimerized PA chains within the SSH model we get Contrary to Landauer, the Kubo approach, more manifestly, describes the response to a specific external perturbation.Further, the current vertices, or velocities, Eq. ( 16), have direct knowledge of the internal structure of the molecule, e.g., the dimerization and the buckling angle , as the electromagnetic vector potential couples to atomic positions.( 15) T

Green's functions and transmission paths in polyacetylene
Processes that contribute to the amplitude of the Landauer and the Kubo transmission are illustrated in Fig. 2a and b, respectively: The Landauer transmission, Eq. ( 9), is completely determined by the nonlocal Green's function between the AO connected to the L and R reservoirs.The Kubo response, instead, consists of many contributions.Each contribution is built from nearest-neighbor hopping processes t i,i±1 and spectral functions ℑG ij .In particular, some contri- butions involve local spectral functions ℑG ii (i.e., projected on a single AO i) which are highlighted by colorful lines in Fig. 2b.
Owing to the particle-hole symmetry of the SSH Hamiltonian, the Fermi energy falls in the middle of the HOMO-LUMO gap (and hence will be set to E F = 0 from here on) and the Green's function components are either even or odd functions of frequency, see Fig. 3. Specifically, the following relations hold [66] This in particular implies ℑG 1,i+2n+1 (0) = 0 , so that for all index pairs (i, j) with opposite parity (i.e., belonging to different sublattices) the corresponding contributions to the Kubo response vanish.Interestingly, at the Fermi energy, the contribution of ℑG r is absent in the Kubo response.As is clearly visible in Fig. 3, the electron transport is dominated by the even components of the Green's function-ℑG ii in the case of Kubo, and ℜG r for Landauer.Crucially, only the former displays a notable dependence on the scattering rate.Indeed, ℑG ii ∝ Γ .As we shall see in the following, this is the most relevant difference between Kubo and Landauer.

Numerical results
We now compare numerical data obtained within the Landauer and the Kubo formalism (see Figs. 4 and 5) and we show that there exist distinct transport regimes, depending on the length of the PA chain N (unit cells) and the scattering rate Γ .Throughout, we set the molecule-to-lead coupling to Γ L = Γ R = 2.5 × 10 −3 eV, which is representative of a weak chemical bond between PA and metallic (e.g., Au) electrodes through thiol (SH) or amino (NH 2 ) anchoring groups.We first discuss the transmission function, which, generally, displays resonances for poles of the Green's function, corresponding to the MO of the scattering region.In the following, we will focus on the transmission T(E F ) directly at the Fermi level E F .The transmission, as well as the conduct- ance, G = e 2 ∕h × T(E F ) , defined in terms of the electron charge e and Planck's constant h, is mainly controlled by the value of Γ L .While this is evident from the Landauer formula (10), it is less obvious, yet equally true, in the case of the Kubo formula (15), as we show in Sect.5.1.1.

Variation with chain length N
Within Landauer, we observe an exponential decay in conductance where L = Nd − a 1 is the length of the PA chain, the pre-factor Landauer is representative of the resist- ance at the molecule-lead interface, and is the attenuation factor.The numerical data are compatible with Landauer ∼ Γ 2 L and = 0.216Å −1 .The attenuation factor is therewith in fair agreement with the empirical law [81] (21) G = e 2 h Landauer exp(− L), Further, the Landauer conductance is insensitive to Γ .Indeed, for Γ ≪ Γ L , the electron transport is phase-coherent and the Landauer is identical to the Kubo response for short chains, see Fig. 4b  and d.For longer chains, scattering processes within the Kubo formalism become dominant and there exists a crossover length scale L c = N c d − a 1 at which the Kubo response displays a gradual transition between an exponential to a hyperbolic regime, for which we find as expected for non-directional diffusion [82].The filled circles in Fig. 4d are a guide-to-the-eye to identify the crossover length between the exponential and the hyperbolic decay.The crossover scale L c decreases as Γ increases.The expo- nential scaling also disappears for Γ ≳ Γ L when transport is dominated by scattering at all length scales.

Temperature dependence
The phenomenological ansatz for finite lifetimes, Eq. ( 6) essentially describes electron-electron scattering between Landau quasi-particles.Within this approximation, the electron Green's function of the scattering region and, hence, the transmission function, do not depend on temperature.Still, there is an effect of temperature on the transport properties, if, instead, we consider the electron current, defined as where f ( ) is the Fermi-Dirac distribution of the lead, which governs thermal activation.The resistance is then obtained as the inverse of the differential conductance Within the linear response regime eV b ≪ Δ , it is safe to assume a weak dependence of the transmission function on the external bias, i.e., that T( , V b ) ≈ T( ) .Moreover, in this regime the current is linear in the external bias, I ∝ V b , and the resistance is bias-independent.In Fig. 5 we show the resistance R as function of temperature for PA chains of different lengths N. At high temperatures, we observe an activated behavior corresponding to transitions across the HOMO-LUMO gap Δ .This regime is bounded from below by a characteristic temperature scale, T * ∝ Δ∕k B that is also controlled by the gap (see the discus- sion below).The saturation resistance below, R(T < T * ), is dominated by the transmission at the Fermi energy T(E F ) and thus follows either the exponential (Landauer) or the hyperbolic (Kubo) decay with the chain length.Noteworthy, as Γ → 0 , the Kubo resistance converges towards the Lan- dauer result, cf.Fig. 5a and c.

Discussion
We will now discuss the above numerical results for the finite PA molecules.First, using analytical arguments, we will rationalize the evidenced dependencies on the molecules' length N and scattering rate Γ .Then, we will put the observed saturation of the molecules' resistance into relation to recent theories for transport in correlated narrow-gap bulk semiconductors.

Landauer versus Kubo: analytical insights
The differences between the Landauer and the Kubo transmission function T(E F ) simulated for low energies and tem- peratures shown in Fig. 4 can be understood from analytical arguments: We symbolically evaluate T(E F ) from Eqs. (7)  and ( 15) and investigate the leading behaviour vis-à-vis the scattering rate Γ and the molecule-lead couplings Γ L = Γ R .

The small scattering limit
For a molecular wire with N PA unit cells one finds for the Landauer response, Eq. ( 10), to lowest order in Γ L and Γ with n ≥ 2, n + m ≥ 4 .For the Kubo response, Eq. ( 15), the leading order comes out as exactly the same, albeit with p ≥ 0 , p + q ≥ 2 .Terms beyond the leading order will be discussed in Sect.5.1.2.The above expansions readily explain why the simulated Landauer and Kubo transmissions (Fig. 4) coincide in the limit of vanishing Γ , irrespective of the chain length N. In fact, in both cases, the lowest order term in the transmission is proportional to Γ 2 L and independent of Γ .For the typical case Γ L ≫ Γ , the above leading term is a good approximation (lower bound) for the transmission T(E F ) in the zero tempera- ture limit, where the scattering rate Γ is expected to become small.Crucially, a finite coupling Γ L leads to a residual trans- mission and, hence, a finite conductance G(T → 0) > 0 , which will be discussed in Sect.5.2 below.

Sizable scattering 0
Above a crossover value N c , we observe a strong deviation between the Landauer and the Kubo transmission, with the Kubo one decaying much more slowly for N → ∞ .The Γ -dependence of this phenomenon (see Fig. 4a, c) provides a definite clue as to its microscopic origin.Indeed, while the leading order terms of Landauer and Kubo (discussed above) coincide, terms beyond the low-scattering limit are different, see the different powers in the corrections indicated in Eqs. ( 25), (26).
By construction, the lowest order in scattering within Landauer is Γ 2 L , i.e., there are no pure-Γ terms and no mixed terms ∝ Γ L Γ .Within Kubo, on the other hand, there are both pure-Γ terms (the lowest order is Γ 2 ) and a mixed term Γ L Γ .Plausibly assuming Γ L ≫ Γ , this mixed term will dominate the correc- tions to the leading ∝ (Γ L ) 2 term in the Kubo transmission Eq. ( 26).Expanding to this next order, we find with, for even N, merely decays hyperbolically.For small Γ and N, the transmission is dominated by c 1 (N)Γ 2 L leading to T Kubo ≈ T Landauer .Owing to the different scaling with system size, however, the c 2 (N)Γ L Γ term (absent in Lan- dauer) is gaining importance for longer molecules, all the more rapid the larger Γ .Indeed, for a given Γ , there must be a crossover length N c d , above which the c 2 (N)Γ L term surpasses c 1 (N)Γ 2 L .For N > N c the Γ L Γ term then dominates the transmission, turning the molecule's scattering rate Γ into a relevant energy scale and causing a strong deviation from Landauer's transmission.
While the Γ L Γ term is absent in the plain Landauer for- malism, the vertex corrections to the lead-molecule coupling introduce such a term.However, as is evident from Eq. ( 14) the vertex corrections only renormalize the pre-factor, not the structure and length-dependence of the tunneling processes.Therefore, the vertex-corrected Landauer transmission still decays exponentially with system size, see Fig. 4.
The reason for the very slow decay of the Kubo conductance with the chain length is the incoherence (finite lifetime) of charge carriers inside the molecule.Indeed, as discussed in Sect.3.3, the ingredients entering the transmission function are significantly more Γ-dependent within the Kubo formalism than within Landauer's.Intuitively, the finite lifetimes of the HOMO and LUMO states lead to an energy broadening, which smoothes the gap edges and spills spectral weight into the gap [83,84].These incoherent carriers can then diffuse through the molecule, which is described by the sequential hopping processes that make up the Kubo response (visualized in Fig. 2).In this sense, the length scale L c = N c d − a 1 (see Sect. 4.1) separates regimes in which molecular transport dominantly occurs through tunneling ( L ≲ L c ) and diffusion ( L ≳ L c ), respectively.

Resistance saturation in semi-conducting molecules
At elevated temperatures, the resistances in Fig. 5 exhibit, as expected, an Arrhenius regime, in which conduction through the semi-conducting molecular chain is activated.The largest resistance is then obtained for the shortest molecules, owing to their larger gaps.Cooling below a crossover temperature T * , however, the resistance from, both, Landauer and Kubo cedes to increase exponentially and, instead, levels off to saturate towards T → 0 .This saturation phenomenon has first been discussed [23] for correlated narrow-gap semiconductors [25], such as Kondo insulators.In these periodic bulk systems, the low-temperature regime is dominated by intra-band transitions [24]: For a band 0 k endowed with a quasi-particle lifetime ℏ∕(2Γ) and weight Z, the Kubo conductivity (without vertex corrections) can be expressed analytically as [85] with the unit cell volume V, the intra-band group velocity with k in the Brillouin zone and derivatives of the digamma function where k = Z 0 k .The characteristic temperature T * that delimits the resistivity saturation regime encoded in the Kubo Eq. ( 31) can be crudely estimated as [23] where Δ is the renormalized charge gap. 3 Note that the crossover from an activated behaviour to a low-temperature conductance regime with weak temperature dependence was recently also suggested based on the Meir-Wingreen formula applied to a single molecular level [86].
The digamma function accounts for thermal ( k B T ) and lifetime ( = ℏ∕(2Γ) ) broadening on an equal footing.For infinitely long-lived charge carriers, Γ → 0 , Eq. ( 31) reduces-to leading 1∕Γ order-to the Boltzmann conduc- tivity, since [85] In that semi-classical limit, the scattering rate Γ becomes a mere prefactor and the response is activated for all temperatures.In that sense, the resistance saturation found, here, in molecules within both the Kubo and the Landauer formalism is a quantum effect.
The above formulas in Eqs.(31,32) describe charge transport from transitions taking place within the same band k that disperses owing to unit cell-to-unit cell hopping.Already for periodic lattice models with multiple atoms per unit cell, one has to extend this setting, allowing in particular for inter-band and intra-unit cell transitions.For the Fermi velocities, this is achieved in the generalized Peierls approach [79], in which (31) where (H k ) nm is the Hamiltonian expressed in a local basis, with n = (i, l) indexing, both, the hosting atomic site i and orbital l. 4 The above formula has the virtue of "interpolating" between the momentum-space description of a periodic system (first term: inter-unit cell transitions) and the large real-space unit cells with open boundary conditions of finite molecules (second term: intra-unit cell transitions).Indeed, Eq. ( 35) assures that transport observables for a periodic solid with a primitive one-atomic unit cell can be equivalently described by non-primitive unit cell that has been, say, doubled in the x direction.
In finite systems (i.e., with open boundary conditions), such as the SSH model Eq. ( 1), the real-space formulation manifestly only involves inter-orbital inter-atomic transitions and the Fermi velocity reduces to the second term, and with (H k ) nm = −t nm for k = (0, 0, 0) , to Eq. ( 16).Likewise, Eq. ( 31) has to be replaced with the (lengthy) expression for inter-band transitions, which can be found in Eqs.(3,11,29) of Ref. [37].Using the same procedure as above, we can crudely estimate the dependencies of the resistance saturation regime, finding A comparison to Eq. ( 33) reveals that the saturation regime, T < T * , is roughly larger by a factor of two for gapped extended systems than in periodic semiconductors with the same gap Δ .According to Eq. ( 36), the dominant control parameter for T * is the charge gap Δ , explaining why shorter chains (with their larger Δ , see Fig. exhibit a basically flat resistance up to far beyond room temperature.The scattering rate Γ only has a sub-leading effect on T * , in congru- ence with the numerical data in Fig. 5c.From the arguments presented in Sect.5.1.2for the transmission function we further understand, that the Kubo resistance converges the Landauer result in the limit Γ → 0.

Conclusion
Using realistic parameters to model semi-conducting polyacetylene molecular wires, our results suggest that a resistance saturation regime, extending up to at least room-temperature, should be commonly observed-provided that extrinsic factors (disorder, doping) do not destroy the charge gap.The residual conduction is provided by incoherent spectral weight (36) inter-band: spilling into the HOMO-LUMO gap.We find that the extent of the temperature regime, 0 ≤ T ≤ T * , that exhibits satura- tion is directly controlled by the size of the gap Δ .Thus, for a given temperature, recovering a more conventional activated (Arrhenius) type of conduction should be possible by straining the molecular wire [86,89].This strain tuning of the crossover temperature T * separating both regimes is analogous to the pressure tuning of transport in, e.g., Kondo insulators [23,90,91].The only conceptual difference between resistivity saturation in these 3D bulk systems and resistance saturation in the extended 1D wire is conduction being dominated, respectively, by intra-and inter-band transitions that live on different energy scales, namely ∼ Δ∕2 and ∼ Δ .Thus, for a given gap, the resistance plateau in extended systems is much wider in temperature than for periodic systems.We further stress that the resistance plateau is sizable even for extremely small scattering rates.It is thus not inconceivable that already the inclusion of electron-phonon scattering from zero-point motion could be sufficient to create it.While results from the Landauer approach and Kubo linear response theory oftentimes appear to be similar, we find qualitative differences to abound when the scattering rate is non-negligible and the molecule is long.In particular, the Landauer (Kubo) conductance decays exponentially (hyperbolically) with system size.At low temperatures, the Kubo response is far more sensitive to the charge carriers' incoherence than the Landauer approach.For a given scattering rate, even if the crossover temperature T * is comparable, the satu- ration value of the Kubo resistance can be orders of magnitude smaller-especially for long molecules.Experimentally observing the length-induced crossover from ballistic to a diffusive transmission in semiconducting polyacetylene wires might be challenging, because the effect of electronic scattering manifests only for chains of considerable length.
On the level of methodology, a clear hierarchy of transport approaches for carriers with finite lifetimes emerges from our discussion.In order of decreasing lifetimes = ℏ∕(2Γ) , the methodology required to fully capture signatures of electronic incoherence escalates as follows: Boltzmann → Landauer → Kubo.Indeed, for a diminishing scattering rate Γ , results from Kubo converge to those of Landauer.For infinitely long-lived quasi-particles, the coefficient of the leading 1∕Γ order reduces to the Boltzmann expression.In all, our model calculations provide guidance for the choice of transport methodology in future ab initio simulations for semi-conducting molecular systems.

Appendix A The macroscopic limit of the linear undimerized chain ( t 1 = t 2 )
We find it instructive to derive relevant elements of the Green's function for the linear undimerized chain that is obtained from the SSH model for t 1 = t 2 ≡ t in the limit of long chains.Note that the latter (in contrast to the SSH model) will be metallic.For an undimerized chain of N sites the Hamiltonian has the form of a symmetric tridiagonal Toeplitz matrix matrix notation the corresponding Green's function then is (G [N] ) −1 (z) = zI [N] − H [N] .While the general inversion of (G [N] ) −1 in terms of Chebyshev polynomials of the second kind is well-known [92,93], we will here follow two different paths.
In summary, the results are: In the macroscopic limit, N → ∞ , the density of states in the middle of the molecule, unsurprisingly, converges to that of the gapless periodic 1D bulk system.The DOS at the extremities of the very long but finite and non-periodic one-dimensional chain becomes semi-circular.Intriguingly, the same DOS is obtained for the infinite dimensional Bethe lattice [45,94].

A.1 The ends of the long linear chain
For Green's function elements at (or connecting) the borders, we follow arguments from Schur's complement.If with matrices F 11 and F 22 that are of square shape, and, fur- ther, F 11 and B = F 22 − F 21 F −1 11 F 12 that are invertible, one can show that Setting up F [N+1] = (G [N+1] ) −1 , the inverse Green's function for N + 1 sites, from F 11 = zI [N] − H [N] = F [N] ∈ ℂ N×N , F 22 = z , a n d F 21 = (0, 0, ⋯ , 0, −t) ∈ ℝ 1×N , w i t h F 12 = (F 21 ) T ∈ ℝ N×1 , one easily finds the recursion relations and (A.1) For the diagonal component, Eq. (A.4), one easily solves the recursion in the N → ∞ limit and, interestingly, finds for the infinite chain with open boundary conditions.With z = + i , to obtain the non-interacting retarded Green's function, one then sees that the DOS at the borders of the infinitely-long one-dimensional chain is semi-circular, as is also the case for the infinite dimensional Bethe lattice [45,94].

A.2 The centre of the long linear chain
For information on the centre of the linear molecule, we devise a different, physics-oriented, strategy: Using a chain with an odd number of atoms, N = 2m + 1 with m ∈ ℕ , we consecutively integrate out the outermost sites to obtain an effective Green's function for the center three sites: Using the path-integral formalism, it is quite elementary to show, for N ≥ 5 , the recursion formula z N = z − t 2 ∕z N−2 .In the limit of large N, the solution is Then, for the centre of the molecule whence which is the same as the local DOS of the dispersion, (k) = −2t cos(ka) , of a periodic one-dimensional system with nearest-neighbour hopping t. (A.5) (G [N+1] ) 1,N+1 = − t(G [N] ) 1,N z − t 2 (G [N] ) N,N N≥2 = (G [2] ) 1,2 These results for the (metallic) undimerized chain could be generalized to the (gapped) dimerized chain.Basically, there will be an alternation of hoppings in the continued fraction recursions, that will lead to the gapping and results very close to the numerical data shown for N = 32 in Fig. 1f.

Fig. 1 a
Fig. 1 a Representations of a two-terminal junction, with a PA chain bridging Au electrodes through thiol (SH) anchoring groups.b Mapping to an SSH model, single (C−C) and double (C= C) bonds correspond to lattice spacing a 1 and a 2 , respectively (with d = a 1 + a 2 the length of the unit cell) while Γ L and Γ R describe the left and right molecule-lead coupling, respectively.c Evolution of the HOMO-LUMO gap as a function of the PA length N (unit cells) converging

Fig. 2 Fig. 3 Fig. 4 a
Fig. 2 Schematic representation of the electron transport across the PA wire within the (a) Landauer and (b) Kubo formalism.The transmission amplitude, Tr(MM † ) , consists of processes M and M † that can be associated with paths in real-space (vertical axis).The Landauer transmission is given by a single path (i.e., from to r).The Kubo

Fig.
Fig. Resistance R(T) within the (a) Landauer without (solid line) and with (dashed line) vertex corrections, and (b, c) Kubo formalism.a, b The saturation value below the length-dependent crossover temperature R(T < T * ) is dominated by T(E F ) .The filled circles in (a, b)