Thomas L. C. Jansen

A transferable electrostatic map for solvation effects on amide I vibrations and its application to linear and two-dimensional spectroscopy

A method for modeling infrared solvent shifts using the electrostatic field generated by the solvent is presented. The method is applied to the amide I vibration of N-methyl acetamide. Using ab initio calculations the fundamental frequency, anharmonicity, and the transition dipoles between the three lowest vibrational states are parametrized in terms of the electrostatic field. The generated map, which takes into account the electric field and its gradients at four molecular positions, is tested in a number of common solvents. Agreement of solvent shift and linewidths with experimental Fourier transform infrared (FTIR) data is found to within seven and four wave numbers, respectively, for polar solvents. This shows that in these solvents electrostatic contributions dominate solvation effects and the map is transferable between these types of solvents. The effect of motional narrowing arising from the fast solvent fluctuations is found to be substantial for the FTIR spectra. Also the two-dimensional infrared (2DIR) spectra, simulated using the constructed map, reproduce experimental results very well. The effect of anharmonicity fluctuations on the 2DIR spectra was found to be negligible.

From atomistic modeling to excitation transfer and two-dimensional spectra of the FMO light-harvesting complex

The experimental observation of long-lived quantum coherences in the Fenna-Matthews-Olson (FMO) light-harvesting complex at low temperatures has challenged general intuition in the field of complex molecular systems and provoked considerable theoretical effort in search of explanations. Here we report on room-temperature calculations of the excited-state dynamics in FMO using a combination of molecular dynamics simulations and electronic structure calculations. Thus we obtain trajectories for the Hamiltonian of this system which contains time-dependent vertical excitation energies of the individual bacteriochlorophyll molecules and their mutual electronic couplings. The distribution of energies and couplings is analyzed together with possible spatial correlations. It is found that in contrast to frequent assumptions the site energy distribution is non-Gaussian. In a subsequent step, averaged wave packet dynamics is used to determine the exciton dynamics in the system. Finally, with the time-dependent Hamiltonian, linear and two-dimensional spectra are determined. The thus-obtained linear absorption line shape agrees well with experimental observation and is largely determined by the non-Gaussian site energy distribution. The two-dimensional spectra are in line with what one would expect by extrapolation of the experimental observations at lower temperatures and indicate almost total loss of long-lived coherences.

Two-dimensional spectroscopy of a molecular dimer unveils the effects of vibronic coupling on exciton coherences

The observation of persistent oscillatory signals in multidimensional spectra of protein-pigment complexes has spurred a debate on the role of coherence-assisted electronic energy transfer as a key operating principle in photosynthesis. Vibronic coupling has recently been proposed as an explanation for the long lifetime of the observed spectral beatings. However, photosynthetic systems are inherently complicated, and tractable studies on simple molecular compounds are needed to fully understand the underlying physics. In this work, we present measurements and calculations on a solvated molecular homodimer with clearly resolvable oscillations in the corresponding two-dimensional spectra. Through analysis of the various contributions to the nonlinear response, we succeed in isolating the signal due to inter-exciton coherence. We find that although calculations predict a prolongation of this coherence due to vibronic coupling, the combination of dynamic disorder and vibrational relaxation leads to a coherence decay on a timescale comparable to the electronic dephasing time.

Two-dimensional infrared spectroscopy and ultrafast anisotropy decay of water

We introduce a sparse-matrix algorithm that allows for the simulation of two-dimensional infrared (2DIR) spectra in systems with many coupled chromophores. We apply the method to bulk water, and our results are based on the recently developed ab initio maps for the vibrational Hamiltonian. Qualitative agreement between theory and experiment is found for the 2DIR spectra without the use of any fitting or scaling parameters in the Hamiltonian. The calculated spectra for bulk water are not so different from those for HOD in D(2)O, which we can understand by considering the spectral diffusion time-correlation functions in both cases. We also calculate the ultrafast anisotropy decay, which is dominated by population transfer, finding very good agreement with experiment. Finally, we determine the vibrational excitation diffusion rate, which is more than two orders of magnitude faster than the diffusion of the water molecules themselves.

Waiting Time Dynamics in Two-Dimensional Infrared Spectroscopy

We review recent work on the waiting time dynamics of coherent two-dimensional infrared (2DIR) spectroscopy. This dynamics can reveal chemical and physical processes that take place on the femto- and picosecond time scale, which is faster than the time scale that may be probed by, for example, nuclear magnetic resonance spectroscopy. A large number of chemically relevant processes take place on this time scale. Such processes range from forming and breaking hydrogen bonds and proton transfer to solvent exchange and vibrational population transfer. In typical 2DIR spectra, multiple processes contribute to the waiting time dynamics and the spectra are often congested. This makes the spectra challenging to interpret, and the aid of theoretical models and simulations is often needed. To be useful, such models need to account for all dynamical processes in the sample simultaneously. The numerical integration of the Schrodinger equation (NISE) method has proven to allow for a very general treatment of the dynamical processes. It accounts for both the motional narrowing resulting from solvent-induced frequency fluctuations and population transfer between coupled vibrations. At the same time, frequency shifts arising from chemical-exchange reactions and changes of the transition dipoles because of either non-Condon effects or molecular reorientation are included in the treatment. This method therefore allows for the disentanglement of all of these processes. The NISE method has thus far been successfully applied to study chemical-exchange processes. It was demonstrated that 2DIR is not only sensitive to reaction kinetics but also to the more detailed reaction dynamics. NISE has also been applied to the study of population transfer within the amide I band (CO stretch) and between the amide I and amide II bands (CN stretch and NH bend) in polypeptides. From the amide I studies, it was found that the population transfer can be used to enhance cross-peaks that act as structural markers for beta-sheet structure in proteins. From the amide I/II investigation, it was found that the amide II band and the hydrogen-bond stretch vibration are important parts of the relaxation pathway for the amide I vibration. With the development of simple approximations, it becomes possible to apply the NISE method even to very big systems, such as the OH stretch of bulk water, which can only be described well when large numbers of coupled vibrations are taken into account.

The third- and fifth-order nonlinear Raman response of liquid CS2 calculated using a finite field nonequilibrium molecular dynamics method

A finite field molecular dynamics (MD) method has been developed to calculate the off-resonant Raman response of liquids. The method has been used to calculate the third- and fifth-order optical responses of CS2. From the third-order response, the intensity of third-order cascading processes has been estimated. The calculated ratio between the fifth-order intensity and the intensity of the third-order cascading processes supports experimental observations, claiming that two-dimensional Raman spectra are dominated by third-order cascading processes.

Stochastic Liouville equation simulation of multidimensional vibrational line shapes of trialanine.

The line shapes detected in coherent femtosecond vibrational spectroscopies contain direct signatures of peptide conformational fluctuations through their effect on vibrational frequencies and intermode couplings. These effects are simulated in trialanine using a Greens function solution of a stochastic Liouville equation constructed for four collective bath coordinates (two Ramachandran angles affecting the mode couplings and two diagonal energies). We find that fluctuations of the Ramachandran angles which hardly affect the linear absorption can be effectively probed by two-dimensional spectra. The signal generated at k(1)+k(2)-k(3) is particularly sensitive to such fluctuations.

Interaction induced effects in the nonlinear Raman response of liquid CS2: A finite field nonequilibrium molecular dynamics approach

The third- and fifth-order time-domain Raman responses of liquid carbon disulfide have been calculated, taking local field effects into account through the dipole-induced dipole approximation to the polarizability. The third-order response is shown to be in excellent agreement with experimental data. The calculated two-dimensional shape of the fifth-order response is compared with recently reported experimental observations of what is claimed to be pure fifth-order response. Considerable discrepancies are observed which might be explained by contamination of the experimental results with sequential and especially parallel third-order cascaded Raman response. A new choice of polarization conditions is proposed, which increases the discrimination against these unwanted cascading effects, as compared to the previously discussed fully polarized and magic angle conditions.

Stochastic Liouville equations for hydrogen-bonding fluctuations and their signatures in two-dimensional vibrational spectroscopy of water

The effects of hydrogen-bond forming and breaking kinetics on the linear and coherent third-order infrared spectra of the OH stretch of HOD in D2O are described by Markovian, not necessarily Gaussian, fluctuations and simulated using the stochastic Liouville equations. Slow (0.5 ps) fluctuations are represented by a collective electrostatic coordinate, whereas fast (<100 fs) frequency fluctuations are described using either a second collective electrostatic coordinate or a four-state jump (FSJ) model for hydrogen-bonding configurations. Parameters for both models were obtained using a 1-ns molecular-dynamics trajectory calculated using the TIP4P force field combined with an electrostatic ab initio map. The asymmetry of the photon-echo spectra (larger linewidth on the blue side than on the red side) predicted by the FSJ is in better agreement with recent experiments.

Dissimilar Dynamics of Coupled Water Vibrations

Dissimilar dynamics of coupled stretch vibrations of a water molecule are revealed by two-dimensional IR correlation spectroscopy. These are caused by essentially non-Gaussian fluctuations of the electric field exerted by the environment on the individual OH stretch vibrations. Non-Gaussian statistics of the individual site frequency fluctuations results in distinctively different dephasing of the symmetric and asymmetric eigenmodes. This phenomenon can only be described if the assumption of Gaussian dynamics in the traditional theories is abandoned.