Heat conduction in polymer chains with controlled end-to-end distance

Molecular Electronic Junctions Achieved High Thermal Switch Ratios in Atomistic Simulations

The development of devices that improve thermal energy management requires thermal regulation with efficiency comparable to the ratios R ∼ 105 in electric regulation. Unfortunately, current materials and devices in thermal regulators have only been reported to achieve R ∼ 10. We use atomistic simulations to demonstrate that Ferrocenyl (Fc) molecules under applied external electric fields can alter charge states and achieve high thermal switch ratios R = Gq/G0, where Gq and G0 are the high and low limiting conductances. When an electric field is applied, Fc molecules are positively charged, and the SAM-Au interfacial interaction is strong, leading to high heat conductance Gq. On the other hand, with no electric field, the Fc molecules are charge neutral and the SAM-Au interfacial interaction is weak, leading to low heat conductance G0. We optimized various design parameters for the device performance, including the Au-to-Au gap distance L, the system operation temperature T, the net charge on Fc molecules q, the Au surface charge number Z, and the SAM number N. We find that Gq can be very large and increases with increasing q, Z, or N, while G0 is near 0 at L > 3.0 nm. As a result, R > 100 was achieved for selected parameter ranges reported here.

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.

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 Transfer Enhancement in Tree-Structured Polymer Linked Gold Nanoparticle Networks

Human brains use a tree-like neuron network for information processing at high efficiency and low energy consumption. Tree-like structures have also been engineered to enhance mass and heat transfer in various applications. In this work, we reveal the heat transfer mechanism in tree-structured polymer linked gold nanoparticle (AuNP) networks using atomistic simulations. We report both upward and downward heat fluxes between root and leaf nodes in tree-structured polyethylene (PE) and poly(p-phenylene) (PPP) linked AuNP networks at tree levels from 1 to 5. We found that the heat conductance increases with an increasing polymer tree level. The heat transfer enhancement is due to the resulting increase in the low-frequency vibrational modes. This and other thermal properties are affected by the location of the AuNPs in the tree. Moreover, complex tree structures with at least five levels were found to be robust in the sense that disabling half of the leaves did not change the overall heat conductance.

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.

Phonon transport along long polymer chains with varying configurations: Effects of phonon scattering

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.

Thermal Transport through Polymer-Linked Gold Nanoparticles

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.

Heat conduction in polymer chains: Effect of substrate on the thermal conductance

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.

Electric Fields Influence Intramolecular Vibrational Energy Relaxation and Line Widths

Intramolecular vibrational energy relaxation (IVR) is fundamentally important to chemical dynamics. We show that externally applied electric fields affect IVR and vibrational line widths by changing the anharmonic couplings and frequency detunings between modes. We demonstrate this effect in benzonitrile for which prior experimental results show a decrease in vibrational line width as a function of applied electric field. We identify three major channels for IVR that depend on electric field. In the dominant channel, the electric field affects the frequency detuning, while in the other two channels, variation of anharmonic couplings as a function of field is the underlying mechanism. Consistent with experimental results, we show that the combination of all channels gives rise to reduced line widths with increasing electric field in benzonitrile. Our results are relevant for controlling IVR with external or internal fields and for gaining a more complete interpretation of line widths of vibrational Stark probes.

Role of the phonon confinement effect and boundary scattering in reducing the thermal conductivity of argon nanowire

The thermal conductivity of model argon nanowires over a wide range of temperatures from 20 K to 70 K has been calculated using the formula obtained by solving the Boltzmann equation and independently by molecular dynamic (MD) simulations. The theoretical predictions for thermal conductivity take into account the effect of phonon confinement and boundary scattering. Two known theoretical approaches were used. The first approach is based on the solution of the Boltzmann equation with given boundary conditions and uses bulk acoustic phonon dispersion and neglects the phonon confinement effect. The second approach includes also the modification of acoustic phonon dispersion due to spatial confinement. In simulations, the square and circular shapes of wire with the transverse size of nanowires from 4.3 nm to 42.9 nm have been considered. It was found that MD simulation results match the theoretical predictions reasonably well. The obtained results showed that the phonon confinement effect influences the thermal conductivity of nanowires, but the dominant factor decreasing the thermal conductivity with the thickness of nanowires is boundary scattering. Moreover, the values of the interface specular parameter indicate that the specular phonon-boundary scattering prevails over diffuse phonon-boundary scattering.

Critical Shear Rate of Polymer-Enhanced Hydraulic Fluids

Many application-relevant fluids exhibit shear thinning, where viscosity decreases with shear rate above some critical shear rate. For hydraulic fluids formulated with polymeric additives, the critical shear rate is a function of the molecular weight and concentration of the polymers. Here we present a model for predicting the critical shear rate and Newtonian viscosity of fluids, with the goal of identifying a fluid that shear thins in a specific range relevant to hydraulic pumps. The model is applied to predict the properties of fluids comprising polyisobutene polymer and polyalphaolefin base oil. The theoretical predictions are validated by comparison to viscosities obtained from experimental measurements and molecular dynamics simulations across many decades of shear rates. Results demonstrate that the molecular weight of the polymer plays a key role in determining the critical shear rate, whereas the concentration of polymer primarily affects the Newtonian viscosity. The simulations are further used to show the molecular origins of shear thinning and critical shear rate. The atomistic simulations and simple model developed in this work can ultimately be used to formulate polymer-enhanced fluids with ideal shear thinning profiles that maximize the efficiency of hydraulic systems.