Ballistic and diffusive vibrational energy transport in molecules

Ballistic Energy Transport via Long Alkyl Chains: A New Initiation Mechanism

In an effort to increase the speed and efficiency of ballistic energy transport via oligomeric chains, we performed measurements of the transport in compounds featuring long alkyl chains of up to 37 methylene units. Compounds of the N3-(CH2)n-COOMe type (denoted as aznME) were synthesized with n = 5, 10, 15, 19, 28, 37 and studied using relaxation-assisted two-dimensional infrared spectroscopy. The speed of the ballistic transport, initiated by the N3 tag excitation, increased ca. 3-fold for the longer chains (n = 19-37) compared to the shorter chains, from 14.7 to 48 Å/ps, in line with an earlier prediction (Nawagamuwage et al. 2021, J. Phys. Chem. B, 125, 7546). Modeling, based on solving numerically the Liouville equation, was capable of reproducing the experimental data only if three wavepackets are included, involving CH2 twisting (Tw), wagging (W), and rocking (Ro) chain bands. The approaches for designing molecular systems featuring a higher speed and efficiency of energy transport are discussed.

Two stage decoherence of optical phonons in long oligomers

Molecular vibrations are generally responsible for chemical energy transport and dissipation in molecular systems. This transport is fast and efficient if energy is transferred by optical phonons in periodic oligomers, but its efficiency is limited by decoherence emerging due to anharmonic interactions with acoustic phonons. Using a general theoretical model, we show that in the most common case of the optical phonon band being narrower than the acoustic bands, decoherence takes place in two stages. The faster stage involves optical phonon multiple forward scattering due to absorption and emission of transverse acoustic phonons, i.e., collective bending modes with a quadratic spectrum; the transport remains ballistic and the speed can be altered. The subsequent slower stage involves phonon backscattering in multiphonon processes involving two or more acoustic phonons resulting in a switch to diffusive transport. If the initially excited optical phonon possesses a relatively small group velocity, then it is accelerated in the first stage due to its transitions to states propagating faster. This theoretical expectation is consistent with the recent measurements of optical phonon transport velocity in alkane chains, increasing with increasing the chain length.

Electron transfer rate modulation with mid-IR in butadiyne-bridged donor-bridge-acceptor compounds

Controlling electron transfer (ET) processes in donor-bridge-acceptor (DBA) compounds by mid-IR excitation can enhance our understanding of the ET dynamics and may find practical applications in molecular sensing and molecular-scale electronics. Alkyne moieties are attractive to serve as ET bridges, as they offer the possibility of fast ET and present convenient vibrational modes to perturb the ET dynamics. Yet, these bridges introduce complexity because of the strong torsion angle dependence of the ET rates and transition dipoles among electronic states and a shallow torsion barrier. In this study, we implemented ultrafast 3-pulse laser spectroscopy to investigate how the ET from the dimethyl aniline (D) electron donor to the N-isopropyl-1,8-napthalimide (NAP) electron acceptor can be altered by exciting the CC stretching mode (νCC) of the butadiyne bridge linking the donor and acceptor. The electron transfer was initiated by electronically exciting the acceptor moiety at 400 nm, followed by vibrational excitation of the alkyne, νCC, and detecting the changes in the absorption spectrum in the visible spectral region. The experiments were performed at different delay times t1 and t2, which are the delays between UV-mid-IR and mid-IR-Vis pulses, respectively. Two sets of torsion-angle conformers were identified, one featuring a very fast mean ET time of 0.63 ps (group A) and another featuring a slower mean ET time of 4.3 ps (group B), in the absence of the mid-IR excitation. TD-DFT calculations were performed to determine key torsion angle dependent molecular parameters, including the electronic and vibrational transition dipoles, transition frequencies, and electronic couplings. To describe the 3-pulse data, we developed a kinetic model that includes a locally excited, acceptor-based S2 state, a charge separated S1 state, and their vibrationally excited counterparts, with either excited νCC (denoted as S1Atr, S1Btr, S2Atr, and S2Btr, where tr stands for the excited triplet bond, νCC) or excited daughter modes of the νCC relaxation (S1Ah, S1Bh, S2Ah, and S2Bh, where h stands for vibrationally hot species). The kinetic model was solved analytically, and the species-associated spectra (SAS) were determined numerically using a matrix approach, treating first the experiments with longer t1 delays and then using the already determined SAS for modeling the experiments with shorter t1 delays. Strong vibronic coupling of νCC and of vibrationally hot states makes the analysis complicated. Nevertheless, the SAS were identified and the ET rates of the vibrationally excited species, S2Atr, S2Btr and S2Bh, were determined. The results show that the ET rate for the S2A species is ca. 1.2-fold slower when the νCC mode is excited. The ET rate for species S2B is slower by ca. 1.3-fold if the compound is vibrationally hot and is essentially unchanged when the νCC mode is excited. The SAS determined for the tr and h species resemble the SAS for their respective precursor species in the 2-pulse transient absorption experiments, which validates the procedure used and the results.

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.

Mapping protein dynamics at high spatial resolution with temperature-jump X-ray crystallography

Understanding and controlling protein motion at atomic resolution is a hallmark challenge for structural biologists and protein engineers because conformational dynamics are essential for complex functions such as enzyme catalysis and allosteric regulation. Time-resolved crystallography offers a window into protein motions, yet without a universal perturbation to initiate conformational changes the method has been limited in scope. Here we couple a solvent-based temperature jump with time-resolved crystallography to visualize structural motions in lysozyme, a dynamic enzyme. We observed widespread atomic vibrations on the nanosecond timescale, which evolve on the submillisecond timescale into localized structural fluctuations that are coupled to the active site. An orthogonal perturbation to the enzyme, inhibitor binding, altered these dynamics by blocking key motions that allow energy to dissipate from vibrations into functional movements linked to the catalytic cycle. Because temperature jump is a universal method for perturbing molecular motion, the method demonstrated here is broadly applicable for studying protein dynamics.

Thermal rectification in mass-asymmetric one-dimensional anharmonic oscillator lattices with and without a ballistic spacer

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.

Combining temperature perturbations with X-ray crystallography to study dynamic macromolecules: A thorough discussion of experimental methods

Temperature is an important state variable that governs the behavior of microscopic systems, yet crystallographers rarely exploit temperature changes to study the structure and dynamics of biological macromolecules. In fact, approximately 90% of crystal structures in the Protein Data Bank were determined under cryogenic conditions, because sample cryocooling makes crystals robust to X-ray radiation damage and facilitates data collection. On the other hand, cryocooling can introduce artifacts into macromolecular structures, and can suppress conformational dynamics that are critical for function. Fortunately, recent advances in X-ray detector technology, X-ray sources, and computational data processing algorithms make non-cryogenic X-ray crystallography easier and more broadly applicable than ever before. Without the reliance on cryocooling, high-resolution crystallography can be combined with various temperature perturbations to gain deep insight into the conformational landscapes of macromolecules. This Chapter reviews the historical reasons for the prevalence of cryocooling in macromolecular crystallography, and discusses its potential drawbacks. Next, the Chapter summarizes technological developments and methodologies that facilitate non-cryogenic crystallography experiments. Finally, the chapter discusses the theoretical underpinnings and practical aspects of multi-temperature and temperature-jump crystallography experiments, which are powerful tools for understanding the relationship between the structure, dynamics, and function of proteins and other biological macromolecules.

Evidence for a short-lived resonance state in enzyme catalysis via rate-equation convolution

At the cellular level, all biological function relies on enzymes to provide catalytic acceleration of essential biochemical processes driving cellular metabolism. The enzyme is presumed to lower the activation energy barrier separating reactants from products, but the precise mechanism remains unresolved. Here we examine the temperature dependence of the enzyme-catalyzed dissociation of p-nitrophenyl-α-D-glucopyranoside (pNPG), a chromogenic analog for maltose, isomaltose, and sucrose disaccharide sugars, into p-nitrophenol (pNP) and glucose (monosaccharide). The enzymes of interest are the wild type and mutant forms of glucosidase MalL produced by the probiotic bacterium Bacillus subtilis. The per-enzyme production rates k(T) for the pNPG→ glucose reaction all show a characteristic temperature profile with an Arrhenius-like (approximately exponential) slow acceleration at low temperatures, rising through a point of inflexion to reach a maximum, then turning over to decline steeply towards zero production at high temperatures. This asymmetric profile is found to be well fitted by convolving an exponential growth function f(T) with a Gaussian temperature distribution g(T) to produce an exponentially modified Gaussian function h(T). To give a physical interpretation of the convolution components, we make the temperature mapping Θ≡T_{ref}-T where T_{ref} marks the temperature at which a given mutant becomes fully denatured (unfolded) and therefore inactive, then convert the convolution components to probability density functions which obey the convolution theorem of statistics. Working in Θ space, we identify f(Θ) as the density function for an Arrhenius-like transition from ground-state A to metastable-state B, and g(Θ) as the Gaussian distribution of offset-temperature fluctuations for the metastable state. By mapping the standard thermodynamic relations for temperature and energy fluctuations to the enzyme frame of reference, we are able to derive an expression for the lifetime for the metastable B state. For the 15 enzyme experiments, we obtain a mean value 〈Δt〉≳(29.0±1.3)×10^{-15}s, in remarkably good agreement with the ∼30-fs estimate for the period of glycosidic bond oscillations extracted from published infrared spectroscopy. We suggest that the metastable B state provides a low-energy target that has the effect of lowering the activation energy barrier by presenting an alternative axis for the reaction coordinate.

Brownian Motion in Optical Tweezers, a Comparison between MD Simulations and Experimental Data in the Ballistic Regime

The four most popular water models in molecular dynamics were studied in large-scale simulations of Brownian motion of colloidal particles in optical tweezers and then compared with experimental measurements in the same time scale. We present the most direct comparison of colloidal polystyrene particle diffusion in molecular dynamics simulations and experimental data on the same time scales in the ballistic regime. The four most popular water models, all of which take into account electrostatic interactions, are tested and compared based on yielded results and resources required. Three different conditions were simulated: a freely moving particle and one in a potential force field with two different strengths based on 1 pN/nm and 10 pN/nm. In all cases, the diameter of the colloidal particle was 50 nm. The acquired data were compared with experimental measurements performed using optical tweezers with position capture rates as high as 125 MHz. The experiments were performed in pure water on polystyrene particles with a 1 μm diameter in special microchannel cells.

Time-resolved spectroscopic mapping of vibrational energy flow in proteins: Understanding thermal diffusion at the nanoscale

Vibrational energy exchange between various degrees of freedom is critical to barrier-crossing processes in proteins. Hemeproteins are well suited for studying vibrational energy exchange in proteins because the heme group is an efficient photothermal converter. The released energy by heme following photoexcitation shows migration in a protein moiety on a picosecond timescale, which is observed using time-resolved ultraviolet resonance Raman spectroscopy. The anti-Stokes ultraviolet resonance Raman intensity of a tryptophan residue is an excellent probe for the vibrational energy in proteins, allowing the mapping of energy flow with the spatial resolution of a single amino acid residue. This Perspective provides an overview of studies on vibrational energy flow in proteins, including future perspectives for both methodologies and applications.

Simulations of heat transport in single-molecule junctions: Investigations of the thermal diode effect

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.

Energy Transport and Its Function in Heptahelical Transmembrane Proteins

Photoproteins such as bacteriorhodopsin (bR) and rhodopsin (Rho) need to effectively dissipate photoinduced excess energy to prevent themselves from damage. Another well-studied seven transmembrane (TM) helices protein is the β₂ adrenergic receptor (β₂AR), a G protein-coupled receptor for which energy dissipation paths have been linked with allosteric communication. To study the vibrational energy transport in the active and inactive states of these proteins, a master equation approach [J. Chem. Phys.2020, 152, 045103] is employed, which uses scaling rules that allow us to calculate energy transport rates solely based on the protein structure. Despite their overall structural similarity, the three 7TM proteins reveal quite different strategies to redistribute excess energy. While bR quickly removes the energy using the TM7 helix as a "lightning rod", Rho exhibits a rather poor energy dissipation, which might eventually require the hydrolysis of the Schiff base between the protein and the retinal chromophore to prevent overheating. Heating the ligand adrenaline of β₂AR, the resulting energy transport network of the protein is found to change significantly upon switching from the active state to the inactive state. While the energy flow may highlight aspects of the inter-residue couplings of β₂AR, it seems not particularly suited to explain allosteric phenomena.

Tracking Energy Transfer across a Platinum Center

Rigid, conjugated alkyne bridges serve as important components in various transition-metal complexes used for energy conversion, charge separation, sensing, and molecular electronics. Alkyne stretching modes have potential for modulating charge separation in donor–bridge–acceptor compounds. Understanding the rules of energy relaxation and energy transfer across the metal center in such compounds can help optimize their electron transfer switching properties. We used relaxation-assisted two-dimensional infrared spectroscopy to track energy transfer across metal centers in platinum complexes featuring a triazole-terminated alkyne ligand of two or six carbons, a perfluorophenyl ligand, and two tri(p-tolyl)phosphine ligands. Comprehensive analyses of waiting-time dynamics for numerous cross and diagonal peaks were performed, focusing on coherent oscillation, energy transfer, and cooling parameters. These observables augmented with density functional theory computations of vibrational frequencies and anharmonic force constants enabled identification of different functional groups of the compounds. Computations of vibrational relaxation pathways and mode couplings were performed, and two regimes of intramolecular energy redistribution are described. One involves energy transfer between ligands via high-frequency modes; the transfer is efficient only if the modes involved are delocalized over both ligands. The energy transport pathways between the ligands are identified. Another regime involves redistribution via low-frequency delocalized modes, which does not lead to interligand energy transport.

Versatile Vibrational Energy Sensors for Proteins

Vibrational energy transfer (VET) is emerging as key mechanism for protein functions, possibly playing an important role for energy dissipation, allosteric regulation, and enzyme catalysis. A deep understanding of VET is required to elucidate its role in such processes. Ultrafast VIS-pump/IR-probe spectroscopy can detect pathways of VET in proteins. However, the requirement of having a VET donor and a VET sensor installed simultaneously limits the possible target proteins and sites; to increase their number we compare six IR labels regarding their utility as VET sensors. We compare these labels in terms of their FTIR, and VET signature in VET donor-sensor dipeptides in different solvents. Furthermore, we incorporated four of these labels in PDZ3 to assess their capabilities in more complex systems. Our results show that different IR labels can be used interchangeably, allowing for free choice of the right label depending on the system under investigation and the methods available.

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.

Competition of Several Energy-Transport Initiation Mechanisms Defines the Ballistic Transport Speed

The ballistic regime of vibrational energy transport in oligomeric molecular chains occurs with a constant, often high, transport speed and high efficiency. Such a transport regime can be initiated by exciting a chain end group with a mid-infrared (IR) photon. To better understand the wavepacket formation process, two chemically identical end groups, azido groups with normal, ¹⁴N₃-, and isotopically substituted, ¹⁵N₃-, nitrogen atoms, were tested for wavepacket initiation in compounds with alkyl chains of n = 5, 10, and 15 methylene units terminated with a carboxylic acid (-a) group, denoted as ¹⁴N₃Cn-a and ¹⁵N₃Cn-a. The transport was initiated by exciting the azido moiety stretching mode, the ν_N≡N tag, at 2100 cm^–1 (¹⁴N₃Cn-a) or 2031 cm^–1 (¹⁵N₃Cn-a). Opposite to the expectation, the ballistic transport speed was found to decrease upon ¹⁴N₃ → ¹⁵N₃ isotope editing. Three mechanisms of the transport initiation of a vibrational wavepacket are described and analyzed. The first mechanism involves the direct formation of a wavepacket via excitation with IR photons of several strong Fermi resonances of the tag mode with the ν_N=N + ν_N–C combination state while each of the combination state components is mixed with delocalized chain states. The second mechanism relies on the vibrational relaxation of an end-group-localized tag into a mostly localized end-group state that is strongly coupled to multiple delocalized states of a chain band. Harmonic mixing of ν_N=N of the azido group with CH₂ wagging states of the chain permits a wavepacket formation within a portion of the wagging band, suggesting a fast transport speed. The third mechanism involves the vibrational relaxation of an end-group-localized mode into chain states. Two such pathways were found for the ν_N≡N initiation: The ν_N=N mode relaxes efficiently into the twisting band states and low-frequency acoustic modes, and the ν_N–C mode relaxes into the rocking band states and low-frequency acoustic modes. The contributions of the three initiation mechanisms in the ballistic energy transport initiated by ν_N≡N tag are quantitatively evaluated and related to the experiment. We conclude that the third mechanism dominates the transport in alkane chains of 5–15 methylene units initiated with the ν_N≡N tag and the wavepacket generated predominantly at the CH₂ twisting band. The isotope effect of the transport speed is attributed to a larger contribution of the faster wavepackets for ¹⁴N₃Cn-a or to the different breadth of the wavepacket within the twisting band. The study offers a systematic description of different transport initiation mechanisms and discusses the requirements and features of each mechanism. Such analysis will be useful for designing novel materials for energy management.

Through bonds or contacts? Mapping protein vibrational energy transfer using non-canonical amino acids

Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.

Harmonic chains and the thermal diode effect

Harmonic oscillator chains connecting two harmonic reservoirs at different constant temperatures cannot act as thermal diodes, irrespective of structural asymmetry. However, here we prove that perfectly harmonic junctions can rectify heat once the reservoirs (described by white Langevin noise) are placed under temperature gradients, which are asymmetric at the two sides, an effect that we term "temperature-gradient harmonic oscillator diodes." This nonlinear diode effect results from the additional constraint-the imposed thermal gradient at the boundaries. We demonstrate the rectification behavior based on the exact analytical formulation of steady-state heat transport in harmonic systems coupled to Langevin baths, which can describe quantum and classical transport, both regimes realizing the diode effect under the involved boundary conditions. Our study shows that asymmetric harmonic systems, such as room-temperature hydrocarbon molecules with varying side groups and end groups, or a linear lattice of trapped ions may rectify heat by going beyond simple boundary conditions.

Unidirectional coherent energy transport via conjugated oligo(p-phenylene) chains

We discovered a way to funnel high-frequency vibrational quanta rapidly and unidirectionally over large distances using oligo(p-phenylene) chains. After mid-IR photon photoexcitation of a -COOH end group, the excess energy is injected efficiently into the chain, forming vibrational wavepackets that propagate freely along the chain. The transport delivers high-energy vibrational quanta with a range of transport speeds reaching 8.6 km/s, which exceeds the speed of sound in common metals (∼5 km/s) and polymers (∼2 km/s). Efficiencies of energy injection into the chain and transport along the chain are found to be very high and dependent on the extent of conjugation across the structure. By tuning the degree of conjugation via electronic doping of the chain, the transport speed and efficiency can be controlled. The study opens avenues for developing materials with controllable energy transport properties for heat management, schemes with efficient energy delivery to hard-to-reach regions, including transport against thermal gradients, and ways for initiating chemical reactions remotely.

On the definitions and simulations of vibrational heat transport in nanojunctions

Thermal transport through nanosystems is central to numerous processes in chemistry, material sciences, and electrical and mechanical engineering, with classical molecular dynamics as the key simulation tool. Here, we focus on thermal junctions with a molecule bridging two solids that are maintained at different temperatures. The classical steady state heat current in this system can be simulated in different ways, either at the interfaces with the solids, which are represented by thermostats, or between atoms within the conducting molecule. We show that while the latter, intramolecular definition feasibly converges to the correct limit, the molecule-thermostat interface definition is more challenging to converge to the correct result. The problem with the interface definition is demonstrated by simulating heat transport in harmonic and anharmonic one-dimensional chains illustrating unphysical effects such as thermal rectification in harmonic junctions.

Mimicking Natural Photosynthesis: Designing Ultrafast Photosensitized Electron Transfer into Multiheme Cytochrome Protein Nanowires

Efficient nanomaterials for artificial photosynthesis require fast and robust unidirectional electron transfer (ET) from photosensitizers through charge-separation and accumulation units to redox-active catalytic sites. We explored the ultrafast time-scale limits of photo-induced charge transfer between a Ru(II)tris(bipyridine) derivative photosensitizer and PpcA, a 3-heme c-type cytochrome serving as a nanoscale biological wire. Four covalent attachment sites (K28C, K29C, K52C, and G53C) were engineered in PpcA enabling site-specific covalent labeling with expected donor-acceptor (DA) distances of 4-8 Å. X-ray scattering results demonstrated that mutations and chemical labeling did not disrupt the structure of the proteins. Time-resolved spectroscopy revealed three orders of magnitude difference in charge transfer rates for the systems with otherwise similar DA distances and the same number of covalent bonds separating donors and acceptors. All-atom molecular dynamics simulations provided additional insight into the structure-function requirements for ultrafast charge transfer and the requirement of van der Waals contact between aromatic atoms of photosensitizers and hemes in order to observe sub-nanosecond ET. This work demonstrates opportunities to utilize multi-heme c-cytochromes as frameworks for designing ultrafast light-driven ET into charge-accumulating biohybrid model systems, and ultimately for mimicking the photosynthetic paradigm of efficiently coupling ultrafast, light-driven electron transfer chemistry to multi-step catalysis within small, experimentally versatile photosynthetic biohybrid assemblies.

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.

Crucial effect of transverse vibrations on the transport through polymer chains

The low temperature transport of electron, or vibrational or electronic exciton toward polymer chains, turns out to be dramatically sensitive to its interaction with transverse acoustic vibrations. We show that this interaction leads to a substantial polaron effect and decoherence, which are generally stronger than those associated with longitudinal vibrations. For site-dependent interactions, transverse phonons form subohmic bath leading to the quantum phase transition accompanied by full suppression of the transport at zero temperature and fast decoherence characterized by temperature dependent rate k₂ ∝ T^3/4 at low temperature, while k₂ ∝ T² for site-independent interactions. The latter dependence was used to interpret recent measurements of temperature dependent vibrational energy transport in polyethylene glycol oligomers.

Low-Temperature Vibrational Energy Transport via PEG Chains

We used relaxation-assisted two-dimensional infrared spectroscopy to study the temperature dependence (10-295 K) of end-to-end energy transport across end-decorated PEG oligomers of various chain lengths. The excess energy was introduced by exciting the azido end-group stretching mode at 2100 cm^-1 (tag); the transport was recorded by observing the asymmetric C═O stretching mode of the succinimide ester end group at 1740 cm^-1. The overall transport involves diffusive steps at the end groups and a ballistic step through the PEG chain. We found that at lower temperatures the through-chain energy transport became faster, while the end-group diffusive transport time and the tag lifetime increase. The modeling of the transport using a quantum Liouville equation linked the observations to the reduction of decoherence rate and an increase of the mean-free-path for the vibrational wavepacket. The energy transport at the end groups slowed down at low temperatures due to the decreased number and efficiency of the anharmonic energy redistribution pathways.

Local Atomic Heat Currents and Classical Interference in Single-Molecule Heat Conduction

We consider interference effects in vibrational heat conduction across single-molecule junctions. Previous theoretical descriptions of such effects have relied on the quantum Landauer-type expression for heat transport by harmonic molecules, and such observations are sometimes termed "quantum interference". Here we demonstrate via classical atomistic simulations of heat conduction in benzenedithiol single-molecule junctions that the room-temperature effect is essentially classical. In fact, classical simulations and quantum evaluation of room-temperature heat conduction in these systems yield similar results. Simulations of para-, meta-, and ortho-connected benzenedithiols between gold substrates yield conductions in the order para > ortho > meta, which is similar to trends found in the electronic conduction of these structures. The (essentially classical) interference origin of this ordering is indicated by the similarity of the quantum and classical behaviors of these systems.

Direct observation of nuclear reorganization driven by ultrafast spin transitions

One of the most basic molecular photophysical processes is that of spin transitions and intersystem crossing between excited states surfaces. The change in spin states affects the spatial distribution of electron density through the spin orbit coupling interaction. The subsequent nuclear reorganization reports on the full extent of the spin induced change in electron distribution, which can be treated similarly to intramolecular charge transfer with effective reaction coordinates depicting the spin transition. Here, single-crystal [Fe^II(bpy)₃](PF₆)₂, a prototypical system for spin crossover (SCO) dynamics, is studied using ultrafast electron diffraction in the single-photon excitation regime. The photoinduced SCO dynamics are resolved, revealing two distinct processes with a (450 ± 20)-fs fast component and a (2.4 ± 0.4)-ps slow component. Using principal component analysis, we uncover the key structural modes, ultrafast Fe–N bond elongations coupled with ligand motions, that define the effective reaction coordinate to fully capture the relevant molecular reorganization.

Electron spin is a fundamental property of molecules, and changes in spin state affect both molecular structure and dynamics. Here, the authors resolve, by ultrafast electron diffraction, the nuclear reorganization stabilizing spin transitions in a [Fe^II(bpy)₃](PF₆)₂ crystal.

Intramolecular vibrational energy redistribution and the quantum ergodicity transition: a phase space perspective

The onset of facile intramolecular vibrational energy flow can be related to features in the connected network of anharmonic resonances in the classical phase space.

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.

Topology, landscapes, and biomolecular energy transport

While ubiquitous, energy redistribution remains a poorly understood facet of the nonequilibrium thermodynamics of biomolecules. At the molecular level, finite-size effects, pronounced nonlinearities, and ballistic processes produce behavior that diverges from the macroscale. Here, we show that transient thermal transport reflects macromolecular energy landscape architecture through the topological characteristics of molecular contacts and the nonlinear processes that mediate dynamics. While the former determines transport pathways via pairwise interactions, the latter reflects frustration within the landscape for local conformational rearrangements. Unlike transport through small-molecule systems, such as alkanes, nonlinearity dominates over coherent processes at even quite short time- and length-scales. Our exhaustive all-atom simulations and novel local-in-time and space analysis, applicable to both theory and experiment, permit dissection of energy migration in biomolecules. The approach demonstrates that vibrational energy transport can probe otherwise inaccessible aspects of macromolecular dynamics and interactions that underly biological function.

Understanding vibrational energy transfer in macromolecules has been challenging to both theory and experiment. Here the authors use non-equilibrium molecular dynamics to reveal the relationship between heat transport in a model peptide, emergent nonlinearity, and the underlying free energy landscape.

Energy Transport Pathways in Proteins: A Non-equilibrium Molecular Dynamics Simulation Study

To facilitate the observation of biomolecular energy transport in real time and with single-residue resolution, recent experiments by Baumann et al. ( Angew. Chem. Int. Ed. 2019 , 58 , 2899 , DOI: 10.1002/anie.201812995 ) have used unnatural amino acids β-(1-azulenyl)alanine (Azu) and azidohomoalanine (Aha) to site-specifically inject and probe vibrational energy in proteins. To aid the interpretation of such experiments, non-equilibrium molecular dynamics simulations of the anisotropic energy flow in proteins TrpZip2 and PDZ3 domains are presented. On this account, an efficient simulation protocol is established that accurately mimics the excitation and probing steps of Azu and Aha. The simulations quantitatively reproduce the experimentally found cooling times of the solvated proteins at room temperature and predict that the cooling slows by a factor 2 below the glass temperature of water. In PDZ3, vibrational energy is shown to travel from the initially excited peptide ligand via a complex network of inter-residue contacts and backbone transport to distal regions of the protein. The supposed connection of these energy transport pathways with pathways of allosteric communication is discussed.

Ultrafast Vibrational Cooling Inside of a Molecular Container

Vibrational cooling of azulene encapsulated in a hemicarcerand molecular container was studied by pump-probe spectroscopy. Within 1.5 ps of excitation of azulene to the S₁ state, rapid internal conversion through a conical intersection leads to the formation of a vibrationally hot (∼1080 K) ground state, the subsequent cooling of which can be monitored by tracking the evolution of the red-shifted hot band at the edge of the ground-state absorption. It was found that the cooling of the hot S₀ state of azulene in the host-guest complex (hemicarceplex) is 2-4 times faster than that in common organic solvents. Such large acceleration points to a high density of matching vibrational modes and efficient mechanical coupling between the guest and the host. The experimental observations were fully corroborated by the results of molecular dynamics simulations.