1
|
Chen R, Gibson T, Craven GT. Molecular heat transport across a time-periodic temperature gradient. J Chem Phys 2024; 160:194305. [PMID: 38767255 DOI: 10.1063/5.0204819] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Accepted: 04/23/2024] [Indexed: 05/22/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Renai Chen
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Tammie Gibson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Galen T Craven
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
2
|
Li J, Wang JJ, Segal D. Thermal transport in fullerene-based molecular junctions: molecular dynamics simulations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2024; 36:325901. [PMID: 38688291 DOI: 10.1088/1361-648x/ad459b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Accepted: 04/30/2024] [Indexed: 05/02/2024]
Abstract
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.
Collapse
Affiliation(s)
- Joanna Li
- Department of Physics, University of Toronto, 60 Saint George St., Toronto, Ontario M5S 1A7, Canada
- Division of Engineering Science, University of Toronto, 42 Saint George St., Toronto, Ontario M5S 2E4, Canada
| | - Jonathan J Wang
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 Saint George St., Toronto, Ontario M5S 3H6, Canada
| | - Dvira Segal
- Department of Physics, University of Toronto, 60 Saint George St., Toronto, Ontario M5S 1A7, Canada
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 Saint George St., Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
3
|
Wang JJ, Gerry M, Segal D. Challenges in molecular dynamics simulations of heat exchange statistics. J Chem Phys 2024; 160:074111. [PMID: 38380748 DOI: 10.1063/5.0187357] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Accepted: 01/29/2024] [Indexed: 02/22/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Jonathan J Wang
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 Saint George St., Toronto, Ontario M5S 3H6, Canada
| | - Matthew Gerry
- Department of Physics, University of Toronto, 60 Saint George St., Toronto, Ontario M5S 1A7, Canada
| | - Dvira Segal
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 Saint George St., Toronto, Ontario M5S 3H6, Canada
- Department of Physics, University of Toronto, 60 Saint George St., Toronto, Ontario M5S 1A7, Canada
| |
Collapse
|
4
|
Abraham E, Dinpajooh M, Climent C, Nitzan A. Heat transport with a twist. J Chem Phys 2023; 159:174904. [PMID: 37916592 DOI: 10.1063/5.0171680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2023] [Accepted: 10/06/2023] [Indexed: 11/03/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Ethan Abraham
- Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Mohammadhasan Dinpajooh
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Clàudia Climent
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
5
|
Romero-Bastida M, Armando Martínez-Torres B. Thermal rectification in mass-asymmetric one-dimensional anharmonic oscillator lattices with and without a ballistic spacer. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 36:025302. [PMID: 37783211 DOI: 10.1088/1361-648x/acff32] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Accepted: 10/02/2023] [Indexed: 10/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- M Romero-Bastida
- SEPI ESIME-Culhuacán, Instituto Politécnico Nacional, Av. Santa Ana No.1000, San Francisco Culhuacán, Culhuacán CTM V, Coyoacán, CDMX 04440, Mexico
| | - Brandon Armando Martínez-Torres
- SEPI ESIME-Culhuacán, Instituto Politécnico Nacional, Av. Santa Ana No.1000, San Francisco Culhuacán, Culhuacán CTM V, Coyoacán, CDMX 04440, Mexico
| |
Collapse
|
6
|
Chen R, Dinpajooh M, Nitzan A. Quantum bath augmented stochastic nonequilibrium atomistic simulations for molecular heat conduction. J Chem Phys 2023; 159:134110. [PMID: 37800644 DOI: 10.1063/5.0168117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Accepted: 09/15/2023] [Indexed: 10/07/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Renai Chen
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- Theoretical Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Mohammadhasan Dinpajooh
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- School of Chemistry, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
7
|
Zimbovskaya NA, Nitzan A. Phonon transport along long polymer chains with varying configurations: Effects of phonon scattering. J Chem Phys 2023; 158:234903. [PMID: 37326160 DOI: 10.1063/5.0155486] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2023] [Accepted: 05/30/2023] [Indexed: 06/17/2023] Open
Abstract
Following recent molecular dynamic simulations [M. Dinpajooh and A. Nitzan, J. Chem. Phys. 153, 164903 (2020)], we theoretically analyze how the phonon heat transport along a single polymer chain may be affected by varying the chain configuration. We suggest that phonon scattering controls the phonon heat conduction in strongly compressed (and tangled) chain when multiple random bends act as scattering centers for vibrational phonon modes, which results in the diffusive character of heat transport. As the chain is straightening up, the number of scatterers decreases, and the heat transport acquires nearly ballistic character. To analyze these effects, we introduce a model of a long atomic chain made out of identical atoms where some atoms are put in contact with scatterers and treat the phonon heat transfer through such a system as a multichannel scattering problem. We simulate the changes in the chain configurations by varying the number of the scatterers and mimic a gradual straightening of the chain by a gradual reducing of the number of scatterers attached to the chain atoms. It is demonstrated, in agreement with recently published simulation results, that the phonon thermal conductance shows a threshold-like transition from the limit where nearly all atoms are attached to the scatterers to the opposite limit where the scatterers vanish, which corresponds to a transition from the diffusive to the ballistic phonon transport.
Collapse
Affiliation(s)
- Natalya A Zimbovskaya
- Department of Physics and Electronics, University of Puerto Rico-Humacao, CUH Station, Humacao, Puerto Rico 00791, USA
| | - Abraham Nitzan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
- School of Chemistry, Tel Aviv University, Tel Aviv, Israel
| |
Collapse
|
8
|
Wang JJ, Gong J, McGaughey AJH, Segal D. Simulations of heat transport in single-molecule junctions: Investigations of the thermal diode effect. J Chem Phys 2022; 157:174105. [PMID: 36347668 DOI: 10.1063/5.0125714] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Jonathan J Wang
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 Saint George St., Toronto, Ontario M5S 3H6, Canada
| | - Jie Gong
- Department of Mechanical Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Alan J H McGaughey
- Department of Mechanical Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Dvira Segal
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, 80 Saint George St., Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
9
|
Wei X, Harazinska E, Zhao Y, Zhuang Y, Hernandez R. Thermal Transport through Polymer-Linked Gold Nanoparticles. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2022; 126:18511-18519. [PMID: 36366755 PMCID: PMC9639611 DOI: 10.1021/acs.jpcc.2c05816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2022] [Revised: 09/28/2022] [Indexed: 06/16/2023]
Abstract
Polymer-nanoparticle networks have potential applications in molecular electronics and nanophononics. In this work, we use all-atom molecular dynamics to reveal the fundamental mechanisms of thermal transport in polymer-linked gold nanoparticle (AuNP) dimers at the molecular level. Attachment of the polymers to AuNPs of varying sizes allows the determination of effects from the flexibility of the chains when their ends are not held fixed. We report heat conductance (G) values for six polymers-viz. polyethylene, poly(p-phenylene), polyacene, polyacetylene, polythiophene, and poly(3,4-ethylenedioxythiophene)-that represent a broad range of stiffness. We address the multimode effects of polymer type, AuNP size, polymer chain length, polymer conformation, system temperature, and number of linking polymers on G. The combination of the mechanisms for phonon boundary scattering and intrinsic phonon scattering has a strong effect on G. We find that the values of G are larger for conjugated polymers because of the stiffness in their backbones. They are also larger in the low-temperature region for all polymers owing to the quenching of segmental rotations at low temperature. Our simulations also suggest that the total G is additive as the number of linking polymers in the AuNP dimer increases from 1 to 2 to 3.
Collapse
Affiliation(s)
- Xingfei Wei
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland21218, United States
| | - Ewa Harazinska
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland21218, United States
| | - Yinong Zhao
- Department
of Chemical & Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland21218, United States
| | - Yi Zhuang
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland21218, United States
| | - Rigoberto Hernandez
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland21218, United States
- Department
of Chemical & Biomolecular Engineering, Johns Hopkins University, Baltimore, Maryland21218, United States
- Department
of Materials Science and Engineering, Johns
Hopkins University, Baltimore, Maryland21218, United States
| |
Collapse
|