## Electrothermal Transistor Effect and Cyclic Electronic Currents in Multithermal Charge Transfer Networks

Galen T. Craven 1 and Abraham Nitzan 1, 2

1 Department of Chemistry, University of Pennsylvania, Philadelphia, PA 19104, USA 2 School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel

A theory is developed to describe the coupled transport of energy and charge in networks of electron donor-acceptor sites which are seated in a thermally heterogeneous environment, where the transfer kinetics are dominated by Marcus-type hopping rates. It is found that the coupling of heat and charge transfer in such systems gives rise to exotic transport phenomena which are absent in thermally homogeneous systems and cannot be described by standard thermoelectric relations. Specifically, the directionality and extent of thermal transistor amplification and cyclical electronic currents in a given network can be controlled by tuning the underlying temperature gradient in the system. The application of these findings toward optimal control of multithermal currents is illustrated on a paradigmatic nanostructure.

Cite as: G. T. Craven and A. Nitzan, Phys. Rev. Lett. 118 , 207201 (2017) DOI:

10.1103/PhysRevLett.118.207201

The coupling between electron transfer (ET) and transport and the underlying thermal environment is a long studied subject [1]. Its manifestations in recent studies of transport in molecular junctions mostly focus on the weak electron-phonon coupling regime. Similarly, thermoelectric phenomena in such junctions [2], where molecules connect between electrodes of different temperatures, are usually treated (with a few exceptions, e.g., Refs. [3-12]) with electron-vibration interaction disregarded. This stands in contrast to electron transfer reactions in condensed molecular systems that are usually dominated by hopping between thermally equilibrated polaron-like states as described by Marcus theory [1316] (analogous kinetics in junction transport is known, mostly in so called redox molecular junctions [17, 18]). We have recently considered the latter type of electronic transport in thermally heterogeneous systems, where an electron hops between two sites of different local temperatures [19]. This study was motivated by recent advances in measurement and control of temperature differences on length scales comparable to those involved in molecular electron transfer processes [20-25]. The corresponding ET rate was obtained as a modified multidimensional Marcus expression that depends on the local temperatures of the two sites. Furthermore, electron hopping was shown to to be accompanied by heat transfer between the sites whose magnitude depends on the temperature difference and on the reorganization (polaron formation) energies at the two sites.

Thermal inhomogeneity can develop spontaneously in driven nonequilibrium systems [26-30] or can be externally controlled as in a thermoelectric device. Considering such models, several recent theoretical studies have discovered interesting thermal transistor effects, whereupon the flux between two sites can be controlled by the temperature on a third site [31-33].

FIG. 1. Network graphs for R L , , and K topologies. Each node represents a donor-acceptor site. As shown explicitly in the R graph, each site is in contact with an independent thermal bath. The center panel shows a schematic of the energy surfaces E a (solid curve) and E b (dashed curve) and Q obt and Q rel for the state transition a → b .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00499v1-figure-cea1bd6f-4c0f-49d0-90d2-065846b92034.jpg)

In this Letter, we develop a theory to describe ET in complex networks of donor-acceptor sites, where each site is associated with normal modes that are in contact with an independent heat bath at the local site temperature. This is an idealization of the standard phenomenology of molecular electron transfer where the electronic processes is most strongly coupled to vibrations that are localized on the donor and acceptor sites. In the limit of strong electron-phonon coupling, electron transport is dominated by hopping-type events and Marcus theory gives the conceptual basis, and sometimes quantitative understanding, of ET reactions [13-16]. The multithermal nature of the examined systems arises because the donor and acceptor sites can be at different temperatures. Purely vibrational heat transfer between sites

is disregarded for simplicity, so in what follows we focus on electron transfer and the associated heat transfer [19, 34]. The development of complex network theory [35] has significantly increased our understanding of the flow of charge, energy, and information in diverse types of systems [36-45] and the present study makes it possible to consider electron and heat transport within such a framework and to study the consequence of their interdependence.

The donor-acceptor networks we consider consist of S sites, where each site s , which has local temperature T s , is associated with N s harmonic modes that are equilibrated with the thermal environment about that site. Specifically, we consider transitions between electronic states for which an excess electron is localized on different sites of this network. In the Marcus picture [13] of ET, the localized occupation of electron density associated with electronic state a of the network is described by the energy surface

$$E _ { a } ( x _ { 1 }, \dots, x _ { N } ) = E _ { a } ^ { ( 0 ) } + \sum _ { s } ^ { S } \sum _ { j \in \mathcal { M } ^ { ( s ) } } \frac { 1 } { 2 } k _ { j } \left ( x _ { j } - \lambda _ { j } ^ { ( a ) } \right ) ^ { 2 }, \quad \begin{array} { c } P _ { a, b } \mathfrak { i k } \\ \text{with} \\ \text{surfac} \\ \text{to the} \\ \text{where} \end{array}$$

where E (0) a is an electronic energy origin of state a , x j is the coordinate of the j th vibrational mode, and λ ( a ) j is a shift in the equilibrium position of the j th mode. Both E (0) a and λ ( a ) j are measured relative to some reference state for which E (0) ref and λ (ref) j vanish (properties of this state do not affect the final results). The total number of modes in the system is denoted by N = ∑ S s N s and the group of modes associated with the s th site is M ( s ) . In a likely special case λ ( a ) j = 0 unless mode j is localized about the site on which the electron is placed in state a (that is, unless j ∈ M ( s ) and a is the electronic state that corresponds to the electron occupying site s ). However, the form (1) can represent more general situations where the modes localized about site s respond to the electronic occupation on a different site. A transition between states a and b is associated with a reorganization energy (assumed temperature independent) in the j th mode, E ( a,b ) R j , and a total reorganization energy, E ( a,b ) R , which are given by:

$$E _ { \text{R} j } ^ { ( a, b ) } = \frac { 1 } { 2 } k _ { j } \left ( \lambda _ { j } ^ { ( a ) } - \lambda _ { j } ^ { ( b ) } \right ) ^ { 2 } \text{ \ and \ } E _ { \text{R} } ^ { ( a, b ) } = \sum _ { j } ^ { N } E _ { \text{R} j } ^ { ( a, b ) }. \\ \text{The shift $\lambda$ and home to r o r a n i z a t i o n o n o r o r o }$$

The shifts λ , and hence the reorganization energies, of a system depend on the intersite distances and therefore the topology of the underlying connectivity network. Shown in Fig. 1 are representative topologies for three typical connectivities: ring R , linear L , and complete K . Associated with each network is an adjacency matrix A where A ab = 1 if sites a and b are connected (and thus the electron can tunnel between sites) and is zero oth- erwise [35]. The topology of A determines which modes are responsive to electron localization on a particular site. We next show that temperature differences between sites give rise to emergent and sometimes exotic thermal and electronic transport properties.

The electron transfer rate between any two sites in the network and the heat transfer rate between the thermal baths involved in this transition can be derived by adopting the formalism put forth in Ref. 19 for bithermal electron hopping between two sites. 'Involved bath' implies that the harmonic modes that are thermalized by this bath are sensitive to the electronic population of at least one of the sites. Under standard transition state theory assumptions the rate of the a → b ≡ a, b transition can be expressed as

$$k _ { a, b } = \frac { 1 } { 2 } \Big < \mathcal { T } _ { a, b } \, \dot { x } _ { \perp } \Big > P _ { a, b }, \quad \quad ( 3 ) \\ \intertext { is the tunnelino} \ n r o b a b i l i t v \, \text{between states}.$$

where T a,b is the tunneling probability between states, P a,b is the probability density about a transition surface (TS) separating the states, both evaluated on the E a surface, and ˙ x ⊥ is the velocity in the direction normal to the TS [46, 47]. The normal velocity ˙ x ⊥ = ˙ x · ˆ u ⊥ ( x ) where ˆ u ⊥ ( x ) is a unit vector normal to the TS at position x on this hypersurface. The factor 〈T a,b ˙ x ⊥ 〉 is a (multi)thermal average that depends on the temperature of each bath involved in the transition. The TS is determined by the requirement that a transition can only take place at nuclear configurations where electronic energy is conserved. For the N -dimensional paraboloid energy surface defined by (1), the TS separating states is an ( N -1)dimensional transition state hypersurface which is the locus of mode configurations where E a = E b , defined by g c ( x , . . . , x 1 N ) = E b ( x , . . . , x 1 N ) -E a ( x , . . . , x 1 N ).

In the adiabatic limit of Marcus rate theory T a,b = 1 and the pre-exponential factor is proportional to 〈 ˙ x ⊥ 〉 [48]. In the nonadiabatic limit [49], T a,b can be approximated by the corresponding limit of the Landau-Zener expression [15, 50, 51] which is evaluated in the direction normal to the TS [52]. In this limit, T a,b ∝ 1 /x ˙ ⊥ , and the expectation value 〈T a,b ˙ x ⊥ 〉 does not depend on the normal velocity [53]. For the a → b transition, the multithermal probability to be on the TS is:

$$\L ^ { \quad } \quad & \text{multithermal probability to be on the $1$} \, \text{is} \, \text{:} \\ \Phi _ { b } \, \Phi _ { \quad } \quad & P _ { a, b } = \int _ { \mathbb { R } ^ { N } } | \nabla g _ { c } | \, \delta \left ( g _ { c } ( x _ { 1 }, \dots, x _ { N } ) \right ) \\ \left ( 2 \right ) \quad & \times \prod _ { s } ^ { S } \prod _ { j \in \mathcal { M } ^ { ( s ) } } \exp \left [ - \beta _ { s } \frac { k _ { j } } { 2 } \left ( x _ { j } - \lambda _ { j } ^ { ( a ) } \right ) ^ { 2 } \right ] d x _ { j } \quad \left ( 4 \right ) \\ \text{ere-} \\ \text{ork.} \quad & \int _ { \mathbb { R } ^ { N } } \prod _ { s } ^ { S } \prod _ { j \in \mathcal { M } ^ { ( s ) } } \exp \left [ - \beta _ { s } \frac { k _ { j } } { 2 } \left ( x _ { j } - \lambda _ { j } ^ { ( a ) } \right ) ^ { 2 } \right ] d x _ { j }, \\ \L ^ { \quad } \, \mathcal { K }. \quad & \text{where} \, \beta _ { s } = 1 / k _ { B } T _ { s } \text{ with } k _ { B } \text{ being the Boltzmann con-}$$

where β s = 1 /k B T s with k B being the Boltzmann constant. The δ -function in (4) constrains the integration over the vibrational subspace in which E a = E b and the

gradient magnitude |∇ g c | = √ ∑ N j 2 k E j ( a,b ) R j gives a precise definition to this constraint [54]. Evaluating the integrals in Eq. (4) we obtain

$$P _ { a, b } & = \sqrt { \sum _ { j } ^ { N } k _ { j } E _ { \text{R} j } ^ { ( a, b ) } \Big / 2 \pi k _ { B } \sum _ { s } ^ { S } T _ { s } \sum _ { j \in \mathcal { M } ^ { ( s ) } } E _ { \text{R} j } ^ { ( a, b ) } } \\ & \times \exp \left [ - \left ( E _ { b a } + E _ { \text{R} } ^ { ( a, b ) } \right ) ^ { 2 } \Big / 4 k _ { B } \sum _ { s } ^ { S } T _ { s } \sum _ { j \in \mathcal { M } ^ { ( s ) } } E _ { \text{R} j } ^ { ( a, b ) } \right ], \quad \text{fig.} \\ \text{which expresses the probability density on the reach} \quad \text{syster}$$

which expresses the probability density on the reaction path in terms of the temperature of each bath, and the reorganization energy and reaction free energy E ba = -E ab = E (0) b -E (0) a of the respective transition. In the uniform temperature limit, combining Eqs. (3) and (5) recovers the Marcus rate [13-16].

With the multidimensional-multithermal transition rate k a,b now derived, the kinetic equations for the occupation probability P of each state can be constructed. For a reaction network (see Fig. 1) with adjacency matrix A these equations take the form

$$\frac { d \mathcal { P } _ { a } } { d t } = \sum _ { b } \Lambda _ { b a } k _ { b, a } \mathcal { P } _ { b } ( t ) - \Lambda _ { a b } k _ { a, b } \mathcal { P } _ { a } ( t ), \quad ( 6 ) \quad \Big \langle \zeta$$

for each state a . At steady-state d P a /dt = 0 ∀ a , which implies that the net electronic flux between sites vanishes.

Next consider the heat transfer associated with a given electron transfer step [19, 55]. In Marcus theory, the nuclear motion leading to the a → b transition proceeds through a point x TS on the TS, and the corresponding heat transferred into a specific bath during the transition is Q ( a,b ) ( x TS ) = -Q ( a ) obt ( x TS )+ Q ( b ) rel ( x TS ) , where the first term is the heat obtained from the bath during the ascent to x TS on the E a surface and the second term is the heat released to the bath during the descent to equilibrium on the E b surface (see Fig. 1). Using the energy surfaces of Eq. (1), the contribution of mode j to the net heat exchange with the bath associated with it during the a → b transition is given by

$$\mathcal { Q } _ { j } ^ { ( a, b ) } = - \frac { 1 } { 2 } k _ { j } \left [ x _ { j } ^ { \text{TS} } - \lambda _ { j } ^ { ( a ) } \right ] ^ { 2 } + \frac { 1 } { 2 } k _ { j } \left [ x _ { j } ^ { \text{TS} } - \lambda _ { j } ^ { ( b ) } \right ] ^ { 2 }. \ ( 7 ) \quad \text{dropo} \\ \dots \. \. \ \dots \. \ \. \ \dots \ \ \text{Th}$$

This quantity should be averaged over the probability to reach the particular configuration on the TS when coming from the a side:

$$\left \langle \mathcal { Q } _ { j } ^ { ( a, b ) } \right \rangle = \int _ { \mathbb { R } ^ { N } } \mathcal { Q } _ { j } ^ { ( a, b ) } \, P _ { a, b } ^ { \ddagger } ( x _ { 1 }, \dots, x _ { N } ) \, d x _ { 1 } \dots d x _ { N }, \ ( 8 ) \quad \text{Note} _ { \mathcal { Q } _ { s } ^ { ( a, b ) } }$$

where P ‡ a,b is the probability density on the TS when the

FIG. 2. Network graph for a three-site ring ( R 3 ) network.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00499v1-figure-70e4e1d7-de97-41bf-bb1f-7d4b3fe73976.jpg)

system is in state a [48]:

$$\text{rec} _ { \text{$\kappa$} } \quad \kappa \cdots \quad \kappa \cdots \quad \kappa \, \\ \text{$\text{$\text{$\text{$array}$} } } _ { \text{$\text{$\text{$\text{$array}$} } } P _ { a, b } ^ { \ddagger } = \delta \left ( g _ { c } \right ) \prod _ { s } ^ { S } \prod _ { j \in \mathcal { M } ^ { ( s ) } } \exp \left [ - \beta _ { s } \frac { 1 } { 2 } k _ { j } \left ( x _ { j } - \lambda _ { j } ^ { ( a ) } \right ) ^ { 2 } \right ] \\ \text{$\text{$\text{$array}$} } \quad \text{$\int} _ { \mathbb { R } ^ { N } } \delta \left ( g _ { c } \right ) \prod _ { s } ^ { S } \prod _ { j \in \mathcal { M } ^ { ( s ) } } \exp \left [ - \beta _ { s } \frac { 1 } { 2 } k _ { j } \left ( x _ { j } - \lambda _ { j } ^ { ( a ) } \right ) ^ { 2 } \right ] d x _ { j }. \\ \text{$\text{$\text{$occ}$} } \quad \text{$\text{$\text{$totalting the expectation integral gives} } }$$

Evaluating the expectation integral gives

/negationslash

$$\begin{array} { c } \colon & - \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ \text{rix} \\ ( 6 ) & \left \langle \mathcal { Q } _ { j } ^ { ( a, b ) } \right \rangle = \frac { E _ { \text{R} j } ^ { ( a, b ) } \left [ E _ { a b } T _ { j } + \sum _ { k \neq j } ^ { N } E _ { \text{R} k } ^ { ( a, b ) } ( T _ { k } - T _ { j } ) \right ] } { \sum _ { k } ^ { N } T _ { k } E _ { \text{R} k } ^ { ( a, b ) } } & ( 1 0 ) \\ \text{ich} \\ \text{one} & \text{ where } T _ { i } = T _ { s } \text{ if } j \in \mathcal { M } ^ { ( s ) }. \text{ The total heat transferred} \end{array}$$

where T j = T s if j ∈ M ( s ) . The total heat transferred to the thermal environment of site s during the a → b transition is the sum of contributions over all modes associated with that site:

/negationslash

$$\text{sads} \\ \text{len} \\ \text{ext} \\ \text{len} \\ \text{eq} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \text{s} \\ \end{cases} ( 1 1 ) \text{$2$} \text{$3$} \sum _ { q \neq s } E _ { R j } ^ { ( a, b ) } \left ( T _ { q } - T _ { s } \right ) \right ] \\ \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{s} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text` \\ \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text{$2$} \text`$$

Here, the term proportional to E ab expresses the heat released to/taken from the baths from the free energy difference between these electronic states, while the term proportional to ( T q -T s ) is the actual heat transfer between baths q and s associated with the a → b transition.

The heat current into the thermal environment of site s is

$$\frac { d \mathcal { Q } _ { s } } { d t } = \sum _ { a, b } \mathbb { A } _ { a b } k _ { a, b } \mathcal { P } _ { a } ( t ) \left \langle \mathcal { Q } _ { s } ^ { ( a, b ) } \right \rangle. \quad \ \ ( 1 2 )$$

Note that while P a ( ) t ≥ 0 ∀ t , the expectation value of Q ( a,b ) s and the heat current ˙ Q s can be either positive or negative implying that the corresponding bath absorbs

FIG. 3. Cyclic flux J c (top) and heat current d Q 1 /dt (bottom) in the nonadiabatic limit as functions of T 3 for a threesite R 3 network. Each curve corresponds to the respective value of E (0) 3 marked in the legend. In each panel: E (0) 1 = -4, E (0) 2 = 0, and T 2 = 3 2 with (a)-(b) / T 1 = 3/2 and (c)-(d) T 1 = 3 4. / The circular markers denote unithermal points. All reorganization energies are E R j = 1 2 for each mode in-/ volved in a particular a transition and zero otherwise. All quantities are shown in reduced units [53, 56].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00499v1-figure-d673f8cd-84e8-4a47-8b79-2aba548a9d66.jpg)

/negationslash or releases energy, respectively. In unithermal systems at steady-state the heat currents vanish ( ˙ Q s = 0 ∀ s ), while in contrast, for multithermal systems, 〈Q 〉 /negationslash s = 0 and thus ˙ Q s = 0. Note that even as the occupation probabilities approach electronic quasi-equilibrium where the net electronic currents are zero, the net flow of heat does not vanish. This phenomenon is not a standard thermoelectric effect and reveals a novel pathway for energy transport in multithermal charge transfer networks.

Another interesting behavior that is observed in multithermal networks with closed loops is the persistence of steady-state net electronic bond currents, J a,b -J b,a = k a,b P a -k b,a P b , i.e., the breaking of detailed balance, and the formation of cyclical current loops that obviously vanish in full equilibrium where detailed balance is maintained. An example is seen in the three-site ring ( R 3 ) shown in Fig. 2. We emphasize this simple system because of its experimental realizability and its general applications in molecular electronics and devices [33, 57, 58], but note that other more complex networks can also be analyzed using the developed theory [59]. In a multithermal R 3 network, the direction and magnitude of the cyclic flux J c = J 1 2 , -J 2 1 , = J 2 3 , -J 3 2 , = J 3 1 , -J 1 3 , can be altered by tuning the temperature T 3 . This is illustrated in Fig. 3(a) where the donor and acceptor sites 1,2 are at the same temperature while variation of the

FIG. 4. (a) Amplification factor α and (b) heat current derivatives ∂ ˙ Q s /∂T 3 as functions of T 3 for a multithermal R 3 system with T 1 = 3 4, / T 2 = 3 2, and / E R j = 1 2 / ∀ j over each state transition. In both panels E (0) 1 = -4 and E (0) 2 = E (0) 3 = 0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00499v1-figure-d34ad80a-b1fc-4fb1-903c-9a1c54d3b1f3.jpg)

temperature on another site, T 3 , determines the direction and magnitude of the cyclical current. We find that this thermally-induced current persists, except in the case that the electronic occupation energy E (0) a is the same for all sites involved in the cycle. Note that at the unithermal point in Figs. 3(a)-(b), J a,b = J b,a for every connection, so J c = 0. This is simply a statement that with no temperature gradient there is no heat current or cyclic electron flux. Similar trends are also observed in Figs. 3(c)-(d) for networks where the local temperature of each site is different. By comparing the temperatures at which J c = 0 in Fig. 3(c) with the heat currents at the same temperature in Fig. 3(d) it can be seen that even when the electron flux vanishes there is still a net heat flow between baths. At full equilibrium, which is achieved in the uniform temperature limit, all electronic and heat fluxes vanish.

The theory of coupled electron and heat transfer developed in this Letter can also be applied to elucidate the electronic contribution to transport phenomena in thermal transistors [31-33], which are recently studied model devices that can be used to control and amplify heat flow. Following Ref. 31, we quantify the magnitude of amplification in thermal transistors by the factor α s = ∂ ˙ Q s /∂ ˙ Q 3 : s ∈ { 1 2 , } , which measures the effect of pumping heat into site 3 on the heat current between sites 1 and 2. If | α s | &gt; 1 a transistor effect is present. In the three-site R 3 system, variation of the heat current ˙ Q 3 by alteration of T 3 can give rise to significant ampli-

fications, as shown in Fig. 4(a). The reason for the electrothermal transistor effect is shown in Fig. 4(b), where ∂ ˙ Q 3 /∂T 3 → 0 at T 3 ≈ 0 4 while . | ∂ ˙ Q s /∂T 3 | /greatermuch ∂ ˙ Q 3 /∂T 3 for s ∈ { 1 2 , } as this limit is approached. This aligns with the region of amplification shown in Fig. 4(a).

We have shown that in nanoscale systems where localized modes are in contact with environments at different temperatures, complex nonlinearities in the thermal gradient can induce currents which are characterized by multiple temperatures. A theory has been developed to unify the description of heat and charge transfer in these systems with multithermal temperature gradients between donor-acceptor states. This work provides a bridge connecting theories of electron transfer, heat transport, and thermoelectricity in systems where electron transport is dominated by intersite hopping, and will be useful in the design of electronic and thermoelectric devices that operate in this limit.

In regimes where the magnitude of heat conduction due to electron transport dominates over the magnitude of conduction from phonons, the developed theory will be directly applicable. Otherwise, a complete picture of the conduction process will require a theory that considers both electrothermal and phononic heat transport, and their coupling. Examination of thermopower, efficiency, and their relation to electrothermal transport in molecular junctions (and other complex donor/acceptor networks with molecule-metal and molecule-semiconductor interfaces) in the phonon-assisted hopping limit of electronic transport is a potential application of the present theory [55, 60, 61]. Generalizations that go beyond the semiclassical Marcus treatment are obviously needed and will be taken on in future work.

AN is supported by the Israel Science Foundation, the US-Israel Binational Science Foundation and the University of Pennsylvania.

- [1] M. Galperin, M. A. Ratner, and A. Nitzan, J. Phys.: Condens. Matter 19 , 103201 (2007).
- [2] Y. Dubi and M. Di Ventra, Rev. Mod. Phys. 83 , 131 (2011), doi:10.1103/RevModPhys.83.131.
- [3] K. Walczak, Physica B 392 , 173 (2007), ISSN 0921-4526, doi:10.1016/j.physb.2006.11.013.
- [4] M. Galperin, A. Nitzan, and M. A. Ratner, Mol. Phys. 106 , 397 (2008), doi:10.1080/00268970701837784.
- [5] Y. Wang, J. Zhou, and R. Yang, J. Phys. Chem. C 115 , 24418 (2011), doi:10.1021/jp208490q.
- [6] J. Ren, J.-X. Zhu, J. E. Gubernatis, C. Wang, and B. Li, Phys. Rev. B 85 , 155443 (2012), doi:10.1103/PhysRevB.85.155443.
- [7] M. B. Tagani and H. R. Soleimani, J. Appl. Phys. 112 , 103719 (2012), doi:10.1063/1.4767376.
- [8] M. B. Tagani and H. R. Soleimani, Physica B 413 , 86 (2013), doi:10.1016/j.physb.2013.01.007.
- [9] C. A. Perroni, D. Ninno, and V. Cataudella, Phys. Rev.
- B 90 , 125421 (2014), doi:10.1103/PhysRevB.90.125421.
- [10] C. A. Perroni, D. Ninno, and V. Cataudella, New J. Phys. 17 , 083050 (2015), http://iopscience.iop.org/article/10.1088/13672630/17/8/083050.
- [11] T. Koch, J. Loos, and H. Fehske, Phys. Rev. B 89 , 155133 (2014), doi:10.1103/PhysRevB.89.155133.
- [12] N. A. Zimbovskaya, J. Phys.: Condens. Matter 28 , 183002 (2016), http://stacks.iop.org/09538984/28/i=18/a=183002.
- [13] R. A. Marcus, J. Chem. Phys. 24 , 966 (1956), doi:10.1063/1.1742723.
- [14] R. A. Marcus, Annu. Rev. Phys. Chem. 15 , 155 (1964), doi:10.1146/annurev.pc.15.100164.001103.
- [15] A. Nitzan, Chemical Dynamics in Condensed Phases: Relaxation, Transfer, and Reactions in Condensed Molecular Systems (Oxford University Press, 2006).
- [16] B. Peters, J. Phys. Chem. B 119 , 6349 (2015), doi:10.1021/acs.jpcb.5b02547.
- [17] J. Zhang, A. M. Kuznetsov, I. G. Medvedev, Q. Chi, T. Albrecht, P. S. Jensen, and J. Ulstrup, Chem. Rev. 108 , 2737 (2008), doi:10.1021/cr068073+.
- [18] A. Migliore and A. Nitzan, J. Am. Chem. Soc. 135 , 9420 (2013), doi:10.1021/ja401336u.
- [19] G. T. Craven and A. Nitzan, Proc. Natl. Acad. Sci. 113 , 9421 (2016), doi:10.1073/pnas.1609141113.
- [20] C. W. Chang, D. Okawa, A. Majumdar, and A. Zettl, Science 314 , 1121 (2006), doi:10.1126/science.1132898.
- [21] S. Sadat, A. Tan, Y. J. Chua, and P. Reddy, Nano Lett. 10 , 2613 (2010), doi:10.1021/nl101354e.
- [22] J. A. Malen, S. K. Yee, A. Majumdar, and R. A. Segalman, Chem. Phys. Lett. 491 , 109 (2010), doi:10.1016/j.cplett.2010.03.028.
- [23] A. Tan, J. Balachandran, S. Sadat, V. Gavini, B. D. Dunietz, S.-Y. Jang, and P. Reddy, J. Am. Chem. Soc. 133 , 8838 (2011), doi:10.1021/ja202178k.
- [24] W. Lee, K. Kim, W. Jeong, L. A. Zotti, F. Pauly, J. C. Cuevas, and P. Reddy, Nature 498 , 209 (2013), doi:10.1038/nature12183.
- [25] Y. Kim, W. Jeong, K. Kim, W. Lee, and P. Reddy, Nature Nanotech. 9 , 881 (2014), doi:10.1038/nnano.2014.209.
- [26] A. V. Popov and R. Hernandez, J. Chem. Phys. 126 , 244506 (2007), doi:10.1063/1.2743032.
- [27] C. Jarzynski, Nature Phys. 11 , 105 (2015), doi:10.1038/nphys3229.
- [28] G. T. Craven and R. Hernandez, Phys. Rev. Lett. 115 , 148301 (2015), doi:10.1103/PhysRevLett.115.148301.
- [29] Q. Li, I. Duchemin, S. Xiong, G. C. Solomon, and D. Donadio, J. Phys. Chem. C 119 , 24636 (2015), doi:10.1021/acs.jpcc.5b07429.
- [30] M. Einax and A. Nitzan, J. Chem. Phys. 145 , 014108 (2016), doi:10.1063/1.4955160.
- [31] B. Li, L. Wang, and G. Casati, Appl. Phys. Lett. 88 , 143501 (2006), doi:10.1063/1.2191730.
- [32] P. Ben-Abdallah and S.-A. Biehs, Phys. Rev. Lett. 112 , 044301 (2014), doi:10.1103/PhysRevLett.112.044301.
- [33] K. Joulain, J. Drevillon, Y. Ezzahri, and J. OrdonezMiranda, Phys. Rev. Lett. 116 , 200601 (2016), doi:10.1103/PhysRevLett.116.200601.
- [34] D. V. Matyushov, Proc. Natl. Acad. Sci. 113 , 9401 (2016), doi:10.1073/pnas.1610542113.
- [35] J. Schnakenberg, Rev. Mod. Phys. 48 , 571 (1976), doi:10.1103/RevModPhys.48.571.

- [36] S. H. Strogatz, Nature 410 , 268 (2001), doi:10.1038/35065725.
- [37] R. Albert and A.-L. Barab´si, Rev. Mod. Phys. a 74 , 47 (2002), doi:10.1103/RevModPhys.74.47.
- [38] E. Estrada and N. Hatano, Chem. Phys. Lett. 439 , 247 (2007), doi:10.1016/j.cplett.2007.03.098.
- [39] J. Gao, S. V. Buldyrev, S. Havlin, and H. E. Stanley, Phys. Rev. Lett. 107 , 195701 (2011), doi:10.1103/PhysRevLett.107.195701.
- [40] G. C. Solomon, C. Herrmann, T. Hansen, V. Mujica, and M. A. Ratner, Nature Chem. 2 , 223 (2010), doi:10/1038/nchem.546.
- [41] D. Rai, O. Hod, and A. Nitzan, J. Phys. Chem. Lett. 2 , 2118 (2011), doi:10.1021/jz200862r.
- [42] D. Rai, O. Hod, and A. Nitzan, Phys. Rev. B 85 , 155440 (2012), doi:10.1103/PhysRevB.85.155440.
- [43] B. M. Savoie, K. L. Kohlstedt, N. E. Jackson, L. X. Chen, M. Olvera de la Cruz, G. C. Schatz, T. J. Marks, and M. A. Ratner, Proc. Natl. Acad. Sci. 111 , 10055 (2014), doi:10.1073/pnas.1409514111.
- [44] M. Einax and A. Nitzan, J. Phys. Chem. C 118 , 27226 (2014), doi:10.1021/jp5084373.
- [45] T. Hansen and G. C. Solomon, J. Phys. Chem. C 120 , 6295 (2016), doi:10.1021/acs.jpcc.5b11211.
- [46] S. Hammes-Schiffer and J. C. Tully, J. Chem. Phys. 103 , 8528 (1995), doi:10.1063/1.470162.
- [47] G. H. J´hannesson and H. J´nsson, J. Chem. Phys. o o 115 , 9644 (2001), doi:10.1063/1.1415499.
- [48] See the Supplemental Material for detailed derivations of each component in the rate and heat transfer expressions.
- [49] For experimental setups in which sites of appreciably different temperatures are relatively far from each other, implying relatively weak electronic coupling, the nonadiabatic limit will be the most pertinent.
- [50] D. L. Thompson, Modern Methods for Multidimensional
- Dynamics Computations in Chemistry (World Scientific, 1998).
- [51] A. O. Lykhin, D. S. Kaliakin, G. E. dePolo, A. A. Kuzubov, and S. A. Varganov, Int. J. Quantum Chem. 116 , 750 (2016), doi:10.1002/qua.25124.
- [52] A. J. Marks and D. L. Thompson, J. Chem. Phys. 96 , 1911 (1992), doi:10.1063/1.462092.
- [53] In the nonadiabatic limit a coupling constant between diabatic surfaces: H a,b = 0 01 eV has been used in all . calculations.
- [54] C. Hartmann, J. C. Latorre, and G. Ciccotti, Eur. Phys. J. Spec. Top. 200 , 73 (2011), doi:10.1140/epjst/e201101519-7.
- [55] G. T. Craven and A. Nitzan, J. Chem. Phys. 146 , 092305 (2017), doi:10.1063/1.4971293.
- [56] Characteristic units are: /epsilon1 = 0 02 eV for . energy, τ = 0 5 ps . for time, and σ = 1nm for length. Using these units, a reduced temperature of k B T//epsilon1 ≈ 1 3 corresponds . to T = 300 K.
- [57] N. Li, J. Ren, L. Wang, G. Zhang, P. H¨ anggi, and B. Li, Rev. Mod. Phys. 84 , 1045 (2012), doi:10.1103/RevModPhys.84.1045.
- [58] A. E. Allahverdyan, S. G. Babajanyan, N. H. Martirosyan, and A. V. Melkikh, Phys. Rev. Lett. 117 , 030601 (2016), doi:10.1103/PhysRevLett.117.030601.
- [59] Analysis of electric and heat currents for other multithermal networks is given in the Supplemental Material.
- [60] C. Jia, A. Migliore, N. Xin, S. Huang, J. Wang, Q. Yang, S. Wang, H. Chen, D. Wang, B. Feng, et al., Science 352 , 1443 (2016), doi:10.1126/science.aaf6298.
- [61] A. G. Gagorik, B. Savoie, N. Jackson, A. Agrawal, A. Choudhary, M. A. Ratner, G. C. Schatz, and K. L. Kohlstedt, J. Phys. Chem. Lett. 8 , 415 (2017), doi:10.1021/acs.jpclett.6b02921.
Galen T. Craven Abraham Nitzan
Electrothermal Transistor Effect and Cyclic Electronic Currents in Multithermal Charge Transfer Networks
A theory is developed to describe the coupled transport of energy and charge
in networks of electron donor-acceptor sites which are seated in a thermally
heterogeneous environment, where the transfer kinetics are dominated by
Marcus-type hopping rates. It is found that the coupling of heat and charge
transfer in such systems gives rise to exotic transport phenomena which are
absent in thermally homogeneous systems and cannot be described by standard
thermoelectric relations. Specifically, the directionality and extent of
thermal transistor amplification and cyclical electronic currents in a given
network can be controlled by tuning the underlying temperature gradient in the
system. The application of these findings toward optimal control of
multithermal currents is illustrated on a paradigmatic nanostructure.