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.

Self-consistent microscopic treatment of the effects of self-motion of the probe on ionic and dipolar solvation dynamics

A simple but self-consistent microscopic theory for the time dependent solvation energy of both ions and dipoles is presented which includes, for the first time, the details of the self-motion of the probe on its own solvation dynamics. The theory leads to several interesting predictions. The most important of them is that, for dipolar solvation, both the rotational and the translational motions of the dipolar solute probe can significantly accelerate the rate of solvation. In addition, the rotational self-motion of the solute can also give rise to an additional mechanism of nonexponentiality in solvation time correlation functions in otherwise slow liquids. A comparison between the present theoretical predictions and the recent experimental studies of Maroncelli et al. on solvation dynamics of aniline in l-propanol seems to indicate that the said experiments have missed the initial solvent response up to about 45 ps. After mapping the experimental results on the redefined time scale, the theoretical results can explain the experimental results for solvation of aniline in 1-propanol very well. For ionic solvation, the translational motion is significant for light solutes only. For example, for Li+ in water, translational motion speeds up the solvation by about 20%. The present theory demonstrates that in dipolar solvation the partial quenching of the self-motion due to the presence of specific solute-solvent interactions (such as H-bonding) may lead to a much slower solvation than that when the self-motion is present. This point has been discussed. In addition, we present the theoretical results for solvation of aniline in propylene carbonate, Here, the solvation is predicted to be complete within 15-20 ps.

Solvation dynamics in slow, viscous liquids: Application to amides

As the viscosity of a liquid increases rapidly in the supercooled regime, the nature of molecular relaxation can exhibit dynamics rather different from the fast dynamics observed in the normal regime. In this article, we present theoretical studies of solvation dynamics and orientational relaxation in slow liquids. As the local short-range correlations are important in the slow liquids, we have extended our previous theory to take into account the shea-range pair correlations between the polar solute and the dipolar solvent molecules. Application of the generalized theory To the study of solvation dynamics of amide systems gives nice agreement with the experimental results of Maroncelli and co-workers (J. Phys. Chem. 1990, 94, 4929). The theory also provides valuable insight into the orientational relaxation precesses in the viscous liquids.

Diffusion Anomaly in Silicalite and VPI-5 from Molecular Dynamics Simulations

Molecular dynamics (MD) simulations on rigid and flexible framework models of silicalite and a rigid framework model of the aluminophosphate VPI-5 for different sorbate diameters are reported. The sorbate-host interactions are modeled in terms of simple atom-atom Lennard-Jones interactions. The results suggest that the diffusion coefficient exhibits an anomaly as gamma approaches unity. The MD results confirm the existence of a linear regime for sorbate diameters significantly smaller than the channel diameter and an anomalous regime observed for sorbate diameters comparable to the channel diameter. The power spectra obtained by Fourier transformation of the velocity autocorrelation function indicate that there is an increase in the intensity of the low-frequency component for the velocity component parallel to the direction of motion for the sorbate diameter in the anomalous regime. The present results suggest that the diffusion anomaly is observed irrespective of (1) the geometry and topology of the pore structure and (2) the nature of the host material. The results are compared with the work of Derouane and co-workers, who have suggested the existence of ''floating molecules'' on the basis of earlier theoretical and computational approaches.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.

Effect of Alkyl Chain Arrangement on Conformation and Dynamics in a Surfactant Intercalated Layered Double Hydroxide: Spectroscopic Measurements and MD Simulations

The anionic surfactant dodecyl sulfate (DDS) has been intercalated in an Mg-Al layered double hydroxide (LDH). Monolayer and bilayer arrangements of the alkyl chains of the intercalated surfactant can be engineered by tuning the Al/Mg ratio of the LDH. In both arrangements the anionic headgroup of the surfactant is tethered to the LDH sheets, and consequently translational mobility of the chains is absent. The degrees of freedom of the confined alkyl chains are restricted to changes in conformation. The effects of the arrangement of the intercalated surfactant chains on conformational order and dynamics have been,investigated by spectroscopic measurements and molecular dynamics simulations. Infrared, Raman, and C-13 NMR spectroscopies were used to investigate conformation of the alkyl chains in the monolayer and bilayer arrangements and variable contact time cross-polarization magic angle spinning (VCT CPMAS) NMR measurements to probe molecular motion. The alkyl chains in the monolayer arrangement of the intercalated DDS chains showed considerably greater conformational disorder and faster dynamics as compared to chains in the bilayer arrangement, in spite of the fact that the volume available per chain in the monolayer is smaller than that in the bilayer. Atomistic MD simulations of the two arrangements of the intercalated surfactant were carried out using an isothermal-isobaric ensemble. The simulations are able to reproduce the essential results of the experiment-greater conformational disorder and faster dynamics for the alkyl chains in the monolayer arrangement of the intercalated surfactant. The MD simulations show that these results are a consequence of the fact that the nature of conformational disorder in the two arrangements is different. In the monolayer arrangement the alkyl chains can sustain isolated gauche defects, whereas in the bilayer arrangement gauche conformers occur only as part of a kink a gauche(+) trans gauche(-) sequence.

Dynamics of Fractional Brownian Walks

We investigate the dynamics of polymers whose solution configurations are represented by fractional Brownian walks. The calculation of the two dynamical quantities considered here, the longest relaxation time tau(r) and the intrinsic viscosity [eta], is formulated in terms of Langevin equations and is carried out within the continuum approach developed in an earlier paper. Our results for tau(r) and [eta] reproduce known scaling relations and provide reasonable numerical estimates of scaling amplitudes. The possible relevance of the work to the study of globular proteins and other compact polymeric phases is discussed.

Lyapunov exponents and predictability of the tropical coupled ocean-atmosphere system

It has recently been proposed that the broad spectrum of interannual variability in the tropics with a peak around four years results from an interaction between the linear low-frequency oscillatory mode of the coupled system and the nonlinear higher-frequency modes of the system. In this study we determine the Lyapunov exponents of the conceptual model consisting of a nonlinear low-order model coupled to a linear oscillator for various values of the coupling constants.

Molecular-dynamics studies of systems of confined dumbbell molecules

We present the results of molecular-dynamics simulations of systems of dumbbell molecules confined by parallel molecular walls. We have carried out systematic studies of three cases: freezing, steady flows, and stick-slip friction. We find that the molecular orientational degrees of freedom cause the surface layers to deviate from a planar configuration. Nevertheless, steady flows, in a channel as narrow as 15 molecular sizes, display continuum behavior. A range of mechanisms in the dynamics of the freezing of a confined fluid is found, as a function of the wall-fluid interactions and the bond length of the dumbbell molecules. The simple order-disorder transition associated with stick-slip motion in the presence of a layer of monoatomic lubricant molecules is supplanted by more complex behavior due to rotational degrees of freedom of the diatomic molecules.