Polariton induced conical intersection and berry phase.

We investigate the Polariton induced conical intersection (PICI) created from coupling a diatomic molecule with the quantized photon mode inside an optical cavity, and the corresponding Berry Phase effects. We use the rigorous Pauli-Fierz Hamiltonian to describe the quantum light-matter interactions between a LiF molecule and the cavity, and use the exact quantum propagation to investigate the polariton quantum dynamics. The molecular rotations relative to the cavity polarization direction play a role as the tuning mode of the PICI, resulting in an effective CI even within a diatomic molecule. To clearly demonstrate the dynamical effects of the Berry phase, we construct two additional models that have the same Born-Oppenheimer surface, but the effects of the geometric phase are removed. We find that when the initial wavefunction is placed in the lower polaritonic surface, the Berry phase causes a π phase-shift in the wavefunction after the encirclement around the CI, indicated from the nuclear probability distribution. On the other hand, when the initial wavefunction is placed in the upper polaritonic surface, the geometric phase significantly influences the couplings between polaritonic states and therefore, the population dynamics between them. These BP effects are further demonstrated through the photo-fragment angular distribution. PICI created from the quantized radiation field has the promise to open up new possibilities to modulate photochemical reactivities.

Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model.

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.

The origin of absorptive features in the two-dimensional electronic spectra of rhodopsin.

In rhodopsin, the absorption of a photon causes the isomerization of the 11-cis isomer of the retinal chromophore to its all-trans isomer. This isomerization is known to occur through a conical intersection (CI) and the internal conversion through the CI is known to be vibrationally coherent. Recently measured two-dimensional electronic spectra (2DES) showed dramatic absorptive spectral features at early waiting times associated with the transition through the CI. The common two-state two-mode model Hamiltonian was unable to elucidate the origin of these features. To rationalize the source of these features, we employ a three-state three-mode model Hamiltonian where the hydrogen out-of plane (HOOP) mode and a higher-lying electronic state are included. The 2DES of the retinal chromophore in rhodopsin are calculated and compared with the experiment. Our analysis shows that the source of the observed features in the measured 2DES is the excited state absorption to a higher-lying electronic state and not the HOOP mode.

Spectral line shapes in linear absorption and two-dimensional spectroscopy with skewed frequency distributions.

The effect of Gaussian dynamics on the line shapes in linear absorption and two-dimensional correlation spectroscopy is well understood as the second-order cumulant expansion provides exact spectra. Gaussian solvent dynamics can be well analyzed using slope line analysis of two-dimensional correlation spectra as a function of the waiting time between pump and probe fields. Non-Gaussian effects are not as well understood, even though these effects are common in nature. The interpretation of the spectra, thus far, relies on complex case to case analysis. We investigate spectra resulting from two physical mechanisms for non-Gaussian dynamics, one relying on the anharmonicity of the bath and the other on non-linear couplings between bath coordinates. These results are compared with outcomes from a simpler log-normal dynamics model. We find that the skewed spectral line shapes in all cases can be analyzed in terms of the log-normal model, with a minimal number of free parameters. The effect of log-normal dynamics on the spectral line shapes is analyzed in terms of frequency correlation functions, maxline slope analysis, and anti-diagonal linewidths. A triangular line shape is a telltale signature of the skewness induced by log-normal dynamics. We find that maxline slope analysis, as for Gaussian dynamics, is a good measure of the solvent dynamics for log-normal dynamics.

The Primary Photochemistry of Vision Occurs at the Molecular Speed Limit.

Ultrafast photochemical reactions are initiated by vibronic transitions from the reactant ground state to the excited potential energy surface, directly populating excited-state vibrational modes. The primary photochemical reaction of vision, the isomerization of retinal in the protein rhodopsin, is known to be a vibrationally coherent reaction, but the Franck-Condon factors responsible for initiating the process have been difficult to resolve with conventional time-resolved spectroscopies. Here we employ experimental and theoretical 2D photon echo spectroscopy to directly resolve for the first time the Franck-Condon factors that initiate isomerization on the excited potential energy surface and track the reaction dynamics. The spectral dynamics reveal vibrationally coherent isomerization occurring on the fastest possible time scale, that of a single period of the local torsional reaction coordinate. We successfully model this process as coherent wavepacket motion through a conical intersection on a ∼30 fs time scale, confirming the reaction coordinate as a local torsional coordinate with a frequency of ∼570 cm-1. As a result of spectral features being spread out along two frequency coordinates, we unambiguously assign reactant and product states following passage through the conical intersection, which reveal the key vibronic transitions that initiate the vibrationally coherent photochemistry of vision.

Probing the Interstate Coupling near a Conical Intersection by Optical Spectroscopy.

Conical intersections are points where adiabatic potential energy surfaces cross. The interstate coupling between the potential energy surfaces plays a crucial role in many processes associated with conical intersections. Still no method exists to measure this coupling driving the chemical reactions between the potential energy surfaces involved. In this Letter, using a generic model for photoisomerization, we propose a novel experimental approach to estimate the coupling that mixes the electronic states near a conical intersection. The approach is based on analyzing the vibrational wavepacket of the reactant in the adiabatic ground and excited electronic states. The nuclear wavepacket dynamics are extracted from linear absorption and two-dimensional electronic spectroscopy. Comparing the frequencies of the coupling mode in the adiabatic ground and excited states from models with and without coupling between the potential energy surfaces suggests an experimental tool to determine the interstate coupling.

Hydration Effect on Amide I Infrared Bands in Water: An Interpretation Based on an Interaction Energy Decomposition Scheme.

The sensitivity of some infrared bands to the local environment can be exploited to shed light on the structure and the dynamics of biological systems. In particular, the amide I band, which is specifically related to vibrations within the peptide bonds, can give information on the ternary structure of proteins, and can be used as a probe of energy transfer. In this work, we propose a model to quantitatively interpret the frequency shift on the amide I band of a model peptide induced by the formation of hydrogen bonds in the first solvation shell. This method allows us to analyze to what extent the electrostatic interaction, electronic polarization and charge transfer affect the position of the amide I band. The impact of the anharmoniticy of the pontential energy surface on the hydration induced shift is elucidated as well.

Vibrational energy relaxation of the amide I mode of N-methylacetamide in D2O studied through Born-Oppenheimer molecular dynamics.

The vibrational relaxation of the amide I mode of deuterated N-methylacetamide in D2O solution is studied through nonequilibrium simulations using the semiempirical Born-Oppenheimer molecular dynamics (SEBOMD) approach to describe the whole solute-solvent system. Relaxation pathways and lifetimes are determined using the instantaneous normal mode (INM) analysis. The relaxation of the amide I mode is characterized by three different time scales; most of the excess energy (80%) is redistributed through intramolecular vibrational energy redistribution processes, with a smaller contribution (20%) of intermolecular energy flowing into the solvent. The amide II mode is found to contribute modestly (7%) to the relaxation mechanism. The amide I mode and the total vibrational energy decay curves obtained using SEBOMD and INM are in satisfactory agreement with the experimental measurements.

Instantaneous normal mode analysis of the vibrational relaxation of the amide I mode of alanine dipeptide in water.

Nonequilibrium Molecular Dynamics (MD) simulations coupled to instantaneous normal modes (INMs) analysis are used to study the vibrational relaxation of the acetyl and amino-end amide I modes of the alanine dipeptide (AlaD) molecule dissolved in water (D2O). The INMs are assigned in terms of the equilibrium normal modes using the Effective Atomic Min-Cost algorithm as adapted to make use of the outputs of standard MD packages, a method which is well suited for the description of flexible molecules. The relaxation energy curves of both amide I modes show multiexponential decays, in good agreement with the experimental findings. It is found that ~85%-90% of the energy relaxes through intramolecular vibrational redistribution. The main relaxation pathways are also identified. The rate at which energy is transferred into the solvent is similar for the acetyl-end and amino-end amide I modes. The conformational changes occurring during relaxation are investigated, showing that the populations of the alpha and beta region conformers are altered by energy transfer in such a way that it takes 15 ps for the equilibrium conformational populations to be recovered after the initial excitation of the AlaD molecule.

A method for analyzing the vibrational energy flow in biomolecules in solution.

A method is proposed to analyze the intra- and intermolecular vibrational energy flow occurring in biomolecules in solution during relaxation processes. It is based on the assumption that the total energy exchanged between the vibrational modes is minimal and the global process is essentially statistical. This statistical minimum flow method is shown to provide very useful information about the amount and the rate at which energy is transferred between the individual vibrations of the molecule. To demonstrate the performance of the method, an application is made to the relaxation of the amide I mode of N-methylacetamide-d in aqueous D(2)O solution which yields a detailed quantitative description of the process.