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

One-Dimensional van der Waals Heterojunction Comprising Carbon Nanotube Half-Wrapped in Boron Nitride Nanotube: Deep Investigation of Thermal Rectification

One-dimensional van der Waals (vdWs) heterostructures are celebrated for their exceptional thermal management capabilities, garnering significant research interest. Consequently, our research focused on the one-dimensional vdWs heterojunction comprising carbon nanotube half-wrapped in boron nitride nanotube (BNCNT), specifically their thermal rectification (TR) properties. We employed non-equilibrium molecular dynamics to explore the TR mechanism and assess the impacts of temperature, strain, and coupling strength on heat flux and TR ratio. Our findings reveal that the backward heat flux demonstrates greater atomic vibration instability, as indicated by mean square displacement (MSD), compared to forward heat flux. This instability leads to a higher concentration of localized phonons, thereby diminishing the backward heat flux and enhancing TR. Additionally, we utilized MSD to shed light on the negative differential thermal resistance phenomenon and the influence of stress on forward and backward heat fluxes. Remarkably, TR ratios reached 344% at 3% strain and 400% at -1% strain. Calculations of phonon density of states revealed a competitive mechanism between in-plane and out-of-plane phonons coupling in the inner carbon nanotube and an overlap degree of out-of-plane phonon spectra between the inner carbon nanotube and outer boron nitride nanotube. This accounts for the differing trends in forward and backward heat fluxes as coupling strength χ increases, with TR ratios exceeding 1000% at χ = 7.5. This study provides vital insights for advancing one-dimensional vdWs thermal rectifiers.

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.

Significant thermal rectification induced by phonon mismatch of functional groups in a single-molecule junction

Single-molecule junctions (SMJs) may bring exotic physical effects. In this work, a significant thermal rectification effect is observed in a cross-dimensional system, comprising a diamond, a single-molecule junction, and a carbon nanotube (CNT). The molecular dynamics simulations indicate that the interfacial thermal resistance varies with the direction of heat flow, the orientation of the crystal planes of the diamond, and the length of the CNT. We find that the thermal rectification ratio escalates with the length of the CNT, achieving a peak value of 730% with the CNT length of 200 nm. A detailed analysis of phonon vibrations suggests that the primary cause of thermal rectification is the mismatched vibrations between the biphenyl and carbonyl groups. This discovery may offer theoretical insights for both the experimental exploration and practical application of SMJs in efficient thermal management strategy for high power and highly integrated chips.

Nonequilibrium thermal rectification at the junction of harmonic chains

A thermal diode or rectifier is a system that transmits heat or energy in one direction better than in the opposite direction. We investigate the influence of the distribution of energy among wave numbers on the diode effect for the junction of two dissimilar harmonic chains. An analytical expression for the diode coefficient, characterizing the difference between heat fluxes through the junction in two directions, is derived. It is shown that the diode coefficient depends on the distribution of energy among wave numbers. For an equilibrium energy distribution, the diode effect is absent, while for non-equilibrium energy distributions the diode effect is observed even though the system is harmonic. We show that the diode effect can be maximized by varying the energy distribution and relative position of spectra of the two harmonic chains. Conditions are formulated under which the system acts as an ideal thermal rectifier, i.e., transmits heat only in one direction. The results obtained are important for understanding the heat transfer in heterogeneous low-dimensional nanomaterials.

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.