An official website of the United States government
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.
more »
« less
Energy transfer and thermoelectricity in molecular junctions in non-equilibrated solvents
We consider a molecular junction immersed in a solvent where the electron transfer is dominated by Marcus-type steps. However, the successive nature of the charge transfer through the junction does not imply that the solvent reaches thermal equilibrium throughout the transport. In our previous work [Kirchberg et al., J. Phys. Chem. Lett. 11, 1729 (2020)], we have determined the nonequilibrium distribution of the solvent where its dynamics, expressed by a friction, is considered in two limiting regimes of fast and slow solvent relaxation. In dependence of the nonequilibrium solvent dynamics, we investigate now the electrical, thermal, and thermoelectric properties of the molecular junction. We show that by suitable tuning of the friction, we can reduce the heat dissipation into the solvent and enhance the heat transfer between the electrodes. Interestingly, we find that the Seebeck coefficient grows significantly by adapting the solvent friction in both regimes.
more »
« less
- Award ID(s):
- 1953701
- PAR ID:
- 10363563
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 156
- Issue:
- 9
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- Article No. 094306
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Reverse nonequilibrium molecular dynamics simulations were used to study heat transport in solvated gold interfaces which have been functionalized with a low-molecular weight thiolated polyethylene glycol (PEG). The gold interfaces studied included (111), (110), and (100) facets as well as spherical nanoparticles with radii of 10 and 20 Å. The embedded atom model (EAM) and the polarizable density-readjusted embedded atom model (DR-EAM) were implemented to determine the effect of metal polarizability on heat transport properties. We find that the interfacial thermal conductance values for thiolated PEG-capped interfaces are higher than those for pristine gold interfaces. Hydrogen bonding between the thiolated PEG and solvent differs between planar facets and the nanospheres, suggesting one mechanism for enhanced transfer of energy, while the covalent gold sulfur bond appears to create the largest barrier to thermal conduction. Through analysis of vibrational power spectra, we find an enhanced population at low-frequency heat-carrying modes for the nanospheres, which may also explain the higher mean interfacial thermal conductance (G) value.more » « less
-
Vertically stacked van der Waals (vdW) heterostructures exhibit unique electronic, optical, and thermal properties that can be manipulated by twist-angle engineering. However, the weak phononic coupling at a bilayer interface imposes a fundamental thermal bottleneck for future two-dimensional devices. Using ultrafast electron diffraction, we directly investigated photoinduced nonequilibrium phonon dynamics in MoS2/WS2at 4° twist angle and WSe2/MoSe2heterobilayers with twist angles of 7°, 16°, and 25°. We identified an interlayer heat transfer channel with a characteristic timescale of ~20 picoseconds, about one order of magnitude faster than molecular dynamics simulations assuming initial intralayer thermalization. Atomistic calculations involving phonon-phonon scattering suggest that this process originates from the nonthermal phonon population following the initial interlayer charge transfer and scattering. Our findings present an avenue for thermal management in vdW heterostructures by tailoring nonequilibrium phonon populations.more » « less
-
The realization of single-molecule thermal conductance measurements has driven the need for theoretical tools to describe conduction processes that occur over atomistic length scales. In macroscale systems, the principle that is typically used to understand thermal conductivity is Fourier’s law. At molecular length scales, however, deviations from Fourier’s law are common in part because microscale thermal transport properties typically depend on the complex interplay between multiple heat conduction mechanisms. Here, the thermal transport properties that arise from electron transfer across a thermal gradient in a molecular conduction junction are examined theoretically. We illustrate how transport in a model junction is affected by varying the electronic structure and length of the molecular bridge in the junction as well as the strength of the coupling between the bridge and its surrounding environment. Three findings are of note: First, the transport properties can vary significantly depending on the characteristics of the molecular bridge and its environment; second, the system’s thermal conductance commonly deviates from Fourier’s law; and third, in properly engineered systems, the magnitude of electron hopping thermal conductance is similar to what has been measured in single-molecule devices.more » « less
-
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.more » « less
