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
Electron hopping heat transport in molecules
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
- Award ID(s):
- 1953701
- PAR ID:
- 10593597
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 17
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
NA (Ed.)Gold-dithiol molecular junctions have been studied both experimentally and theoretically. However, the nature of the gold-thiolate bond as it relates to the solvent has been seldom investigated. It is known that solvents can impact the electronic structure of single molecule junctions, but the correlation between the solvent and dithiol-linked single-molecule junction conductance is not well understood. We study molecular junctions formed with thiol terminated phenylenes from both 1-chloronaphthalene and 1-bromonaphthalene solutions. We find that the most probable conductance and the distribution of conductances are both affected by the solvent. First-principles calculations show that junction conductance depends on the binding configurations (adatom, atop, bridge) of the thiolate on the Au surface as has been shown previously. More importantly, we find that brominated solvents can restrict the binding of thiols to specific Au sites. This mechanism offers new insight into the effects of the solvent environment on covalent bonding in molecular junctions.more » « less
-
Abstract In the present work we theoretically analyze thermoelectric transport in single-molecule junctions (SMJ) characterized by strong interactions between electrons on the molecular linkers and phonons in their nuclear environments where electron hopping between the electrodes and the molecular bridge states predominates in the steady state electron transport. The analysis is based on the modified Marcus theory accounting for the lifetime broadening of the bridge’s energy levels. We show that the reorganization processes in the environment accompanying electron transport may significantly affect SMJ thermoelectric properties both within and beyond linear transport regime. Specifically, we study the effect of environmental phonons on the electron conductance, the thermopower and charge current induced by the temperature gradient applied across the system.more » « less
-
We demonstrate that imidazole based π–π stacked dimers form strong and efficient conductance pathways in single-molecule junctions using the scanning-tunneling microscope-break junction (STM-BJ) technique and density functional theory-based calculations. We first characterize an imidazole-gold contact by measuring the conductance of imidazolyl-terminated alkanes ( im-N-im , N = 3–6). We show that the conductance of these alkanes decays exponentially with increasing length, indicating that the mechanism for electron transport is through tunneling or super-exchange. We also reveal that π–π stacked dimers can be formed between imidazoles and have better coupling than through-bond tunneling. These experimental results are rationalized by calculations of molecular junction transmission using non-equilibrium Green's function formalism. This study verifies the capability of imidazole as a Au-binding ligand to form stable single- and π-stacked molecule junctions at room temperature.more » « less
-
The correlation of electron transfer with molecular conductance ( g : electron transport through single molecules) by Nitzan and others has contributed to a fundamental understanding of single-molecule electronic materials. When an unsymmetric, dipolar molecule spans two electrodes, the possibility exists for different conductance values at equal, but opposite electrode biases. In the device configuration, these molecules serve as rectifiers of the current and the efficiency of the device is given by the rectification ratio (RR = g forward / g reverse ). Experimental determination of the RR is challenging since the orientation of the rectifying molecule with respect to the electrodes and with respect to the electrode bias direction is difficult to establish. Thus, while two different values of g can be measured and a RR calculated, one cannot easily assign each conductance value as being aligned with or opposed to the molecular dipole, and calculations are often required to resolve the uncertainty. Herein, we describe the properties of two isomeric, triplet ground state biradical molecules that serve as constant-bias analogs of single-molecule electronic devices. Through established theoretical relationships between g and electronic coupling, H 2 , and between H 2 and magnetic exchange coupling, J ( g ∝ H 2 ∝ J ), we use the ratio of experimental J -values for our two isomers to calculate a RR for an unsymmetric bridge molecule with known geometry relative to the two radical fragments of the molecule and at a spectroscopically-defined potential bias. Our experimental results are compared with device transport calculations.more » « less
