Phononic heat transport in molecular junctions: Quantum effects and vibrational mismatch

Challenges in molecular dynamics simulations of heat exchange statistics

We study heat exchange in temperature-biased metal-molecule-metal molecular junctions by employing the molecular dynamics simulator LAMMPS. Generating the nonequilibrium steady state with Langevin thermostats at the boundaries of the junction, we show that the average heat current across a gold-alkanedithiol-gold nanojunction behaves physically, with the thermal conductance value matching the literature. In contrast, the full probability distribution function for heat exchange, as generated by the simulator, violates the fundamental fluctuation symmetry for entropy production. We trace this failure back to the implementation of the thermostats and the expression used to calculate the heat exchange. To rectify this issue and produce the correct statistics, we introduce single-atom thermostats as an alternative to conventional many-atom thermostats. Once averaging heat exchange over the hot and cold thermostats, this approach successfully generates the correct probability distribution function, which we use to study the behavior of both the average heat current and its noise. We further examine the thermodynamic uncertainty relation in the molecular junction and show that it holds, albeit demonstrating nontrivial trends. Our study points to the need to carefully implement nonequilibrium molecular dynamics solvers in atomistic simulation software tools for future investigations of noise phenomena in thermal transport.

Simulations of heat transport in single-molecule junctions: Investigations of the thermal diode effect

With the objective of understanding microscopic principles governing thermal energy flow in nanojunctions, we study phononic heat transport through metal-molecule-metal junctions using classical molecular dynamics (MD) simulations. Considering a single-molecule gold-alkanedithiol-gold junction, we first focus on aspects of method development and compare two techniques for calculating thermal conductance: (i) The Reverse Nonequilibrium MD (RNEMD) method, where heat is inputted and extracted at a constant rate from opposite metals. In this case, the thermal conductance is calculated from the nonequilibrium temperature profile that is created at the junction. (ii) The Approach-to-Equilibrium MD (AEMD) method, with the thermal conductance of the junction obtained from the equilibration dynamics of the metals. In both methods, simulations of alkane chains of a growing size display an approximate length-independence of the thermal conductance, with calculated values matching computational and experimental studies. The RNEMD and AEMD methods offer different insights, and we discuss their benefits and shortcomings. Assessing the potential application of molecular junctions as thermal diodes, alkane junctions are made spatially asymmetric by modifying their contact regions with the bulk, either by using distinct endgroups or by replacing one of the Au contacts with Ag. Anharmonicity is built into the system within the molecular force-field. We find that, while the temperature profile strongly varies (compared with the gold-alkanedithiol-gold junctions) due to these structural modifications, the thermal diode effect is inconsequential in these systems-unless one goes to very large thermal biases. This finding suggests that one should seek molecules with considerable internal anharmonic effects for developing nonlinear thermal devices.

Long-range charge transport in homogeneous and alternating-rigidity chains

We study the interplay of intrinsic-electronic and environmental factors in long-range charge transport across molecular chains with up to N ∼ 80 monomers. We describe the molecular electronic structure of the chain with a tight-binding Hamiltonian. Thermal effects in the form of electron decoherence and inelastic scattering are incorporated with the Landauer–Büttiker probe method. In short chains of up to ten units, we observe the crossover between coherent (tunneling, ballistic) motion and thermally-assisted conduction, with thermal effects enhancing the current beyond the quantum coherent limit. We further show that unconventional (nonmonotonic with size) transport behavior emerges when monomer-to-monomer electronic coupling is made large. In long chains, we identify a different behavior, with thermal effects suppressing the conductance below the coherent-ballistic limit. With the goal to identify a minimal model for molecular chains displaying unconventional and effective long-range transport, we simulate a modular polymer with alternating regions of high and low rigidity. Simulations show that, surprisingly, while charge correlations are significantly affected by structuring environmental conditions, reflecting charge delocalization, the electrical resistance displays an averaging effect, and it is not sensitive to this patterning. We conclude by arguing that efficient long-range charge transport requires engineering both internal electronic parameters and environmental conditions.

Heat conduction in polymer chains: Effect of substrate on the thermal conductance

In standard molecular junctions, a molecular structure is placed between and connected to metal leads. Understanding how mechanical tuning in such molecular junctions can change heat conductance has interesting applications in nanoscale energy transport. In this work, we use nonequilibrium molecular dynamics simulations to address the effect of stretching on the phononic contribution to the heat conduction of molecular junctions consisting of single long-chain alkanes and various metal leads, such as Ag, Au, Cu, Ni, and Pt. The thermal conductance of such junctions is found to be much smaller than the intrinsic thermal conductance of the polymer and significantly depends on the nature of metal leads as expressed by the metal–molecule coupling and metal vibrational density of states. This behavior is expected and reflects the mismatch of phonon spectra at the metal molecule interfaces. As a function of stretching, we find a behavior similar to what was observed earlier [M. Dinpajooh and A. Nitzan, J. Chem. Phys. 153, 164903 (2020)] for pure polymeric structures. At relatively short electrode distances, where the polyethylene chains are compressed, it is found that the thermal conductances of the molecular junctions remain almost constant as one stretches the polymer chains. At critical electrode distances, the thermal conductances start to increase, reaching the values of the fully extended molecular junctions. Similar behaviors are observed for junctions in which several long-chain alkanes are sandwiched between various metal leads. These findings indicate that this behavior under stretching is an intrinsic property of the polymer chain and not significantly associated with the interfacial structures.

Environment-dependent vibrational heat transport in molecular junctions: Rectification, quantum effects, vibrational mismatch

Vibrational heat transport in molecular junctions is a central issue in different contemporary research areas such as chemistry, materials science, mechanical engineering, thermoelectrics, and power generation. Our model system consists of a chain of molecules which are sandwiched between two solids that are maintained at different temperatures. We employ a quantum self-consistent reservoir model, which is built on a generalized quantum Langevin equation, to investigate quantum effects and far from equilibrium conditions on thermal conduction at nanoscale. The present self-consistent reservoir model can easily mimic the phonon-phonon scattering mechanisms. Different thermal environments are modeled as (i) Ohmic, (ii) sub-Ohmic, and (iii) super-Ohmic environments, and their effects are demonstrated for the thermal rectification properties of the system with spring graded or mass graded features. The behavior of heat current across molecular junctions as a function of chain length, temperature gradient, and phonon scattering rates are studied. Further, our analysis reveals the effects of vibrational mismatch between the solids phonon spectra on heat transfer characteristics in molecular junctions for different thermal environments.

Efficient Steady-State Solver for the Hierarchical Equations of Motion Approach: Formulation and Application to Charge Transport through Nanosystems

An iterative approach is introduced, which allows the efficient solution of the hierarchical equations of motion (HEOM) for the steady-state of open quantum systems. The approach combines the method of matrix equations with an efficient preconditioning technique to reduce the numerical effort of solving the HEOM. Illustrative applications to simulate nonequilibrium charge transport in single-molecule junctions demonstrate the performance of the method.

Harmonic chains and the thermal diode effect

Harmonic oscillator chains connecting two harmonic reservoirs at different constant temperatures cannot act as thermal diodes, irrespective of structural asymmetry. However, here we prove that perfectly harmonic junctions can rectify heat once the reservoirs (described by white Langevin noise) are placed under temperature gradients, which are asymmetric at the two sides, an effect that we term "temperature-gradient harmonic oscillator diodes." This nonlinear diode effect results from the additional constraint-the imposed thermal gradient at the boundaries. We demonstrate the rectification behavior based on the exact analytical formulation of steady-state heat transport in harmonic systems coupled to Langevin baths, which can describe quantum and classical transport, both regimes realizing the diode effect under the involved boundary conditions. Our study shows that asymmetric harmonic systems, such as room-temperature hydrocarbon molecules with varying side groups and end groups, or a linear lattice of trapped ions may rectify heat by going beyond simple boundary conditions.

On the definitions and simulations of vibrational heat transport in nanojunctions

Thermal transport through nanosystems is central to numerous processes in chemistry, material sciences, and electrical and mechanical engineering, with classical molecular dynamics as the key simulation tool. Here, we focus on thermal junctions with a molecule bridging two solids that are maintained at different temperatures. The classical steady state heat current in this system can be simulated in different ways, either at the interfaces with the solids, which are represented by thermostats, or between atoms within the conducting molecule. We show that while the latter, intramolecular definition feasibly converges to the correct limit, the molecule-thermostat interface definition is more challenging to converge to the correct result. The problem with the interface definition is demonstrated by simulating heat transport in harmonic and anharmonic one-dimensional chains illustrating unphysical effects such as thermal rectification in harmonic junctions.

Local Atomic Heat Currents and Classical Interference in Single-Molecule Heat Conduction

We consider interference effects in vibrational heat conduction across single-molecule junctions. Previous theoretical descriptions of such effects have relied on the quantum Landauer-type expression for heat transport by harmonic molecules, and such observations are sometimes termed "quantum interference". Here we demonstrate via classical atomistic simulations of heat conduction in benzenedithiol single-molecule junctions that the room-temperature effect is essentially classical. In fact, classical simulations and quantum evaluation of room-temperature heat conduction in these systems yield similar results. Simulations of para-, meta-, and ortho-connected benzenedithiols between gold substrates yield conductions in the order para > ortho > meta, which is similar to trends found in the electronic conduction of these structures. The (essentially classical) interference origin of this ordering is indicated by the similarity of the quantum and classical behaviors of these systems.

Strong Coupling Thermodynamics of Open Quantum Systems

A general thermodynamic framework is presented for open quantum systems in fixed contact with a thermal reservoir. The first and second law are obtained for arbitrary system-reservoir coupling strengths, and including both factorized and correlated initial conditions. The thermodynamic properties are adapted to the generally strong coupling regime, approaching the ones defined for equilibrium, and their standard weak-coupling counterparts as appropriate limits. Moreover, they can be inferred from measurements involving only system observables. Finally, a thermodynamic signature of non-Markovianity is formulated in the form of a negative entropy production rate.