Quantum and non-Markovian effects in the electron transfer reaction dynamics in the Marcus inverted region

Electron transfer reactions in large molecules may often be coupled to both the polar solvent modes and the intramolecular vibrational modes of the molecule. This can give rise to a complex dynamics which may in some systems, like betaine, be controlled more by vibrational rather than by solvent effects. Additionally, a significant contribution from an ultrafast relaxation component in the solvation dynamics may enhance the complexity. To explain the wide range of behavior that has been observed experimentally, Barbara et al. recently proposed that a model of an electron transfer reaction should minimally consist of a low-frequency classical solvent mode (X), a low-frequency vibrational mode (Q), and a high-frequency quantum mode (q) (J. Phys. Chem. 1991, 96, 3728). In the present work, a theoretical study of this model is described. This study generalizes earlier work by including the biphasic solvent response and the dynamics of the low-frequency vibrational mode in the presence of a delocalized, extended reaction zone. A novel Green's function technique has been developed which allowed us to study the non-Markovian dynamics on a multidimensional surface. The contributions from the high-frequency vibrational mode and the ultrafast component in the non-Markovian solvent dynamics are found to be primarily responsible for the dramatic increase in charge transfer rate over the prediction of the classical theories that neglect both these factors. These, along with a large coupling between the reactant and the product states, may combine to render the electron transfer rate both very large and constant over a wide range of solvent relaxation rates. A study on the free energy gap dependence of the electron transfer rate reveals that the rates are sensitive to changes in the quantum frequency particularly when the free energy gap is very large.

Biological water: Femtosecond dynamics of macromolecular hydration

The unique features of a macromolecule and water as a solvent make the issue of solvation unconventional, with questions about the static versus dynamic nature of hydration and the, physics of orientational and translational diffusion at the boundary. For proteins, the hydration shell that covers the surface is critical to the stability of its structure and function. Dynamically speaking, the residence time of water at the surface is a signature of its mobility and binding. With femtosecond time resolution it is possible to unravel the shortest residence times which are key for the description of the hydration layer, static or dynamic. In this article we review these issues guided by experimental studies, from this laboratory, of polar hydration dynamics at the surfaces of two proteins (Subtilisin Carlsberg (SC) and Monellin). The natural probe tryptophan amino acid was used for the interrogation of the dynamics, and for direct comparison we also studied the behavior in bulk water - a complete hydration in 1 ps. We develop a theoretical description of solvation and relate the theory to the experimental observations. In this - theoretical approach, we consider the dynamical equilibrium in the hydration shell, defining the rate processes for breaking and making the transient hydrogen bonds, and the effective friction in the layer which is defined by the translational and orientational motions of water molecules. The relationship between the residence time of water molecules and the observed slow component in solvation dynamics is a direct one. For the two proteins studied, we observed a "bimodal decay" for the hydration correlation function, with two primary relaxation times: ultrafast, typically 1 ps or less, and longer, typically 15-40 ps, and both are related to the residence time at the protein surface, depending on the binding energies. We end by making extensions to studies of the denatured state of the protein, random coils, and the biomimetic micelles, and conclude with our thoughts on the relevance of the dynamics of native structures to their functions.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Molecular dynamics and simulations study on the vibrational and electronic solvatochromism of benzophenone

Solvent plays a key role in diverse physico-chemical and biological processes. Therefore, understanding solute-solvent interactions at the molecular level of detail is of utmost importance. A comprehensive solvatochromic analysis of benzophenone (Bzp) was carried out in various solvents using Raman and electronic spectroscopy, in conjunction with Density Functional Theory (DFT) calculations of supramolecular solute-solvent clusters generated using classical Molecular Dynamics Simulations (c-MDSs). The >C=O stretching frequency undergoes a bathochromic shift with solvent polarity. Interestingly, in protic solvents this peak appears as a doublet: c-MDS and ad hoc explicit solvent ab initio calculations suggest that the lower and higher frequency peaks are associated with the hydrogen bonded and dangling carbonyl group of Bzp, respectively. Additionally, the dangling carbonyl in methanol (MeOH) solvent is 4 cm(-1) blue-shifted relative to acetonitrile solvent, despite their similar dipolarity/polarizability. This suggests that the cybotactic region of the dangling carbonyl group in MeOH is very different from its bulk solvent structure. Therefore, we propose that this blue-shift of the dangling carbonyl originates in the hydrophobic solvation shell around it resulting from extended hydrogen bonding network of the protic solvents. Furthermore, the 1(1)n pi* (band I) and 1(1)pi pi* (band II) electronic transitions show a hypsochromic and bathochromic shift, respectively. In particular, these shifts in protic solvents are due to differences in their excited state-hydrogen bonding mechanisms. Additionally, a linear relationship is obtained for band I and the >C=O stretching frequency (cm(-1)), which suggests that the different excitation wavelengths in band I correspond to different solvation states. Therefore, we hypothesize that the variation in excitation wavelengths in band I could arise from different solvation states leading to varying solvation dynamics. This will have implications for ultrafast processes associated with electron-transfer, charge transfer, and also the photophysical aspects of excited states. (C) 2016 AIP Publishing LLC.

The energy landscape for solvent dynamics in electron transfer reactions: A minimalist model

Energy fluctuations of a solute molecule embedded in a polar solvent are investigated to depict the energy landscape for solvation dynamics. The system is modeled by a charged molecule surrounded by two layers of solvent dipolar molecules with simple rotational dynamics. Individual solvent molecules are treated as simple dipoles that can point toward or away from the central charge (Ising spins). Single-spin-flip Monte Carlo kinetics simulations are carried out in a two-dimensional lattice for different central charges, radii of outer shell, and temperatures. By analyzing the density of states as a function of energy and temperatures, we have determined the existence of multiple freezing transitions. Each of them can be associated with the freezing of a different layer of the solvent. (C) 2002 American Institute of Physics.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.

Electron transfer reaction in the marcus inverted region: role of high frequency vibrational modes

A theoretical study of the dynamics of photo-electron transfer reactions in the Marcus inverted regime is presented. This study is motivated partly by the recent proposal of Barbara et al. (J. Phys. Chem. 96, 3728, 1991) that a minimal model of an electron transfer reaction should consist of a polar solvent mode (X), a low frequency vibrational mode (Q) and one high frequency mode (q). Interplay between these modes may be responsible for the crossover observed in the dynamics from a solvent controlled to a vibrational controlled electron transfer. The following results have been obtained. (i) In the case of slowly relaxing solvents, the proximity of the point of excitation to an effective sink on the excited surface is critical in determining the decay of the reactant population. This is because the Franck-Condon overlap between the reactant ground and the product excited states decreases rapidly with increase in the quantum number of the product vibrational state. (ii) Non-exponential solvation dynamics has an important effect in determining the rates of electron transfer. Especially, a biphasic solvation and a large coupling between the reactant and the product states both may be needed to explain the experimental results. ©1996 American Institute of Physics

Orientational relaxation in dipolar systems: How much do we understand the role of correlations

Dipolar systems, both liquids and solids, constitute a class of naturally abundant systems that are important in all branches of natural science. The study of orientational relaxation provides a powerful method to understand the microscopic properties of these systems and, fortunately, there are many experimental tools to study orientational relaxation in the condensed phases. However, even after many years of intense research, our understanding of orientational relaxation in dipolar systems has remained largely imperfect. A major hurdle towards achieving a comprehensive understanding is the long range and complex nature of dipolar interactions which also made reliable theoretical study extremely difficult. These difficulties have led to the development of continuum model based theories, which although they provide simple, elegant expressions for quantities of interest, are mostly unsatisfactory as they totally neglect the molecularity of inter-molecular interactions. The situation has improved in recent years because of renewed studies, led by computer simulations. In this review, we shall address some of the recent advances, with emphasis on the work done in our laboratory at Bangalore. The reasons for the failure of the continuum model, as revealed by the recent Brownian dynamics simulations of the dipolar lattice, are discussed. The main reason is that the continuum model predicts too fast a decay of the torque-torque correlation function. On the other hand, a perturbative calculation, based on Zwanzig's projection operator technique, provides a fairly satisfactory description of the single particle orientational dynamics for not too strongly polar dipolar systems. A recently developed molecular hydrodynamic theory that properly includes the effects of intermolecular orientational pair correlations provides an even better description of the single-particle orientational dynamics. We also discuss the rank dependence of the dielectric friction. The other topics reviewed here includes dielectric relaxation and solvation dynamics, as they are intimately connected with orientational relaxation. Recent molecular dynamics simulations of the dipolar lattice are also discussed. The main theme of the present review is to understand the effects of intermolecular interactions on orientational relaxation. The presence of strong orientational pair correlation leads to a strong coupling between the single particle and the collective dynamics. This coupling can lead to rich dynamical properties, some of which are detailed here, while a major part remains yet unexplored.

Role of Specific Cations and Water Entropy on the Stability of Branched DNA Motif Structures

DNA three-way junctions (TWJs) are important intermediates in various cellular processes and are the simplest of a family of branched nucleic acids being considered as scaffolds for biomolecular nanotechnology. Branched nucleic acids are stabilized by divalent cations such as Mg2+, presumably due to condensation and neutralization of the negatively charged DNA backbone. However, electrostatic screening effects point to more complex solvation dynamics and a large role of interfacial waters in thermodynamic stability. Here, we report extensive computer simulations in explicit water and salt on a model TWJ and use free energy calculations to quantify the role of ionic character and strength on stability. We find that enthalpic stabilization of the first and second hydration shells by Mg2+ accounts for 1/3 and all of the free energy gain in 50% and pure MgCl2 solutions, respectively. The more distorted DNA molecule is actually destabilized in pure MgCl2 compared to pure NaCl. Notably, the first shell, interfacial waters have very low translational and rotational entropy (i.e., mobility) compared to the bulk, an entropic loss that is overcompensated by increased enthalpy from additional electrostatic interactions with Mg2+. In contrast, the second hydration shell has anomalously high entropy as it is trapped between an immobile and bulklike layer. The nonmonotonic entropic signature and long-range perturbations of the hydration shells to Mg2+ may have implications in the molecular recognition of these motifs. For example, we find that low salt stabilizes the parallel configuration of the three-way junction, whereas at normal salt we find antiparallel configurations deduced from the NMR. We use the 2PT analysis to follow the thermodynamics of this transition and find that the free energy barrier is dominated by entropic effects that result from the decreased surface area of the antiparallel form which has a smaller number of low entropy waters in the first monolayer.

Effects of phase-breaking on long-range charge transfer in DNA: Partially-coherent-tunneling model study

The mechanism of hole charge transfer in DNA of various lengths and sequences is investigated based on a partially coherent tunneling theory (Zhang et al., J Chem Phys 117:4578, 2002), where the effects of phase-breaking in adenine-thymine and guanine-cytosine base pairs are treated on equal foot. This work aims at providing a self-consistent microscopic interpretation for rate experiments on various DNA systems. We will also clarify the condition under which the simple superexchange-mediated-hopping picture is valid, and make some comments on the further development of present theory.

Scattering matrix approach to electronic dephasing in long-range electron transfer

Based on the Buttiker dephasing model, we propose an analytical scattering matrix approach to the long-range electron transfer phenomena. The present efficient scheme smoothly interpolates between the superexchange and the sequential hopping mechanisms. Various properties such as the drastic dephasing-assisted enhancement and turnover behaviors are demonstrated in good agreement with those obtained via the dynamical reduced density-matrix methods. These properties are further elucidated as results of the interplay among the dephasing strength, the tunneling parameter, and the bridge length of the electron transfer system. (C) 2001 American Institute of Physics.

Influência de parâmetros moleculares em funções de correlação temporal na dinâmica de solvatação mecânica

No presente trabalho descrevemos nossos resultados relativos à investigação da dinâmica de solvatação mecânica por meio de simulações por dinâmica molecular, respeitando o regime da resposta linear, em sistemas-modelo de argônio líquido com um soluto monoatômico ou diatômico dissolvido. Estudamos sistematicamente a influência dos parâmetros moleculares dos solutos (tamanho, polarizabilidade) e da densidade frente a vários modelos de solvatação. Funções de Correlação Temporal da Energia de Solvatação foram calculadas com relação à correlações de n-corpos (n = 2; 3) distinguindo interações repulsivas e atrativas para ambos os sistemas líquidos. Também obtivemos segundas derivadas temporais dessas funções referindo-se à parcelas translacionais, rotacionais e roto-translacionais na solução do diatômico. Encontramos que funções de correlação temporal coletivas podem ser razoavelmente bem aproximadas por correlações binárias a densidades baixas e, a densidades altas, correlações ternárias tornam-se mais importantes produzindo um descorrelacionamento mais rápido das funções coletivas devido a efeitos de cancelamento parciais. As funções de correlação para interações repulsivas e atrativas exibem comportamentos dinâmicos independentes do modelo de solvatação devido a fatores de escalonamento linear que afetam apenas as amplitudes das dessas funções de correlação temporal. Em geral, os sistemas com grau de liberdade rotacional apresentam tempos de correlação mais curtos para a dinâmica coletiva e tempos de correlação mais longos para as funções binárias e ternárias. Finalmente, esse estudo mostra que os sistemas contendo o diatômico relaxam-se predominantemente por mecanismos translacionais binários em modelos de solvatação envolvendo alterações apenas na polarizabilidade do soluto, e por mecanismos rotacionais atrativos binários em modelos envolvendo alterações no comprimento de ligação.

Solvatationsdynamik: lokales Relaxationsverhalten glasbildender Systeme

In dieser Arbeit wurde das lokale Relaxationsverhalten niedermolekularer glasbildender Materialien mit der Methode der Solvatationsdynamik untersucht. Bei dieser Methode werden phosphoreszente Farbstoffe als molekulare Sonden für die lokale Dynamik im Glasbildner eingesetzt. Je nach verwendeter Sonde sind unterschiedliche Eigenschaften zugänglich:Mit Farbstoffen, die ihr Dipolmoment bei optischer Anregung nicht ändern, wird eine mechanische Solvatation gemessen, die die lokale viskoelastische Antwort des Materials widerspiegelt. Mit Farbstoffen, die eine Änderung des Dipolmoments bei optischer Anregung aufweisen,sind zusätzlich lokale dielektrische Eigenschaften des Materials zugänglich, die bei endlicher Polarität des Lösungsmittels dominieren. Das heterogene Relaxationsverhalten des Glasbildners 2-Methyltetrahydrofuran wurde quantitativ untersucht.Auf einem Temperaturbereich, auf dem die mittlere Relaxationszeit um mehr als 4 Dekaden variiert, konnte gezeigt werden, daß die Relaxation lokal rein exponentiell verläuft. Zur Abschätzung von Reichweiteneffekten wurden Solvatationsmessungen in einschränkenden Geometrien in Form von porösen Sol-Gel Gläsern mit Porendurchmessern zwischen 7,5-2,5 nm durchgeführt. Der Einfluß der Beschaffenheit der Porenoberfläche wurde durch Vergleichzwischen der Solvatation in nativen und in silanisierten Gläsern untersucht. Es wurde gezeigt, daß auch in den kleinsten verwendeten silanisierten Poren alle Charakteristika des jeweiligen Bulk-Systems erhalten bleiben.In den nativen porösen Gläsern dagegen führt der Einfluß der stark polaren Oberfläche zu einer Änderung der Dynamik.

Dielektrische Relaxation und Retardation lokaler Prozesse in glasbildenden Materialien

Lokale sekundäre Prozesse in glasbildenden Materialien sind wegen ihrer Wechselwirkung mit der primären Relaxation von besonderem Interesse. Für D-Sorbitol wurde die beta-Relaxation mit drei verschiedenen dielektrischen Meßtechniken untersucht und die Ergebnisse miteinander verglichen. Im Gegensatz zu konventionellen Messungen der dielektrischen Retardation, detektiert eine echte dielektrische Relaxationsmessung den Zerfall des elektrischen Feldes unter der Bedingung einer konstanten dielektrischen Verschiebung. Eine weitere dielektrische Relaxationsmethode ist die Solvatationsdynamik. Sie detektiert die dielektrische Relaxation in der direkten Umgebung eines Farbstoffmoleküls. Die Übereinstimmung der lokal ermittelten und der makroskopisch gemittelten Ergebnisse für den Glaszustand weisen darauf hin, daß die Sekundärrelaxation in D-Sorbitol eine räumlich homogene Eigenschaft ist. Im Gegensatz zu beta-Prozessen anderer Materialien, zeigt nur der abgeschreckte Glaszustand von ortho-Terphenyl einen sekundären Prozeß. Es wurde beobachtet, daß die beta-Amplitude beim Ausheilen langsam abnimmt und im Gleichgewichtszustand der Flüssigkeit oberhalb der Glasübergangstemperatur völlig verschwindet. Viele glasbildende Materialien, wie z. B. Salol zeigen keinen dielektrischen beta-Prozeß. Im Gegensatz zu Messungen mit Standardkühlraten, zeigt eine stark abgeschreckte Salol-Probe eine symmetrische dielektrische Sekundärrelaxation. Diese neuartige Eigenschaft von Salol verschwindet irreversibel, wenn sich die Temperatur der Glasübergangstemperatur nähert.

Criteria for quantum coherent transfer of excitations between chromophores in a polar solvent

We show that the quantum decoherence of Forster resonant energy transfer between two optically active molecules can be described by a spin-boson model. This allows us to give quantitative criteria that are necessary for coherent quantum oscillations of excitations between the chromophores. Experimental tests of our results should be possible with flourescent resonant energy transfer (FRET) spectroscopy. Although we focus on the case of protein-pigment complexes our results are also relevant to quantum dots and organic molecules in a dielectric medium. (c) 2006 Elsevier B.V. All rights reserved.