Stochastic simulation of nonequilibrium heat conduction in extended molecular junctions

The effect of temperature oscillations on energy storage rectification in harmonic systems

Rectification, the preferential transport of a current in one direction through a system, has garnered significant attention in molecules because of its importance for controlling thermal and electronic currents at the nanoscale. Here, we report the presence of energy storage rectification effects in a molecular chain. This phenomenon is generated by subjecting a harmonic molecular chain to an oscillating temperature gradient and showing that the energy absorption rate of the system depends on the direction of the gradient. We examine how the energy storage rectification ratios in the chain are affected by the oscillating gradient, asymmetry in the chain, and the system parameters. We find that energy storage rectification can be observed in harmonic lattice structures with time-dependent temperatures and that, correspondingly, anharmonicity is not required to generate this rectification mechanism in such systems.

Molecular heat transport across a time-periodic temperature gradient

The time-periodic modulation of a temperature gradient can alter the heat transport properties of a physical system. Oscillating thermal gradients give rise to behaviors such as modified thermal conductivity and controllable time-delayed energy storage that are not present in a system with static temperatures. Here, we examine how the heat transport properties of a molecular lattice model are affected by an oscillating temperature gradient. We use analytical analysis and molecular dynamics simulations to investigate the vibrational heat flow in a molecular lattice system consisting of a chain of particles connected to two heat baths at different temperatures, where the temperature difference between baths is oscillating in time. We derive expressions for heat currents in this system using a stochastic energetics framework and a nonequilibrium Green's function approach that is modified to treat the nonstationary average energy fluxes. We find that emergent energy storage, energy release, and thermal conductance mechanisms induced by the temperature oscillations can be controlled by varying the frequency, waveform, and amplitude of the oscillating gradient. The developed theoretical approach provides a general framework to describe how vibrational heat transmission through a molecular lattice is affected by temperature gradient oscillations.

Thermal transport in fullerene-based molecular junctions: molecular dynamics simulations

We investigate phonon thermal transport of fullerene-based single-molecule junctions by employing classical molecular dynamics (MD) simulations. We compute the thermal conductances of C60fullerene monomers, dimers, and trimers utilizing three distinct MD methods. We observe the equilibration dynamics in one approach, and employ two other nonequilibrium steady state simulation methods. We discuss technical aspects of each simulation technique, and show that their predictions for the thermal conductance agree. Our simulations reveal that while the thermal conductance of fullerene monomer and dimer junctions remains similar, that of trimer junctions experiences a significant reduction. This study could assist in the design of high-performing thermoelectric junctions, where low thermal conductance is desired.

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.

Heat transport with a twist

Despite the desirability of polymers for use in many products due to their flexibility, light weight, and durability, their status as thermal insulators has precluded their use in applications where thermal conductors are required. However, recent results suggest that the thermal conductance of polymers can be enhanced and that their heat transport behaviors may be highly sensitive to nanoscale control. Here we use non-equilibrium molecular dynamics simulations to study the effect of mechanical twist on the steady-state thermal conductance across multi-stranded polyethylene wires. We find that a highly twisted double-helical polyethylene wire can display a thermal conductance up to three times that of its untwisted form, an effect which can be attributed to a structural transition in the strands of the double helix. We also find that in thicker wires composed of many parallel strands, adding just one twist can increase its thermal conductance by over 30%. However, we find that unlike stretching a polymer wire, which causes a monotonic increase in thermal conductance, the effect of twist is highly non-monotonic, and certain amounts of twist can actually decrease the thermal conductance. Finally, we apply the Continuous Chirality Measure (CCM) in an attempt to explore the correlation between heat conductance and chirality. The CCM is found to correlate with twist as expected, but we attribute the observed heat transport behaviors to structural factors other than chirality.

Thermal rectification in mass-asymmetric one-dimensional anharmonic oscillator lattices with and without a ballistic spacer

In this work we perform a systematic analysis of various structural parameters that have influence on the thermal rectification effect, i.e. asymmetrical heat flow, and the negative differential thermal resistance -reduction of the heat flux as the applied thermal bias is increased- present in a one-dimensional, segmented mass-graded system consisting of a coupled nearest-neighbor harmonic oscillator lattice (ballistic spacer) and two diffusive leads (modeled by a substrate potential) attached to the lattice at both boundaries. At variance with previous works, we consider the size of the spacer as smaller than that of the leads. Also considered is the case where the leads are connected along the whole length of the oscillator lattice; that is, in the absence of the ballistic spacer. Upon variation of the system's parameters it was determined that the performance of the device, as quantified by the spectral properties, is largely enhanced in the absence of the ballistic spacer for the small system-size limit herein considered.

Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction

Classical molecular dynamics (MD) has been shown to be effective in simulating heat conduction in certain molecular junctions since it inherently takes into account some essential methodological components which are lacking in the quantum Landauer-type transport model, such as many-body full force-field interactions, anharmonicity effects and nonlinear responses for large temperature biases. However, the classical MD reaches its limit in the environments where the quantum effects are significant (e.g. with low-temperatures substrates, presence of extremely high frequency molecular modes). Here, we present an atomistic simulation methodology for molecular heat conduction that incorporates the quantum Bose-Einstein statistics into an "effective temperature" in the form of a modified Langevin equation. We show that the results from such a quasi-classical effective temperature MD method deviates drastically when the baths temperature approaches zero from classical MD simulations and the results converge to the classical ones when the bath approaches the high-temperature limit, which makes the method suitable for full temperature range. In addition, we show that our quasi-classical thermal transport method can be used to model the conducting substrate layout and molecular composition (e.g. anharmonicities, high-frequency modes). Anharmonic models are explicitly simulated via the Morse potential and compared to pure harmonic interactions to show the effects of anharmonicities under quantum colored bath setups. Finally, the chain length dependence of heat conduction is examined for one-dimensional polymer chains placed in between quantum augmented baths.

Energy transport between heat baths with oscillating temperatures

Energy transport is a fundamental physical process that plays a prominent role in the function and performance of myriad systems and technologies. Recent experimental measurements have shown that subjecting a macroscale system to a time-periodic temperature gradient can increase thermal conductivity in comparison to a static temperature gradient. Here, we theoretically examine this mechanism in a nanoscale model by applying a stochastic Langevin framework to describe the energy transport properties of a particle connecting two heat baths with different temperatures, where the temperature difference between baths is oscillating in time. Analytical expressions for the energy flux of each heat bath and for the system itself are derived for the case of a free particle and a particle in a harmonic potential. We find that dynamical effects in the energy flux induced by temperature oscillations give rise to complex energy transport hysteresis effects. The presented results suggest that applying time-periodic temperature modulations is a potential route to control energy storage and release in molecular devices and nanosystems.

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.

Brownian bridges for stochastic chemical processes-An approximation method based on the asymptotic behavior of the backward Fokker-Planck equation

A Brownian bridge is a continuous random walk conditioned to end in a given region by adding an effective drift to guide paths toward the desired region of phase space. This idea has many applications in chemical science where one wants to control the endpoint of a stochastic process-e.g., polymer physics, chemical reaction pathways, heat/mass transfer, and Brownian dynamics simulations. Despite its broad applicability, the biggest limitation of the Brownian bridge technique is that it is often difficult to determine the effective drift as it comes from a solution of a Backward Fokker-Planck (BFP) equation that is infeasible to compute for complex or high-dimensional systems. This paper introduces a fast approximation method to generate a Brownian bridge process without solving the BFP equation explicitly. Specifically, this paper uses the asymptotic properties of the BFP equation to generate an approximate drift and determine ways to correct (i.e., re-weight) any errors incurred from this approximation. Because such a procedure avoids the solution of the BFP equation, we show that it drastically accelerates the generation of conditioned random walks. We also show that this approach offers reasonable improvement compared to other sampling approaches using simple bias potentials.

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.