Recent progress in transient infrared spectroscopy has made it possible to monitor the transient flow of vibrational energy along a peptide helix [V. Botan et al., Proc. Natl. Acad. Sci. U.S.A. 104, 12749 (2007)]. To provide a theoretical description of these experiments, extensive nonequilibrium molecular dynamics simulations of the photoinduced energy transport in a photoswitchable Aib peptide are performed. By calculating the response of the molecule caused by its excitation via optical and infrared pulses as well as temperature jump and stationary heating, it is shown that these methods are equivalent in that they provide approximately the same molecular energy transfer times. The resulting thermal diffusivity of 10 A(2) ps(-1) qualitatively agrees with the results of previous normal mode calculations for proteins and with experimental bulk values (e.g., 14 A(2) ps(-1) for water). To compare to experiment, a new way of approximating the measured signals is suggested which leads to an improved agreement with the experimental results and explains previous discrepancies. To elucidate the mechanism of energy transfer, modifications to the molecular dynamics force field are introduced, which reveal that the energy transfer occurs mainly through the peptide backbone and depends surprisingly little on the force field parametrization. Employing a harmonic model, quantum-mechanical effects are estimated to moderately (about a factor of 2) speed up the energy transport along the peptide.

A systematic molecular dynamics (MD) simulation study of the photoinduced heat transfer from the model peptide N-methylacetamide (NMA) to various solvents is presented, which considers four types of solvent (water, dimethyl sulfoxide, chloroform, and carbon tetrachloride), and in total 24 different force field models for these solvents. To initiate nonstationary energy flow, an initial temperature jump of NMA is assumed and nonequilibrium MD simulations are performed. As expected from simple theoretical models of heat transfer, the cooling process is proportional to the heat capacity C(V) and--to some extent--to the viscosity eta of the solvent. The complex interplay of Coulomb and Lennard-Jones interactions is studied by scaling these interaction energies. The study reveals that realistic changes (< or approximately 10%) of the Lennard-Jones and Coulomb parameters do not change the cooling time considerably. Including polarizibility, on the other hand, appears to enhance the energy dissipation. Moreover, the solvent's internal degrees of freedom may significantly participate in the heat transfer. This is less so for water, which possesses only three high-frequency vibrational modes, but certainly so for the larger solvent molecules dimethyl sulfoxide and chloroform, which possess several low-frequency vibrational modes. For water, the simulated cooling rate is in excellent agreement with experiment, while only qualitative agreement (up to a factor of 2) is found for the other considered solvents. The importance of the force field model and quantum-mechanical effects to correctly describe the cooling process is discussed in some detail.