Rectification of heat transfer in nanodevices can be realized by combining the system inherent anharmonicity with structural asymmetry. We analyze this phenomenon within the simplest anharmonic system-a spin-boson nanojunction model. We consider two variants of the model that yield, for the first time, analytical solutions: a linear separable model in which the heat reservoirs contribute additively, and a nonseparable model suitable for a stronger system-bath interaction. Both models show asymmetric (rectifying) heat conduction when the couplings to the heat reservoirs are different.

We consider phononic heat transport through molecular chains connecting two thermal reservoirs. For relatively short molecules at normal temperatures we find, using classical stochastic simulations, that heat conduction is dominated by the harmonic part of the molecular force-field. We develop a general theory for the heat conduction through harmonic chains in 3-dimensions. Our approach uses the standard formalism that leads to the generalized Langevin equation for a system coupled to a harmonic heat bath, however the driving and relaxation terms are considered separately in a way that leads directly to the steady state response and the heat current under non-equilibrium driving. A Landauer-type expression for the heat conduction is obtained, in agreement with other recent studies. We used this general formalism to study the heat conduction properties of alkane. We find that for relatively short (1-30 carbon molecules) the length and temperature dependence of the molecular heat conduction results from the balance of three factors: (i) The molecular frequency spectrum in relation to the frequency cutoff of the thermal reservoirs, (ii) the degree of localization of the molecular normal modes and (iii) the molecule-heat reservoirs coupling. The fact that molecular modes at different frequency regimes have different localization properties gives rise to intricate dependence of the heat conduction on molecular length at different temperature. For example, the heat conduction increases with molecular length for short molecular chains at low temperatures. Isotopically substituted disordered chains are also studied and their behavior can be traced to the above factors together with the increased mode localization in disordered chain and the increase in the density of low frequency modes associated with heavier mass substitution. Finally, we compare the heat conduction obtained from this microscopic calculation to that estimated by considering the molecule as a cylinder characterized by a macroscopic heat conduction typical to 2 organic solids. We find that this classical model overestimates the heat conduction of single alkane molecules by about an order of magnitude at room temperature.Implications of the present study to the problem of heating in electrically conducting molecular junctions are pointed out.3

The effect of dephasing and relaxation on electron transfer in bridged molecular systems is investigated using a simple molecular model. The interaction between the molecular system and the thermal environment is described on the level of the Redfield theory, modified when needed for the description of steady-state situations. Noting that transient as well as steady-state measurements are possible in such system, we discuss the relationship between the rates obtained from these different types of experiments and, in particular, the conditions under which these rates are the same. Also, a formal relation between the steady-state rate for electron transfer across a molecular bridge and the conductance of this bridge when placed between two metal contacts is established. The effect of dephasing and relaxation on the electron transfer is investigated, and new observations are made with regard to the transition from the superexchange to the thermal (hopping through bridge) regime of the transfer process. In particular, the rate is temperature-independent in the superexchange regime, and its dependence on the bridge length (N) is exponential, exp(-N). The rate behaves like (R 1 + R 2 N) -1 exp(-∆E/k B T) beyond a crossover value of N, where ∆E is the energy gap between the donor/acceptor and the bridge levels, and where R 1 and R 2 are characteristic times for activation onto the bridge and diffusion in the bridge, respectively. We find that, in typical cases, R 1 . R 2 , and therefore, a region of very weak N dependence is expected before the Ohmic behavior, N -1 , is established for large enough N. In addition, a relatively weak exponential dependence, exp(-RN), is expected for long bridges if competing processes capture electrons away from the bridge sites. Finally, we consider ways to distinguish experimentally between the thermal and the tunneling routes.

Photothermal catalysis is an emerging sub-discipline of heterogeneous catalysis that exploits broad absorption of the solar spectrum to stimulate a combination of thermochemical and photochemical processes, which contribute synergistically to driving catalytic reactions. In particular, it is proving an effective and promising strategy for converting CO2 to synthetic fuels.

We investigate the heat conduction properties of molecular junctions comprising anharmonic interactions. We find that nonlinear interactions can lead to novel phenomena: negative differential thermal conductance and heat rectification. Based on analytically solvable models we derive an expression for the heat current that clearly reflects the interplay between anharmonic interactions, strengths of coupling to the thermal reservoirs, and junction asymmetry. This expression indicates that negative differential thermal conductance shows up when the molecule is strongly coupled to the thermal baths, even in the absence of internal molecular nonlinearities. In contrast, diode like behavior is expected for a highly anharmonic molecule with an inherent structural asymmetry.Anharmonic interactions are also necessary for manifesting Fourier type transport. We briefly present an extension of our model system that can lead to this behavior.

We develop an iterative, numerically exact approach for the treatment of nonequilibrium quantum transport and dissipation problems that avoids the real-time sign problem associated with standard Monte Carlo techniques. The method requires a well-defined decorrelation time of the non-local influence functional for proper convergence to the exact limit. Since finite decorrelation times may arise either from temperature or from a voltage drop at zero temperature, the approach is well suited for the description of the real-time dynamics of single-molecule devices and quantum dots driven to a steady-state via interaction with two or more electron leads. We numerically investigate two non-trivial models: the evolution of the nonequilibrium population of a two-level system coupled to two electronic reservoirs, and quantum transport in the nonequilibrium Anderson model. For the latter case, two distinct formulations are described. Results are compared to those obtained by other techniques.

We focus on the non-equilibrium two-bath spin-boson model, a toy model for examining quantum thermal transport in many-body open systems. Describing the dynamics within the noninteracting-blip approximation equations, applicable, e.g., in the strong system-bath coupling limit and/or at high temperatures, we derive expressions for the cumulant generating function in both the Markovian and non-Markovian limits by energy-resolving the quantum master equation of the subsystem. For a Markovian bath, we readily demonstrate the validity of a steady-state heat exchange fluctuation theorem. In the non-Markovian limit a "weaker" symmetry relation generally holds, a general outcome of microreversibility. We discuss the reduction of this symmetry relation to the universal steady-state fluctuation theorem. Using the cumulant generating function, an analytic expression for the heat current is obtained. Our results establish the validity of the steady-state heat exchange fluctuation theorem in quantum systems with strong system-bath interactions. From the practical point of view, this study provides tools for exploring transport characteristics of the two-bath spin-boson model, a prototype for a nonlinear thermal conductor.

We propose a novel molecular device that pumps heat against a thermal gradient. The system consists of a molecular element connecting two thermal reservoirs that are characterized by different spectral properties. The pumping action is achieved by applying an external force that periodically modulates molecular levels. This modulation affects periodic oscillations of the internal temperature of the molecule and the strength of its coupling to each reservoir resulting in a net heat flow in the desired direction. The heat flow is examined in the slow and fast modulation limits and for different modulation waveforms, thus making it possible to optimize the device performance.