Vibrational energy transport in molecules and the statistical properties of vibrational modes

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.

Inhibition of vibrational energy flow within an aromatic scaffold via heavy atom effect

The regulation of intramolecular vibrational energy redistribution (IVR) to influence energy flow within molecular scaffolds provides a way to steer fundamental processes of chemistry, such as chemical reactivity in proteins and design of molecular diodes. Using two-dimensional infrared (2D IR) spectroscopy, changes in the intensity of vibrational cross-peaks are often used to evaluate different energy transfer pathways present in small molecules. Previous 2D IR studies of para-azidobenzonitrile (PAB) demonstrated that several possible energy pathways from the N3 to the cyano-vibrational reporters were modulated by Fermi resonance, followed by energy relaxation into the solvent [Schmitz et al., J. Phys. Chem. A 123, 10571 (2019)]. In this work, the mechanisms of IVR were hindered via the introduction of a heavy atom, selenium, into the molecular scaffold. This effectively eliminated the energy transfer pathway and resulted in the dissipation of the energy into the bath and direct dipole-dipole coupling between the two vibrational reporters. Several structural variations of the aforementioned molecular scaffold were employed to assess how each interrupted the energy transfer pathways, and the evolution of 2D IR cross-peaks was measured to assess the changes in the energy flow. By eliminating the energy transfer pathways through isolation of specific vibrational transitions, through-space vibrational coupling between an azido (N3) and a selenocyanato (SeCN) probe is facilitated and observed for the first time. Thus, the rectification of this molecular circuitry is accomplished through the inhibition of energy flow using heavy atoms to suppress the anharmonic coupling and, instead, favor a vibrational coupling pathway.

Stochastic simulation of nonequilibrium heat conduction in extended molecular junctions

Understanding phononic heat transport processes in molecular junctions is a central issue in the developing field of nanoscale heat conduction. Here, we present a Langevin dynamics simulation framework to investigate heat transport processes in molecular junctions at and beyond the linear response regime and apply it to saturated and unsaturated linear hydrocarbon chains connecting two gold substrates. Thermal boundary conditions represented by Markovian noise and damping are filtered through several (up to four) gold layers to provide a realistic and controllable bath spectral density. Classical simulations using the full universal force field are compared with quantum calculations that use only the harmonic part of this field. The close agreement found at about room temperature between these very different calculations suggests that heat transport at such temperatures is dominated by lower frequency vibrations whose dynamics is described well by classical mechanics. The results obtained for alkanedithiol molecules connecting gold substrates agree with previous quantum calculations based on the Landauer formula and match recent experimental measurements [e.g., thermal conductance around 20 pW/K for alkanedithiols in single-molecule junctions (SMJs)]. Heat conductance simulations on polyynes of different lengths illuminate the effects of molecular conjugation on thermal transport. The difference between alkanes and polyynes is not large but correlates with the larger rigidity and stronger mode localization that characterize the polyyne structure. This computational approach has been recently used [R. Chen, I. Sharony, and A. Nitzan, J. Phys. Chem. Lett. 11, 4261-4268 (2020)] to unveil local atomic heat currents and phononic interference effect in aromatic-ring based SMJs.

Tuning Molecular Vibrational Energy Flow within an Aromatic Scaffold via Anharmonic Coupling

From guiding chemical reactivity in synthesis or protein folding to the design of energy diodes, intramolecular vibrational energy redistribution harnesses the power to influence the underlying fundamental principles of chemistry. To evaluate the ability to steer these processes, the mechanism and time scales of intramolecular vibrational energy redistribution through aromatic molecular scaffolds have been assessed by utilizing two-dimensional infrared (2D IR) spectroscopy. 2D IR cross peaks reveal energy relaxation through an aromatic scaffold from the azido- to the cyano-vibrational reporters in para-azidobenzonitrile (PAB) and para-(azidomethyl)benzonitrile (PAMB) prior to energy relaxation into the solvent. The rates of energy transfer are modulated by Fermi resonances, which are apparent by the coupling cross peaks identified within the 2D IR spectrum. Theoretical vibrational mode analysis allowed the determination of the origins of the energy flow, the transfer pathway, and a direct comparison of the associated transfer rates, which were in good agreement with the experimental results. Large variations in energy-transfer rates, approximately 1.9 ps for PAB and 23 ps for PAMB, illustrate the importance of strong anharmonic coupling, i.e., Fermi resonance, on the transfer pathways. In particular, vibrational energy rectification is altered by Fermi resonances of the cyano- and azido-modes allowing control of the propensity for energy flow.

Molecules and the Eigenstate Thermalization Hypothesis

We review a theory that predicts the onset of thermalization in a quantum mechanical coupled non-linear oscillator system, which models the vibrational degrees of freedom of a molecule. A system of N non-linear oscillators perturbed by cubic anharmonic interactions exhibits a many-body localization (MBL) transition in the vibrational state space (VSS) of the molecule. This transition can occur at rather high energy in a sizable molecule because the density of states coupled by cubic anharmonic terms scales as N³, in marked contrast to the total density of states, which scales as exp(aN), where a is a constant. The emergence of a MBL transition in the VSS is seen by analysis of a random matrix ensemble that captures the locality of coupling in the VSS, referred to as local random matrix theory (LRMT). Upon introducing higher order anharmonicity, the location of the MBL transition of even a sizable molecule, such as an organic molecule with tens of atoms, still lies at an energy that may exceed the energy to surmount a barrier to reaction, such as a barrier to conformational change. Illustrative calculations are provided, and some recent work on the influence of thermalization on thermal conduction in molecular junctions is also discussed.

Thermalization and Thermal Transport in Molecules

The nature and rate of thermal transport through molecular junctions depend on the length over which thermalization occurs. For junctions formed by alkane chains, in which thermalization occurs only slowly, measurements reveal that thermal resistance is controlled by bonding with the substrates, whereas fluorination can introduce thermal resistance within the molecules themselves, although the mechanism remains unclear. Here we present results of quantum-mechanical calculations of elastic and inelastic scattering rates, the length over which thermalization occurs, and thermal conductance in alkane and perfluoroalkane junctions. The contribution to thermalization of quantum effects that give rise to many-body localization (MBL) in isolated molecules is examined. While MBL does not occur due to dephasing, thermalization is typically too slow to establish local temperature if the same molecule in isolation exhibits MBL. The results indicate limitations on the applicability of classical molecular simulations in modeling thermal transport in molecular junctions.