Dynamical solvation effects

Electron and charge transfer reactions play an important role in many chemical and biochemical processes. Dynamic solvation effects, among other factors, can largely contribute to determine the reaction rate of these processes and can be studied either by quantum mechanical or simulation methods. 

All types of time evolution are present in dynamical solvation effects. It is difficult, and perhaps not convenient, to define a general formulation of the Hamiltonian which can be used to treat all the possible cases. It is better to treat separately more homogeneous families of phenomena. The usual classification into three main types adiabatic, impulsive and oscillatory, may be of some help. The time dependence of the phenomenon may remain in the solute, and this will be the main case in our survey, but also in the solvent in both cases the coupling will oblige us to consider the dynamic behaviour of the whole system. We shall limit ourselves here to a selection of phenomena which will be considered in the following contributions for which extensions of the basic equilibrium QM approach are used, mainly phenomena related to spectroscopy. Other phenomena will be considered in the next section. 

The second approach to the approximate description of the dynamic solvation effects is based on the semiempirical account for the time-dependent electrical polarization of the medium in the field of the solute molecule. In this case, the statistical averaging over the solute-solvent intennolecular distances and configurations is presumed before the solution of the SchrOdinger equation for the solute and correspondingly, the solvent is described as a polarizable dielectric continuum. The respective electrostatic solvation energy of a solute molecule is given by the following equation[13]... 

The effective matrix elements Hfj describe only the intramolecular terms associated with the chemical bonding but do not take into account long range and intermolecular interactions. For instance, the dipolar interaction between a solute and the molecules of a polar solvent are not accounted by the plain EHT matrix elements. Since semiempirical methods are much faster, the limitations imposed by the use of a continuum dielectric model for the solvent, which do not provide a good approximation for the immediate solvation shells in the vicinity of the solute or near the solid surface, can be overcome by atomistic quantum mechanical models for the solvent. Dynamic solvation effects can also be included through the semiempirical models. The hybrid QM/MM methods are also a valuable alternative to describe the dynamic effects of solvents on the quantum dynamics of the solute. The dipoles can be either intrinsic or induced. In the case of polar solvents, the electronic part of the dipole moment produced by the kth solvent molecule is f k f),... 

Some properties of excited merocyanine dyes are compiled in Table 36.6. The most frequently studied is merocyanine 540 (Ml). Ml is readily soluble in alcohols and Hpid membranes, making it attractive for many technical applications and biological activity studies.The solvent polarity has a significant influence on Tf and the related quantum yields. - An intrinsic barrier height of Ej, = 3 kj mol" for isomerization in Sj was determined from dynamical solvation effects on the rate of Ml in polar solvents. Environmental effects on radiative and nonradiative transitions have been studied for mero-cyanines in homogeneous and microheterogeneous systems. ... 

Onganer, Y, Yin, M., Bessire, D.R., and Quitevis, E.L., Dynamical solvation effects on the cis-trans isomerization reaction photoisomerization of merocyanine 540 in polar solvents, /. Phys. 

To enable an atomic interpretation of the AFM experiments, we have developed a molecular dynamics technique to simulate these experiments [49], Prom such force simulations rupture models at atomic resolution were derived and checked by comparisons of the computed rupture forces with the experimental ones. In order to facilitate such checks, the simulations have been set up to resemble the AFM experiment in as many details as possible (Fig. 4, bottom) the protein-ligand complex was simulated in atomic detail starting from the crystal structure, water solvent was included within the simulation system to account for solvation effects, the protein was held in place by keeping its center of mass fixed (so that internal motions were not hindered), the cantilever was simulated by use of a harmonic spring potential and, finally, the simulated cantilever was connected to the particular atom of the ligand, to which in the AFM experiment the linker molecule was connected. 

D. Beglov and B. Roux. Dominant solvations effects from the primary shell of hydration Approximation for molecular dynamics simulations. Biopolymers, 35 171-178, 1994. 

W. M. Kwok and D. L. Phillips, Solvation effects and short-time photodissociation dynamics of CH2I2 in solution from resonance Raman spectroscopy. Chem. Phys. Lett. 235(3-4), 260-267 (1995). 

Even at this level of dynamical theory, one is not restricted to considering equilibrium solvation of the gas-phase saddle point or of configurations along the gas-phase reaction path [109, 338-344], and to the extent that the solvent is allowed to affect the choice of the reaction path itself, dynamic (i.e., nonequilibrium) solvation effects begin to appear in the theory. 

The reader is referred to review articles concerned with dynamic solvent effects for further discussion of the interesting issues involved in applying continuum and explicit solvation models to dynamical situations [333,381-385],... 

One should take careful note of the fact that in the nonadiabatic solvation, or frozen solvent" limit, it is the absence of solvation dynamics that is important. But is just this lack that is responsible for the deviation from equilibrium solvation, which instead assumes the dynamics are effective in always maintaining equilibrium. 

At the beginning of this decade, Zewail and coworkers reported a fundamental work of solvation effect on a proton transfer reaction [195]. a-naphthol and n-ammonia molecules were studied in real-time for the reaction dynamics on the number of solvent molecules involved in the proton transfer reaction from alcohol towards the ammonia base. Nanosecond dynamics was observed for n=l and 2, while no evidence for proton transfer was found. For n=3 and 4, proton transfer reaction was measured at pisosecond time scale. The nanosecond dynamics appears to be related to the global cluster behavior. The idea of a critical solvation number required to onset proton transfer... 

Solvation of carbonium ions has also been taken into account. The dynamics of CH5+ solvated by H2 and other molecules have been studied experimentally, and in the case of solvation by H2 through classical AIMD studies (92-94). The solvation effect on the 3c2e bonds leads to freezing of the scrambling of protons within CH5 that could be considered due to hydrogen bonding between o bonds (Figure 25). This solvation leads to a more resolved IR spectrum (93). 

TvaroSka, KoS r and Hricovini in this book). One way to account for the effect of solvent on conforxnation might be to represent the molecule without environmental influences, and then explicitly include the solvent or other environmental molecules in the calculation. While avoiding built-in influences of environment is a satisfying concept, it is difficult to obtain by experiment parameters that lack those influences. Several methods have been used to study solvation effects, including continuum descriptions (24) and the explicit treatment of solvent molecules in Monte Carlo and molecular dynamics simulation. 

The primary area where classical PB equations find application is to biomolecules, whose size for the most part precludes application of quantum chemical methods. The dynamics of such macromolecules in solution is often of particular interest, and considerable work has gone into including PB solvation effects in the dynamics equations (see, for instance, Lu and Luo 2003). Typically, force-field atomic partial charges are used for the primary solute charge distribution. 

The quantity 17(f) is the time-dependent friction kernel. It characterizes the dissipation effects of the solvent motion along the reaction coordinate. The dynamic solute-solvent interactions in the case of charge transfer are analogous to the transient solvation effects manifested in C(t) (see Section II). We assume that the underlying dynamics of the dielectric function for BA and other molecules are similar to the dynamics for the coumarins. Thus we quantify t](t) from the experimental C(t) values using the relationship discussed elsewhere [139], The solution to the GLE is in the form of p(z, t), the probability distribution function. 

The reader is also referred to the innovative nonphotochemical electron transfer studies of Weaver et al. [147], These authors have been exploring dynamical solvent effects on ground state self-exchange kinetics for or-ganometallic compounds. This work has explored many aspects of solvent control on intermediate barrier electron transfer reactions, including the effect on a distribution of solvation times. The experimental C(t) data on various solvents have been incorporated into the theoretical modeling of the ground state electron transfer reactions studied by Weaver et al. [147]. 

The second explanation for the solvent isotope effect arises from the dynamic medium effect . At 25 °C the rotational and translational diffusion of DjO molecules in D20 is some 20% slower than H20 molecules in H20 (Albery, 1975a) the viscosity of D20 is also 20% greater than H20. Hence any reaction which is diffusion controlled will be 20% slower in D20 than in H20. This effect would certainly apply to transition state D in Fig. 3 where in the transition state the leaving group is diffusing away. A similar effect may also apply to the classical SN1 and SN2 transition states, if the rotational diffusion of water molecules to form the solvation shell is part of the motion along the reaction co-ordinate in the transition state. Robertson (Laughton and Robertson, 1959 Heppolette and Robertson, 1961) has indeed correlated solvent isotope effects for both SN1 and SN2 reactions with the relative fluidities of H20 and D20. However, while the correlation shows that this is a possible explanation, it may also be that the temperature variation of the solvent isotope effect and of the relative fluidities just happen to be very similar (see below). 

---