Solvent reorganization entropy of electron transfer in polar solvents

We report the results of molecular dynamics simulations of the solvent reorganization energy of intramolecular electron transfer in a charge-transfer molecule dissolved in water and acetonitrile at varying temperatures. The simulations confirm the prediction of microscopic solvation theories of a positive reorganization entropy in polar solvents. The results of simulations are analyzed in terms of the splitting of the reorganization entropy into the contributions from the solute-solvent interaction and from the alteration of the solvent structure induced by the solute. These two contributions mutually cancel each other, resulting in the reorganization entropy amounting to only a fraction of each component.     

Energetics of electron transfer reactions in polar solvents

The free energy change of an electron transfer reaction in a polar solvent is rigorously analyzed within the framework of the dielectric continuum model. An appropriate expression for the electrostatic energy between the two product ions separated by R is derived. The present result does not support a recent claim by Suppan that, if R is close to the contact distance, the electrostatic energy should be much larger in magnitude than estimated from the usual expression −e²/ε_sR.     

Proton transfer in nanoconfined polar solvents. 1. Free energies and solute position

The reaction free energy curves for a model phenol-amine proton-transfer system in a confined CH3Cl solvent have been calculated by Monte Carlo simulations. The free energy curves, as a function of a collective solvent coordinate, have been obtained for several fixed reaction complex radial positions (based on the center-of-mass). A smooth, hydrophobic spherical cavity was used to confine the solvent, and radii of 10 and 15 A have been considered. Quantum effects associated with the transferring proton have been included by adding the proton zero-point energy to the classical free energy. The results indicate the reaction complex position can be an important component of the reaction coordinate for proton-transfer reactions in nanoconfined solvents.     

Ligand reorganization and activation energies in nonadiabatic electron transfer reactions

The activation energy and ligand reorganization energy for nonadiabatic electron transfer reactions in chemical and biological systems are investigated in this paper. The free energy surfaces and the activation energy are derived exactly in the general case in which the ligand vibration frequencies are not equal. The activation energy is derived by free energy minimization at the transition state. Our formulation leads to the Marcus-Hush [J. Chem. Phys. 24, 979 (1956); 98, 7170 (1994); 28, 962 (1958)] results in the equal-frequency limit and also generalizes the Marcus-Sumi [J. Chem. Phys. 84, 4894 (1986)] model in the context of studying the solvent dynamic effect on electron transfer reactions. It is found that when the ligand vibration frequencies are different, the activation energy derived from the Marcus-Hush formula deviates by 5%-10% from the exact value. If the reduced reorganization energy approximation is introduced in the Marcus-Hush formula, the result is almost exact.     

Anomalous reorganization free energies in optical electron transfer in solution

Autoionization bands are observed in the photoelectron emission spectroscopy of aqueous solutions of cyanometalate complexes (Mn, Fe, W, Mo), anions (NO₃, CIO⁻₄ and cations (Ag⁺ TI⁺. Reorganization free energies for autoionization bands are anomalously low in absolute value (by≈1 eV) in comparison with direct transitions to the continuum. Interpretation is based on potential energy profiles and model calculations for the reorganization free energy.     

Solvent reorganization energy of electron-transfer reactions in polar solvents

A microscopic theory of solvent reorganization energy in polar molecular solvents is developed. The theory represents the solvent response as a combination of the density and polarization fluctuations of the solvent given in terms of the density and polarization structure factors. A fully analytical formulation of the theory is provided for a solute of arbitrary shape with an arbitrary distribution of charge. A good agreement between the analytical procedure and the results of Monte Carlo simulations of model systems is achieved. The reorganization energy splits into the contributions from density fluctuations and polarization fluctuations. The polarization part is dominated by longitudinal polarization response. The density part is inversely proportional to temperature. The dependence of the solvent reorganization energy on the solvent dipole moment and refractive index is discussed.     

p-Dimethylaminobenzonitrile in polar solution. L Time-dependent rate in intramolecular electron transfer reaction

An experimental study in the pico- and nano-second range of the fluorescence of DMABN in propanol solution revealed for a large range of viscosities (20?3×10³cP) a non-exponential decay of the planar excited state disappearing by a twisted charge transfer state formation. These results have been analyzed by a theoretical model of electronic relaxation in the absence of a potential barrier, making evident a time dependence of the reaction rate which is confirmed by the analysis of the appearance of the CT state.     

The reorganization energy of electron transfer in nonpolar solvents: molecular level treatment of the solvent

The intermolecular electron transfer in a solute pair consisting of pyrene and dimethylaniline is investigated in a nonpolar solvent, n-hexane. The earlier elaborated approach [M. Tachiya, J. Phys Chem. 97, 5911 (1993)] is used; this method provides a physically relevant background for separating inertial and inertialess polarization responses for both nonpolarizable and polarizable molecular level simulations. The molecular-dynamics technique was implemented for obtaining the equilibrium ensemble of solvent configurations. The nonpolar solvent, n-hexane, was treated in terms of OPLS-AA parametrization. Solute Lennard-Jones parameters were taken from the same parametrization. Solute charge distributions of the initial and final states were determined using ab initio level [HF/6-31G(d,p)] quantum-chemical calculations. Configuration analysis was performed explicitly taking into account the anisotropic polarizability of n-hexane. It is shown that the Gaussian law well describes calculated distribution functions of the solvent coordinate, therefore, the rate constant of the ET reaction can be characterized by the reorganization energy. Evaluated values of the reorganization energies are in a range of 0.03-0.11 eV and significant contribution (more then 40% of magnitude) comes from anisotropic polarizability. Investigation of the reorganization energy lambda dependence on the solute pair separation distance d revealed unexpected behavior. The dependence has a very sharp peak at the distance d=7 A where solvent molecules are able to penetrate into the intermediate space between the solute pair. The reason for such behavior is clarified. This new effect has a purely molecular origin and cannot be described within conventional continuum solvent models.     

Theory and computation of electron transfer reorganization energies with continuum and molecular solvent models

Contemporary continuum-based models of solvation in polar media are surveyed and assessed, with special focus on non-equilibrium solvation. A new hybrid approach combining molecular-level treatment of inertial solvent response, and inclusion of inertialess solvent response at the continuum level, is presented and illustrated in terms of calculated equilibrium solvation free energies for small molecular ions and reorganization free energies for model dumbbell systems.     

Solvent reorganization of electron transitions in viscous solvents

We develop a model of electron transfer reactions at conditions of nonergodicity when the time of solvent relaxation crosses the observation time window set up by the reaction rate. Solvent reorganization energy of intramolecular electron transfer in a charge-transfer molecule dissolved in water and acetonitrile is studied by molecular dynamics simulations at varying temperatures. We observe a sharp decrease of the reorganization energy at a temperature identified as the temperature of structural arrest due to cage effect, as discussed by the mode-coupling theory. This temperature also marks the onset of the enhancement of translational diffusion relative to rotational relaxation signaling the breakdown of the Stokes-Einstein relation. The change in the reorganization energy at the transition temperature reflects the dynamical arrest of the slow, collective relaxation of the solvent related to the relaxation of the solvent dipolar polarization. An analytical theory proposed to describe this effect agrees well with both the simulations and experimental Stokes shift data. The theory is applied to the analysis of charge-transfer kinetics in a low-temperature glass former. We show that the reorganization energy is substantially lower than its equilibrium value for the low-temperature portion of the data. The theory predicts the possibility of discontinuous changes in the dependence of the electron transfer rate on the free energy gap when the reaction switches between ergodic and nonergodic regimes.     

Molecular level approaches for investigation of electron transfer in nonpolar solvents

The authors extend their previous work published in Leontyev and TachiyaJ. Chem. Phys. 123, 224502 (2005) and study not only forward but also reverse electron transfer between pyrene and dimethylaniline in a nonpolar solvent, n-hexane. The distribution function methodology and molecular dynamics technique adopted in their previous work are used. Two algorithms (I and II) are formulated for obtaining the reorganization energy and the solvation free energy difference in the linear response approximation. The two algorithms are combined with different cutoff schemes and tested for polarizable and nonpolarizable solvent models. Agreement between the results obtained by the two algorithms was achieved only for simulations employing the particle mesh Ewald treatment. It is concluded that algorithm I provides a reliable scheme for evaluation of the reorganization energy and the solvation free energy difference. Moreover, a new algorithm referred to as the G-function algorithm is formulated which does not assume the linear response approximation, and is tested on evaluation of the solvation free energy difference. Agreement between the results from the G-function algorithm and those from algorithms I and II is fairly good, although it depends on the degree of statistical consistency of the simulations. In the case of nonpolar solvents the G-function method has practical importance because, unlike the conventional thermodynamic integration approach, it requires equilibrium molecular configuration ensembles only for the initial and final states of the system.     

Proton transfer in nanoconfined polar solvents. II. Adiabatic proton transfer dynamics

The reaction dynamics for a model phenol-amine proton transfer system in a confined methyl chloride solvent have been simulated by mixed quantum-classical molecular dynamics. In this approach, the proton vibration is treated quantum mechanically (and adiabatically), while the rest of the system is described classically. Nonequilibrium trajectories are used to determine the proton transfer reaction rate constant. The reaction complex and methyl chloride solvent are confined in a smooth, hydrophobic spherical cavity, and radii of 10, 12, and 15 A have been considered. The effects of the cavity radius and the heavy atom (hydrogen bond) distance on the reaction dynamics are considered, and the mechanism of the proton transfer is examined in detail by analysis of the trajectories.     

Polarizability response in polar solvents: molecular-dynamics simulations of acetonitrile and chloroform

The relaxation of the many-body polarizability in liquid acetonitrile and chloroform at room temperature was studied by molecular-dynamics simulations. The collective polarizability induced by intermolecular interactions was included using first- and all-orders dipole-induced-dipole models and calculated considering both molecule-centered and distributed site polarizabilities. The anisotropic response was analyzed using a separation scheme that allows a decomposition of the total response in terms of orientational and collision-induced effects. We found the method effective in approximately separating the contributions of these relaxation mechanisms, although the orientational-collision-induced interference makes a non-negligible contribution to the total response. In both liquids the main contribution to the anisotropic response is due to orientational dynamics, but intermolecular collision-induced (or translational) effects are important, especially at short times. We found that higher-order interaction-induced effects were essentially negligible for both liquids. Larger differences were found between the center-center and site-site models, with the latter showing faster polarizability relaxation and better agreement with experiment. Isotropic and anisotropic spectra were computed from the corresponding time correlation functions. The lowest-frequency contributions are largely suppressed in the isotropic spectra and their overall shape is similar to the purely collision-induced contribution to the anisotropic spectra, but with an amplitude which is smaller by a factor of approximately 5 in acetonitrile and approximately 3 in chloroform.     

The appearance of reorganization of intramolecular low-frequency modes in the kinetics of electron transfer from the second excited state in zinc-porphyrin derivatives

An analysis of theoretical modeling results of ultrafast kinetics of photoinduced intramolecular charge separation from the second excited singlet state in the dyad Zn-tetraphenylporphyrin-aminonaphthalenediimide (Zn-TPP-ANDI) in a solution of toluene is presented. The calculations are performed within the framework of the stochastic multi-channel model, which includes four electron states (the ground, first and second excited singlet states, the state with charge separation), as well as their vibration sublevels corresponding to the excitation of highfrequency intramolecular vibration modes. A bimodal kinetic curve of population of the state with charge separation observed in experiments is quantitatively reproduced. The absolute yield values of the state with charge separation are determined. The results of the modeling show that intramolecular modes make a significant contribution to the reorganization of low-frequency modes. Quantum chemical calculations were performed, determining the degrees of freedom related to the intramolecular slow motion of nuclei of high amplitude in the dyad Zn-TPPANDI on going from the ground state to the state with charge separation.     

Molecular dynamics simulations of porphyrin-dendrimer systems: toward modeling electron transfer in solution

We have performed computational simulations of porphyrin-dendrimer systems--a cationic porphyrin electrostatically associated to a negatively charged dendrimer--using the method of classical molecular dynamics (MD) with an atomistic force field. Previous experimental studies have shown a strong quenching effect of the porphyrin fluorescence that was assigned to electron transfer (ET) from the dendrimer's tertiary amines (Paulo, P. M. R.; Costa, S. M. B. J. Phys. Chem. B 2005, 109, 13928). In the present contribution, we evaluate computationally the role of the porphyrin-dendrimer conformation in the development of a statistical distribution of ET rates through its dependence on the donor-acceptor distance. We started from simulations without explicit solvent to obtain trajectories of the donor-acceptor distance and the respective time-averaged distributions for two dendrimer sizes and different initial configurations of the porphyrin-dendrimer pair. By introducing explicit solvent (water) in our simulations, we were able to estimate the reorganization energy of the medium for the systems with the dendrimer of smaller size. The values obtained are in the range 0.6-1.5 eV and show a linear dependence with the inverse of the donor-acceptor distance, which can be explained by a two-phase dielectric continuum model taking into account the medium heterogeneity provided by the dendrimer organic core. Dielectric relaxation accompanying ET was evaluated from the simulations with explicit solvent showing fast decay times of some tens of femtoseconds and slow decay times in the range of hundreds of femtoseconds to a few picoseconds. The variations of the slow relaxation times reflect the heterogeneity of the dendrimer donor sites which add to the complexity of ET kinetics as inferred from the experimental fluorescence decays.